fixed typos

MintonGroup · Oct 30, 2023 · 7b19426 · 7b19426
1 parent 80d2cfd
commit 7b19426
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/src/swiftest/swiftest_sph.f90 b/src/swiftest/swiftest_sph.f90
@@ -78,7 +78,7 @@ module subroutine swiftest_sph_g_acc_one(GMcb, r_0, phi, theta, rh, c_lm, g_sph,
                                                               - ccss * (dplm * sin(theta) / (rh(3) * sin(phi)) &
                                                                         + plm * (l + 1) * rh(2) / r_mag**2)) ! g_y
                 g_sph(3) = g_sph(3) + GMcb * r_0**l / r_mag**(l + 1) * ccss * (dplm * sin(theta) / sqrt(r_mag**2 - rh(3)**2) &
-                                                                     + plm * (l + 1) * rh(3) / r_mag**2) ! g_z
+                                                                     + plm * (l + 1) * rh(3) / r_mag**2)     ! g_z
             end do
         end do
 

diff --git a/src/whm/whm_kick.f90 b/src/whm/whm_kick.f90
@@ -237,6 +237,7 @@ module subroutine whm_kick_vh_pl(self, nbody_system, param, t, dt, lbeg)
       logical,                      intent(in)    :: lbeg   !! Logical flag indicating whether this is the beginning of the half step or not. 
       ! Internals
       integer(I4B) :: i, npl
+      real(DP) :: tmp ! to check the acceleration kick due to the CB
 
       associate(pl => self, cb => nbody_system%cb)
          npl = self%nbody
@@ -260,6 +261,7 @@ module subroutine whm_kick_vh_pl(self, nbody_system, param, t, dt, lbeg)
          do concurrent(i = 1:npl, pl%lmask(i))
 #endif
             pl%vh(:, i) = pl%vh(:, i) + pl%ah(:, i) * dt
+            ! tmp = .mag.(pl%ah(:, i) - pl%aobl(:, i))
          end do
       end associate