j back to i symba_frag_pos

MintonGroup · May 12, 2021 · 9d967c9 · 9d967c9
1 parent 18642ce
commit 9d967c9
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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