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

Merge branch 'master' into reconnection

Browse Source
This commit is contained in:
Grzegorz Kowal 2024-07-01 20:12:15 -03:00
commit 1c63d15b32
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
sources/evolution.F90
View File

@ -76,6 +76,7 @@ module evolution
real(kind=8) , save :: dte = 0.0d+00
real(kind=8) , save :: dtp = huge(dt) ! dt for precise snapshots
real(kind=8) , save :: dtsafe = 1.0d-02
real(kind=8) , save :: facjmp = 1.05d+00
! the absolute and relative tolerances, limiting factors, the maximum error,
! the maximum number of passes for the adaptive step,
@ -197,6 +198,7 @@ module evolution
call get_parameter("minimum_factor" , facmin)
call get_parameter("maximum_factor" , facmax)
call get_parameter("dt_safe_factor" , dtsafe)
call get_parameter("dt_jump_factor" , facjmp)
! select the integration method, check the correctness of the integration
! parameters and adjust the CFL coefficient if necessary
@ -1046,6 +1048,7 @@ module evolution
call print_parameter(verbose, "time advance", name_int)
call print_parameter(verbose, "no. of registers", registers)
call print_parameter(verbose, "CFL coefficient", cfl)
call print_parameter(verbose, "timestep jump factor", facjmp)
if (magnetized) then
call print_parameter(verbose, "GLM alpha coefficient", glm_alpha)
end if
@ -1459,7 +1462,7 @@ module evolution
dth = dx_min / max(cmax &
+ 2.0d+00 * max(viscosity, resistivity) / dx_min, eps)
dt = cfl * dth
dt = min(facjmp * dt, cfl * dth)
dt = min(dt, dtp)
if (error_control) dt = min(dt, dte)