diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index fc4d50a08..0e882abd6 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -426,10 +426,8 @@ function symba_frag_pos_fragment_velocity(nfrag, m_frag, x_frag, v_r_unit, Lambd f_vec_const(1:3) = Gam(:) + tau(:) f_vec_const(4) = Beta - Lambda - ! Our initial guess for the first 4 fragments and the values of will be based on an equipartition of KE with some random variation + ! The secant method requires two guesses, so we will use small values to start it off v_r_mag_01(:) = 0.0_DP - ! The secant method requires two guesses, so we will use a small random variate to update the initial guesses - !call random_number(v_r_mag_02(:)) v_r_mag_02(:) = 1e-2_DP * sum(v_r_mag(:)) / (nfrag - 4) f_vec_01(:) = f_vec_const(:)