MPI initialized with 2 MPI processes MPI initialized with thread support level 0 OMP initialized with 4 OMP threads AMReX (22.01-26-g65f906216a5e) initialized Calling Setup() Calling ReadParameters() reading extern runtime parameters ... Calling VariableSetup() Calling BCSetup() Calling BaseStateGeometry::Init() Calling Init() Calling InitData() initdata model_File = cutoff densities: low density cutoff (for mapping the model) = 1e-10 buoyancy cutoff density (for zeroing rho - rho_0, centrifugal term) = 5e-10 anelastic cutoff = 1e-10 Maximum HSE Error = 4.547473509e-13 (after putting initial model into base state arrays, and for density < base_cutoff_density) Writing plotfile RT_pltInitData after InitData Time to write plotfile: 0.06770944595 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002477764 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile RT_pltafter_InitProj after InitProj Time to write plotfile: 0.05981254578 Call to firstdt for level 0 gives dt_lev = 0.0006089835859 Multiplying dt_lev by init_shrink; dt_lev = 6.089835859e-05 Minimum firstdt over all levels = 6.089835859e-05 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.00235743 Iter = 0 Bottom = 0 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 0.08660444773 Minimum estdt over all levels = 0.08660444773 Call to estdt at end of istep_divu_iter = 1 gives dt = 0.08660444773 Multiplying dt by init_shrink; dt = 0.008660444773 Ignoring this new dt since it's larger than the previous dt = 6.089835859e-05 Writing plotfile RT_pltafter_DivuIter after final DivuIter Time to write plotfile: 0.06118154526 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 6.089835859e-05 Cell Count: Level 0, 131072 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.002421164472 MLMG: Initial residual (resid0) = 0.002421164472 MLMG: Final Iter. 9 resid, resid/bnorm = 1.751907767e-13, 7.235806522e-11 MLMG: Timers: Solve = 0.034296523 Iter = 0.032333428 Bottom = 0.002866098 <<< 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 = 0.002421205691 MLMG: Initial residual (resid0) = 2.151449931e-06 MLMG: Final Iter. 7 resid, resid/bnorm = 9.927909435e-14, 4.100399017e-11 MLMG: Timers: Solve = 0.024769678 Iter = 0.023187752 Bottom = 0.002111428 <<< 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 = 58.12657035 MLMG: Initial residual (resid0) = 58.12657035 MLMG: Final Iter. 8 resid, resid/bnorm = 2.15578666e-11, 3.708780077e-13 MLMG: Timers: Solve = 0.061342755 Iter = 0.059976914 Bottom = 0.002760517 Done calling nodal solver Timestep 0 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05 Timing summary: Advection :0.07319092751 seconds MAC Proj :0.06539869308 seconds Nodal Proj :0.06742429733 seconds Reactions :0.00680565834 seconds Misc :0.02654147148 seconds Base State :0.0004138946533 seconds Time to advance time step: 0.2423553467 Writing plotfile 0 after all initialization Time to write plotfile: 0.04795432091 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 6.089835859e-05 Cell Count: Level 0, 131072 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 2.394238647e-05 MLMG: Initial residual (resid0) = 2.394238647e-05 MLMG: Final Iter. 9 resid, resid/bnorm = 3.181792119e-16, 1.328936914e-11 MLMG: Timers: Solve = 0.031739591 Iter = 0.030197333 Bottom = 0.002944348 <<< 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 = 2.394218844e-05 MLMG: Initial residual (resid0) = 1.928423574e-06 MLMG: Final Iter. 9 resid, resid/bnorm = 6.036515824e-16, 2.521288244e-11 MLMG: Timers: Solve = 0.031440923 Iter = 0.029884873 Bottom = 0.002581847 <<< 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 = 0.003539812449 MLMG: Initial residual (resid0) = 0.003539812449 MLMG: Final Iter. 8 resid, resid/bnorm = 1.336170757e-15, 3.77469365e-13 MLMG: Timers: Solve = 0.067933833 Iter = 0.066434158 Bottom = 0.004614481 Done calling nodal solver Timestep 1 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05 Timing summary: Advection :0.06690216064 seconds MAC Proj :0.06840872765 seconds Nodal Proj :0.07419800758 seconds Reactions :0.006360530853 seconds Misc :0.01514554024 seconds Base State :0.0004138946533 seconds Time to advance time step: 0.2349522114 Call to estdt for level 0 gives dt_lev = 0.3128878797 Minimum estdt over all levels = 0.3128878797 Call to estdt at beginning of step 2 gives dt =0.3128878797 dt_growth factor limits the new dt = 6.698819445e-05 Timestep 2 starts with TIME = 6.089835859e-05 DT = 6.698819445e-05 Cell Count: Level 0, 131072 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 3.050043585e-05 MLMG: Initial residual (resid0) = 3.050043585e-05 MLMG: Final Iter. 9 resid, resid/bnorm = 2.284813777e-15, 7.491085662e-11 MLMG: Timers: Solve = 0.034915275 Iter = 0.033579851 Bottom = 0.004712444 <<< 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 = 3.059063412e-05 MLMG: Initial residual (resid0) = 5.66083759e-06 MLMG: Final Iter. 10 resid, resid/bnorm = 4.892619004e-16, 1.59938463e-11 MLMG: Timers: Solve = 0.038114502 Iter = 0.036823244 Bottom = 0.00486832 <<< 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 = 0.003894013591 MLMG: Initial residual (resid0) = 0.003894013591 MLMG: Final Iter. 8 resid, resid/bnorm = 1.456300358e-15, 3.739844056e-13 MLMG: Timers: Solve = 0.074243826 Iter = 0.072876818 Bottom = 0.005167488 Done calling nodal solver Timestep 2 ends with TIME = 0.000127886553 DT = 6.698819445e-05 Timing summary: Advection :0.06616067886 seconds MAC Proj :0.07856583595 seconds Nodal Proj :0.08031654358 seconds Reactions :0.006130456924 seconds Misc :0.01497530937 seconds Base State :0.000417470932 seconds Time to advance time step: 0.2465054989 Call to estdt for level 0 gives dt_lev = 0.3130286362 Minimum estdt over all levels = 0.3130286362 Call to estdt at beginning of step 3 gives dt =0.3130286362 dt_growth factor limits the new dt = 7.36870139e-05 Timestep 3 starts with TIME = 0.000127886553 DT = 7.36870139e-05 Cell Count: Level 0, 131072 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 3.168800511e-05 MLMG: Initial residual (resid0) = 3.168800511e-05 MLMG: Final Iter. 10 resid, resid/bnorm = 5.650717727e-16, 1.783235552e-11 MLMG: Timers: Solve = 0.041888481 Iter = 0.038468962 Bottom = 0.005821868 <<< 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 = 3.178722337e-05 MLMG: Initial residual (resid0) = 1.012935749e-05 MLMG: Final Iter. 10 resid, resid/bnorm = 6.83363312e-16, 2.149804983e-11 MLMG: Timers: Solve = 0.039157758 Iter = 0.037752541 Bottom = 0.005521156 <<< 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 = 0.004283414244 MLMG: Initial residual (resid0) = 0.004283414244 MLMG: Final Iter. 8 resid, resid/bnorm = 1.625435897e-15, 3.794720297e-13 MLMG: Timers: Solve = 0.065251609 Iter = 0.063714101 Bottom = 0.005214195 Done calling nodal solver Timestep 3 ends with TIME = 0.0002015735669 DT = 7.36870139e-05 Timing summary: Advection :0.0868537426 seconds MAC Proj :0.09165811539 seconds Nodal Proj :0.0718357563 seconds Reactions :0.007849216461 seconds Misc :0.01528453827 seconds Base State :0.0007870197296 seconds Time to advance time step: 0.2740221024 Call to estdt for level 0 gives dt_lev = 0.3130285358 Minimum estdt over all levels = 0.3130285358 Call to estdt at beginning of step 4 gives dt =0.3130285358 dt_growth factor limits the new dt = 8.105571529e-05 Timestep 4 starts with TIME = 0.0002015735669 DT = 8.105571529e-05 Cell Count: Level 0, 131072 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 3.570538736e-05 MLMG: Initial residual (resid0) = 3.570538736e-05 MLMG: Final Iter. 10 resid, resid/bnorm = 7.571110825e-16, 2.120439347e-11 MLMG: Timers: Solve = 0.039287217 Iter = 0.037882292 Bottom = 0.005603278 <<< 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 = 3.581452755e-05 MLMG: Initial residual (resid0) = 1.504387741e-05 MLMG: Final Iter. 10 resid, resid/bnorm = 8.870129024e-16, 2.476684639e-11 MLMG: Timers: Solve = 0.039301175 Iter = 0.037932418 Bottom = 0.005387071 <<< 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 = 0.004711755906 MLMG: Initial residual (resid0) = 0.004711755906 MLMG: Final Iter. 8 resid, resid/bnorm = 1.773321073e-15, 3.763609806e-13 MLMG: Timers: Solve = 0.065113442 Iter = 0.063643017 Bottom = 0.005069529 Done calling nodal solver Timestep 4 ends with TIME = 0.0002826292822 DT = 8.105571529e-05 Timing summary: Advection :0.06682038307 seconds MAC Proj :0.08416295052 seconds Nodal Proj :0.0710952282 seconds Reactions :0.00653219223 seconds Misc :0.02336883545 seconds Base State :0.0003838539124 seconds Time to advance time step: 0.2523562908 Call to estdt for level 0 gives dt_lev = 0.313028505 Minimum estdt over all levels = 0.313028505 Call to estdt at beginning of step 5 gives dt =0.313028505 dt_growth factor limits the new dt = 8.916128682e-05 Timestep 5 starts with TIME = 0.0002826292822 DT = 8.916128682e-05 Cell Count: Level 0, 131072 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 4.013180326e-05 MLMG: Initial residual (resid0) = 4.013180326e-05 MLMG: Final Iter. 10 resid, resid/bnorm = 9.681417181e-16, 2.412405224e-11 MLMG: Timers: Solve = 0.03933305 Iter = 0.0379415 Bottom = 0.005598281 <<< 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 = 4.025185748e-05 MLMG: Initial residual (resid0) = 2.044985444e-05 MLMG: Final Iter. 10 resid, resid/bnorm = 1.111034482e-15, 2.760206737e-11 MLMG: Timers: Solve = 0.039196466 Iter = 0.037839917 Bottom = 0.00537932 <<< 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 = 0.005182931901 MLMG: Initial residual (resid0) = 0.005182931901 MLMG: Final Iter. 8 resid, resid/bnorm = 1.960671209e-15, 3.782938395e-13 MLMG: Timers: Solve = 0.064987567 Iter = 0.063479433 Bottom = 0.005110862 Done calling nodal solver Timestep 5 ends with TIME = 0.000371790569 DT = 8.916128682e-05 Timing summary: Advection :0.0669593811 seconds MAC Proj :0.08438777924 seconds Nodal Proj :0.07089829445 seconds Reactions :0.005714893341 seconds Misc :0.01448607445 seconds Base State :0.0003845691681 seconds Time to advance time step: 0.2427909374 Writing plotfile 5 Time to write plotfile: 0.04955267906 Total Time: 1.83842063 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2] [TOP]::amr.check_file(nvals = 1) :: [RT_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] AMReX (22.01-26-g65f906216a5e) finalized