Initializing AMReX (26.02-52-g55c25f77198e)... MPI initialized with 2 MPI processes MPI initialized with thread support level 0 OMP initialized with 4 OMP threads AMReX (26.02-52-g55c25f77198e) 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.062682887 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.003210878 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile RT_pltafter_InitProj after InitProj Time to write plotfile: 0.049356495 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.003147885 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.047352031 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.031760484 Iter = 0.030270053 Bottom = 0.003724739 <<< 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.024704805 Iter = 0.023189153 Bottom = 0.002729355 <<< 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.051896368 Iter = 0.050411895 Bottom = 0.002514335 Done calling nodal solver Timestep 0 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05 Timing summary: Advection :0.079395452 seconds MAC Proj :0.063002985 seconds Nodal Proj :0.05810737 seconds Reactions :0.003943561 seconds Misc :0.02482555 seconds Base State :0.0004079629998 seconds Time to advance time step: 0.232668968 Writing plotfile 0 after all initialization Time to write plotfile: 0.037941527 Beginning main evolution  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.181792236e-16, 1.328936963e-11 MLMG: Timers: Solve = 0.031162983 Iter = 0.029683708 Bottom = 0.00382822 <<< 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.036511589e-16, 2.521286475e-11 MLMG: Timers: Solve = 0.0304018 Iter = 0.028922975 Bottom = 0.003283462 <<< 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.329231863e-15, 3.755091216e-13 MLMG: Timers: Solve = 0.050065361 Iter = 0.048563702 Bottom = 0.002462405 Done calling nodal solver Timestep 1 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05 Timing summary: Advection :0.071637847 seconds MAC Proj :0.068325037 seconds Nodal Proj :0.055877509 seconds Reactions :0.003512118 seconds Misc :0.011882045 seconds Base State :0.0003851600002 seconds Time to advance time step: 0.214610655 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.284814624e-15, 7.491088439e-11 MLMG: Timers: Solve = 0.030705015 Iter = 0.029287808 Bottom = 0.003651754 <<< 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.892627475e-16, 1.599387399e-11 MLMG: Timers: Solve = 0.033728944 Iter = 0.032252589 Bottom = 0.003705996 <<< 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.490994828e-15, 3.828940996e-13 MLMG: Timers: Solve = 0.05061445 Iter = 0.049123367 Bottom = 0.002522092 Done calling nodal solver Timestep 2 ends with TIME = 0.000127886553 DT = 6.698819445e-05 Timing summary: Advection :0.071924837 seconds MAC Proj :0.071245365 seconds Nodal Proj :0.056391955 seconds Reactions :0.003785043 seconds Misc :0.012174744 seconds Base State :0.000395953 seconds Time to advance time step: 0.215836169 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.650700787e-16, 1.783230205e-11 MLMG: Timers: Solve = 0.034098464 Iter = 0.032678965 Bottom = 0.003944531 <<< 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.833616179e-16, 2.149799653e-11 MLMG: Timers: Solve = 0.050202178 Iter = 0.048754009 Bottom = 0.00360405 <<< 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.614810716e-15, 3.769914894e-13 MLMG: Timers: Solve = 0.050881511 Iter = 0.049281739 Bottom = 0.002360465 Done calling nodal solver Timestep 3 ends with TIME = 0.0002015735669 DT = 7.36870139e-05 Timing summary: Advection :0.083373889 seconds MAC Proj :0.092032447 seconds Nodal Proj :0.058917159 seconds Reactions :0.005962846 seconds Misc :0.012536834 seconds Base State :0.000407952 seconds Time to advance time step: 0.253262855 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.033094028 Iter = 0.031546692 Bottom = 0.003661529 <<< 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.036640763 Iter = 0.035068295 Bottom = 0.003532969 <<< 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.768550584e-15, 3.753485153e-13 MLMG: Timers: Solve = 0.050115965 Iter = 0.048595613 Bottom = 0.00236566 Done calling nodal solver Timestep 4 ends with TIME = 0.0002826292822 DT = 8.105571529e-05 Timing summary: Advection :0.073659647 seconds MAC Proj :0.076713378 seconds Nodal Proj :0.056355765 seconds Reactions :0.003841595 seconds Misc :0.016970211 seconds Base State :0.0008547179999 seconds Time to advance time step: 0.227858871 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.681425651e-16, 2.412407334e-11 MLMG: Timers: Solve = 0.034693279 Iter = 0.033160551 Bottom = 0.003624142 <<< 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.111036176e-15, 2.760210946e-11 MLMG: Timers: Solve = 0.033446959 Iter = 0.032021701 Bottom = 0.00355802 <<< 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.942022931e-15, 3.746958224e-13 MLMG: Timers: Solve = 0.05015381 Iter = 0.048670572 Bottom = 0.002303504 Done calling nodal solver Timestep 5 ends with TIME = 0.000371790569 DT = 8.916128682e-05 Timing summary: Advection :0.069482833 seconds MAC Proj :0.074633195 seconds Nodal Proj :0.055856517 seconds Reactions :0.002657368 seconds Misc :0.010810554 seconds Base State :0.000363119 seconds Time to advance time step: 0.213756357 Writing plotfile 5 Time to write plotfile: 0.03767344 Total Time: 1.656627401 Unused ParmParse Variables: [TOP]::amr.check_file(nvals = 1) :: [RT_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] [TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2] AMReX (26.02-52-g55c25f77198e) finalized