diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index e32130998..92d177456 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -445,7 +445,7 @@ function symba_frag_pos_fragment_velocity(nfrag, m_frag, x_frag, v_r_unit, Lambd v_r_mag(i) = v_r_mag_02(i) - f_vec_02(i) * (v_r_mag_02(i) - v_r_mag_01(i)) / (f_vec_02(i) - f_vec_01(i)) end do - if (norm2(f_vec_2(:)) < TOL) exit + if (norm2(f_vec_02(:)) < TOL) exit write(*,*) 'v_r_mag: ',v_r_mag(1:4) write(*,*) 'KE_after: ',KE_after