MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.01-44-gfd6d6494bb75) 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.881750528
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.02431267316, 9.066345114e-11
MLMG: Timers: Solve = 0.144857305 Iter = 0.137226678 Bottom = 0.000143457
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.877546541
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.088460624, 1.078792884e-08
MLMG: Timers: Solve = 0.07499677 Iter = 0.068906388 Bottom = 6.8961e-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.76832059
Minimum estdt over all levels = 8.76832059
Call to estdt at end of istep_divu_iter = 1 gives dt = 8.76832059
Multiplying dt by init_shrink; dt = 0.876832059
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.01948056556, 3.943341443e-10
MLMG: Timers: Solve = 0.075527977 Iter = 0.069270994 Bottom = 6.884e-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.962174028
Minimum estdt over all levels = 8.962174028
Call to estdt at end of istep_divu_iter = 2 gives dt = 8.962174028
Multiplying dt by init_shrink; dt = 0.8962174028
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.001457821783, 5.023185039e-11
MLMG: Timers: Solve = 0.110922603 Iter = 0.104812594 Bottom = 0.000100277
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.055396171
Minimum estdt over all levels = 9.055396171
Call to estdt at end of istep_divu_iter = 3 gives dt = 9.055396171
Multiplying dt by init_shrink; dt = 0.9055396171
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.83745515
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 = 264654289.7
MLMG: Initial residual (resid0) = 264654289.7
MLMG: Final Iter. 19 resid, resid/bnorm = 0.2254616693, 8.51910126e-10
MLMG: Timers: Solve = 0.211104938 Iter = 0.20264882 Bottom = 0.000844546
<<< 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 = 264654290.1
MLMG: Initial residual (resid0) = 4556.639188
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1724912524, 6.517606511e-10
MLMG: Timers: Solve = 0.029240205 Iter = 0.022844139 Bottom = 9.8973e-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 = 14853677.88
MLMG: Initial residual (resid0) = 14853677.88
MLMG: Final Iter. 12 resid, resid/bnorm = 0.001004018429, 6.759392772e-11
MLMG: Timers: Solve = 0.435119522 Iter = 0.428780417 Bottom = 0.000394454
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.241643394 seconds
MAC Proj :0.274538712 seconds
Nodal Proj :0.494707239 seconds
Reactions :1.069547922 seconds
Misc :0.252007107 seconds
Base State :0.007140020001 seconds
Time to advance time step: 3.363214391
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.847290728
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 = 264654226.6
MLMG: Initial residual (resid0) = 264654226.6
MLMG: Final Iter. 19 resid, resid/bnorm = 0.2255878225, 8.52387001e-10
MLMG: Timers: Solve = 0.209229452 Iter = 0.20281365 Bottom = 0.000894099
<<< 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 = 264654226.4
MLMG: Initial residual (resid0) = 6.414984538
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1746363845, 6.598662218e-10
MLMG: Timers: Solve = 0.029414525 Iter = 0.022986754 Bottom = 0.000124762
<<< 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. 5 resid, resid/bnorm = 0.001143789077, 5.808184232e-11
MLMG: Timers: Solve = 0.180894702 Iter = 0.174541293 Bottom = 0.000166876
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.23860362 seconds
MAC Proj :0.27214083 seconds
Nodal Proj :0.239978772 seconds
Reactions :1.056636834 seconds
Misc :0.251198895 seconds
Base State :0.006488501 seconds
Time to advance time step: 3.101491929
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.849930494
Call to estdt for level 0 gives dt_lev = 3.005531336
Call to estdt for level 1 gives dt_lev = 0.6541868857
Minimum estdt over all levels = 0.6541868857
Call to estdt at beginning of step 2 gives dt =0.6541868857
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 = 263576951.5
MLMG: Initial residual (resid0) = 263576951.5
MLMG: Final Iter. 5 resid, resid/bnorm = 0.177768942, 6.74447978e-10
MLMG: Timers: Solve = 0.060864662 Iter = 0.054266736 Bottom = 0.000219491
<<< 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 = 263576951.3
MLMG: Initial residual (resid0) = 0.288285995
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01677477685, 6.364280628e-11
MLMG: Timers: Solve = 0.018634383 Iter = 0.012164062 Bottom = 4.3491e-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.18
MLMG: Initial residual (resid0) = 14829531.18
MLMG: Final Iter. 12 resid, resid/bnorm = 0.001177727245, 7.941769909e-11
MLMG: Timers: Solve = 0.410141676 Iter = 0.403836967 Bottom = 0.000365015
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection :1.251115447 seconds
MAC Proj :0.113110421 seconds
Nodal Proj :0.469115451 seconds
Reactions :1.060551368 seconds
Misc :0.252411083 seconds
Base State :0.006607003 seconds
Time to advance time step: 3.146827478
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.852299779
Time to regrid: 0.134313703
Call to estdt for level 0 gives dt_lev = 4.022374534
Call to estdt for level 1 gives dt_lev = 0.8210384495
Minimum estdt over all levels = 0.8210384495
Call to estdt at beginning of step 3 gives dt =0.8210384495
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 = 263106335.3
MLMG: Initial residual (resid0) = 263106335.3
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1847731462, 7.022755495e-10
MLMG: Timers: Solve = 0.06061011 Iter = 0.054191882 Bottom = 0.000221781
<<< 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 = 263106335.1
MLMG: Initial residual (resid0) = 0.2530228249
MLMG: No iterations needed
MLMG: Timers: Solve = 0.006374312 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.0005713972444, 4.640447523e-11
MLMG: Timers: Solve = 0.214101027 Iter = 0.207796913 Bottom = 0.000188485
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection :1.256412011 seconds
MAC Proj :0.100795735 seconds
Nodal Proj :0.273128713 seconds
Reactions :1.057653575 seconds
Misc :0.252885777 seconds
Base State :0.006447390001 seconds
Time to advance time step: 2.941411619
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.856890989
Total Time: 22.57313397
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-44-gfd6d6494bb75) finalized