From ae69ef59c34270e408fc21333b8e7fe5b6f81521 Mon Sep 17 00:00:00 2001 From: anand43 Date: Thu, 16 Nov 2023 11:18:57 -0500 Subject: [PATCH] fixed typo --- src/swiftest/swiftest_sph.f90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/swiftest/swiftest_sph.f90 b/src/swiftest/swiftest_sph.f90 index bce58bfcb..7eff7e16a 100644 --- a/src/swiftest/swiftest_sph.f90 +++ b/src/swiftest/swiftest_sph.f90 @@ -72,8 +72,10 @@ module subroutine swiftest_sph_g_acc_one(GMcb, r_0, phi, theta, rh, c_lm, g_sph, if(abs(cos_theta) < epsilon(0.0_DP)) then cos_theta = 0.0_DP + end if if(abs(sin_theta) < epsilon(0.0_DP)) then sin_theta = 0.0_DP + end if ! call PlmBar_d1(p, p_deriv, l_max, cos_theta) ! Associated Legendre Polynomials and the 1st Derivative call PlmBar(p, l_max, cos_theta) @@ -145,7 +147,7 @@ module subroutine swiftest_sph_g_acc_one(GMcb, r_0, phi, theta, rh, c_lm, g_sph, else fac1 = m * plm / sin_theta end if - + fac2 = plm * (l + m + 1) * sin_theta + plm1 * cos_theta fac3 = fac2 - fac1 r_fac = -GMcb * r_0**l / r_mag**(l + 2)