MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.02-12-g1d566b101fa3) 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.823731417 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 = 268163993.9 MLMG: Initial residual (resid0) = 268163993.9 MLMG: Final Iter. 5 resid, resid/bnorm = 0.02495675907, 9.306528704e-11 MLMG: Timers: Solve = 0.175667953 Iter = 0.169997812 Bottom = 0.00016948 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.835786419 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 = 100896162.8 MLMG: Initial residual (resid0) = 100896162.8 MLMG: Final Iter. 2 resid, resid/bnorm = 1.090743365, 1.081055349e-08 MLMG: Timers: Solve = 0.07424643 Iter = 0.068831607 Bottom = 6.4674e-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 = 8.768321343 Minimum estdt over all levels = 8.768321343 Call to estdt at end of istep_divu_iter = 1 gives dt = 8.768321343 Multiplying dt by init_shrink; dt = 0.8768321343 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 = 49401163.56 MLMG: Initial residual (resid0) = 49401163.56 MLMG: Final Iter. 2 resid, resid/bnorm = 0.01984390663, 4.016890535e-10 MLMG: Timers: Solve = 0.074834129 Iter = 0.069705578 Bottom = 6.3321e-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 = 8.962174815 Minimum estdt over all levels = 8.962174815 Call to estdt at end of istep_divu_iter = 2 gives dt = 8.962174815 Multiplying dt by init_shrink; dt = 0.8962174815 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 = 29021861.06 MLMG: Initial residual (resid0) = 29021861.06 MLMG: Final Iter. 3 resid, resid/bnorm = 0.001461790608, 5.036860335e-11 MLMG: Timers: Solve = 0.110379204 Iter = 0.105355122 Bottom = 9.2692e-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 = 9.055396973 Minimum estdt over all levels = 9.055396973 Call to estdt at end of istep_divu_iter = 3 gives dt = 9.055396973 Multiplying dt by init_shrink; dt = 0.9055396973 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.843682603 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 = 264654307 MLMG: Initial residual (resid0) = 264654307 MLMG: Final Iter. 18 resid, resid/bnorm = 0.2569767833, 9.709903692e-10 MLMG: Timers: Solve = 0.204342935 Iter = 0.195778241 Bottom = 0.000797353 <<< 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 = 264654307.4 MLMG: Initial residual (resid0) = 3790.869879 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1915032453, 7.235976892e-10 MLMG: Timers: Solve = 0.035696447 Iter = 0.029173664 Bottom = 9.5378e-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 = 14853631.31 MLMG: Initial residual (resid0) = 14853631.31 MLMG: Final Iter. 12 resid, resid/bnorm = 0.001000724396, 6.737237348e-11 MLMG: Timers: Solve = 0.425903374 Iter = 0.419752256 Bottom = 0.000372369 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :1.246406775 seconds MAC Proj :0.27229667 seconds Nodal Proj :0.480346008 seconds Reactions :1.080439481 seconds Misc :0.239676364 seconds Base State :0.005408241001 seconds Time to advance time step: 3.349909922 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.82594299 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 = 264654243.8 MLMG: Initial residual (resid0) = 264654243.8 MLMG: Final Iter. 18 resid, resid/bnorm = 0.2570389137, 9.712253615e-10 MLMG: Timers: Solve = 0.202571302 Iter = 0.195847765 Bottom = 0.00087257 <<< 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 = 264654243.7 MLMG: Initial residual (resid0) = 5.816711846 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1906003778, 7.201863653e-10 MLMG: Timers: Solve = 0.034274796 Iter = 0.027263995 Bottom = 0.000105368 <<< 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 = 19692713.43 MLMG: Initial residual (resid0) = 19692713.43 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001091005139, 5.540146324e-11 MLMG: Timers: Solve = 0.182395355 Iter = 0.176108379 Bottom = 0.000177516 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :1.197853931 seconds MAC Proj :0.271725287 seconds Nodal Proj :0.237795334 seconds Reactions :1.051937846 seconds Misc :0.242743702 seconds Base State :0.005354905999 seconds Time to advance time step: 3.033606666 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.838703094 Call to estdt for level 0 gives dt_lev = 3.005530282 Call to estdt for level 1 gives dt_lev = 0.6541869355 Minimum estdt over all levels = 0.6541869355 Call to estdt at beginning of step 2 gives dt =0.6541869355 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 = 263576960.2 MLMG: Initial residual (resid0) = 263576960.2 MLMG: Final Iter. 5 resid, resid/bnorm = 0.1716586081, 6.512656038e-10 MLMG: Timers: Solve = 0.061230781 Iter = 0.054780681 Bottom = 0.000230839 <<< 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 = 263576960 MLMG: Initial residual (resid0) = 2.087379048 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01388106373, 5.266417721e-11 MLMG: Timers: Solve = 0.018543286 Iter = 0.012169621 Bottom = 5.157e-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 = 14829531.16 MLMG: Initial residual (resid0) = 14829531.16 MLMG: Final Iter. 16 resid, resid/bnorm = 0.001309256069, 8.828708439e-11 MLMG: Timers: Solve = 0.548928378 Iter = 0.542698561 Bottom = 0.000536054 Done calling nodal solver Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974 Timing summary: Advection :1.197390854 seconds MAC Proj :0.111684761 seconds Nodal Proj :0.603950855 seconds Reactions :1.05589247 seconds Misc :0.238787329 seconds Base State :0.004957174 seconds Time to advance time step: 3.208186618 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.845580834 Time to regrid: 0.134559405 Call to estdt for level 0 gives dt_lev = 4.022405073 Call to estdt for level 1 gives dt_lev = 0.8210388141 Minimum estdt over all levels = 0.8210388141 Call to estdt at beginning of step 3 gives dt =0.8210388141 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 = 263106329.4 MLMG: Initial residual (resid0) = 263106329.4 MLMG: Final Iter. 5 resid, resid/bnorm = 0.1879902836, 7.145030834e-10 MLMG: Timers: Solve = 0.060637808 Iter = 0.054591318 Bottom = 0.000242422 <<< 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 = 263106329.2 MLMG: Initial residual (resid0) = 0.5969235472 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01681184224, 6.389752119e-11 MLMG: Timers: Solve = 0.0178635 Iter = 0.012014287 Bottom = 4.3628e-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 = 12313408.17 MLMG: Initial residual (resid0) = 12313408.17 MLMG: Final Iter. 8 resid, resid/bnorm = 0.001010330394, 8.205123878e-11 MLMG: Timers: Solve = 0.277052417 Iter = 0.271028261 Bottom = 0.000252612 Done calling nodal solver Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571 Timing summary: Advection :1.165542023 seconds MAC Proj :0.109011774 seconds Nodal Proj :0.332212394 seconds Reactions :1.055261713 seconds Misc :0.242049051 seconds Base State :0.005352681999 seconds Time to advance time step: 2.90452837 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.840565894 Total Time: 21.91705742 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 (22.02-12-g1d566b101fa3) finalized