diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index cde64d629..ac3f5bd7c 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -410,8 +410,7 @@ function symba_frag_pos_ke_objective_function(v_r_mag, m_frag, v_r_unit, L_lin_t real(DP) :: fnorm !! The objective function result: norm of the vector composed of the tangential momentum and energy !! Minimizing this brings us closer to our objective ! Internals - integer(I4B), parameter :: Neqs = 4 - real(DP), dimension(:), allocatable :: f_vec + real(DP), dimension(4) :: f_vec f_vec(1) = sum(m_frag(:) * v_r_mag(:) * v_r_unit(1,:)) f_vec(2) = sum(m_frag(:) * v_r_mag(:) * v_r_unit(2,:))