MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.03-17-g7b4d10daf9c4) 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.5226889 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. 7 resid, resid/bnorm = 0.006740517914, 8.164084659e-11 MLMG: Timers: Solve = 0.123037874 Iter = 0.117782206 Bottom = 0.000134804 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.510125654 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.809456168, 1.1247993e-07 MLMG: Timers: Solve = 0.037530484 Iter = 0.033751452 Bottom = 4.2747e-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.01186252502, 6.390142429e-10 MLMG: Timers: Solve = 0.054227966 Iter = 0.050425936 Bottom = 5.5614e-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.76 MLMG: Initial residual (resid0) = 11712254.76 MLMG: Final Iter. 4 resid, resid/bnorm = 0.000591928605, 5.053925289e-11 MLMG: Timers: Solve = 0.071263231 Iter = 0.067478533 Bottom = 7.1274e-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.497702303 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 = 97092871.11 MLMG: Initial residual (resid0) = 97092871.11 MLMG: Final Iter. 35 resid, resid/bnorm = 0.007416933775, 7.63900963e-11 MLMG: Timers: Solve = 0.163755887 Iter = 0.160774182 Bottom = 0.000921468 <<< 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 = 97092863.31 MLMG: Initial residual (resid0) = 1047.097046 MLMG: Final Iter. 12 resid, resid/bnorm = 0.008365057409, 8.615522422e-11 MLMG: Timers: Solve = 0.057056638 Iter = 0.055061965 Bottom = 0.000324366 <<< 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 = 140047068.4 MLMG: Initial residual (resid0) = 140047068.4 MLMG: Final Iter. 9 resid, resid/bnorm = 0.0008406338748, 6.002509616e-12 MLMG: Timers: Solve = 0.153823548 Iter = 0.150006768 Bottom = 0.000158313 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.43298483 seconds MAC Proj :0.230018949 seconds Nodal Proj :0.177180955 seconds Reactions :0.714363741 seconds Misc :0.122830504 seconds Base State :0.001942264 seconds Time to advance time step: 1.70189322 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.494079922 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 = 97078002.75 MLMG: Initial residual (resid0) = 97078002.75 MLMG: Final Iter. 34 resid, resid/bnorm = 0.009658396244, 9.949108933e-11 MLMG: Timers: Solve = 0.157162948 Iter = 0.155154572 Bottom = 0.000878896 <<< 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 = 97078002.63 MLMG: Initial residual (resid0) = 116.4032 MLMG: Final Iter. 12 resid, resid/bnorm = 0.008642971516, 8.903120461e-11 MLMG: Timers: Solve = 0.057078992 Iter = 0.055085756 Bottom = 0.000326314 <<< 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 = 8056354.934 MLMG: Initial residual (resid0) = 8056354.934 MLMG: Final Iter. 7 resid, resid/bnorm = 2.363556996e-05, 2.933779625e-12 MLMG: Timers: Solve = 0.121271443 Iter = 0.117447722 Bottom = 0.000120651 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.414679733 seconds MAC Proj :0.222554241 seconds Nodal Proj :0.144448549 seconds Reactions :0.710845639 seconds Misc :0.122090569 seconds Base State :0.001855944 seconds Time to advance time step: 1.639241435 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.505515711 Call to estdt for level 0 gives dt_lev = 1.410017548 Minimum estdt over all levels = 1.410017548 Call to estdt at beginning of step 2 gives dt =1.410017548 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 = 92795856.71 MLMG: Initial residual (resid0) = 92795856.71 MLMG: Final Iter. 9 resid, resid/bnorm = 0.007149539888, 7.704589559e-11 MLMG: Timers: Solve = 0.043490098 Iter = 0.041436711 Bottom = 0.000224313 <<< 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 = 92795856.26 MLMG: Initial residual (resid0) = 1015.497526 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005852021277, 6.306339003e-11 MLMG: Timers: Solve = 0.020583041 Iter = 0.018545253 Bottom = 0.000109437 <<< 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.084 MLMG: Initial residual (resid0) = 5904880.084 MLMG: Final Iter. 29 resid, resid/bnorm = 4.671188071e-05, 7.910724697e-12 MLMG: Timers: Solve = 0.479141138 Iter = 0.475308215 Bottom = 0.000397467 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.42974437 seconds MAC Proj :0.07294394 seconds Nodal Proj :0.502838122 seconds Reactions :0.713282928 seconds Misc :0.12284712 seconds Base State :0.00206476 seconds Time to advance time step: 1.842129802 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.510112034 Call to estdt for level 0 gives dt_lev = 1.702937995 Minimum estdt over all levels = 1.702937995 Call to estdt at beginning of step 3 gives dt =1.702937995 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 = 91006740.53 MLMG: Initial residual (resid0) = 91006740.53 MLMG: Final Iter. 8 resid, resid/bnorm = 0.007962588221, 8.749448859e-11 MLMG: Timers: Solve = 0.039531482 Iter = 0.037474169 Bottom = 0.000224422 <<< 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 = 91006740.25 MLMG: Initial residual (resid0) = 2.846421008 MLMG: Final Iter. 3 resid, resid/bnorm = 0.008897874504, 9.777159889e-11 MLMG: Timers: Solve = 0.015853767 Iter = 0.013826496 Bottom = 8.5433e-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.324 MLMG: Initial residual (resid0) = 4615830.324 MLMG: Final Iter. 24 resid, resid/bnorm = 3.47928144e-05, 7.537715202e-12 MLMG: Timers: Solve = 0.405379553 Iter = 0.401410532 Bottom = 0.000455294 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.433221499 seconds MAC Proj :0.064309714 seconds Nodal Proj :0.428962266 seconds Reactions :0.732885909 seconds Misc :0.124340502 seconds Base State :0.001961544 seconds Time to advance time step: 1.784048086 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.515164103 Total Time: 13.11756252 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.03-17-g7b4d10daf9c4) finalized