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

EVOLUTION: Rename RK coefficients in evolve_rk34().

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2014-08-21 14:01:27 -03:00
parent f4236e791e
commit b75ad37011
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/evolution.F90
View File

@ -868,8 +868,8 @@ module evolution
! local integration parameters
!
real(kind=8), parameter :: fh = 1.0d+00 / 2.0d+00, fo = 1.0d+00 / 3.0d+00
real(kind=8), parameter :: ft = 2.0d+00 / 3.0d+00, fs = 1.0d+00 / 6.0d+00
real(kind=8), parameter :: b1 = 1.0d+00 / 2.0d+00, b3 = 1.0d+00 / 6.0d+00
real(kind=8), parameter :: a31 = 2.0d+00 / 3.0d+00, a33 = 1.0d+00 / 3.0d+00
!
!-------------------------------------------------------------------------------
!
@ -877,7 +877,7 @@ module evolution
!
! calculate fractional time step
!
ds = fh * dt
ds = b1 * dt
! update fluxes for the first step of the RK2 integration
!
@ -961,7 +961,7 @@ module evolution
!
! calculate fractional time step
!
ds = fs * dt
ds = b3 * dt
! update fluxes for the first step of the RK2 integration
!
@ -982,8 +982,8 @@ module evolution
! update the solution for the fluid variables
!
pblock%u1(1:nv,:,:,:) = ft * pblock%u0(1:nv,:,:,:) &
+ fo * pblock%u1(1:nv,:,:,:) + ds * du(1:nv,:,:,:)
pblock%u1(1:nv,:,:,:) = a31 * pblock%u0(1:nv,:,:,:) &
+ a33 * pblock%u1(1:nv,:,:,:) + ds * du(1:nv,:,:,:)
! assign pointer to the next block
!
@ -1003,7 +1003,7 @@ module evolution
!
! calculate fractional time step
!
ds = fh * dt
ds = b1 * dt
! update fluxes for the second step of the RK2 integration
!