Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
fixed typo
Browse files Browse the repository at this point in the history
  • Loading branch information
anand43 committed Nov 16, 2023
1 parent d8e23dc commit ae69ef5
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion src/swiftest/swiftest_sph.f90
Original file line number Diff line number Diff line change
Expand Up @@ -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)
Expand Down Expand Up @@ -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)
Expand Down

0 comments on commit ae69ef5

Please sign in to comment.