MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.04-28-g56661522b4cd) 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.81723721 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. 6 resid, resid/bnorm = 0.01291736215, 4.816963667e-11 MLMG: Timers: Solve = 0.194441015 Iter = 0.188641416 Bottom = 0.000180214 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.806846386 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.089379756, 1.079703852e-08 MLMG: Timers: Solve = 0.068988055 Iter = 0.063954865 Bottom = 6.3476e-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.768321628 Minimum estdt over all levels = 8.768321628 Call to estdt at end of istep_divu_iter = 1 gives dt = 8.768321628 Multiplying dt by init_shrink; dt = 0.8768321628 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.01966028614, 3.97972127e-10 MLMG: Timers: Solve = 0.069486908 Iter = 0.064406131 Bottom = 6.3697e-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.962175113 Minimum estdt over all levels = 8.962175113 Call to estdt at end of istep_divu_iter = 2 gives dt = 8.962175113 Multiplying dt by init_shrink; dt = 0.8962175113 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.001470201781, 5.065842531e-11 MLMG: Timers: Solve = 0.102749012 Iter = 0.097743767 Bottom = 9.2039e-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.055397278 Minimum estdt over all levels = 9.055397278 Call to estdt at end of istep_divu_iter = 3 gives dt = 9.055397278 Multiplying dt by init_shrink; dt = 0.9055397278 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.804310068 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 = 264654309.1 MLMG: Initial residual (resid0) = 264654309.1 MLMG: Final Iter. 21 resid, resid/bnorm = 0.2200590931, 8.31496354e-10 MLMG: Timers: Solve = 0.219125846 Iter = 0.213082825 Bottom = 0.000972454 <<< 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 = 264654309.5 MLMG: Initial residual (resid0) = 7827.673425 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1555062775, 5.875826385e-10 MLMG: Timers: Solve = 0.025713726 Iter = 0.021139038 Bottom = 9.1077e-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 = 14853790.41 MLMG: Initial residual (resid0) = 14853790.41 MLMG: Final Iter. 12 resid, resid/bnorm = 0.0009969010987, 6.711425642e-11 MLMG: Timers: Solve = 0.401840781 Iter = 0.39673401 Bottom = 0.000406135 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :0.978573836 seconds MAC Proj :0.270129168 seconds Nodal Proj :0.457072059 seconds Reactions :1.046855049 seconds Misc :0.223208745 seconds Base State :0.004714247 seconds Time to advance time step: 3.00018951 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.815710312 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 = 264654246 MLMG: Initial residual (resid0) = 264654246 MLMG: Final Iter. 21 resid, resid/bnorm = 0.2202278152, 8.321340713e-10 MLMG: Timers: Solve = 0.2142853 Iter = 0.210110022 Bottom = 0.00106785 <<< 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 = 264654245.8 MLMG: Initial residual (resid0) = 731.8634252 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1751463152, 6.617929544e-10 MLMG: Timers: Solve = 0.025042567 Iter = 0.020876442 Bottom = 9.5308e-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 = 19692713.47 MLMG: Initial residual (resid0) = 19692713.47 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001111145892, 5.642421468e-11 MLMG: Timers: Solve = 0.166400238 Iter = 0.161425512 Bottom = 0.000176785 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :0.945316649 seconds MAC Proj :0.262547112 seconds Nodal Proj :0.215492727 seconds Reactions :1.045666429 seconds Misc :0.228411405 seconds Base State :0.004832062 seconds Time to advance time step: 2.721766868 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.813786868 Call to estdt for level 0 gives dt_lev = 3.00553387 Call to estdt for level 1 gives dt_lev = 0.6541868959 Minimum estdt over all levels = 0.6541868959 Call to estdt at beginning of step 2 gives dt =0.6541868959 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 = 263576961.7 MLMG: Initial residual (resid0) = 263576961.7 MLMG: Final Iter. 5 resid, resid/bnorm = 0.1873782054, 7.109050963e-10 MLMG: Timers: Solve = 0.05644047 Iter = 0.052234915 Bottom = 0.000245201 <<< 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 = 263576961.5 MLMG: Initial residual (resid0) = 2.879265447 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01544932099, 5.861407954e-11 MLMG: Timers: Solve = 0.015131232 Iter = 0.010938437 Bottom = 4.6124e-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.21 MLMG: Initial residual (resid0) = 14829531.21 MLMG: Final Iter. 16 resid, resid/bnorm = 0.001333530992, 8.992401537e-11 MLMG: Timers: Solve = 0.503716666 Iter = 0.498692582 Bottom = 0.000557197 Done calling nodal solver Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974 Timing summary: Advection :0.959132296 seconds MAC Proj :0.094726725 seconds Nodal Proj :0.555876075 seconds Reactions :1.045725062 seconds Misc :0.221213895 seconds Base State :0.004763611 seconds Time to advance time step: 2.877674868 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.822101361 Time to regrid: 0.129779988 Call to estdt for level 0 gives dt_lev = 4.022409037 Call to estdt for level 1 gives dt_lev = 0.8210388353 Minimum estdt over all levels = 0.8210388353 Call to estdt at beginning of step 3 gives dt =0.8210388353 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 = 263106328.9 MLMG: Initial residual (resid0) = 263106328.9 MLMG: Final Iter. 5 resid, resid/bnorm = 0.1797550907, 6.832032186e-10 MLMG: Timers: Solve = 0.05505373 Iter = 0.05078178 Bottom = 0.000218806 <<< 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 = 263106328.7 MLMG: Initial residual (resid0) = 120.9403621 MLMG: Final Iter. 1 resid, resid/bnorm = 0.1138474956, 4.327052723e-10 MLMG: Timers: Solve = 0.015359336 Iter = 0.010977835 Bottom = 4.2804e-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.18 MLMG: Initial residual (resid0) = 12313408.18 MLMG: Final Iter. 8 resid, resid/bnorm = 0.001084003132, 8.803437005e-11 MLMG: Timers: Solve = 0.258726921 Iter = 0.253686352 Bottom = 0.00025553 Done calling nodal solver Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571 Timing summary: Advection :0.966087498 seconds MAC Proj :0.093700714 seconds Nodal Proj :0.307329063 seconds Reactions :1.044120086 seconds Misc :0.223569049 seconds Base State :0.004701741 seconds Time to advance time step: 2.635234595 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.825198222 Total Time: 20.70893304 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.04-28-g56661522b4cd) finalized