From 27f24ab300a54ff775b112d679ba8e06e4ca0e29 Mon Sep 17 00:00:00 2001 From: anand43 Date: Wed, 1 Nov 2023 09:41:31 -0400 Subject: [PATCH] Fixed sign typos --- src/swiftest/swiftest_sph.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/swiftest/swiftest_sph.f90 b/src/swiftest/swiftest_sph.f90 index 02412d207..aed6d4f0a 100644 --- a/src/swiftest/swiftest_sph.f90 +++ b/src/swiftest/swiftest_sph.f90 @@ -72,10 +72,10 @@ module subroutine swiftest_sph_g_acc_one(GMcb, r_0, phi, theta, rh, c_lm, g_sph, ! m > 0 g_sph(1) = g_sph(1) + GMcb * r_0**l / r_mag**(l + 2) * (cssc * m * plm * sin(phi) / sin(theta) & - - ccss * sin(theta) * cos(phi) & + + ccss * sin(theta) * cos(phi) & * (dplm * cos(theta) + plm * (l + 1))) ! g_x - g_sph(2) = g_sph(2) - GMcb * r_0**l / r_mag**(l + 2) * (cssc * m * plm * cos(phi) / sin(theta) & - - ccss * sin(theta) * sin(phi) & + g_sph(2) = g_sph(2) + GMcb * r_0**l / r_mag**(l + 2) * (-1.0_DP * cssc * m * plm * cos(phi) / sin(theta) & + + ccss * sin(theta) * sin(phi) & * (dplm * cos(theta) + plm * (l + 1))) ! g_y g_sph(3) = g_sph(3) + GMcb * r_0**l / r_mag**(l + 2) * ccss * (-1.0_DP * dplm * sin(theta)**2 & + plm * (l + 1) * cos(theta)) ! g_z