MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.01-7-g7a313e260ccf) 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.569198425 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. 9 resid, resid/bnorm = 0.0036513675, 8.763184089e-11 MLMG: Timers: Solve = 0.152188021 Iter = 0.148263833 Bottom = 0.000166103 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.534345701 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.589968046, 2.042875028e-07 MLMG: Timers: Solve = 0.036413735 Iter = 0.032534383 Bottom = 4.1321e-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.79 MLMG: Initial residual (resid0) = 13466729.79 MLMG: Final Iter. 3 resid, resid/bnorm = 0.01553362957, 1.153481938e-09 MLMG: Timers: Solve = 0.052980415 Iter = 0.04906937 Bottom = 5.2624e-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.458 MLMG: Initial residual (resid0) = 9200997.458 MLMG: Final Iter. 4 resid, resid/bnorm = 0.0007020710036, 7.63037928e-11 MLMG: Timers: Solve = 0.068692666 Iter = 0.064861759 Bottom = 6.9073e-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.527070529 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 = 156168945.8 MLMG: Initial residual (resid0) = 156168945.8 MLMG: Final Iter. 35 resid, resid/bnorm = 0.01235002279, 7.908116898e-11 MLMG: Timers: Solve = 0.169083125 Iter = 0.165621701 Bottom = 0.00085474 <<< 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 = 156168960.3 MLMG: Initial residual (resid0) = 3426.279444 MLMG: Final Iter. 5 resid, resid/bnorm = 0.01363201439, 8.729016552e-11 MLMG: Timers: Solve = 0.025768915 Iter = 0.023818886 Bottom = 0.000140231 <<< 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 = 140238249.3 MLMG: Initial residual (resid0) = 140238249.3 MLMG: Final Iter. 7 resid, resid/bnorm = 0.0003455430269, 2.463971339e-12 MLMG: Timers: Solve = 0.117299807 Iter = 0.113419489 Bottom = 0.000112925 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.4526733 seconds MAC Proj :0.204971963 seconds Nodal Proj :0.140091697 seconds Reactions :0.802493127 seconds Misc :0.127881416 seconds Base State :0.001825498 seconds Time to advance time step: 1.755354238 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.539129287 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.1 MLMG: Initial residual (resid0) = 156143863.1 MLMG: Final Iter. 34 resid, resid/bnorm = 0.01555861533, 9.964282307e-11 MLMG: Timers: Solve = 0.162133598 Iter = 0.160187001 Bottom = 0.000836509 <<< 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 = 156143862.9 MLMG: Initial residual (resid0) = 37.32845885 MLMG: Final Iter. 8 resid, resid/bnorm = 0.01373890042, 8.798873149e-11 MLMG: Timers: Solve = 0.039565204 Iter = 0.037639654 Bottom = 0.000202206 <<< 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.776 MLMG: Initial residual (resid0) = 7408896.776 MLMG: Final Iter. 5 resid, resid/bnorm = 9.003269952e-06, 1.21519711e-12 MLMG: Timers: Solve = 0.084302475 Iter = 0.080348791 Bottom = 8.279e-05 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.423406267 seconds MAC Proj :0.210541358 seconds Nodal Proj :0.10703819 seconds Reactions :0.784867838 seconds Misc :0.124785537 seconds Base State :0.001864026001 seconds Time to advance time step: 1.69037053 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.534350093 Call to estdt for level 0 gives dt_lev = 1.338765343 Minimum estdt over all levels = 1.338765343 Call to estdt at beginning of step 2 gives dt =1.338765343 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 = 148819014 MLMG: Initial residual (resid0) = 148819014 MLMG: Final Iter. 8 resid, resid/bnorm = 0.009800747037, 6.585682012e-11 MLMG: Timers: Solve = 0.039874063 Iter = 0.037936868 Bottom = 0.000195899 <<< 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 = 148819013.1 MLMG: Initial residual (resid0) = 56.34133842 MLMG: Final Iter. 2 resid, resid/bnorm = 0.01291100925, 8.67564499e-11 MLMG: Timers: Solve = 0.011564767 Iter = 0.009629555 Bottom = 5.8485e-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.54 MLMG: Initial residual (resid0) = 6145962.54 MLMG: Final Iter. 25 resid, resid/bnorm = 5.331914872e-05, 8.675475709e-12 MLMG: Timers: Solve = 0.408376945 Iter = 0.404387687 Bottom = 0.000348419 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.425252176 seconds MAC Proj :0.060326362 seconds Nodal Proj :0.431303843 seconds Reactions :0.787596588 seconds Misc :0.123550632 seconds Base State :0.001843722 seconds Time to advance time step: 1.82835159 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.537338276 Call to estdt for level 0 gives dt_lev = 1.641626298 Minimum estdt over all levels = 1.641626298 Call to estdt at beginning of step 3 gives dt =1.641626298 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.4 MLMG: Initial residual (resid0) = 146032253.4 MLMG: Final Iter. 8 resid, resid/bnorm = 0.01270103455, 8.697417354e-11 MLMG: Timers: Solve = 0.040078234 Iter = 0.038144538 Bottom = 0.000196717 <<< 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.5 MLMG: Initial residual (resid0) = 216.7870271 MLMG: Final Iter. 3 resid, resid/bnorm = 0.003212187439, 2.199642465e-11 MLMG: Timers: Solve = 0.01619976 Iter = 0.014267172 Bottom = 8.1528e-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.970926628e-05, 7.630277541e-12 MLMG: Timers: Solve = 0.309817174 Iter = 0.305813843 Bottom = 0.000262534 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.428511313 seconds MAC Proj :0.065254518 seconds Nodal Proj :0.332731187 seconds Reactions :0.789036077 seconds Misc :0.124075296 seconds Base State :0.001837725 seconds Time to advance time step: 1.739924598 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.540500448 Total Time: 13.71649136 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-7-g7a313e260ccf) finalized