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

GRAVITY: Reset acc in gacc_none().

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-03-09 16:28:30 -03:00
parent 322b8a141b
commit 31bf3006b7
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/gravity.F90
View File

@ -216,11 +216,11 @@ module gravity
!
! t, dt - time and the time increment;
! x, y, z - rectangular coordinates;
! gacc - vector of the gravitational acceleration;
! acc - vector of the gravitational acceleration;
!
!===============================================================================
!
subroutine gacc_none(t, dt, x, y, z, gacc)
subroutine gacc_none(t, dt, x, y, z, acc)
! local variables are not implicit by default
!
@ -230,7 +230,7 @@ module gravity
!
real(kind=8) , intent(in) :: t, dt
real(kind=8) , intent(in) :: x, y, z
real(kind=8), dimension(3), intent(out) :: gacc
real(kind=8), dimension(3), intent(out) :: acc
!
!-------------------------------------------------------------------------------
!
@ -240,6 +240,10 @@ module gravity
call start_timer(imc)
#endif /* PROFILE */
! reset gravitational acceleration
!
acc(:) = 0.0d+00
#ifdef PROFILE
! stop accounting time for the gravitational acceleration calculation
!