Initializing AMReX (25.05-26-gfc469a32d34d)...
MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (25.05-26-gfc469a32d34d) 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.504934637
inner sponge: r_sp , r_tp : 152343750, 183593750
outer sponge: r_sp_outer, r_tp_outer: 183593750, 214843750
[1m[32mDoing initial projection[0m
Calling nodal solver
MLMG: Initial rhs = 41667132.21
MLMG: Initial residual (resid0) = 41667132.21
MLMG: Final Iter. 20 resid, resid/bnorm = 0.003808366135, 9.139976602e-11
MLMG: Timers: Solve = 0.271492656 Iter = 0.26714173 Bottom = 0.000319716
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.500071913
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
[1m[32mDoing initial divu iteration #1[0m
Calling nodal solver
MLMG: Initial rhs = 22468178.34
MLMG: Initial residual (resid0) = 22468178.34
MLMG: Final Iter. 2 resid, resid/bnorm = 4.603428559, 2.048865952e-07
MLMG: Timers: Solve = 0.031092258 Iter = 0.027034318 Bottom = 3.7553e-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
[1m[32mDoing initial divu iteration #2[0m
Calling nodal solver
MLMG: Initial rhs = 13466729.78
MLMG: Initial residual (resid0) = 13466729.78
MLMG: Final Iter. 3 resid, resid/bnorm = 0.009654234862, 7.168952684e-10
MLMG: Timers: Solve = 0.050109176 Iter = 0.043705198 Bottom = 5.909e-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
[1m[32mDoing initial divu iteration #3[0m
Calling nodal solver
MLMG: Initial rhs = 9200997.461
MLMG: Initial residual (resid0) = 9200997.461
MLMG: Final Iter. 7 resid, resid/bnorm = 0.0007183607668, 7.807422725e-11
MLMG: Timers: Solve = 0.098584223 Iter = 0.093628183 Bottom = 0.000113754
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.497018699
[1m[32mDoing initial pressure iteration #1[0m
[1m[32m
Timestep 0 starts with TIME = 0 DT = 0.0005784714117[0m
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 = 156168946.3
MLMG: Initial residual (resid0) = 156168946.3
MLMG: Final Iter. 8 resid, resid/bnorm = 0.009963668883, 6.380057699e-11
MLMG: Timers: Solve = 0.042642875 Iter = 0.038985823 Bottom = 0.000193492
<<< 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 = 156168960.8
MLMG: Initial residual (resid0) = 555.4440325
MLMG: Final Iter. 3 resid, resid/bnorm = 0.009710792452, 6.218132211e-11
MLMG: Timers: Solve = 0.018288495 Iter = 0.014750727 Bottom = 7.1865e-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 = 140238141.7
MLMG: Initial residual (resid0) = 140238141.7
MLMG: Final Iter. 27 resid, resid/bnorm = 0.00126494345, 9.019967287e-12
MLMG: Timers: Solve = 0.36505252 Iter = 0.36080936 Bottom = 0.000419878
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.442825927 seconds
MAC Proj :0.07190437 seconds
Nodal Proj :0.391475241 seconds
Reactions :0.772814145 seconds
Misc :0.11995013 seconds
Base State :0.002026786 seconds
Time to advance time step: 1.81505678
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.500113794
[1m[32mBeginning main evolution[0m
[1m[32m
Timestep 1 starts with TIME = 0 DT = 0.0005784714117[0m
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 = 156143863.7
MLMG: Initial residual (resid0) = 156143863.7
MLMG: Final Iter. 8 resid, resid/bnorm = 0.009885139763, 6.330789781e-11
MLMG: Timers: Solve = 0.042526359 Iter = 0.039132261 Bottom = 0.000195118
<<< 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 = 156143863.4
MLMG: Initial residual (resid0) = 1.447776845
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01058171352, 6.776900027e-11
MLMG: Timers: Solve = 0.013304401 Iter = 0.009884395 Bottom = 4.8041e-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 = 7408896.825
MLMG: Initial residual (resid0) = 7408896.825
MLMG: Final Iter. 16 resid, resid/bnorm = 5.175825208e-05, 6.985959355e-12
MLMG: Timers: Solve = 0.21762635 Iter = 0.213472353 Bottom = 0.000256598
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.435967302 seconds
MAC Proj :0.066842155 seconds
Nodal Proj :0.244260545 seconds
Reactions :0.774339046 seconds
Misc :0.119670427 seconds
Base State :0.001980711 seconds
Time to advance time step: 1.660707045
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.503725744
Call to estdt for level 0 gives dt_lev = 1.338762569
Minimum estdt over all levels = 1.338762569
Call to estdt at beginning of step 2 gives dt =1.338762569
dt_growth factor limits the new dt = 0.0006363185529
[1m[32m
Timestep 2 starts with TIME = 0.0005784714117 DT = 0.0006363185529[0m
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 = 148819014.4
MLMG: Initial residual (resid0) = 148819014.4
MLMG: Final Iter. 8 resid, resid/bnorm = 0.01129304618, 7.588443065e-11
MLMG: Timers: Solve = 0.042655823 Iter = 0.039111114 Bottom = 0.000193157
<<< 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 = 148819013.5
MLMG: Initial residual (resid0) = 1.432991188
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01314416307, 8.832314337e-11
MLMG: Timers: Solve = 0.013518002 Iter = 0.009998522 Bottom = 5.1661e-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 = 6145962.493
MLMG: Initial residual (resid0) = 6145962.493
MLMG: Final Iter. 20 resid, resid/bnorm = 4.935404286e-05, 8.030319566e-12
MLMG: Timers: Solve = 0.277505404 Iter = 0.273263825 Bottom = 0.000328522
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.441673753 seconds
MAC Proj :0.067604408 seconds
Nodal Proj :0.303995726 seconds
Reactions :0.76974667 seconds
Misc :0.119781855 seconds
Base State :0.00191052 seconds
Time to advance time step: 1.704950492
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.504779986
Call to estdt for level 0 gives dt_lev = 1.641614762
Minimum estdt over all levels = 1.641614762
Call to estdt at beginning of step 3 gives dt =1.641614762
dt_growth factor limits the new dt = 0.0006999504082
[1m[32m
Timestep 3 starts with TIME = 0.001214789965 DT = 0.0006999504082[0m
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 = 146032253.1
MLMG: Initial residual (resid0) = 146032253.1
MLMG: Final Iter. 8 resid, resid/bnorm = 0.01191277802, 8.157634883e-11
MLMG: Timers: Solve = 0.042347265 Iter = 0.038920222 Bottom = 0.000197996
<<< 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 = 146032252.2
MLMG: Initial residual (resid0) = 1.363398294
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01211341214, 8.295025213e-11
MLMG: Timers: Solve = 0.013418228 Iter = 0.00992447 Bottom = 4.7672e-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 = 5204170.627
MLMG: Initial residual (resid0) = 5204170.627
MLMG: Final Iter. 21 resid, resid/bnorm = 3.837421536e-05, 7.373742737e-12
MLMG: Timers: Solve = 0.28428461 Iter = 0.280113692 Bottom = 0.00033437
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.458083291 seconds
MAC Proj :0.066935225 seconds
Nodal Proj :0.311218512 seconds
Reactions :0.773830048 seconds
Misc :0.120417365 seconds
Base State :0.004194504 seconds
Time to advance time step: 1.73263975
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.508684775
Total Time: 13.47363363
Unused ParmParse Variables:
[TOP]::amr.check_file(nvals = 1) :: [wdconvect_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
[TOP]::amr.ref_ratio(nvals = 6) :: [2, 2, 2, 2, 2, 2]
AMReX (25.05-26-gfc469a32d34d) finalized