gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

GRAVITY, SOURCES: Pass time and dt to gravitational acceleration module.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-03-07 15:47:22 -03:00
parent d63d6bbe6f
commit 42dcce439f
2 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/gravity.F90
View File

@ -212,12 +212,13 @@ module gravity
!
! Arguments:
!
! t, dt - time and the time increment;
! x, y, z - rectangular coordinates;
! gacc - vector of the gravitational acceleration;
!
!===============================================================================
!
subroutine gacc_none(x, y, z, gacc)
subroutine gacc_none(t, dt, x, y, z, gacc)
! local variables are not implicit by default
!
@ -225,6 +226,7 @@ module gravity
! subroutine arguments
!
real(kind=8) , intent(in) :: t, dt
real(kind=8) , intent(in) :: x, y, z
real(kind=8), dimension(3), intent(out) :: gacc
!
@ -256,12 +258,13 @@ module gravity
!
! Arguments:
!
! t, dt - time and the time increment;
! x, y, z - rectangular coordinates;
! gacc - vector of the gravitational acceleration;
!
!===============================================================================
!
subroutine gacc_rayleigh_taylor(x, y, z, gacc)
subroutine gacc_rayleigh_taylor(t, dt, x, y, z, gacc)
! include external procedures and variables
!
@ -273,6 +276,7 @@ module gravity
! subroutine arguments
!
real(kind=8) , intent(in) :: t, dt
real(kind=8) , intent(in) :: x, y, z
real(kind=8), dimension(3), intent(out) :: gacc

2
src/sources.F90
View File

@ -293,7 +293,7 @@ module sources
! get gravitational acceleration components
!
call gravitational_acceleration(x(i), y(j), z(k), ga(1:3))
call gravitational_acceleration(t, dt, x(i), y(j), z(k), ga(1:3))
! calculate the gravitational source terms
!