EVOLUTION: Fix intermediate times in evolve_ssprk3_m().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2019-02-08 13:33:11 -02:00
parent 3c649ec71f
commit 603622886b

View File

@ -1880,7 +1880,7 @@ module evolution
! local variables
!
integer :: i
real(kind=8) :: ds, dh
real(kind=8) :: ds
real(kind=8) :: tm, dtm
!
!-------------------------------------------------------------------------------
@ -2034,8 +2034,7 @@ module evolution
!
! prepare times
!
dh = 0.5d+00 * dt
tm = time + dh
tm = time + dt
dtm = ds
! update fluxes
@ -2085,7 +2084,7 @@ module evolution
! prepare times
!
tm = time + dh + i * ds
tm = time + (i - n) * ds
dtm = ds
! update fluxes