diff --git a/src/gravity.F90 b/src/gravity.F90 index e7d0ab8..a5a505a 100644 --- a/src/gravity.F90 +++ b/src/gravity.F90 @@ -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 diff --git a/src/sources.F90 b/src/sources.F90 index baf9b7f..763acf3 100644 --- a/src/sources.F90 +++ b/src/sources.F90 @@ -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 !