Initializing AMReX (25.04-16-gae7024c40673)...
MPI initialized with 2 MPI processes
MPI initialized with thread support level 0
OMP initialized with 4 OMP threads
AMReX (25.04-16-gae7024c40673) initialized
Calling Setup()
Calling ReadParameters()
reading extern runtime parameters ...
Calling VariableSetup()
Calling BCSetup()
Calling BaseStateGeometry::Init()
Calling Init()
Calling InitData()
initdata model_File =
cutoff densities:
low density cutoff (for mapping the model) = 1e-10
buoyancy cutoff density
(for zeroing rho - rho_0, centrifugal term) = 5e-10
anelastic cutoff = 1e-10
Maximum HSE Error = 2.273736754e-13
(after putting initial model into base state arrays, and
for density < base_cutoff_density)
Writing plotfile RT_pltInitData after InitData
Time to write plotfile: 0.083244442
[1m[32mDoing initial projection[0m
Calling nodal solver
MLMG: Initial rhs = 0
MLMG: Initial residual (resid0) = 0
MLMG: No iterations needed
MLMG: Timers: Solve = 0.003972814 Iter = 0 Bottom = 0
Done calling nodal solver
Writing plotfile RT_pltafter_InitProj after InitProj
Time to write plotfile: 0.06006856
Call to firstdt for level 0 gives dt_lev = 0.0006089835859
Multiplying dt_lev by init_shrink; dt_lev = 6.089835859e-05
Minimum firstdt over all levels = 6.089835859e-05
[1m[32mDoing initial divu iteration #1[0m
Calling nodal solver
MLMG: Initial rhs = 0
MLMG: Initial residual (resid0) = 0
MLMG: No iterations needed
MLMG: Timers: Solve = 0.003562443 Iter = 0 Bottom = 0
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 0.08660444773
Minimum estdt over all levels = 0.08660444773
Call to estdt at end of istep_divu_iter = 1 gives dt = 0.08660444773
Multiplying dt by init_shrink; dt = 0.008660444773
Ignoring this new dt since it's larger than the previous dt = 6.089835859e-05
Writing plotfile RT_pltafter_DivuIter after final DivuIter
Time to write plotfile: 0.057387009
[1m[32mDoing initial pressure iteration #1[0m
[1m[32m
Timestep 0 starts with TIME = 0 DT = 6.089835859e-05[0m
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 0.002421164472
MLMG: Initial residual (resid0) = 0.002421164472
MLMG: Final Iter. 9 resid, resid/bnorm = 1.751907767e-13, 7.235806522e-11
MLMG: Timers: Solve = 0.03481192 Iter = 0.032666457 Bottom = 0.002297206
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 0.002421205691
MLMG: Initial residual (resid0) = 2.151449931e-06
MLMG: Final Iter. 7 resid, resid/bnorm = 9.927909435e-14, 4.100399017e-11
MLMG: Timers: Solve = 0.026063978 Iter = 0.024177497 Bottom = 0.001680914
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 58.12657035
MLMG: Initial residual (resid0) = 58.12657035
MLMG: Final Iter. 8 resid, resid/bnorm = 2.15578666e-11, 3.708780077e-13
MLMG: Timers: Solve = 0.050456629 Iter = 0.048769401 Bottom = 0.001524661
Done calling nodal solver
Timestep 0 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05
Timing summary:
Advection :0.093304247 seconds
MAC Proj :0.068168831 seconds
Nodal Proj :0.056805025 seconds
Reactions :0.004790248001 seconds
Misc :0.029638142 seconds
Base State :0.0003981910008 seconds
Time to advance time step: 0.256331436
Writing plotfile 0 after all initialization
Time to write plotfile: 0.045953892
[1m[32mBeginning main evolution[0m
[1m[32m
Timestep 1 starts with TIME = 0 DT = 6.089835859e-05[0m
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 2.394238647e-05
MLMG: Initial residual (resid0) = 2.394238647e-05
MLMG: Final Iter. 9 resid, resid/bnorm = 3.181792236e-16, 1.328936963e-11
MLMG: Timers: Solve = 0.034333351 Iter = 0.032321607 Bottom = 0.002368908
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 2.394218844e-05
MLMG: Initial residual (resid0) = 1.928423574e-06
MLMG: Final Iter. 9 resid, resid/bnorm = 6.036511589e-16, 2.521286475e-11
MLMG: Timers: Solve = 0.03570189 Iter = 0.034012103 Bottom = 0.003061338
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 0.003539812449
MLMG: Initial residual (resid0) = 0.003539812449
MLMG: Final Iter. 8 resid, resid/bnorm = 1.329231863e-15, 3.755091216e-13
MLMG: Timers: Solve = 0.057953519 Iter = 0.056175894 Bottom = 0.002306844
Done calling nodal solver
Timestep 1 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05
Timing summary:
Advection :0.088524238 seconds
MAC Proj :0.077035335 seconds
Nodal Proj :0.064461166 seconds
Reactions :0.004399717 seconds
Misc :0.014560105 seconds
Base State :0.0004181200006 seconds
Time to advance time step: 0.252543779
Call to estdt for level 0 gives dt_lev = 0.3128878797
Minimum estdt over all levels = 0.3128878797
Call to estdt at beginning of step 2 gives dt =0.3128878797
dt_growth factor limits the new dt = 6.698819445e-05
[1m[32m
Timestep 2 starts with TIME = 6.089835859e-05 DT = 6.698819445e-05[0m
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 3.050043585e-05
MLMG: Initial residual (resid0) = 3.050043585e-05
MLMG: Final Iter. 9 resid, resid/bnorm = 2.284814624e-15, 7.491088439e-11
MLMG: Timers: Solve = 0.03623421 Iter = 0.034568654 Bottom = 0.003294952
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 3.059063412e-05
MLMG: Initial residual (resid0) = 5.66083759e-06
MLMG: Final Iter. 10 resid, resid/bnorm = 4.892627475e-16, 1.599387399e-11
MLMG: Timers: Solve = 0.039735216 Iter = 0.038056573 Bottom = 0.003429293
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 0.003894013591
MLMG: Initial residual (resid0) = 0.003894013591
MLMG: Final Iter. 8 resid, resid/bnorm = 1.490994828e-15, 3.828940996e-13
MLMG: Timers: Solve = 0.070218443 Iter = 0.068457342 Bottom = 0.002307096
Done calling nodal solver
Timestep 2 ends with TIME = 0.000127886553 DT = 6.698819445e-05
Timing summary:
Advection :0.0854875 seconds
MAC Proj :0.08284048 seconds
Nodal Proj :0.078208702 seconds
Reactions :0.004286689 seconds
Misc :0.014425914 seconds
Base State :0.0004027890009 seconds
Time to advance time step: 0.26580992
Call to estdt for level 0 gives dt_lev = 0.3130286362
Minimum estdt over all levels = 0.3130286362
Call to estdt at beginning of step 3 gives dt =0.3130286362
dt_growth factor limits the new dt = 7.36870139e-05
[1m[32m
Timestep 3 starts with TIME = 0.000127886553 DT = 7.36870139e-05[0m
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 3.168800511e-05
MLMG: Initial residual (resid0) = 3.168800511e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 5.650700787e-16, 1.783230205e-11
MLMG: Timers: Solve = 0.039110114 Iter = 0.037349521 Bottom = 0.003601979
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 3.178722337e-05
MLMG: Initial residual (resid0) = 1.012935749e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 6.833616179e-16, 2.149799653e-11
MLMG: Timers: Solve = 0.039073137 Iter = 0.037304599 Bottom = 0.003415191
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 0.004283414244
MLMG: Initial residual (resid0) = 0.004283414244
MLMG: Final Iter. 8 resid, resid/bnorm = 1.614810716e-15, 3.769914894e-13
MLMG: Timers: Solve = 0.051428403 Iter = 0.049708426 Bottom = 0.002319407
Done calling nodal solver
Timestep 3 ends with TIME = 0.0002015735669 DT = 7.36870139e-05
Timing summary:
Advection :0.085971559 seconds
MAC Proj :0.085139465 seconds
Nodal Proj :0.057610439 seconds
Reactions :0.004306422 seconds
Misc :0.02581123 seconds
Base State :0.0005359560005 seconds
Time to advance time step: 0.259136225
Call to estdt for level 0 gives dt_lev = 0.3130285358
Minimum estdt over all levels = 0.3130285358
Call to estdt at beginning of step 4 gives dt =0.3130285358
dt_growth factor limits the new dt = 8.105571529e-05
[1m[32m
Timestep 4 starts with TIME = 0.0002015735669 DT = 8.105571529e-05[0m
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 3.570538736e-05
MLMG: Initial residual (resid0) = 3.570538736e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 7.571110825e-16, 2.120439347e-11
MLMG: Timers: Solve = 0.039232196 Iter = 0.037563206 Bottom = 0.003629119
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 3.581452755e-05
MLMG: Initial residual (resid0) = 1.504387741e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 8.870129024e-16, 2.476684639e-11
MLMG: Timers: Solve = 0.038337075 Iter = 0.036599007 Bottom = 0.003414467
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 0.004711755906
MLMG: Initial residual (resid0) = 0.004711755906
MLMG: Final Iter. 8 resid, resid/bnorm = 1.768550584e-15, 3.753485153e-13
MLMG: Timers: Solve = 0.050953387 Iter = 0.049246274 Bottom = 0.002303432
Done calling nodal solver
Timestep 4 ends with TIME = 0.0002826292822 DT = 8.105571529e-05
Timing summary:
Advection :0.084341555 seconds
MAC Proj :0.084407848 seconds
Nodal Proj :0.057134496 seconds
Reactions :0.003991348 seconds
Misc :0.01426713 seconds
Base State :0.0003831880003 seconds
Time to advance time step: 0.244461774
Call to estdt for level 0 gives dt_lev = 0.313028505
Minimum estdt over all levels = 0.313028505
Call to estdt at beginning of step 5 gives dt =0.313028505
dt_growth factor limits the new dt = 8.916128682e-05
[1m[32m
Timestep 5 starts with TIME = 0.0002826292822 DT = 8.916128682e-05[0m
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 4.013180326e-05
MLMG: Initial residual (resid0) = 4.013180326e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 9.681425651e-16, 2.412407334e-11
MLMG: Timers: Solve = 0.038410527 Iter = 0.036688985 Bottom = 0.003554781
<<< STEP 4 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs = 4.025185748e-05
MLMG: Initial residual (resid0) = 2.044985444e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 1.111036176e-15, 2.760210946e-11
MLMG: Timers: Solve = 0.038429229 Iter = 0.036692607 Bottom = 0.00342338
<<< STEP 8 : advect base >>>
: density_advance >>>
: tracer_advance >>>
: enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs = 0.005182931901
MLMG: Initial residual (resid0) = 0.005182931901
MLMG: Final Iter. 8 resid, resid/bnorm = 1.942022931e-15, 3.746958224e-13
MLMG: Timers: Solve = 0.051019613 Iter = 0.049277187 Bottom = 0.002297934
Done calling nodal solver
Timestep 5 ends with TIME = 0.000371790569 DT = 8.916128682e-05
Timing summary:
Advection :0.083638784 seconds
MAC Proj :0.084001347 seconds
Nodal Proj :0.057244892 seconds
Reactions :0.003971067 seconds
Misc :0.014157791 seconds
Base State :0.0004150320005 seconds
Time to advance time step: 0.24330245
Writing plotfile 5
Time to write plotfile: 0.045621271
Total Time: 1.898361098
Unused ParmParse Variables:
[TOP]::amr.check_file(nvals = 1) :: [RT_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
[TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2]
AMReX (25.04-16-gae7024c40673) finalized