MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (23.01-27-g756468d42106) 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.590619094
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 = 41667132.21
MLMG: Initial residual (resid0) = 41667132.21
MLMG: Final Iter. 13 resid, resid/bnorm = 0.003028079867, 7.267310483e-11
MLMG: Timers: Solve = 0.221948732 Iter = 0.216761474 Bottom = 0.000266305
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.57311638
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 = 22468178.34
MLMG: Initial residual (resid0) = 22468178.34
MLMG: Final Iter. 2 resid, resid/bnorm = 4.59930351, 2.04703e-07
MLMG: Timers: Solve = 0.038832296 Iter = 0.034006627 Bottom = 4.8997e-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 = 13466729.78
MLMG: Initial residual (resid0) = 13466729.78
MLMG: Final Iter. 3 resid, resid/bnorm = 0.01356752368, 1.00748466e-09
MLMG: Timers: Solve = 0.056458196 Iter = 0.051798954 Bottom = 6.1731e-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 = 9200997.459
MLMG: Initial residual (resid0) = 9200997.459
MLMG: Final Iter. 5 resid, resid/bnorm = 0.0004817973822, 5.236360344e-11
MLMG: Timers: Solve = 0.088380596 Iter = 0.083528823 Bottom = 0.00012004
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.56762094
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 = 156168946.3
MLMG: Initial residual (resid0) = 156168946.3
MLMG: Final Iter. 27 resid, resid/bnorm = 0.01547729969, 9.910612869e-11
MLMG: Timers: Solve = 0.139905451 Iter = 0.135191917 Bottom = 0.00068489
<<< 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 = 156168961.1
MLMG: Initial residual (resid0) = 995.1267602
MLMG: Final Iter. 13 resid, resid/bnorm = 0.0151938647, 9.729119406e-11
MLMG: Timers: Solve = 0.067019369 Iter = 0.064058164 Bottom = 0.000395856
<<< 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 = 140238774.2
MLMG: Initial residual (resid0) = 140238774.2
MLMG: Final Iter. 7 resid, resid/bnorm = 0.0002321262526, 1.655221631e-12
MLMG: Timers: Solve = 0.121408904 Iter = 0.116918041 Bottom = 0.000148079
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.55088016 seconds
MAC Proj :0.220109132 seconds
Nodal Proj :0.147335969 seconds
Reactions :0.843593114 seconds
Misc :0.141946162 seconds
Base State :0.002281216 seconds
Time to advance time step: 1.935748445
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.564430587
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 = 156143863.6
MLMG: Initial residual (resid0) = 156143863.6
MLMG: Final Iter. 27 resid, resid/bnorm = 0.01503571868, 9.629400949e-11
MLMG: Timers: Solve = 0.134605464 Iter = 0.13158309 Bottom = 0.000726269
<<< 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) = 12.90985522
MLMG: Final Iter. 14 resid, resid/bnorm = 0.01083907485, 6.94172324e-11
MLMG: Timers: Solve = 0.071640973 Iter = 0.068627417 Bottom = 0.00032973
<<< 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.529
MLMG: Initial residual (resid0) = 7408896.529
MLMG: Final Iter. 5 resid, resid/bnorm = 6.248708814e-05, 8.434061388e-12
MLMG: Timers: Solve = 0.08814017 Iter = 0.083334955 Bottom = 8.4918e-05
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.540839383 seconds
MAC Proj :0.217989174 seconds
Nodal Proj :0.114204371 seconds
Reactions :0.840742246 seconds
Misc :0.138395293 seconds
Base State :0.002683047 seconds
Time to advance time step: 1.883195341
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.574700393
Call to estdt for level 0 gives dt_lev = 1.338762745
Minimum estdt over all levels = 1.338762745
Call to estdt at beginning of step 2 gives dt =1.338762745
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 = 148819013.9
MLMG: Initial residual (resid0) = 148819013.9
MLMG: Final Iter. 8 resid, resid/bnorm = 0.009842343628, 6.613633146e-11
MLMG: Timers: Solve = 0.042185113 Iter = 0.039164102 Bottom = 0.000190646
<<< 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 = 148819012.7
MLMG: Initial residual (resid0) = 378.8848725
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01253152778, 8.420649728e-11
MLMG: Timers: Solve = 0.012959431 Iter = 0.010001306 Bottom = 5.8953e-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.801
MLMG: Initial residual (resid0) = 6145962.801
MLMG: Final Iter. 25 resid, resid/bnorm = 5.296524614e-05, 8.617892405e-12
MLMG: Timers: Solve = 0.421555701 Iter = 0.416837786 Bottom = 0.000349566
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.525775105 seconds
MAC Proj :0.066942222 seconds
Nodal Proj :0.447690479 seconds
Reactions :0.84111674 seconds
Misc :0.138348566 seconds
Base State :0.002258278 seconds
Time to advance time step: 2.020261095
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.590357118
Call to estdt for level 0 gives dt_lev = 1.641616814
Minimum estdt over all levels = 1.641616814
Call to estdt at beginning of step 3 gives dt =1.641616814
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 = 146032253.6
MLMG: Initial residual (resid0) = 146032253.6
MLMG: Final Iter. 8 resid, resid/bnorm = 0.01268048584, 8.683345996e-11
MLMG: Timers: Solve = 0.042119371 Iter = 0.039087638 Bottom = 0.000194162
<<< 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.7
MLMG: Initial residual (resid0) = 3975.129505
MLMG: Final Iter. 3 resid, resid/bnorm = 0.005708731944, 3.909226789e-11
MLMG: Timers: Solve = 0.018064332 Iter = 0.015050968 Bottom = 8.2018e-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.629
MLMG: Initial residual (resid0) = 5204170.629
MLMG: Final Iter. 19 resid, resid/bnorm = 3.883335739e-05, 7.461968518e-12
MLMG: Timers: Solve = 0.320026335 Iter = 0.31509378 Bottom = 0.000261598
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.527277194 seconds
MAC Proj :0.071924422 seconds
Nodal Proj :0.346021567 seconds
Reactions :0.840060572 seconds
Misc :0.137978064 seconds
Base State :0.002208801 seconds
Time to advance time step: 1.923635037
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.576770901
Total Time: 14.87135386
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 (23.01-27-g756468d42106) finalized