Removed rlim2 calc from non radius version of kick

MintonGroup · Feb 28, 2023 · d04c449 · d04c449
1 parent 150e165
commit d04c449
Showing 1 changed file with 0 additions and 2 deletions.
diff --git a/src/swiftest/swiftest_kick.f90 b/src/swiftest/swiftest_kick.f90
@@ -234,7 +234,6 @@ module subroutine swiftest_kick_getacch_int_all_triangular_pl(npl, nplm, r, Gmas
                   ry = r(2, j) - r(2, i) 
                   rz = r(3, j) - r(3, i) 
                   rji2 = rx**2 + ry**2 + rz**2
-                  rlim2 = (radius(i) + radius(j))**2
                   call swiftest_kick_getacch_int_one_pl(rji2, rx, ry, rz, Gmass(i), Gmass(j), &
                                              ahi(1,i), ahi(2,i), ahi(3,i), ahj(1,j), ahj(2,j), ahj(3,j))
                end do
@@ -252,7 +251,6 @@ module subroutine swiftest_kick_getacch_int_all_triangular_pl(npl, nplm, r, Gmas
                   ry = r(2,j) - r(2,i)
                   rz = r(3,j) - r(3,i)
                   rji2 = rx**2 + ry**2 + rz**2
-                  rlim2 = (radius(i) + radius(j))**2
                   fac = Gmass(j) / (rji2 * sqrt(rji2))
                   acc(1,i) = acc(1,i) + fac * rx
                   acc(2,i) = acc(2,i) + fac * ry