MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.02-10-g76d08651adb9) 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.576841438
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. 9 resid, resid/bnorm = 0.004268283956, 5.169726126e-11
MLMG: Timers: Solve = 0.156057784 Iter = 0.150661658 Bottom = 0.000144928
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.50953821
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.815701978, 1.12664347e-07
MLMG: Timers: Solve = 0.042270593 Iter = 0.03775189 Bottom = 4.1093e-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.14
MLMG: Initial residual (resid0) = 18563788.14
MLMG: Final Iter. 3 resid, resid/bnorm = 0.01394330838, 7.511025379e-10
MLMG: Timers: Solve = 0.054408394 Iter = 0.05035826 Bottom = 5.2268e-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.0005240472965, 4.474350222e-11
MLMG: Timers: Solve = 0.09288904 Iter = 0.088572469 Bottom = 8.6292e-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.506759043
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 = 97092886.47
MLMG: Initial residual (resid0) = 97092886.47
MLMG: Final Iter. 38 resid, resid/bnorm = 0.007781125605, 8.014104728e-11
MLMG: Timers: Solve = 0.180589828 Iter = 0.177468013 Bottom = 0.000994861
<<< 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 = 97092878.36
MLMG: Initial residual (resid0) = 1102.546728
MLMG: Final Iter. 9 resid, resid/bnorm = 0.008318845183, 8.56792519e-11
MLMG: Timers: Solve = 0.044586476 Iter = 0.042335113 Bottom = 0.000247742
<<< 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 = 140046438.7
MLMG: Initial residual (resid0) = 140046438.7
MLMG: Final Iter. 6 resid, resid/bnorm = 0.0009652432054, 6.892308112e-12
MLMG: Timers: Solve = 0.104724847 Iter = 0.100647209 Bottom = 9.5576e-05
Done calling nodal solver
Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.462353455 seconds
MAC Proj :0.234853603 seconds
Nodal Proj :0.128769563 seconds
Reactions :0.729971971 seconds
Misc :0.127251777 seconds
Base State :0.002012366 seconds
Time to advance time step: 1.706467676
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.507953642
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 = 97078018.02
MLMG: Initial residual (resid0) = 97078018.02
MLMG: Final Iter. 38 resid, resid/bnorm = 0.007706178352, 7.938129052e-11
MLMG: Timers: Solve = 0.180072975 Iter = 0.177813402 Bottom = 0.001001617
<<< 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 = 97078017.86
MLMG: Initial residual (resid0) = 27.61699302
MLMG: Final Iter. 9 resid, resid/bnorm = 0.008268721402, 8.517604278e-11
MLMG: Timers: Solve = 0.044584589 Iter = 0.042338139 Bottom = 0.000254928
<<< 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.33
MLMG: Initial residual (resid0) = 8056355.33
MLMG: Final Iter. 5 resid, resid/bnorm = 2.859182314e-05, 3.548977418e-12
MLMG: Timers: Solve = 0.088341942 Iter = 0.084172566 Bottom = 9.0087e-05
Done calling nodal solver
Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117
Timing summary:
Advection :0.447966346 seconds
MAC Proj :0.233543947 seconds
Nodal Proj :0.112619241 seconds
Reactions :0.748492658 seconds
Misc :0.12737785 seconds
Base State :0.0021025 seconds
Time to advance time step: 1.694614784
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.518897871
Call to estdt for level 0 gives dt_lev = 1.410016905
Minimum estdt over all levels = 1.410016905
Call to estdt at beginning of step 2 gives dt =1.410016905
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 = 92795868.89
MLMG: Initial residual (resid0) = 92795868.89
MLMG: Final Iter. 9 resid, resid/bnorm = 0.004433698952, 4.777905531e-11
MLMG: Timers: Solve = 0.044598225 Iter = 0.042356535 Bottom = 0.000220111
<<< 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 = 92795868.73
MLMG: Initial residual (resid0) = 3.445192834
MLMG: Final Iter. 4 resid, resid/bnorm = 0.00454864651, 4.901776956e-11
MLMG: Timers: Solve = 0.021211852 Iter = 0.018955887 Bottom = 0.000103389
<<< 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.513
MLMG: Initial residual (resid0) = 5904880.513
MLMG: Final Iter. 22 resid, resid/bnorm = 5.201017484e-05, 8.807997846e-12
MLMG: Timers: Solve = 0.372884212 Iter = 0.368778646 Bottom = 0.000306439
Done calling nodal solver
Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529
Timing summary:
Advection :0.449168984 seconds
MAC Proj :0.074659022 seconds
Nodal Proj :0.396919446 seconds
Reactions :0.734901461 seconds
Misc :0.126931423 seconds
Base State :0.001917096999 seconds
Time to advance time step: 1.782918395
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.520371138
Call to estdt for level 0 gives dt_lev = 1.702931242
Minimum estdt over all levels = 1.702931242
Call to estdt at beginning of step 3 gives dt =1.702931242
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 = 91006735.7
MLMG: Initial residual (resid0) = 91006735.7
MLMG: Final Iter. 9 resid, resid/bnorm = 0.005086801946, 5.589478523e-11
MLMG: Timers: Solve = 0.04457329 Iter = 0.042336319 Bottom = 0.00021935
<<< 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 = 91006735.41
MLMG: Initial residual (resid0) = 2.889003656
MLMG: Final Iter. 4 resid, resid/bnorm = 0.004876181483, 5.358044612e-11
MLMG: Timers: Solve = 0.021171679 Iter = 0.018921617 Bottom = 0.000104802
<<< 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.313
MLMG: Initial residual (resid0) = 4615830.313
MLMG: Final Iter. 20 resid, resid/bnorm = 4.509603605e-05, 9.769864356e-12
MLMG: Timers: Solve = 0.339801915 Iter = 0.335711895 Bottom = 0.000286445
Done calling nodal solver
Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082
Timing summary:
Advection :0.449477061 seconds
MAC Proj :0.074611591 seconds
Nodal Proj :0.364170898 seconds
Reactions :0.735887317 seconds
Misc :0.127115953 seconds
Base State :0.001913038999 seconds
Time to advance time step: 1.751596849
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.521746698
Total Time: 13.52812147
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.02-10-g76d08651adb9) finalized