Initializing AMReX (26.08-31-g06424932b757)... MPI initialized with 2 MPI processes MPI initialized with thread support level 0 OMP initialized with 4 OMP threads AMReX (26.08-31-g06424932b757) 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.69440655 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.522042274 Iter = 0.504854802 Bottom = 0.002281483 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.675195759 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.135340083 Iter = 0.118136495 Bottom = 0.000508343 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.190663387 Iter = 0.173811761 Bottom = 0.00075611 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.001000557095, 2.930690201e-11 MLMG: Timers: Solve = 0.245500978 Iter = 0.22869801 Bottom = 0.000995965 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.691163508 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.8135733251, 3.30567653e-09 MLMG: Timers: Solve = 0.167673551 Iter = 0.15252328 Bottom = 0.004718545 <<< 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.501256215e-10 MLMG: Timers: Solve = 0.075114766 Iter = 0.062780631 Bottom = 0.001819821 <<< 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.001211438328, 8.730369042e-11 MLMG: Timers: Solve = 1.687290213 Iter = 1.670259857 Bottom = 0.007644948 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :2.493647393 seconds MAC Proj :0.337255933 seconds Nodal Proj :1.825025313 seconds Reactions :2.338985475 seconds Misc :0.4947694 seconds Base State :0.014865332 seconds Time to advance time step: 7.493324913 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.671376494 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.813606672, 3.305812427e-09 MLMG: Timers: Solve = 0.168143278 Iter = 0.15306902 Bottom = 0.004807573 <<< 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.8116332125 MLMG: No iterations needed MLMG: Timers: Solve = 0.014882455 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.002423048136, 9.850877948e-11 MLMG: Timers: Solve = 0.468207199 Iter = 0.451007028 Bottom = 0.002028776 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :2.538506665 seconds MAC Proj :0.279179388 seconds Nodal Proj :0.604655599 seconds Reactions :2.334686826 seconds Misc :0.467897904 seconds Base State :0.014738414 seconds Time to advance time step: 6.229711936 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.689853642 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.7995681735, 3.289753359e-09 MLMG: Timers: Solve = 0.165363859 Iter = 0.153251592 Bottom = 0.004847953 <<< 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.015305283 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.308628804 Iter = 0.292979034 Bottom = 0.001306218 Done calling nodal solver Timestep 2 ends with TIME = 0.0006059090797 DT = 0.0003173809465 Timing summary: Advection :2.552767724 seconds MAC Proj :0.275367997 seconds Nodal Proj :0.442611585 seconds Reactions :2.334668341 seconds Misc :0.471258083 seconds Base State :0.014866687 seconds Time to advance time step: 6.077047829 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.704316572 Time to regrid: 0.241657966 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.8497645864, 3.501608071e-09 MLMG: Timers: Solve = 0.166953527 Iter = 0.154303743 Bottom = 0.004806626 <<< 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.8510290105 MLMG: No iterations needed MLMG: Timers: Solve = 0.0151278 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.372031174 Iter = 0.356194931 Bottom = 0.001515466 Done calling nodal solver Timestep 3 ends with TIME = 0.0009550281208 DT = 0.0003491190411 Timing summary: Advection :2.578070831 seconds MAC Proj :0.277122094 seconds Nodal Proj :0.507712939 seconds Reactions :2.334788306 seconds Misc :0.469641301 seconds Base State :0.014860675 seconds Time to advance time step: 6.167720408 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.705672043 Total Time: 44.36413098 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.08-31-g06424932b757) finalized