diff --git a/src/swiftest/swiftest_sph.f90 b/src/swiftest/swiftest_sph.f90 index ca31611bc..69c57b5be 100644 --- 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 index 9c6a2a248..93b448967 100644 --- 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