MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.09-11-g9aa23c202a13) initialized
Calling Setup()
Calling ReadParameters()
reading extern runtime parameters ...
Calling VariableSetup()
Calling BCSetup()
Calling BaseStateGeometry::Init()
Calling Init()
Calling InitData()
initdata model_File = kepler_new_6.25e8.hybrid.hse.320
model file = kepler_new_6.25e8.hybrid.hse.320
reading initial model
320 points found in the initial model
6 variables found in the initial model
WARNING: variable not found: ash
WARNING: magnesium-24 not provided in inputs file
model file mapping (spherical base state)
dr of MAESTRO base state = 781250.000000
dr of input file data = 1562500.000000
maximum radius (cell-centered) of input model = 499218750.000000
setting r_cutoff to 239
radius at r_cutoff 187109375
Maximum HSE Error = 0.002311
(after putting initial model into base state arrays, and
for density < base_cutoff_density)
Writing plotfile wdconvect-amr_pltInitData after InitData
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.857499532
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Doing initial projection
Calling nodal solver
MLMG: Initial rhs = 268163993.9
MLMG: Initial residual (resid0) = 268163993.9
MLMG: Final Iter. 4 resid, resid/bnorm = 0.02479596774, 9.24656863e-11
MLMG: Timers: Solve = 0.172685012 Iter = 0.166876732 Bottom = 0.000119678
Done calling nodal solver
Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.840397217
Call to firstdt for level 0 gives dt_lev = 0.005789736987
Multiplying dt_lev by init_shrink; dt_lev = 0.0005789736987
Call to firstdt for level 1 gives dt_lev = 0.002885750885
Multiplying dt_lev by init_shrink; dt_lev = 0.0002885750885
Minimum firstdt over all levels = 0.0002885750885
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 100896162.8
MLMG: Initial residual (resid0) = 100896162.8
MLMG: Final Iter. 2 resid, resid/bnorm = 1.093377186, 1.083665777e-08
MLMG: Timers: Solve = 0.071471481 Iter = 0.066469123 Bottom = 5.9495e-05
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 21.60061365
Call to estdt for level 1 gives dt_lev = 8.768321865
Minimum estdt over all levels = 8.768321865
Call to estdt at end of istep_divu_iter = 1 gives dt = 8.768321865
Multiplying dt by init_shrink; dt = 0.8768321865
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
Doing initial divu iteration #2
Calling nodal solver
MLMG: Initial rhs = 49401163.56
MLMG: Initial residual (resid0) = 49401163.56
MLMG: Final Iter. 2 resid, resid/bnorm = 0.02217691438, 4.489148187e-10
MLMG: Timers: Solve = 0.071356122 Iter = 0.066645632 Bottom = 6.0296e-05
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 21.60061365
Call to estdt for level 1 gives dt_lev = 8.96217536
Minimum estdt over all levels = 8.96217536
Call to estdt at end of istep_divu_iter = 2 gives dt = 8.96217536
Multiplying dt by init_shrink; dt = 0.896217536
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
Doing initial divu iteration #3
Calling nodal solver
MLMG: Initial rhs = 29021861.07
MLMG: Initial residual (resid0) = 29021861.07
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001463753683, 5.04362446e-11
MLMG: Timers: Solve = 0.104281093 Iter = 0.099617862 Bottom = 8.5117e-05
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 21.60061365
Call to estdt for level 1 gives dt_lev = 9.05539753
Minimum estdt over all levels = 9.05539753
Call to estdt at end of istep_divu_iter = 3 gives dt = 9.05539753
Multiplying dt by init_shrink; dt = 0.905539753
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.841634506
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 0.0002885750885
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 264654315.6
MLMG: Initial residual (resid0) = 264654315.6
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1712112343, 6.469240221e-10
MLMG: Timers: Solve = 0.058623468 Iter = 0.051647015 Bottom = 0.000221006
<<< 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 = 264654315.9
MLMG: Initial residual (resid0) = 8910.020001
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1666150695, 6.295573487e-10
MLMG: Timers: Solve = 0.025370562 Iter = 0.0211831 Bottom = 8.9711e-05
<<< 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 = 14853693.18
MLMG: Initial residual (resid0) = 14853693.18
MLMG: Final Iter. 31 resid, resid/bnorm = 0.001154070022, 7.76958301e-11
MLMG: Timers: Solve = 0.984405025 Iter = 0.979661715 Bottom = 0.000859224
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :0.975084891 seconds
MAC Proj :0.109038 seconds
Nodal Proj :1.033372516 seconds
Reactions :1.073713493 seconds
Misc :0.228545344 seconds
Base State :0.005005416 seconds
Time to advance time step: 3.450000945
Writing plotfile 0 after all initialization
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.826415865
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 0.0002885750885
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 264654252.4
MLMG: Initial residual (resid0) = 264654252.4
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1711969408, 6.468701684e-10
MLMG: Timers: Solve = 0.055889371 Iter = 0.05170096 Bottom = 0.000222446
<<< 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 = 264654252.2
MLMG: Initial residual (resid0) = 0.2043398281
MLMG: No iterations needed
MLMG: Timers: Solve = 0.004274238 Iter = 0 Bottom = 0
<<< 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 = 19692713.44
MLMG: Initial residual (resid0) = 19692713.44
MLMG: Final Iter. 7 resid, resid/bnorm = 0.001801351085, 9.147297504e-11
MLMG: Timers: Solve = 0.227133765 Iter = 0.222387029 Bottom = 0.000201695
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :0.94575961 seconds
MAC Proj :0.083116599 seconds
Nodal Proj :0.276760487 seconds
Reactions :1.083874582 seconds
Misc :0.228096932 seconds
Base State :0.005042799 seconds
Time to advance time step: 2.646847719
Writing plotfile 1
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.837717447
Call to estdt for level 0 gives dt_lev = 3.005533232
Call to estdt for level 1 gives dt_lev = 0.6541869341
Minimum estdt over all levels = 0.6541869341
Call to estdt at beginning of step 2 gives dt =0.6541869341
dt_growth factor limits the new dt = 0.0003174325974
Timestep 2 starts with TIME = 0.0002885750885 DT = 0.0003174325974
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 263576966.2
MLMG: Initial residual (resid0) = 263576966.2
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1796549195, 6.816032605e-10
MLMG: Timers: Solve = 0.056095052 Iter = 0.051873512 Bottom = 0.000215043
<<< 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 = 263576966
MLMG: Initial residual (resid0) = 0.2137418815
MLMG: No iterations needed
MLMG: Timers: Solve = 0.004228835 Iter = 0 Bottom = 0
<<< 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 = 14829531.18
MLMG: Initial residual (resid0) = 14829531.18
MLMG: Final Iter. 5 resid, resid/bnorm = 0.00126108711, 8.503890609e-11
MLMG: Timers: Solve = 0.168718085 Iter = 0.163987745 Bottom = 0.000144726
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection :0.957186165 seconds
MAC Proj :0.08329996 seconds
Nodal Proj :0.218191148 seconds
Reactions :1.08469 seconds
Misc :0.229130807 seconds
Base State :0.005064395 seconds
Time to advance time step: 2.572908228
Writing plotfile 2
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.842784906
Time to regrid: 0.136012634
Call to estdt for level 0 gives dt_lev = 4.022414647
Call to estdt for level 1 gives dt_lev = 0.8210389182
Minimum estdt over all levels = 0.8210389182
Call to estdt at beginning of step 3 gives dt =0.8210389182
dt_growth factor limits the new dt = 0.0003491758571
Timestep 3 starts with TIME = 0.0006060076859 DT = 0.0003491758571
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 263106326.3
MLMG: Initial residual (resid0) = 263106326.3
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1838002102, 6.985776922e-10
MLMG: Timers: Solve = 0.055261694 Iter = 0.051084934 Bottom = 0.000218961
<<< 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 = 263106326.1
MLMG: Initial residual (resid0) = 1974.351877
MLMG: Final Iter. 2 resid, resid/bnorm = 0.038897662, 1.478400865e-10
MLMG: Timers: Solve = 0.025565692 Iter = 0.021244411 Bottom = 8.9362e-05
<<< 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 = 12313408.17
MLMG: Initial residual (resid0) = 12313408.17
MLMG: Final Iter. 6 resid, resid/bnorm = 0.0005361591016, 4.354270517e-11
MLMG: Timers: Solve = 0.201882293 Iter = 0.197129076 Bottom = 0.000173646
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection :0.960811188 seconds
MAC Proj :0.103693253 seconds
Nodal Proj :0.251402837 seconds
Reactions :1.084526347 seconds
Misc :0.228369742 seconds
Base State :0.005247428001 seconds
Time to advance time step: 2.629229282
Writing plotfile 3
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.846446777
Total Time: 20.78532769
Unused ParmParse Variables:
[TOP]::maestro.v(nvals = 1) :: [1]
[TOP]::amr.check_file(nvals = 1) :: [wdconvect-amr_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
AMReX (22.09-11-g9aa23c202a13) finalized