MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.01-38-g3a2a74fc4b38) 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.82460756
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.02424257126, 9.040203687e-11
MLMG: Timers: Solve = 0.146471862 Iter = 0.140584552 Bottom = 0.000125834
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.813747841
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.088073483, 1.078409181e-08
MLMG: Timers: Solve = 0.075415793 Iter = 0.070286805 Bottom = 6.78e-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.768320188
Minimum estdt over all levels = 8.768320188
Call to estdt at end of istep_divu_iter = 1 gives dt = 8.768320188
Multiplying dt by init_shrink; dt = 0.8768320188
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.01943439664, 3.933995728e-10
MLMG: Timers: Solve = 0.075352811 Iter = 0.07014237 Bottom = 6.67e-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.962173609
Minimum estdt over all levels = 8.962173609
Call to estdt at end of istep_divu_iter = 2 gives dt = 8.962173609
Multiplying dt by init_shrink; dt = 0.8962173609
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.06
MLMG: Initial residual (resid0) = 29021861.06
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001475231654, 5.083173854e-11
MLMG: Timers: Solve = 0.112663515 Iter = 0.107500859 Bottom = 9.8917e-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.055395742
Minimum estdt over all levels = 9.055395742
Call to estdt at end of istep_divu_iter = 3 gives dt = 9.055395742
Multiplying dt by init_shrink; dt = 0.9055395742
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.812756657
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 = 264654280.2
MLMG: Initial residual (resid0) = 264654280.2
MLMG: Final Iter. 18 resid, resid/bnorm = 0.2588133365, 9.779299103e-10
MLMG: Timers: Solve = 0.190237436 Iter = 0.184037593 Bottom = 0.000774319
<<< 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 = 264654280.6
MLMG: Initial residual (resid0) = 3954.864047
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1976329442, 7.46758918e-10
MLMG: Timers: Solve = 0.025755406 Iter = 0.021334784 Bottom = 9.3584e-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 = 14853731.07
MLMG: Initial residual (resid0) = 14853731.07
MLMG: Final Iter. 12 resid, resid/bnorm = 0.001012212453, 6.814533321e-11
MLMG: Timers: Solve = 0.422683048 Iter = 0.417498363 Bottom = 0.000382127
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :0.963444324 seconds
MAC Proj :0.240882459 seconds
Nodal Proj :0.470953716 seconds
Reactions :1.045549281 seconds
Misc :0.219842262 seconds
Base State :0.004868777001 seconds
Time to advance time step: 2.966682743
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.81209059
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 = 264654217.1
MLMG: Initial residual (resid0) = 264654217.1
MLMG: Final Iter. 18 resid, resid/bnorm = 0.2589809522, 9.785634819e-10
MLMG: Timers: Solve = 0.187685118 Iter = 0.183209043 Bottom = 0.000774104
<<< 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 = 264654216.9
MLMG: Initial residual (resid0) = 6.7804943
MLMG: Final Iter. 2 resid, resid/bnorm = 0.2229616903, 8.424641513e-10
MLMG: Timers: Solve = 0.025851494 Iter = 0.02143887 Bottom = 9.1193e-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 = 19692713.46
MLMG: Initial residual (resid0) = 19692713.46
MLMG: Final Iter. 5 resid, resid/bnorm = 0.001142572695, 5.802007416e-11
MLMG: Timers: Solve = 0.181466881 Iter = 0.17618634 Bottom = 0.000161216
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :0.947254283 seconds
MAC Proj :0.237093861 seconds
Nodal Proj :0.230262569 seconds
Reactions :1.046496888 seconds
Misc :0.223184062 seconds
Base State :0.004790637999 seconds
Time to advance time step: 2.710727563
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.823869792
Call to estdt for level 0 gives dt_lev = 3.005530273
Call to estdt for level 1 gives dt_lev = 0.6541869199
Minimum estdt over all levels = 0.6541869199
Call to estdt at beginning of step 2 gives dt =0.6541869199
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 = 263576945.6
MLMG: Initial residual (resid0) = 263576945.6
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1808852218, 6.862710296e-10
MLMG: Timers: Solve = 0.057572164 Iter = 0.052428678 Bottom = 0.000232591
<<< 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 = 263576945.4
MLMG: Initial residual (resid0) = 0.2959210049
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01847116854, 7.007884745e-11
MLMG: Timers: Solve = 0.01647675 Iter = 0.011974873 Bottom = 4.4368e-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 = 14829531.19
MLMG: Initial residual (resid0) = 14829531.19
MLMG: Final Iter. 11 resid, resid/bnorm = 0.001472469419, 9.92930525e-11
MLMG: Timers: Solve = 0.374455025 Iter = 0.369318829 Bottom = 0.000331533
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection :0.964421879 seconds
MAC Proj :0.097563513 seconds
Nodal Proj :0.423165663 seconds
Reactions :1.04480824 seconds
Misc :0.221869675 seconds
Base State :0.004721317 seconds
Time to advance time step: 2.752256435
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.826780957
Time to regrid: 0.129168761
Call to estdt for level 0 gives dt_lev = 4.022358291
Call to estdt for level 1 gives dt_lev = 0.8210381227
Minimum estdt over all levels = 0.8210381227
Call to estdt at beginning of step 3 gives dt =0.8210381227
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 = 263106338.8
MLMG: Initial residual (resid0) = 263106338.8
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1831941996, 6.962743672e-10
MLMG: Timers: Solve = 0.055570479 Iter = 0.051172817 Bottom = 0.000217297
<<< 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 = 263106338.6
MLMG: Initial residual (resid0) = 0.2598261177
MLMG: No iterations needed
MLMG: Timers: Solve = 0.004368241 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 = 12313408.17
MLMG: Initial residual (resid0) = 12313408.17
MLMG: Final Iter. 6 resid, resid/bnorm = 0.0005843852191, 4.745925832e-11
MLMG: Timers: Solve = 0.211770706 Iter = 0.206628332 Bottom = 0.000193716
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection :0.96490228 seconds
MAC Proj :0.083440041 seconds
Nodal Proj :0.26054377 seconds
Reactions :1.044713726 seconds
Misc :0.222639889 seconds
Base State :0.004700678 seconds
Time to advance time step: 2.576661057
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.830765618
Total Time: 20.57043074
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.01-38-g3a2a74fc4b38) finalized