MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.06-39-g2d931f63cb4d) 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.540337862 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. 6 resid, resid/bnorm = 0.007489142939, 9.070815886e-11 MLMG: Timers: Solve = 0.101670784 Iter = 0.096289159 Bottom = 0.000108251 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.512393054 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.61 MLMG: Initial residual (resid0) = 33867874.61 MLMG: Final Iter. 2 resid, resid/bnorm = 3.808735294, 1.124586452e-07 MLMG: Timers: Solve = 0.035222628 Iter = 0.031389145 Bottom = 3.8167e-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.0124410399, 6.701778652e-10 MLMG: Timers: Solve = 0.051250836 Iter = 0.047418001 Bottom = 5.0249e-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. 4 resid, resid/bnorm = 0.0008234889247, 7.0310025e-11 MLMG: Timers: Solve = 0.066021388 Iter = 0.062194227 Bottom = 6.0978e-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.495245776 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 = 97092862.8 MLMG: Initial residual (resid0) = 97092862.8 MLMG: Final Iter. 37 resid, resid/bnorm = 0.008145976812, 8.389882198e-11 MLMG: Timers: Solve = 0.172932641 Iter = 0.169845249 Bottom = 0.000934985 <<< 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 = 97092854.46 MLMG: Initial residual (resid0) = 4268.184198 MLMG: Final Iter. 12 resid, resid/bnorm = 0.008214619011, 8.460580396e-11 MLMG: Timers: Solve = 0.058400498 Iter = 0.055340402 Bottom = 0.000300294 <<< 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 = 140045958 MLMG: Initial residual (resid0) = 140045958 MLMG: Final Iter. 9 resid, resid/bnorm = 0.000321155414, 2.293214446e-12 MLMG: Timers: Solve = 0.146446117 Iter = 0.142386014 Bottom = 0.000128248 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.455647369 seconds MAC Proj :0.241745403 seconds Nodal Proj :0.169794599 seconds Reactions :0.710080342 seconds Misc :0.124292739 seconds Base State :0.002053955 seconds Time to advance time step: 1.725626393 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.511675067 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 = 97077994.47 MLMG: Initial residual (resid0) = 97077994.47 MLMG: Final Iter. 37 resid, resid/bnorm = 0.008088413626, 8.331871368e-11 MLMG: Timers: Solve = 0.172789547 Iter = 0.169739931 Bottom = 0.000941674 <<< 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 = 97077994.27 MLMG: Initial residual (resid0) = 32.42851122 MLMG: Final Iter. 12 resid, resid/bnorm = 0.008653838187, 8.914314982e-11 MLMG: Timers: Solve = 0.058285485 Iter = 0.055214754 Bottom = 0.000302215 <<< 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.636 MLMG: Initial residual (resid0) = 8056355.636 MLMG: Final Iter. 5 resid, resid/bnorm = 3.392901272e-05, 4.211459157e-12 MLMG: Timers: Solve = 0.082834128 Iter = 0.078828786 Bottom = 8.2939e-05 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.452844435 seconds MAC Proj :0.241537632 seconds Nodal Proj :0.107040048 seconds Reactions :0.708793722 seconds Misc :0.124878189 seconds Base State :0.001942925 seconds Time to advance time step: 1.677033904 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.50197831 Call to estdt for level 0 gives dt_lev = 1.410018681 Minimum estdt over all levels = 1.410018681 Call to estdt at beginning of step 2 gives dt =1.410018681 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 = 92795854.21 MLMG: Initial residual (resid0) = 92795854.21 MLMG: Final Iter. 8 resid, resid/bnorm = 0.006667355075, 7.18497085e-11 MLMG: Timers: Solve = 0.040203036 Iter = 0.037125256 Bottom = 0.000195991 <<< 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 = 92795854.28 MLMG: Initial residual (resid0) = 236.1724532 MLMG: Final Iter. 3 resid, resid/bnorm = 0.008556066081, 9.2203107e-11 MLMG: Timers: Solve = 0.017222725 Iter = 0.014142361 Bottom = 7.9307e-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.836 MLMG: Initial residual (resid0) = 5904880.836 MLMG: Final Iter. 28 resid, resid/bnorm = 4.61153686e-05, 7.809703512e-12 MLMG: Timers: Solve = 0.441656321 Iter = 0.437646749 Bottom = 0.000370763 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.456373545 seconds MAC Proj :0.067973008 seconds Nodal Proj :0.465342326 seconds Reactions :0.711336535 seconds Misc :0.124666113 seconds Base State :0.001911826 seconds Time to advance time step: 1.827765678 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.504386184 Call to estdt for level 0 gives dt_lev = 1.702939612 Minimum estdt over all levels = 1.702939612 Call to estdt at beginning of step 3 gives dt =1.702939612 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 = 91006742.91 MLMG: Initial residual (resid0) = 91006742.91 MLMG: Final Iter. 9 resid, resid/bnorm = 0.005843035877, 6.420442805e-11 MLMG: Timers: Solve = 0.044663326 Iter = 0.041565333 Bottom = 0.00021059 <<< 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 = 91006742.62 MLMG: Initial residual (resid0) = 2.918235196 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005682192743, 6.243705224e-11 MLMG: Timers: Solve = 0.021777522 Iter = 0.018679203 Bottom = 0.000100045 <<< 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.308 MLMG: Initial residual (resid0) = 4615830.308 MLMG: Final Iter. 18 resid, resid/bnorm = 3.66307795e-05, 7.935902548e-12 MLMG: Timers: Solve = 0.284467963 Iter = 0.280417638 Bottom = 0.000244827 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.459181111 seconds MAC Proj :0.076985209 seconds Nodal Proj :0.308087798 seconds Reactions :0.710903897 seconds Misc :0.125322849 seconds Base State :0.001938328 seconds Time to advance time step: 1.682583016 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.506277896 Total Time: 13.03750334 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.06-39-g2d931f63cb4d) finalized