MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.02-2-gf79f40d4d0cb) 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.515579199 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. 7 resid, resid/bnorm = 0.003170095384, 3.83960512e-11 MLMG: Timers: Solve = 0.124461373 Iter = 0.119044429 Bottom = 0.000119826 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.51138819 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.810359928, 1.125066149e-07 MLMG: Timers: Solve = 0.038074078 Iter = 0.033967349 Bottom = 3.8447e-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.01206794777, 6.500800205e-10 MLMG: Timers: Solve = 0.054844749 Iter = 0.050755694 Bottom = 5.1975e-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.0006594294682, 5.630252091e-11 MLMG: Timers: Solve = 0.071703857 Iter = 0.067603291 Bottom = 6.7106e-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.502826324 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 = 97092850.58 MLMG: Initial residual (resid0) = 97092850.58 MLMG: Final Iter. 35 resid, resid/bnorm = 0.007982999086, 8.222025658e-11 MLMG: Timers: Solve = 0.163275953 Iter = 0.16014409 Bottom = 0.000788315 <<< 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 = 97092842.37 MLMG: Initial residual (resid0) = 617.9399317 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001313373446, 1.352698525e-11 MLMG: Timers: Solve = 0.025327221 Iter = 0.023139959 Bottom = 0.000157191 <<< 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 = 140046226.7 MLMG: Initial residual (resid0) = 140046226.7 MLMG: Final Iter. 9 resid, resid/bnorm = 0.001075320179, 7.678323107e-12 MLMG: Timers: Solve = 0.154242611 Iter = 0.150148287 Bottom = 0.000133231 Done calling nodal solver Timestep 0 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.451504065 seconds MAC Proj :0.198200899 seconds Nodal Proj :0.178131943 seconds Reactions :0.710967607 seconds Misc :0.128460334 seconds Base State :0.001937486999 seconds Time to advance time step: 1.693529876 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.494165678 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 = 97077982.28 MLMG: Initial residual (resid0) = 97077982.28 MLMG: Final Iter. 35 resid, resid/bnorm = 0.007918104529, 8.156437065e-11 MLMG: Timers: Solve = 0.162047048 Iter = 0.159876531 Bottom = 0.000769294 <<< 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 = 97077982.1 MLMG: Initial residual (resid0) = 57.82689221 MLMG: Final Iter. 5 resid, resid/bnorm = 0.006387531757, 6.579794531e-11 MLMG: Timers: Solve = 0.025498635 Iter = 0.023302916 Bottom = 0.000141065 <<< 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.469 MLMG: Initial residual (resid0) = 8056355.469 MLMG: Final Iter. 6 resid, resid/bnorm = 3.902660683e-05, 4.844201199e-12 MLMG: Timers: Solve = 0.103374346 Iter = 0.099236128 Bottom = 9.605e-05 Done calling nodal solver Timestep 1 ends with TIME = 0.0005784714117 DT = 0.0005784714117 Timing summary: Advection :0.437639673 seconds MAC Proj :0.196285709 seconds Nodal Proj :0.126977894 seconds Reactions :0.705292849 seconds Misc :0.124378133 seconds Base State :0.001859933001 seconds Time to advance time step: 1.614324037 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.502802345 Call to estdt for level 0 gives dt_lev = 1.410020434 Minimum estdt over all levels = 1.410020434 Call to estdt at beginning of step 2 gives dt =1.410020434 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 = 92795847.39 MLMG: Initial residual (resid0) = 92795847.39 MLMG: Final Iter. 8 resid, resid/bnorm = 0.008088860661, 8.716834738e-11 MLMG: Timers: Solve = 0.038928916 Iter = 0.036724412 Bottom = 0.000197185 <<< 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 = 92795847.34 MLMG: Initial residual (resid0) = 122.5299544 MLMG: Final Iter. 3 resid, resid/bnorm = 0.009181700647, 9.894516738e-11 MLMG: Timers: Solve = 0.016049212 Iter = 0.013855001 Bottom = 7.9801e-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.668 MLMG: Initial residual (resid0) = 5904880.668 MLMG: Final Iter. 26 resid, resid/bnorm = 5.470076576e-05, 9.263653042e-12 MLMG: Timers: Solve = 0.435809273 Iter = 0.431547329 Bottom = 0.000363166 Done calling nodal solver Timestep 2 ends with TIME = 0.001214789965 DT = 0.0006363185529 Timing summary: Advection :0.438445263 seconds MAC Proj :0.063720847 seconds Nodal Proj :0.45981078 seconds Reactions :0.709152188 seconds Misc :0.125219246 seconds Base State :0.001890893 seconds Time to advance time step: 1.796713597 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.503407568 Call to estdt for level 0 gives dt_lev = 1.702946695 Minimum estdt over all levels = 1.702946695 Call to estdt at beginning of step 3 gives dt =1.702946695 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 = 91006747.56 MLMG: Initial residual (resid0) = 91006747.56 MLMG: Final Iter. 9 resid, resid/bnorm = 0.005484633148, 6.026622525e-11 MLMG: Timers: Solve = 0.043878838 Iter = 0.041686831 Bottom = 0.000217929 <<< 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 = 91006747.27 MLMG: Initial residual (resid0) = 2.885192442 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005470819771, 6.011444135e-11 MLMG: Timers: Solve = 0.020790837 Iter = 0.018581557 Bottom = 0.000106588 <<< 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.3 MLMG: Initial residual (resid0) = 4615830.3 MLMG: Final Iter. 18 resid, resid/bnorm = 3.25688161e-05, 7.055895469e-12 MLMG: Timers: Solve = 0.302750244 Iter = 0.298560604 Bottom = 0.000253335 Done calling nodal solver Timestep 3 ends with TIME = 0.001914740373 DT = 0.0006999504082 Timing summary: Advection :0.439116146 seconds MAC Proj :0.073447464 seconds Nodal Proj :0.326555962 seconds Reactions :0.712866052 seconds Misc :0.12486748 seconds Base State :0.001889378 seconds Time to advance time step: 1.677188746 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.503853835 Total Time: 13.3919016 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.02-2-gf79f40d4d0cb) finalized