MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.09-18-g3e5cc7780280) 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.527937819
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.30325468 Iter = 0.297681627 Bottom = 0.000275089
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.500946915
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.035285421 Iter = 0.031738484 Bottom = 5.7107e-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.049864926 Iter = 0.046238452 Bottom = 5.245e-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.091099892 Iter = 0.086841115 Bottom = 8.302e-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.503642843
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.040354811 Iter = 0.037006504 Bottom = 0.000202807
<<< 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.020464455 Iter = 0.018691143 Bottom = 0.000107142
<<< 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.46074727 Iter = 0.457007832 Bottom = 0.000396242
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.434198203 seconds
MAC Proj :0.069855361 seconds
Nodal Proj :0.483638625 seconds
Reactions :0.754453601 seconds
Misc :0.123200195 seconds
Base State :0.002016093 seconds
Time to advance time step: 1.891891018
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.49394036
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.038963499 Iter = 0.037202585 Bottom = 0.000201911
<<< 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.015581735 Iter = 0.013811737 Bottom = 8.1521e-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.285970915 Iter = 0.282250379 Bottom = 0.000268544
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.401316632 seconds
MAC Proj :0.062282296 seconds
Nodal Proj :0.30863954 seconds
Reactions :0.719556614 seconds
Misc :0.120135897 seconds
Base State :0.001890516001 seconds
Time to advance time step: 1.63835872
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.504581351
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.043489056 Iter = 0.041684643 Bottom = 0.000217437
<<< 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.015927512 Iter = 0.014174076 Bottom = 8.2483e-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.316672668 Iter = 0.312970592 Bottom = 0.000283718
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.409422023 seconds
MAC Proj :0.067186303 seconds
Nodal Proj :0.339229266 seconds
Reactions :0.723299467 seconds
Misc :0.119408898 seconds
Base State :0.001882386 seconds
Time to advance time step: 1.658836862
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.505887924
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.043499834 Iter = 0.041739867 Bottom = 0.000220717
<<< 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.015672925 Iter = 0.013873228 Bottom = 8.5486e-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.331257362 Iter = 0.32767408 Bottom = 0.000290811
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.406768704 seconds
MAC Proj :0.067047319 seconds
Nodal Proj :0.353816601 seconds
Reactions :0.719793208 seconds
Misc :0.119946606 seconds
Base State :0.002759342 seconds
Time to advance time step: 1.667680699
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.507689697
Total Time: 13.22894467
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-18-g3e5cc7780280) finalized