MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.09) 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.533462815
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. 19 resid, resid/bnorm = 0.007501523942, 9.0858117e-11
MLMG: Timers: Solve = 0.310356671 Iter = 0.304457099 Bottom = 0.000291738
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.523031704
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.807827444, 1.124318395e-07
MLMG: Timers: Solve = 0.036366508 Iter = 0.03260387 Bottom = 3.9675e-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.01720862149, 9.269994554e-10
MLMG: Timers: Solve = 0.05197524 Iter = 0.048259637 Bottom = 5.4554e-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. 5 resid, resid/bnorm = 0.00106503116, 9.093305965e-11
MLMG: Timers: Solve = 0.083421473 Iter = 0.07972416 Bottom = 8.0273e-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.515752403
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 = 97092964.44
MLMG: Initial residual (resid0) = 97092964.44
MLMG: Final Iter. 8 resid, resid/bnorm = 0.00966941379, 9.95892323e-11
MLMG: Timers: Solve = 0.039830272 Iter = 0.036554143 Bottom = 0.000205953
<<< 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 = 97092956.34
MLMG: Initial residual (resid0) = 5219.663759
MLMG: Final Iter. 4 resid, resid/bnorm = 0.005041599274, 5.192548938e-11
MLMG: Timers: Solve = 0.021400418 Iter = 0.01819362 Bottom = 0.000104226
<<< 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 = 140046470.4
MLMG: Initial residual (resid0) = 140046470.4
MLMG: Final Iter. 29 resid, resid/bnorm = 0.001360194641, 9.712452138e-12
MLMG: Timers: Solve = 0.456948687 Iter = 0.453264617 Bottom = 0.000397568
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.439217818 seconds
MAC Proj :0.071543274 seconds
Nodal Proj :0.479444215 seconds
Reactions :0.709665145 seconds
Misc :0.123648081 seconds
Base State :0.001935897 seconds
Time to advance time step: 1.851162589
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.529629551
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 = 97078095.67
MLMG: Initial residual (resid0) = 97078095.67
MLMG: Final Iter. 8 resid, resid/bnorm = 0.009577121586, 9.865378508e-11
MLMG: Timers: Solve = 0.039379612 Iter = 0.036179733 Bottom = 0.000216226
<<< 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 = 97078095.5
MLMG: Initial residual (resid0) = 2.550180789
MLMG: Final Iter. 3 resid, resid/bnorm = 0.008997315541, 9.268121191e-11
MLMG: Timers: Solve = 0.0170243 Iter = 0.013814736 Bottom = 8.5976e-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 = 8056355.308
MLMG: Initial residual (resid0) = 8056355.308
MLMG: Final Iter. 18 resid, resid/bnorm = 6.610807031e-05, 8.205704414e-12
MLMG: Timers: Solve = 0.293566862 Iter = 0.289894185 Bottom = 0.000265284
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.434122526 seconds
MAC Proj :0.066704606 seconds
Nodal Proj :0.316304403 seconds
Reactions :0.713633864 seconds
Misc :0.122115354 seconds
Base State :0.00183433 seconds
Time to advance time step: 1.697360052
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.520152201
Call to estdt for level 0 gives dt_lev = 1.410006144
Minimum estdt over all levels = 1.410006144
Call to estdt at beginning of step 2 gives dt =1.410006144
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 = 92795912.03
MLMG: Initial residual (resid0) = 92795912.03
MLMG: Final Iter. 9 resid, resid/bnorm = 0.004761144519, 5.13076968e-11
MLMG: Timers: Solve = 0.043878214 Iter = 0.040674467 Bottom = 0.00021255
<<< 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 = 92795911.86
MLMG: Initial residual (resid0) = 2.563477293
MLMG: Final Iter. 3 resid, resid/bnorm = 0.009068086743, 9.77207569e-11
MLMG: Timers: Solve = 0.017130305 Iter = 0.013925467 Bottom = 8.4195e-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.498
MLMG: Initial residual (resid0) = 5904880.498
MLMG: Final Iter. 20 resid, resid/bnorm = 5.254801363e-05, 8.899081641e-12
MLMG: Timers: Solve = 0.315978212 Iter = 0.312278858 Bottom = 0.000293144
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.435850284 seconds
MAC Proj :0.071240346 seconds
Nodal Proj :0.338520306 seconds
Reactions :0.711603128 seconds
Misc :0.121618841 seconds
Base State :0.001856012 seconds
Time to advance time step: 1.682157624
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.505977535
Call to estdt for level 0 gives dt_lev = 1.702890075
Minimum estdt over all levels = 1.702890075
Call to estdt at beginning of step 3 gives dt =1.702890075
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 = 91006708.33
MLMG: Initial residual (resid0) = 91006708.33
MLMG: Final Iter. 9 resid, resid/bnorm = 0.004961580038, 5.451883855e-11
MLMG: Timers: Solve = 0.044961211 Iter = 0.041773063 Bottom = 0.000219677
<<< 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 = 91006708.02
MLMG: Initial residual (resid0) = 2.401901934
MLMG: Final Iter. 3 resid, resid/bnorm = 0.00848967582, 9.328626432e-11
MLMG: Timers: Solve = 0.017234238 Iter = 0.01403227 Bottom = 8.4904e-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 = 4615830.289
MLMG: Initial residual (resid0) = 4615830.289
MLMG: Final Iter. 21 resid, resid/bnorm = 4.000449553e-05, 8.666803809e-12
MLMG: Timers: Solve = 0.332993992 Iter = 0.329245315 Bottom = 0.000298484
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.438145042 seconds
MAC Proj :0.073286979 seconds
Nodal Proj :0.356240966 seconds
Reactions :0.714261381 seconds
Misc :0.123411489 seconds
Base State :0.001773831 seconds
Time to advance time step: 1.708719046
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.509734591
Total Time: 13.83563149
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.09) finalized