From edd784a4f9e0232f7f5ab5b6d16aaf30cf6dcefa Mon Sep 17 00:00:00 2001 From: David Minton Date: Thu, 13 May 2021 17:40:57 -0400 Subject: [PATCH] Use fixed size array for the objection function vector --- src/symba/symba_frag_pos.f90 | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) 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,:))