MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.08-13-gbd5f6a9f6a1a) 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.545200993
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.319469029 Iter = 0.312687347 Bottom = 0.000325976
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.58277647
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.037568837 Iter = 0.032210581 Bottom = 4.7114e-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.052736744 Iter = 0.048221047 Bottom = 5.5371e-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.082472464 Iter = 0.07799602 Bottom = 7.8113e-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.519780094
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.040570642 Iter = 0.03638441 Bottom = 0.000217061
<<< 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.023043385 Iter = 0.01854815 Bottom = 0.00010913
<<< 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.465303406 Iter = 0.46102594 Bottom = 0.000412026
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.52365622 seconds
MAC Proj :0.076809588 seconds
Nodal Proj :0.492566704 seconds
Reactions :0.720893523 seconds
Misc :0.132712167 seconds
Base State :0.002279809 seconds
Time to advance time step: 1.978176742
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.586054107
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.041390601 Iter = 0.036688688 Bottom = 0.000203697
<<< 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.019004629 Iter = 0.014289186 Bottom = 8.1748e-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.287548504 Iter = 0.282894471 Bottom = 0.000257706
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.553666135 seconds
MAC Proj :0.074481055 seconds
Nodal Proj :0.313530815 seconds
Reactions :0.763963218 seconds
Misc :0.137906095 seconds
Base State :0.002756389 seconds
Time to advance time step: 1.876935496
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.525905996
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.046330547 Iter = 0.041593896 Bottom = 0.000215636
<<< 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.019206738 Iter = 0.014311653 Bottom = 8.0788e-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.324602693 Iter = 0.320034673 Bottom = 0.000282295
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.564253172 seconds
MAC Proj :0.079803642 seconds
Nodal Proj :0.350945514 seconds
Reactions :0.726009207 seconds
Misc :0.136979147 seconds
Base State :0.002634042 seconds
Time to advance time step: 1.861505633
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.598995631
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.046468736 Iter = 0.041757266 Bottom = 0.000222388
<<< 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.019442425 Iter = 0.014534341 Bottom = 8.262e-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.339116481 Iter = 0.334413824 Bottom = 0.000297741
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.585479344 seconds
MAC Proj :0.08041012 seconds
Nodal Proj :0.366169989 seconds
Reactions :0.732577911 seconds
Misc :0.13717922 seconds
Base State :0.002349724 seconds
Time to advance time step: 1.905311648
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.537325512
Total Time: 15.07366559
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.08-13-gbd5f6a9f6a1a) finalized