From 9d967c9ee387d0cf39d6ec5442075697357f70c6 Mon Sep 17 00:00:00 2001 From: Carlisle Wishard Date: Wed, 12 May 2021 10:12:03 -0400 Subject: [PATCH] j back to i symba_frag_pos --- src/symba/symba_frag_pos.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index 1e150b3ca..02b1c841c 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -432,7 +432,7 @@ subroutine symba_frag_pos_fragment_velocity(m_frag, x_frag, v_r_unit, Lambda, v_ call random_number(v_r_mag_02(:)) v_r_mag_02(:) = sqrt(2 * Lambda / nfrag / m_frag(:)) * (v_r_mag_02(:) + 0.5_DP) - do i = 1, MAXITER + do j = 1, MAXITER do i = 1, 4 f_1_01 = f_1_01 + m_frag(i) * v_r_mag_01(i) * x_frag(1,i) f_1_02 = f_1_02 + m_frag(i) * v_r_mag_02(i) * x_frag(1,i) @@ -452,7 +452,7 @@ subroutine symba_frag_pos_fragment_velocity(m_frag, x_frag, v_r_unit, Lambda, v_ KE_after = 0.0_DP - do j = 1, 4 + do i = 1, 4 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)) KE_after = KE_after + 0.5_DP * m_frag(i) * v_r_mag(i)**2 end do