MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.04-7-gec648a92a89f) 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.51503431 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.007488844916, 9.070454922e-11 MLMG: Timers: Solve = 0.101816735 Iter = 0.096621299 Bottom = 0.000106212 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.51009177 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.807062387, 1.124092501e-07 MLMG: Timers: Solve = 0.035083376 Iter = 0.031445907 Bottom = 4.1726e-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.01152177935, 6.206588477e-10 MLMG: Timers: Solve = 0.050928488 Iter = 0.04721171 Bottom = 5.3569e-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.0004292526282, 3.66498712e-11 MLMG: Timers: Solve = 0.067615104 Iter = 0.063990321 Bottom = 6.9719e-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.500185447 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 = 97092811.98 MLMG: Initial residual (resid0) = 97092811.98 MLMG: Final Iter. 34 resid, resid/bnorm = 0.007426172495, 7.648529632e-11 MLMG: Timers: Solve = 0.156791165 Iter = 0.153931843 Bottom = 0.000819787 <<< 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 = 97092803.82 MLMG: Initial residual (resid0) = 3085.592249 MLMG: Final Iter. 9 resid, resid/bnorm = 0.008550554514, 8.806579043e-11 MLMG: Timers: Solve = 0.042643376 Iter = 0.040892955 Bottom = 0.000231868 <<< 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 = 140046317 MLMG: Initial residual (resid0) = 140046317 MLMG: Final Iter. 10 resid, resid/bnorm = 0.0006782647688, 4.843146061e-12 MLMG: Timers: Solve = 0.159616655 Iter = 0.156026709 Bottom = 0.000148723 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.417453092 seconds MAC Proj :0.208068413 seconds Nodal Proj :0.18190838 seconds Reactions :0.709198714 seconds Misc :0.120244892 seconds Base State :0.001986291 seconds Time to advance time step: 1.661406267 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.492586565 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 = 97077943.83 MLMG: Initial residual (resid0) = 97077943.83 MLMG: Final Iter. 34 resid, resid/bnorm = 0.007367417216, 7.589177238e-11 MLMG: Timers: Solve = 0.156945304 Iter = 0.155182888 Bottom = 0.000820575 <<< 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 = 97077943.66 MLMG: Initial residual (resid0) = 29.70406509 MLMG: Final Iter. 9 resid, resid/bnorm = 0.0079767555, 8.216856682e-11 MLMG: Timers: Solve = 0.042657475 Iter = 0.040791269 Bottom = 0.00023214 <<< 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.406 MLMG: Initial residual (resid0) = 8056355.406 MLMG: Final Iter. 7 resid, resid/bnorm = 2.63066031e-05, 3.265323062e-12 MLMG: Timers: Solve = 0.114171648 Iter = 0.11051136 Bottom = 0.000114112 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.397676336 seconds MAC Proj :0.207208048 seconds Nodal Proj :0.136990474 seconds Reactions :0.707926095 seconds Misc :0.120095973 seconds Base State :0.001856637 seconds Time to advance time step: 1.594160069 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.501197775 Call to estdt for level 0 gives dt_lev = 1.410025614 Minimum estdt over all levels = 1.410025614 Call to estdt at beginning of step 2 gives dt =1.410025614 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 = 92795823.99 MLMG: Initial residual (resid0) = 92795823.99 MLMG: Final Iter. 9 resid, resid/bnorm = 0.004293054342, 4.626344331e-11 MLMG: Timers: Solve = 0.04265536 Iter = 0.040896156 Bottom = 0.000217202 <<< 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 = 92795823.86 MLMG: Initial residual (resid0) = 19.33500622 MLMG: Final Iter. 4 resid, resid/bnorm = 0.004988048226, 5.375293863e-11 MLMG: Timers: Solve = 0.020691885 Iter = 0.018893647 Bottom = 0.000110096 <<< 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.57 MLMG: Initial residual (resid0) = 5904880.57 MLMG: Final Iter. 24 resid, resid/bnorm = 4.399148747e-05, 7.450021546e-12 MLMG: Timers: Solve = 0.385188563 Iter = 0.381253369 Bottom = 0.000333404 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.407885826 seconds MAC Proj :0.071076784 seconds Nodal Proj :0.407988399 seconds Reactions :0.714264922 seconds Misc :0.122209453 seconds Base State :0.001824482 seconds Time to advance time step: 1.7237648 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.508747564 Call to estdt for level 0 gives dt_lev = 1.70296962 Minimum estdt over all levels = 1.70296962 Call to estdt at beginning of step 3 gives dt =1.70296962 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 = 91006761.75 MLMG: Initial residual (resid0) = 91006761.75 MLMG: Final Iter. 9 resid, resid/bnorm = 0.005128346384, 5.635126759e-11 MLMG: Timers: Solve = 0.042518122 Iter = 0.040756323 Bottom = 0.000216067 <<< 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 = 91006761.48 MLMG: Initial residual (resid0) = 3.013281926 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005550988019, 6.099533626e-11 MLMG: Timers: Solve = 0.020617498 Iter = 0.018861116 Bottom = 0.000108887 <<< 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.334 MLMG: Initial residual (resid0) = 4615830.334 MLMG: Final Iter. 20 resid, resid/bnorm = 3.991136327e-05, 8.646627018e-12 MLMG: Timers: Solve = 0.323796894 Iter = 0.320042009 Bottom = 0.000282572 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.408511909 seconds MAC Proj :0.070789731 seconds Nodal Proj :0.346269895 seconds Reactions :0.714193389 seconds Misc :0.121993974 seconds Base State :0.001862721 seconds Time to advance time step: 1.662108082 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.509590679 Total Time: 13.23498373 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.04-7-gec648a92a89f) finalized