SOURCES: Fix gravitational acceleration from a point mass.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
parent
8ac8248279
commit
176076dce3
@ -208,7 +208,7 @@ module sources
|
|||||||
! local variables
|
! local variables
|
||||||
!
|
!
|
||||||
integer :: i, j, k
|
integer :: i, j, k
|
||||||
real(kind=8) :: r2, gc, gx, gy, gz
|
real(kind=8) :: r2, r3, gc, gx, gy, gz
|
||||||
|
|
||||||
! local arrays
|
! local arrays
|
||||||
!
|
!
|
||||||
@ -250,10 +250,11 @@ module sources
|
|||||||
#if NDIMS == 3
|
#if NDIMS == 3
|
||||||
r2 = x(i) * x(i) + y(j) * y(j) + z(k) * z(k)
|
r2 = x(i) * x(i) + y(j) * y(j) + z(k) * z(k)
|
||||||
#endif /* NDIMS == 3 */
|
#endif /* NDIMS == 3 */
|
||||||
|
r3 = r2 * sqrt(r2)
|
||||||
|
|
||||||
! calculate gravitational acceleration factors
|
! calculate gravitational acceleration factors
|
||||||
!
|
!
|
||||||
gc = gpoint * pdata%q(idn,i,j,k) / max(1.0d-16, r2)
|
gc = gpoint * pdata%q(idn,i,j,k) / max(1.0d-16, r3)
|
||||||
gx = gc * x(i)
|
gx = gc * x(i)
|
||||||
gy = gc * y(j)
|
gy = gc * y(j)
|
||||||
#if NDIMS == 3
|
#if NDIMS == 3
|
||||||
|
Loading…
x
Reference in New Issue
Block a user