MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.06-36-g478fd8a4ac98) 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.549625609 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. 8 resid, resid/bnorm = 0.004905575886, 5.941611215e-11 MLMG: Timers: Solve = 0.136638677 Iter = 0.131496878 Bottom = 0.000146399 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.501939863 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.813934696, 1.126121654e-07 MLMG: Timers: Solve = 0.035361791 Iter = 0.03168496 Bottom = 5.8984e-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.01313454029, 7.075355631e-10 MLMG: Timers: Solve = 0.051818997 Iter = 0.048157747 Bottom = 5.2785e-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.001105872449, 9.442011567e-11 MLMG: Timers: Solve = 0.066776503 Iter = 0.063157436 Bottom = 6.5416e-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.501071349 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 = 97092872.47 MLMG: Initial residual (resid0) = 97092872.47 MLMG: Final Iter. 36 resid, resid/bnorm = 0.008748561144, 9.010508105e-11 MLMG: Timers: Solve = 0.171603081 Iter = 0.168558646 Bottom = 0.000597181 <<< 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 = 97092864.37 MLMG: Initial residual (resid0) = 1103.178003 MLMG: Final Iter. 5 resid, resid/bnorm = 0.006639543921, 6.838343851e-11 MLMG: Timers: Solve = 0.024828607 Iter = 0.022920765 Bottom = 0.000133624 <<< 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 = 140046451.7 MLMG: Initial residual (resid0) = 140046451.7 MLMG: Final Iter. 8 resid, resid/bnorm = 0.0005394299515, 3.851793065e-12 MLMG: Timers: Solve = 0.129640082 Iter = 0.126010398 Bottom = 0.00011793 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.439826664 seconds MAC Proj :0.2057797 seconds Nodal Proj :0.152949487 seconds Reactions :0.715551721 seconds Misc :0.123084189 seconds Base State :0.001925195 seconds Time to advance time step: 1.663576646 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.496805267 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 = 97078004.09 MLMG: Initial residual (resid0) = 97078004.09 MLMG: Final Iter. 36 resid, resid/bnorm = 0.008675698191, 8.936832059e-11 MLMG: Timers: Solve = 0.166362725 Iter = 0.164459002 Bottom = 0.000589641 <<< 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 = 97078003.92 MLMG: Initial residual (resid0) = 29.40087473 MLMG: Final Iter. 5 resid, resid/bnorm = 0.005938638002, 6.117387834e-11 MLMG: Timers: Solve = 0.025010916 Iter = 0.023090937 Bottom = 0.000131129 <<< 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.323 MLMG: Initial residual (resid0) = 8056355.323 MLMG: Final Iter. 5 resid, resid/bnorm = 3.149034455e-05, 3.908758153e-12 MLMG: Timers: Solve = 0.08511279 Iter = 0.081377462 Bottom = 8.1275e-05 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.416669117 seconds MAC Proj :0.199380124 seconds Nodal Proj :0.108556421 seconds Reactions :0.71450524 seconds Misc :0.121827626 seconds Base State :0.001869714 seconds Time to advance time step: 1.585348091 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.505834321 Call to estdt for level 0 gives dt_lev = 1.410018253 Minimum estdt over all levels = 1.410018253 Call to estdt at beginning of step 2 gives dt =1.410018253 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 = 92795860.9 MLMG: Initial residual (resid0) = 92795860.9 MLMG: Final Iter. 9 resid, resid/bnorm = 0.004699967802, 5.06484638e-11 MLMG: Timers: Solve = 0.04288244 Iter = 0.040958338 Bottom = 0.000213511 <<< 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 = 92795860.72 MLMG: Initial residual (resid0) = 101.6938298 MLMG: Final Iter. 4 resid, resid/bnorm = 0.00459587574, 4.952673217e-11 MLMG: Timers: Solve = 0.021181976 Iter = 0.0192384 Bottom = 0.00010457 <<< 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.497 MLMG: Initial residual (resid0) = 5904880.497 MLMG: Final Iter. 19 resid, resid/bnorm = 5.43827191e-05, 9.209791651e-12 MLMG: Timers: Solve = 0.306618451 Iter = 0.302871856 Bottom = 0.000260369 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.422009912 seconds MAC Proj :0.072390121 seconds Nodal Proj :0.32992736 seconds Reactions :0.717687251 seconds Misc :0.122264081 seconds Base State :0.001836716999 seconds Time to advance time step: 1.664616871 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.505368024 Call to estdt for level 0 gives dt_lev = 1.702936365 Minimum estdt over all levels = 1.702936365 Call to estdt at beginning of step 3 gives dt =1.702936365 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.02 MLMG: Initial residual (resid0) = 91006740.02 MLMG: Final Iter. 9 resid, resid/bnorm = 0.004960529506, 5.450727611e-11 MLMG: Timers: Solve = 0.043132112 Iter = 0.041203068 Bottom = 0.000215708 <<< 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 = 91006739.73 MLMG: Initial residual (resid0) = 2.824363064 MLMG: Final Iter. 4 resid, resid/bnorm = 0.004936657846, 5.424496978e-11 MLMG: Timers: Solve = 0.020395678 Iter = 0.018473326 Bottom = 0.000103996 <<< 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. 21 resid, resid/bnorm = 3.883824684e-05, 8.414140955e-12 MLMG: Timers: Solve = 0.331654953 Iter = 0.327934395 Bottom = 0.000282168 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.416447125 seconds MAC Proj :0.071536358 seconds Nodal Proj :0.355442973 seconds Reactions :0.70749647 seconds Misc :0.122914236 seconds Base State :0.001952038001 seconds Time to advance time step: 1.674183502 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.508323069 Total Time: 12.74229133 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-36-g478fd8a4ac98) finalized