MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.01-1-g30f2290cf7c3) 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 = 781250.000000 dr of input file data = 1562500.000000 maximum radius (cell-centered) of input model = 499218750.000000 setting r_cutoff to 239 radius at r_cutoff 187109375 Maximum HSE Error = 0.002311 (after putting initial model into base state arrays, and for density < base_cutoff_density) Writing plotfile wdconvect-amr_pltInitData after InitData inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.88905262 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Doing initial projection Calling nodal solver MLMG: Initial rhs = 214086148.6 MLMG: Initial residual (resid0) = 214086148.6 MLMG: Final Iter. 5 resid, resid/bnorm = 0.004865899452, 2.272869816e-11 MLMG: Timers: Solve = 0.177020406 Iter = 0.171514102 Bottom = 0.000146432 Done calling nodal solver Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.87154051 Call to firstdt for level 0 gives dt_lev = 0.005789736987 Multiplying dt_lev by init_shrink; dt_lev = 0.0005789736987 Call to firstdt for level 1 gives dt_lev = 0.002885750885 Multiplying dt_lev by init_shrink; dt_lev = 0.0002885750885 Minimum firstdt over all levels = 0.0002885750885 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 98262295.15 MLMG: Initial residual (resid0) = 98262295.15 MLMG: Final Iter. 2 resid, resid/bnorm = 0.3484706796, 3.546331571e-09 MLMG: Timers: Solve = 0.0734546 Iter = 0.067687322 Bottom = 6.9413e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 10.1029103 Minimum estdt over all levels = 10.1029103 Call to estdt at end of istep_divu_iter = 1 gives dt = 10.1029103 Multiplying dt by init_shrink; dt = 1.01029103 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 57185020 MLMG: Initial residual (resid0) = 57185020 MLMG: Final Iter. 2 resid, resid/bnorm = 0.01790044543, 3.130268282e-10 MLMG: Timers: Solve = 0.074313719 Iter = 0.068682975 Bottom = 6.1266e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 10.47579521 Minimum estdt over all levels = 10.47579521 Call to estdt at end of istep_divu_iter = 2 gives dt = 10.47579521 Multiplying dt by init_shrink; dt = 1.047579521 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 36948393.74 MLMG: Initial residual (resid0) = 36948393.74 MLMG: Final Iter. 3 resid, resid/bnorm = 0.002158799259, 5.842741837e-11 MLMG: Timers: Solve = 0.109240431 Iter = 0.10372768 Bottom = 9.087e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 10.6224833 Minimum estdt over all levels = 10.6224833 Call to estdt at end of istep_divu_iter = 3 gives dt = 10.6224833 Multiplying dt by init_shrink; dt = 1.06224833 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.862936778 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 0.0002885750885 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 300731026 MLMG: Initial residual (resid0) = 300731026 MLMG: Final Iter. 18 resid, resid/bnorm = 0.2997429669, 9.967144758e-10 MLMG: Timers: Solve = 0.190612789 Iter = 0.183556725 Bottom = 0.00075298 <<< 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 = 300731026.5 MLMG: Initial residual (resid0) = 10746.31623 MLMG: Final Iter. 3 resid, resid/bnorm = 0.02984782924, 9.925091397e-11 MLMG: Timers: Solve = 0.036583771 Iter = 0.032226791 Bottom = 0.000141497 <<< 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 = 16091202.17 MLMG: Initial residual (resid0) = 16091202.17 MLMG: Final Iter. 12 resid, resid/bnorm = 0.00100694557, 6.257739845e-11 MLMG: Timers: Solve = 0.405291377 Iter = 0.399957016 Bottom = 0.000372314 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :0.963986295 seconds MAC Proj :0.253985176 seconds Nodal Proj :0.454077894 seconds Reactions :1.168537033 seconds Misc :0.229367985 seconds Base State :0.004646494001 seconds Time to advance time step: 3.098422698 Writing plotfile 0 after all initialization inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.90884348 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 0.0002885750885 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 300732045.6 MLMG: Initial residual (resid0) = 300732045.6 MLMG: Final Iter. 18 resid, resid/bnorm = 0.299878031, 9.971602142e-10 MLMG: Timers: Solve = 0.187477491 Iter = 0.183092124 Bottom = 0.00075072 <<< 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 = 300732045.4 MLMG: Initial residual (resid0) = 17.79234161 MLMG: Final Iter. 3 resid, resid/bnorm = 0.03547971882, 1.179778456e-10 MLMG: Timers: Solve = 0.036091655 Iter = 0.031755709 Bottom = 0.000143485 <<< 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 = 27933872.39 MLMG: Initial residual (resid0) = 27933872.39 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001086700381, 3.89026042e-11 MLMG: Timers: Solve = 0.174628025 Iter = 0.169408432 Bottom = 0.000158156 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :0.920683617 seconds MAC Proj :0.248400955 seconds Nodal Proj :0.222230995 seconds Reactions :1.135564894 seconds Misc :0.220465228 seconds Base State :0.005308410999 seconds Time to advance time step: 2.779127705 Writing plotfile 1 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.853889007 Call to estdt for level 0 gives dt_lev = 2.502703181 Call to estdt for level 1 gives dt_lev = 0.5259215388 Minimum estdt over all levels = 0.5259215388 Call to estdt at beginning of step 2 gives dt =0.5259215388 dt_growth factor limits the new dt = 0.0003174325974 Timestep 2 starts with TIME = 0.0002885750885 DT = 0.0003174325974 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 294777753.7 MLMG: Initial residual (resid0) = 294777753.7 MLMG: Final Iter. 5 resid, resid/bnorm = 0.277498275, 9.41381334e-10 MLMG: Timers: Solve = 0.056150547 Iter = 0.05178629 Bottom = 0.000217028 <<< 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 = 294777753.5 MLMG: Initial residual (resid0) = 0.4180259939 MLMG: Final Iter. 1 resid, resid/bnorm = 0.02582601718, 8.761182577e-11 MLMG: Timers: Solve = 0.015704612 Iter = 0.011340821 Bottom = 4.2783e-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 = 21482439.4 MLMG: Initial residual (resid0) = 21482439.4 MLMG: Final Iter. 13 resid, resid/bnorm = 0.001904040546, 8.863241789e-11 MLMG: Timers: Solve = 0.429689543 Iter = 0.424269659 Bottom = 0.000384305 Done calling nodal solver Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974 Timing summary: Advection :0.95300028 seconds MAC Proj :0.0965213 seconds Nodal Proj :0.478357167 seconds Reactions :1.165713591 seconds Misc :0.227408183 seconds Base State :0.004750892001 seconds Time to advance time step: 2.921415298 Writing plotfile 2 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.874420955 Time to regrid: 0.133815745 Call to estdt for level 0 gives dt_lev = 3.210026718 Call to estdt for level 1 gives dt_lev = 0.6510211778 Minimum estdt over all levels = 0.6510211778 Call to estdt at beginning of step 3 gives dt =0.6510211778 dt_growth factor limits the new dt = 0.0003491758571 Timestep 3 starts with TIME = 0.0006060076859 DT = 0.0003491758571 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 292704981.9 MLMG: Initial residual (resid0) = 292704981.9 MLMG: Final Iter. 5 resid, resid/bnorm = 0.275859639, 9.424494153e-10 MLMG: Timers: Solve = 0.056806667 Iter = 0.052491362 Bottom = 0.000224007 <<< 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 = 292704981.7 MLMG: Initial residual (resid0) = 0.4147914393 MLMG: Final Iter. 1 resid, resid/bnorm = 0.02585147805, 8.831922815e-11 MLMG: Timers: Solve = 0.015549965 Iter = 0.011238834 Bottom = 4.6058e-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 = 16800418.13 MLMG: Initial residual (resid0) = 16800418.13 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001508717412, 8.980237278e-11 MLMG: Timers: Solve = 0.172271928 Iter = 0.166951808 Bottom = 0.000159358 Done calling nodal solver Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571 Timing summary: Advection :0.954387256 seconds MAC Proj :0.097196631 seconds Nodal Proj :0.220826587 seconds Reactions :1.179538575 seconds Misc :0.228144276 seconds Base State :0.004920937999 seconds Time to advance time step: 2.680523176 Writing plotfile 3 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.876751218 Total Time: 21.41684875 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::amr.check_file(nvals = 1) :: [wdconvect-amr_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] AMReX (23.01-1-g30f2290cf7c3) finalized