MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.06-33-g6f72de283c38) 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 = 1562500.000000
dr of input file data = 1562500.000000
maximum radius (cell-centered) of input model = 499218750.000000
setting r_cutoff to 119
radius at r_cutoff 186718750
Maximum HSE Error = 0.001223
(after putting initial model into base state arrays, and
for density < base_cutoff_density)
Writing plotfile wdconvect_pltInitData after InitData
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.514445713
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Doing initial projection
Calling nodal solver
MLMG: Initial rhs = 82563057.54
MLMG: Initial residual (resid0) = 82563057.54
MLMG: Final Iter. 7 resid, resid/bnorm = 0.0026114434, 3.162968376e-11
MLMG: Timers: Solve = 0.114655657 Iter = 0.109550514 Bottom = 0.000128297
Done calling nodal solver
Writing plotfile wdconvect_pltafter_InitProj after InitProj
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.50528478
Call to firstdt for level 0 gives dt_lev = 0.005784714117
Multiplying dt_lev by init_shrink; dt_lev = 0.0005784714117
Minimum firstdt over all levels = 0.0005784714117
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 33867874.62
MLMG: Initial residual (resid0) = 33867874.62
MLMG: Final Iter. 2 resid, resid/bnorm = 3.810671231, 1.125158066e-07
MLMG: Timers: Solve = 0.035013867 Iter = 0.031407819 Bottom = 4.1493e-05
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 7.317922768
Minimum estdt over all levels = 7.317922768
Call to estdt at end of istep_divu_iter = 1 gives dt = 7.317922768
Multiplying dt by init_shrink; dt = 0.7317922768
Ignoring this new dt since it's larger than the previous dt = 0.0005784714117
Doing initial divu iteration #2
Calling nodal solver
MLMG: Initial rhs = 18563788.13
MLMG: Initial residual (resid0) = 18563788.13
MLMG: Final Iter. 3 resid, resid/bnorm = 0.01210942189, 6.523141616e-10
MLMG: Timers: Solve = 0.051700639 Iter = 0.048156517 Bottom = 5.4924e-05
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 7.317922768
Minimum estdt over all levels = 7.317922768
Call to estdt at end of istep_divu_iter = 2 gives dt = 7.317922768
Multiplying dt by init_shrink; dt = 0.7317922768
Ignoring this new dt since it's larger than the previous dt = 0.0005784714117
Doing initial divu iteration #3
Calling nodal solver
MLMG: Initial rhs = 11712254.75
MLMG: Initial residual (resid0) = 11712254.75
MLMG: Final Iter. 4 resid, resid/bnorm = 0.0006582401693, 5.620097778e-11
MLMG: Timers: Solve = 0.066672796 Iter = 0.063064528 Bottom = 6.6825e-05
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 7.317922768
Minimum estdt over all levels = 7.317922768
Call to estdt at end of istep_divu_iter = 3 gives dt = 7.317922768
Multiplying dt by init_shrink; dt = 0.7317922768
Ignoring this new dt since it's larger than the previous dt = 0.0005784714117
Writing plotfile wdconvect_pltafter_DivuIter after final DivuIter
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.499824216
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 0.0005784714117
Cell Count:
Level 0, 262144 cells
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 97092858.68
MLMG: Initial residual (resid0) = 97092858.68
MLMG: Final Iter. 35 resid, resid/bnorm = 0.008220948279, 8.467098807e-11
MLMG: Timers: Solve = 0.16437029 Iter = 0.16131036 Bottom = 0.000855312
<<< 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 = 97092850.36
MLMG: Initial residual (resid0) = 1577.675652
MLMG: Final Iter. 10 resid, resid/bnorm = 0.009494312108, 9.778590362e-11
MLMG: Timers: Solve = 0.048230766 Iter = 0.04648764 Bottom = 0.000257531
<<< 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 = 140045999.2
MLMG: Initial residual (resid0) = 140045999.2
MLMG: Final Iter. 9 resid, resid/bnorm = 0.001376963803, 9.832225199e-12
MLMG: Timers: Solve = 0.143848376 Iter = 0.140357358 Bottom = 0.000134263
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.432030988 seconds
MAC Proj :0.221368908 seconds
Nodal Proj :0.166268658 seconds
Reactions :0.737063416 seconds
Misc :0.126581008 seconds
Base State :0.001931272001 seconds
Time to advance time step: 1.708007595
Writing plotfile 0 after all initialization
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.500300438
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 0.0005784714117
Cell Count:
Level 0, 262144 cells
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 97077990.39
MLMG: Initial residual (resid0) = 97077990.39
MLMG: Final Iter. 35 resid, resid/bnorm = 0.008149342611, 8.394634643e-11
MLMG: Timers: Solve = 0.161100234 Iter = 0.159339427 Bottom = 0.000848041
<<< 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 = 97077990.2
MLMG: Initial residual (resid0) = 31.47476672
MLMG: Final Iter. 11 resid, resid/bnorm = 0.00714299269, 7.357993996e-11
MLMG: Timers: Solve = 0.052400913 Iter = 0.050623393 Bottom = 0.000282953
<<< 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 = 8056355.603
MLMG: Initial residual (resid0) = 8056355.603
MLMG: Final Iter. 5 resid, resid/bnorm = 2.859951928e-05, 3.549932586e-12
MLMG: Timers: Solve = 0.081767556 Iter = 0.078219133 Bottom = 8.3904e-05
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.400790805 seconds
MAC Proj :0.221287677 seconds
Nodal Proj :0.104354836 seconds
Reactions :0.719511052 seconds
Misc :0.120551765 seconds
Base State :0.002848471001 seconds
Time to advance time step: 1.593069899
Writing plotfile 1
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.508369428
Call to estdt for level 0 gives dt_lev = 1.410019887
Minimum estdt over all levels = 1.410019887
Call to estdt at beginning of step 2 gives dt =1.410019887
dt_growth factor limits the new dt = 0.0006363185529
Timestep 2 starts with TIME = 0.0005784714117 DT = 0.0006363185529
Cell Count:
Level 0, 262144 cells
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 92795851.03
MLMG: Initial residual (resid0) = 92795851.03
MLMG: Final Iter. 8 resid, resid/bnorm = 0.006255773827, 6.741436991e-11
MLMG: Timers: Solve = 0.039117216 Iter = 0.037358299 Bottom = 0.000204233
<<< 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 = 92795851.08
MLMG: Initial residual (resid0) = 3.609344613
MLMG: Final Iter. 3 resid, resid/bnorm = 0.008871402591, 9.560128485e-11
MLMG: Timers: Solve = 0.01589941 Iter = 0.014152202 Bottom = 8.3368e-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 = 5904880.805
MLMG: Initial residual (resid0) = 5904880.805
MLMG: Final Iter. 28 resid, resid/bnorm = 4.771212116e-05, 8.08011588e-12
MLMG: Timers: Solve = 0.450554058 Iter = 0.446897873 Bottom = 0.000384886
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.409543053 seconds
MAC Proj :0.062724684 seconds
Nodal Proj :0.473867764 seconds
Reactions :0.734998695 seconds
Misc :0.123232813 seconds
Base State :0.001949354 seconds
Time to advance time step: 1.804684024
Writing plotfile 2
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.524077108
Call to estdt for level 0 gives dt_lev = 1.702946918
Minimum estdt over all levels = 1.702946918
Call to estdt at beginning of step 3 gives dt =1.702946918
dt_growth factor limits the new dt = 0.0006999504082
Timestep 3 starts with TIME = 0.001214789965 DT = 0.0006999504082
Cell Count:
Level 0, 262144 cells
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 91006746.01
MLMG: Initial residual (resid0) = 91006746.01
MLMG: Final Iter. 9 resid, resid/bnorm = 0.005927838385, 6.513625248e-11
MLMG: Timers: Solve = 0.04370709 Iter = 0.041940927 Bottom = 0.000217336
<<< 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 = 91006745.71
MLMG: Initial residual (resid0) = 52.22964847
MLMG: Final Iter. 4 resid, resid/bnorm = 0.005878932774, 6.459886823e-11
MLMG: Timers: Solve = 0.020518534 Iter = 0.018738159 Bottom = 0.000103221
<<< 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 = 4615830.301
MLMG: Initial residual (resid0) = 4615830.301
MLMG: Final Iter. 19 resid, resid/bnorm = 3.268825822e-05, 7.081772094e-12
MLMG: Timers: Solve = 0.306891693 Iter = 0.30318507 Bottom = 0.000260992
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.412794624 seconds
MAC Proj :0.072026232 seconds
Nodal Proj :0.330179964 seconds
Reactions :0.748011198 seconds
Misc :0.124255229 seconds
Base State :0.001744467001 seconds
Time to advance time step: 1.687589114
Writing plotfile 3
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
Time to write plotfile: 0.525225062
Total Time: 13.26327753
Unused ParmParse Variables:
[TOP]::maestro.v(nvals = 1) :: [1]
[TOP]::amr.ref_ratio(nvals = 6) :: [2, 2, 2, 2, 2, 2]
[TOP]::amr.check_file(nvals = 1) :: [wdconvect_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
AMReX (22.06-33-g6f72de283c38) finalized