MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.03-12-g806f2373b68e) 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.536146324 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. 9 resid, resid/bnorm = 0.007611370645, 9.218857528e-11 MLMG: Timers: Solve = 0.152560694 Iter = 0.147277115 Bottom = 0.000158542 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.498141324 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.62 MLMG: Initial residual (resid0) = 33867874.62 MLMG: Final Iter. 2 resid, resid/bnorm = 3.819495839, 1.127763664e-07 MLMG: Timers: Solve = 0.036794592 Iter = 0.032878094 Bottom = 4.1562e-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.01410771406, 7.599587945e-10 MLMG: Timers: Solve = 0.054383516 Iter = 0.050405911 Bottom = 5.6408e-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. 5 resid, resid/bnorm = 0.0005255709402, 4.487359192e-11 MLMG: Timers: Solve = 0.087610485 Iter = 0.08367047 Bottom = 8.2758e-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.502050003 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 = 97092897.53 MLMG: Initial residual (resid0) = 97092897.53 MLMG: Final Iter. 37 resid, resid/bnorm = 0.007875282317, 8.111079716e-11 MLMG: Timers: Solve = 0.174807392 Iter = 0.171656895 Bottom = 0.000916674 <<< 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 = 97092889.64 MLMG: Initial residual (resid0) = 4269.78354 MLMG: Final Iter. 13 resid, resid/bnorm = 0.009295057505, 9.573365815e-11 MLMG: Timers: Solve = 0.062732296 Iter = 0.060520782 Bottom = 0.00033559 <<< 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 = 140046873.2 MLMG: Initial residual (resid0) = 140046873.2 MLMG: Final Iter. 6 resid, resid/bnorm = 0.001332590356, 9.515316735e-12 MLMG: Timers: Solve = 0.106109907 Iter = 0.102168998 Bottom = 0.00010501 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.454536233 seconds MAC Proj :0.247314396 seconds Nodal Proj :0.130283755 seconds Reactions :0.728184928 seconds Misc :0.126830342 seconds Base State :0.002002105 seconds Time to advance time step: 1.709891569 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.503971498 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 = 97078028.98 MLMG: Initial residual (resid0) = 97078028.98 MLMG: Final Iter. 37 resid, resid/bnorm = 0.007793787867, 8.028374648e-11 MLMG: Timers: Solve = 0.174924547 Iter = 0.172702515 Bottom = 0.000909999 <<< 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 = 97078028.84 MLMG: Initial residual (resid0) = 25.54434763 MLMG: Final Iter. 13 resid, resid/bnorm = 0.009465780109, 9.750692533e-11 MLMG: Timers: Solve = 0.062680698 Iter = 0.060466976 Bottom = 0.000333429 <<< 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.067 MLMG: Initial residual (resid0) = 8056355.067 MLMG: Final Iter. 5 resid, resid/bnorm = 1.725084076e-05, 2.141271159e-12 MLMG: Timers: Solve = 0.087398053 Iter = 0.083412406 Bottom = 8.674e-05 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.437277589 seconds MAC Proj :0.246516441 seconds Nodal Proj :0.111940958 seconds Reactions :0.729210805 seconds Misc :0.125574019 seconds Base State :0.001984848 seconds Time to advance time step: 1.676126276 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.508206776 Call to estdt for level 0 gives dt_lev = 1.410016307 Minimum estdt over all levels = 1.410016307 Call to estdt at beginning of step 2 gives dt =1.410016307 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 = 92795874.88 MLMG: Initial residual (resid0) = 92795874.88 MLMG: Final Iter. 9 resid, resid/bnorm = 0.006033092737, 6.50146652e-11 MLMG: Timers: Solve = 0.04494571 Iter = 0.042744553 Bottom = 0.000221944 <<< 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 = 92795874.52 MLMG: Initial residual (resid0) = 3.318257083 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005111910403, 5.508769036e-11 MLMG: Timers: Solve = 0.020970165 Iter = 0.018742431 Bottom = 0.00010333 <<< 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.223 MLMG: Initial residual (resid0) = 5904880.223 MLMG: Final Iter. 28 resid, resid/bnorm = 4.179868847e-05, 7.078668304e-12 MLMG: Timers: Solve = 0.472851968 Iter = 0.468847769 Bottom = 0.000421161 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.438193981 seconds MAC Proj :0.074824083 seconds Nodal Proj :0.496818455 seconds Reactions :0.733380299 seconds Misc :0.12571637 seconds Base State :0.001937869 seconds Time to advance time step: 1.86928233 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.529002002 Call to estdt for level 0 gives dt_lev = 1.702929365 Minimum estdt over all levels = 1.702929365 Call to estdt at beginning of step 3 gives dt =1.702929365 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 = 91006732.6 MLMG: Initial residual (resid0) = 91006732.6 MLMG: Final Iter. 9 resid, resid/bnorm = 0.004234939814, 4.653435732e-11 MLMG: Timers: Solve = 0.044051938 Iter = 0.04182157 Bottom = 0.000252047 <<< 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 = 91006732.33 MLMG: Initial residual (resid0) = 121.6655032 MLMG: Final Iter. 4 resid, resid/bnorm = 0.004215355963, 4.631916622e-11 MLMG: Timers: Solve = 0.022452552 Iter = 0.020261571 Bottom = 0.000103064 <<< 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.337 MLMG: Initial residual (resid0) = 4615830.337 MLMG: Final Iter. 23 resid, resid/bnorm = 3.699073568e-05, 8.013885472e-12 MLMG: Timers: Solve = 0.375289464 Iter = 0.371300444 Bottom = 0.000325092 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.440681058 seconds MAC Proj :0.076421757 seconds Nodal Proj :0.398912767 seconds Reactions :0.73009942 seconds Misc :0.123864312 seconds Base State :0.001848358 seconds Time to advance time step: 1.770313743 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.508631973 Total Time: 13.30371813 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.03-12-g806f2373b68e) finalized