Initializing AMReX (26.05-54-g012f8ddff80e)... MPI initialized with 2 MPI processes MPI initialized with thread support level 0 OMP initialized with 4 OMP threads AMReX (26.05-54-g012f8ddff80e) 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.62215713 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.498325936 Iter = 0.481789226 Bottom = 0.002324423 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.606306317 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.129031901 Iter = 0.112835658 Bottom = 0.000521681 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.181415423 Iter = 0.165316143 Bottom = 0.000779519 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.234339603 Iter = 0.217938355 Bottom = 0.001031801 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.605729982 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.166870404 Iter = 0.154267123 Bottom = 0.004898613 <<< 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.076050707 Iter = 0.063882936 Bottom = 0.00187579 <<< 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.604842216 Iter = 1.590126556 Bottom = 0.007802115 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :2.471606472 seconds MAC Proj :0.336495035 seconds Nodal Proj :1.736545829 seconds Reactions :2.135560045 seconds Misc :0.496679729 seconds Base State :0.014269729 seconds Time to advance time step: 7.179268619 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.593967678 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.167797566 Iter = 0.155080941 Bottom = 0.004979192 <<< 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.012700297 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.445362423 Iter = 0.430596054 Bottom = 0.002113128 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :2.473969204 seconds MAC Proj :0.272592787 seconds Nodal Proj :0.576663971 seconds Reactions :2.131437097 seconds Misc :0.466798551 seconds Base State :0.014338482 seconds Time to advance time step: 5.925473964 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.610892154 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.166682013 Iter = 0.154432236 Bottom = 0.00493507 <<< 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.013153429 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.30836464 Iter = 0.293714226 Bottom = 0.001370512 Done calling nodal solver Timestep 2 ends with TIME = 0.0006059090797 DT = 0.0003173809465 Timing summary: Advection :2.500574457 seconds MAC Proj :0.272725063 seconds Nodal Proj :0.439999816 seconds Reactions :2.132082726 seconds Misc :0.466857856 seconds Base State :0.014707608 seconds Time to advance time step: 5.812601617 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.620680374 Time to regrid: 0.235615469 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.167516673 Iter = 0.154431165 Bottom = 0.004968088 <<< 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.01264296 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.353309594 Iter = 0.338604814 Bottom = 0.001585183 Done calling nodal solver Timestep 3 ends with TIME = 0.0009550281208 DT = 0.0003491190411 Timing summary: Advection :2.510105332 seconds MAC Proj :0.273025908 seconds Nodal Proj :0.484733277 seconds Reactions :2.135033025 seconds Misc :0.469575685 seconds Base State :0.014452812 seconds Time to advance time step: 5.872804454 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.624600348 Total Time: 42.37666144 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.05-54-g012f8ddff80e) finalized