Initializing AMReX (26.10-11-ga39ac3a1ed11)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.10-11-ga39ac3a1ed11) 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.1280 model file = kepler_new_6.25e8.hybrid.hse.1280 reading initial model 1280 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 = 390625.000000 maximum radius (cell-centered) of input model = 499804687.500000 setting r_cutoff to 239 radius at r_cutoff 187109375 Maximum HSE Error = 0.003696 (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: 1.327275801 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 = 204587874.9 MLMG: Initial residual (resid0) = 204587874.9 MLMG: Final Iter. 9 resid, resid/bnorm = 0.01924523525, 9.40683081e-11 MLMG: Timers: Solve = 0.322342695 Iter = 0.312027058 Bottom = 0.000172301 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: 1.319624387 Call to firstdt for level 0 gives dt_lev = 0.01172063987 Multiplying dt_lev by init_shrink; dt_lev = 0.001172063987 Call to firstdt for level 1 gives dt_lev = 0.00578879243 Multiplying dt_lev by init_shrink; dt_lev = 0.000578879243 Call to firstdt for level 2 gives dt_lev = 0.002885281332 Multiplying dt_lev by init_shrink; dt_lev = 0.0002885281332 Minimum firstdt over all levels = 0.0002885281332 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 93695655.05 MLMG: Initial residual (resid0) = 93695655.05 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005666818703, 6.048112583e-11 MLMG: Timers: Solve = 0.151625107 Iter = 0.141658845 Bottom = 7.4019e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 66.92844876 Call to estdt for level 1 gives dt_lev = 30.20205688 Call to estdt for level 2 gives dt_lev = 9.589324705 Minimum estdt over all levels = 9.589324705 Call to estdt at end of istep_divu_iter = 1 gives dt = 9.589324705 Multiplying dt by init_shrink; dt = 0.9589324705 Ignoring this new dt since it's larger than the previous dt = 0.0002885281332 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: 1.318532665 Doing initial pressure iteration #1  Timestep 0 starts with TIME = 0 DT = 0.0002885281332 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 280128357.4 MLMG: Initial residual (resid0) = 280128357.4 MLMG: Final Iter. 5 resid, resid/bnorm = 0.7851953404, 2.802984131e-09 MLMG: Timers: Solve = 0.094860598 Iter = 0.085861987 Bottom = 0.000152333 <<< 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 = 280128358.1 MLMG: Initial residual (resid0) = 4804.307766 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1073008775, 3.830418249e-10 MLMG: Timers: Solve = 0.042952037 Iter = 0.03557614 Bottom = 6.0263e-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 = 13780896.98 MLMG: Initial residual (resid0) = 13780896.98 MLMG: Final Iter. 30 resid, resid/bnorm = 0.001205376582, 8.746720796e-11 MLMG: Timers: Solve = 1.048331136 Iter = 1.038005114 Bottom = 0.00060542 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :1.89641468 seconds MAC Proj :0.191276183 seconds Nodal Proj :1.154298291 seconds Reactions :1.872425254 seconds Misc :0.367626222 seconds Base State :0.009872832 seconds Time to advance time step: 5.485110682 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: 1.344050588 Beginning main evolution  Timestep 1 starts with TIME = 0 DT = 0.0002885281332 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 280128393.8 MLMG: Initial residual (resid0) = 280128393.8 MLMG: Final Iter. 5 resid, resid/bnorm = 0.7852297546, 2.803106619e-09 MLMG: Timers: Solve = 0.094084188 Iter = 0.086610075 Bottom = 0.000153184 <<< 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 = 280128393.4 MLMG: Initial residual (resid0) = 0.8028606176 MLMG: No iterations needed MLMG: Timers: Solve = 0.007464678 Iter = 0 Bottom = 0 <<< 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 = 50383856.33 MLMG: Initial residual (resid0) = 50383856.33 MLMG: Final Iter. 7 resid, resid/bnorm = 0.004018843421, 7.976450621e-11 MLMG: Timers: Solve = 0.256743466 Iter = 0.246690383 Bottom = 0.000153109 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :1.929557917 seconds MAC Proj :0.154794626 seconds Nodal Proj :0.362381148 seconds Reactions :1.899298516 seconds Misc :0.366930668 seconds Base State :0.009912148 seconds Time to advance time step: 4.717005654 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: 1.358508819 Call to estdt for level 0 gives dt_lev = 5.556499851 Call to estdt for level 1 gives dt_lev = 1.360999872 Call to estdt for level 2 gives dt_lev = 0.3374737097 Minimum estdt over all levels = 0.3374737097 Call to estdt at beginning of step 2 gives dt =0.3374737097 dt_growth factor limits the new dt = 0.0003173809465  Timestep 2 starts with TIME = 0.0002885281332 DT = 0.0003173809465 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 255394708.2 MLMG: Initial residual (resid0) = 255394708.2 MLMG: Final Iter. 5 resid, resid/bnorm = 0.6347444654, 2.485346974e-09 MLMG: Timers: Solve = 0.094351571 Iter = 0.086667347 Bottom = 0.000159346 <<< 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 = 255394707.9 MLMG: Initial residual (resid0) = 0.6680596396 MLMG: No iterations needed MLMG: Timers: Solve = 0.007579159 Iter = 0 Bottom = 0 <<< 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 = 34141443.43 MLMG: Initial residual (resid0) = 34141443.43 MLMG: Final Iter. 4 resid, resid/bnorm = 0.002735068622, 8.010992939e-11 MLMG: Timers: Solve = 0.153055163 Iter = 0.143038238 Bottom = 0.00010341 Done calling nodal solver Timestep 2 ends with TIME = 0.0006059090797 DT = 0.0003173809465 Timing summary: Advection :1.944908836 seconds MAC Proj :0.15521607 seconds Nodal Proj :0.257084481 seconds Reactions :1.886763072 seconds Misc :0.365414467 seconds Base State :0.010457963 seconds Time to advance time step: 4.610681805 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: 1.363374768 Time to regrid: 0.178396011 Call to estdt for level 0 gives dt_lev = 8.534960391 Call to estdt for level 1 gives dt_lev = 2.077552712 Call to estdt for level 2 gives dt_lev = 0.4347551008 Minimum estdt over all levels = 0.4347551008 Call to estdt at beginning of step 3 gives dt =0.4347551008 dt_growth factor limits the new dt = 0.0003491190411  Timestep 3 starts with TIME = 0.0006059090797 DT = 0.0003491190411 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 244508152.9 MLMG: Initial residual (resid0) = 244508152.9 MLMG: Final Iter. 5 resid, resid/bnorm = 0.6593789328, 2.696756427e-09 MLMG: Timers: Solve = 0.093040894 Iter = 0.085641737 Bottom = 0.000157183 <<< 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 = 244508152.5 MLMG: Initial residual (resid0) = 0.667964071 MLMG: No iterations needed MLMG: Timers: Solve = 0.007636694 Iter = 0 Bottom = 0 <<< 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 = 24597178.87 MLMG: Initial residual (resid0) = 24597178.87 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001891462023, 7.689751874e-11 MLMG: Timers: Solve = 0.192419897 Iter = 0.182393212 Bottom = 0.000105956 Done calling nodal solver Timestep 3 ends with TIME = 0.0009550281208 DT = 0.0003491190411 Timing summary: Advection :1.966098179 seconds MAC Proj :0.154249658 seconds Nodal Proj :0.297352439 seconds Reactions :1.901089999 seconds Misc :0.367045345 seconds Base State :0.009869556 seconds Time to advance time step: 4.686207372 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: 1.347798792 Total Time: 32.35197098 Unused ParmParse Variables: [TOP]::amr.check_file(nvals = 1) :: [wdconvect_amr_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] [TOP]::amr.plot_file(nvals = 1) :: [wdconvect_amr_plt] AMReX (26.10-11-ga39ac3a1ed11) finalized