MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.04-9-g84a1f7d0715b) 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.534023172 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.003142089117, 3.805684057e-11 MLMG: Timers: Solve = 0.105329207 Iter = 0.100119562 Bottom = 0.000119821 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.516065358 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.803742786, 1.123112339e-07 MLMG: Timers: Solve = 0.036712673 Iter = 0.032706344 Bottom = 3.8583e-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.01096891984, 5.908772374e-10 MLMG: Timers: Solve = 0.052570718 Iter = 0.048605453 Bottom = 5.0691e-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.0002620539162, 2.23743354e-11 MLMG: Timers: Solve = 0.069554082 Iter = 0.065825237 Bottom = 6.3963e-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.549892745 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 = 97092795.74 MLMG: Initial residual (resid0) = 97092795.74 MLMG: Final Iter. 32 resid, resid/bnorm = 0.00855819881, 8.814452962e-11 MLMG: Timers: Solve = 0.158832852 Iter = 0.155926034 Bottom = 0.000776798 <<< 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 = 97092787.63 MLMG: Initial residual (resid0) = 479.9637084 MLMG: Final Iter. 8 resid, resid/bnorm = 0.007515966892, 7.741014627e-11 MLMG: Timers: Solve = 0.04008335 Iter = 0.038052652 Bottom = 0.000211027 <<< 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 = 140046431.8 MLMG: Initial residual (resid0) = 140046431.8 MLMG: Final Iter. 11 resid, resid/bnorm = 0.0003913390683, 2.794352296e-12 MLMG: Timers: Solve = 0.182403632 Iter = 0.178849997 Bottom = 0.000161681 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.428975194 seconds MAC Proj :0.208062167 seconds Nodal Proj :0.20528371 seconds Reactions :0.732708533 seconds Misc :0.12410788 seconds Base State :0.002123484 seconds Time to advance time step: 1.76593129 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.505177183 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 = 97077927.64 MLMG: Initial residual (resid0) = 97077927.64 MLMG: Final Iter. 32 resid, resid/bnorm = 0.008489765227, 8.745309498e-11 MLMG: Timers: Solve = 0.155405567 Iter = 0.153668268 Bottom = 0.000779289 <<< 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 = 97077927.47 MLMG: Initial residual (resid0) = 28.7998304 MLMG: Final Iter. 8 resid, resid/bnorm = 0.007158502936, 7.373975859e-11 MLMG: Timers: Solve = 0.03915951 Iter = 0.037411824 Bottom = 0.000210641 <<< 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.336 MLMG: Initial residual (resid0) = 8056355.336 MLMG: Final Iter. 7 resid, resid/bnorm = 5.602696911e-05, 6.954381574e-12 MLMG: Timers: Solve = 0.116968365 Iter = 0.113324042 Bottom = 0.000115426 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.405834427 seconds MAC Proj :0.20338074 seconds Nodal Proj :0.140822019 seconds Reactions :0.733818678 seconds Misc :0.123838176 seconds Base State :0.001882298 seconds Time to advance time step: 1.633196335 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.513161197 Call to estdt for level 0 gives dt_lev = 1.410027137 Minimum estdt over all levels = 1.410027137 Call to estdt at beginning of step 2 gives dt =1.410027137 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 = 92795813.47 MLMG: Initial residual (resid0) = 92795813.47 MLMG: Final Iter. 9 resid, resid/bnorm = 0.004723660648, 5.090381205e-11 MLMG: Timers: Solve = 0.044421059 Iter = 0.042684006 Bottom = 0.000216554 <<< 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 = 92795813.29 MLMG: Initial residual (resid0) = 3.587058391 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005256935954, 5.665057256e-11 MLMG: Timers: Solve = 0.020510619 Iter = 0.018764851 Bottom = 0.00010535 <<< 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.492 MLMG: Initial residual (resid0) = 5904880.492 MLMG: Final Iter. 19 resid, resid/bnorm = 5.400460213e-05, 9.145757007e-12 MLMG: Timers: Solve = 0.3299458 Iter = 0.326276311 Bottom = 0.000279569 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.411303679 seconds MAC Proj :0.07273143 seconds Nodal Proj :0.353103103 seconds Reactions :0.736315063 seconds Misc :0.122751225 seconds Base State :0.001943102 seconds Time to advance time step: 1.696524079 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.520292232 Call to estdt for level 0 gives dt_lev = 1.702977295 Minimum estdt over all levels = 1.702977295 Call to estdt at beginning of step 3 gives dt =1.702977295 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 = 91006767.28 MLMG: Initial residual (resid0) = 91006767.28 MLMG: Final Iter. 9 resid, resid/bnorm = 0.004906795919, 5.391682471e-11 MLMG: Timers: Solve = 0.043982302 Iter = 0.042255786 Bottom = 0.000221318 <<< 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 = 91006767.02 MLMG: Initial residual (resid0) = 3.010248802 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005607455969, 6.161581334e-11 MLMG: Timers: Solve = 0.020603722 Iter = 0.018840674 Bottom = 0.000105331 <<< 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.343 MLMG: Initial residual (resid0) = 4615830.343 MLMG: Final Iter. 21 resid, resid/bnorm = 3.917491995e-05, 8.487079687e-12 MLMG: Timers: Solve = 0.341873416 Iter = 0.338183606 Bottom = 0.000303813 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.409726686 seconds MAC Proj :0.072362894 seconds Nodal Proj :0.36485522 seconds Reactions :0.737311208 seconds Misc :0.123317251 seconds Base State :0.001951286 seconds Time to advance time step: 1.707906746 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.518289243 Total Time: 13.05116298 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-9-g84a1f7d0715b) finalized