Initializing AMReX (26.06-17-g6a6dcd903f99)... MPI initialized with 2 MPI processes MPI initialized with thread support level 0 OMP initialized with 4 OMP threads AMReX (26.06-17-g6a6dcd903f99) 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_omp_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.627079882 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.01924523572, 9.406831038e-11 MLMG: Timers: Solve = 0.502682805 Iter = 0.485628153 Bottom = 0.002402639 Done calling nodal solver Writing plotfile wdconvect_omp_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.609922431 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. 2 resid, resid/bnorm = 18.06786, 1.928356229e-07 MLMG: Timers: Solve = 0.130572845 Iter = 0.113903508 Bottom = 0.00052651 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 66.92844878 Call to estdt for level 1 gives dt_lev = 30.20205689 Call to estdt for level 2 gives dt_lev = 9.58932479 Minimum estdt over all levels = 9.58932479 Call to estdt at end of istep_divu_iter = 1 gives dt = 9.58932479 Multiplying dt by init_shrink; dt = 0.958932479 Ignoring this new dt since it's larger than the previous dt = 0.0002885281332 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 50386913.72 MLMG: Initial residual (resid0) = 50386913.72 MLMG: Final Iter. 3 resid, resid/bnorm = 0.05789050367, 1.14891942e-09 MLMG: Timers: Solve = 0.183467329 Iter = 0.166920196 Bottom = 0.000783034 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 66.92844878 Call to estdt for level 1 gives dt_lev = 30.20205689 Call to estdt for level 2 gives dt_lev = 9.794370248 Minimum estdt over all levels = 9.794370248 Call to estdt at end of istep_divu_iter = 2 gives dt = 9.794370248 Multiplying dt by init_shrink; dt = 0.9794370248 Ignoring this new dt since it's larger than the previous dt = 0.0002885281332 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 34140664.01 MLMG: Initial residual (resid0) = 34140664.01 MLMG: Final Iter. 4 resid, resid/bnorm = 0.001000558957, 2.930695657e-11 MLMG: Timers: Solve = 0.236367214 Iter = 0.220020238 Bottom = 0.001014172 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 66.92844878 Call to estdt for level 1 gives dt_lev = 30.20205689 Call to estdt for level 2 gives dt_lev = 9.920131412 Minimum estdt over all levels = 9.920131412 Call to estdt at end of istep_divu_iter = 3 gives dt = 9.920131412 Multiplying dt by init_shrink; dt = 0.9920131412 Ignoring this new dt since it's larger than the previous dt = 0.0002885281332 Writing plotfile wdconvect_omp_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.606073609 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 = 246114015.6 MLMG: Initial residual (resid0) = 246114015.6 MLMG: Final Iter. 5 resid, resid/bnorm = 0.8135733195, 3.305676507e-09 MLMG: Timers: Solve = 0.167174739 Iter = 0.153748567 Bottom = 0.004937495 <<< 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 = 246114016.4 MLMG: Initial residual (resid0) = 4972.291553 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1107822246, 4.501256217e-10 MLMG: Timers: Solve = 0.075491412 Iter = 0.063399137 Bottom = 0.001848144 <<< 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 = 13876141.11 MLMG: Initial residual (resid0) = 13876141.11 MLMG: Final Iter. 30 resid, resid/bnorm = 0.001211420866, 8.730243198e-11 MLMG: Timers: Solve = 1.620208024 Iter = 1.605069888 Bottom = 0.007964834 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :2.459081995 seconds MAC Proj :0.334018632 seconds Nodal Proj :1.75339816 seconds Reactions :2.144580968 seconds Misc :0.493011712 seconds Base State :0.014398052 seconds Time to advance time step: 7.188772138 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.600045484 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 = 246113985.6 MLMG: Initial residual (resid0) = 246113985.6 MLMG: Final Iter. 5 resid, resid/bnorm = 0.8136066943, 3.305812518e-09 MLMG: Timers: Solve = 0.166815348 Iter = 0.15420061 Bottom = 0.004893031 <<< 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 = 246113985.2 MLMG: Initial residual (resid0) = 0.8116332367 MLMG: No iterations needed MLMG: Timers: Solve = 0.011971981 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 = 24597281.06 MLMG: Initial residual (resid0) = 24597281.06 MLMG: Final Iter. 8 resid, resid/bnorm = 0.002423047554, 9.850875581e-11 MLMG: Timers: Solve = 0.450002277 Iter = 0.433397324 Bottom = 0.002118766 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :2.472118606 seconds MAC Proj :0.269735089 seconds Nodal Proj :0.584835181 seconds Reactions :2.146579805 seconds Misc :0.469191092 seconds Base State :0.014650951 seconds Time to advance time step: 5.946415015 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.619389498 Call to estdt for level 0 gives dt_lev = 10.32095586 Call to estdt for level 1 gives dt_lev = 2.510028204 Call to estdt for level 2 gives dt_lev = 0.49071817 Minimum estdt over all levels = 0.49071817 Call to estdt at beginning of step 2 gives dt =0.49071817 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 = 243048060.5 MLMG: Initial residual (resid0) = 243048060.5 MLMG: Final Iter. 5 resid, resid/bnorm = 0.7995681865, 3.289753413e-09 MLMG: Timers: Solve = 0.166522197 Iter = 0.154278218 Bottom = 0.004947134 <<< 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 = 243048060.2 MLMG: Initial residual (resid0) = 0.7989585372 MLMG: No iterations needed MLMG: Timers: Solve = 0.011797975 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 = 18863837.5 MLMG: Initial residual (resid0) = 18863837.5 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001626458589, 8.622098178e-11 MLMG: Timers: Solve = 0.302976001 Iter = 0.28782925 Bottom = 0.001363787 Done calling nodal solver Timestep 2 ends with TIME = 0.0006059090797 DT = 0.0003173809465 Timing summary: Advection :2.477260317 seconds MAC Proj :0.269712958 seconds Nodal Proj :0.441964053 seconds Reactions :2.143168746 seconds Misc :0.467245983 seconds Base State :0.014580291 seconds Time to advance time step: 5.799937696 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.635540278 Time to regrid: 0.236664857 Call to estdt for level 0 gives dt_lev = 13.17000951 Call to estdt for level 1 gives dt_lev = 3.169691461 Call to estdt for level 2 gives dt_lev = 0.5959748044 Minimum estdt over all levels = 0.5959748044 Call to estdt at beginning of step 3 gives dt =0.5959748044 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 = 242678383.5 MLMG: Initial residual (resid0) = 242678383.5 MLMG: Final Iter. 5 resid, resid/bnorm = 0.8497645995, 3.501608125e-09 MLMG: Timers: Solve = 0.167040904 Iter = 0.154488051 Bottom = 0.004924631 <<< 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 = 242678383.1 MLMG: Initial residual (resid0) = 0.851029031 MLMG: No iterations needed MLMG: Timers: Solve = 0.012498613 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 = 15333980.08 MLMG: Initial residual (resid0) = 15333980.08 MLMG: Final Iter. 6 resid, resid/bnorm = 0.001013564797, 6.609926396e-11 MLMG: Timers: Solve = 0.357827938 Iter = 0.342726812 Bottom = 0.00162498 Done calling nodal solver Timestep 3 ends with TIME = 0.0009550281208 DT = 0.0003491190411 Timing summary: Advection :2.477216444 seconds MAC Proj :0.271351297 seconds Nodal Proj :0.490983058 seconds Reactions :2.144013249 seconds Misc :0.466706582 seconds Base State :0.014709691 seconds Time to advance time step: 5.851475018 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.634313267 Total Time: 42.43945616 Unused ParmParse Variables: [TOP]::amr.check_file(nvals = 1) :: [wdconvect_omp_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] [TOP]::amr.plot_file(nvals = 1) :: [wdconvect_omp_plt] AMReX (26.06-17-g6a6dcd903f99) finalized