Initializing AMReX (23.11-23-g9b733ec45cd9)... MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.11-23-g9b733ec45cd9) 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) = 3000000 buoyancy cutoff density (for zeroing rho - rho_0, centrifugal term) = 15000000 anelastic cutoff = 3000000 Writing plotfile flame_3d_pltInitData after InitData ... Level 0 create thermal coeffs: Time to write plotfile: 1.156478913 Doing initial projection ... Level 0 create thermal coeffs: Calling nodal solver MLMG: Initial rhs = 26401.29973 MLMG: Initial residual (resid0) = 26401.29973 MLMG: Final Iter. 7 resid, resid/bnorm = 1.286627827e-07, 4.873350327e-12 MLMG: Timers: Solve = 0.248398779 Iter = 0.240678417 Bottom = 0.003650333 Done calling nodal solver Writing plotfile flame_3d_pltafter_InitProj after InitProj ... Level 0 create thermal coeffs: Time to write plotfile: 1.122947587 Call to firstdt for level 0 gives dt_lev = 9.531664523e-08 Multiplying dt_lev by init_shrink; dt_lev = 9.531664523e-09 Minimum firstdt over all levels = 9.531664523e-09 Doing initial divu iteration #1 ... Level 0 create thermal coeffs: Calling nodal solver MLMG: Initial rhs = 1745.991399 MLMG: Initial residual (resid0) = 1745.991399 MLMG: Final Iter. 5 resid, resid/bnorm = 4.073751367e-05, 2.333202426e-08 MLMG: Timers: Solve = 0.189376756 Iter = 0.181350296 Bottom = 0.002939456 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 9.529696764e-08 Minimum estdt over all levels = 9.529696764e-08 Call to estdt at end of istep_divu_iter = 1 gives dt = 9.529696764e-08 Multiplying dt by init_shrink; dt = 9.529696764e-09 Doing initial divu iteration #2 ... Level 0 create thermal coeffs: Calling nodal solver MLMG: Initial rhs = 8.63846235e-05 MLMG: Initial residual (resid0) = 8.63846235e-05 MLMG: Final Iter. 6 resid, resid/bnorm = 2.068701113e-14, 2.394756184e-10 MLMG: Timers: Solve = 0.218334205 Iter = 0.210551697 Bottom = 0.003392796 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 9.529696764e-08 Minimum estdt over all levels = 9.529696764e-08 Call to estdt at end of istep_divu_iter = 2 gives dt = 9.529696764e-08 Multiplying dt by init_shrink; dt = 9.529696764e-09 Ignoring this new dt since it's larger than the previous dt = 9.529696764e-09 Doing initial divu iteration #3 ... Level 0 create thermal coeffs: Calling nodal solver MLMG: Initial rhs = 2.038288652e-06 MLMG: Initial residual (resid0) = 2.038288652e-06 MLMG: Final Iter. 6 resid, resid/bnorm = 6.91178885e-19, 3.390976466e-13 MLMG: Timers: Solve = 0.223735251 Iter = 0.215952134 Bottom = 0.003561004 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 9.529696764e-08 Minimum estdt over all levels = 9.529696764e-08 Call to estdt at end of istep_divu_iter = 3 gives dt = 9.529696764e-08 Multiplying dt by init_shrink; dt = 9.529696764e-09 Ignoring this new dt since it's larger than the previous dt = 9.529696764e-09 Writing plotfile flame_3d_pltafter_DivuIter after final DivuIter ... Level 0 create thermal coeffs: Time to write plotfile: 1.231222487 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 9.529696764e-09 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 3970.663627 MLMG: Initial residual (resid0) = 3970.663627 MLMG: Final Iter. 8 resid, resid/bnorm = 3.293567445e-07, 8.294753106e-11 MLMG: Timers: Solve = 0.076065926 Iter = 0.069480445 Bottom = 0.002859665 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 7.394289443e+21 MLMG: Final Iter. 3 resid, resid/bnorm = 3.423088935e+12, 1.309096751e-13 MLMG: Timers: Solve = 0.027911296 Iter = 0.025651392 Bottom = 0.000158195 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 3731.731767 MLMG: Initial residual (resid0) = 237.2771339 MLMG: Final Iter. 8 resid, resid/bnorm = 2.840775981e-07, 7.612487066e-11 MLMG: Timers: Solve = 0.074169897 Iter = 0.068689282 Bottom = 0.003005721 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 7.380668913e+21 MLMG: Final Iter. 3 resid, resid/bnorm = 3.384434229e+12, 1.294313978e-13 MLMG: Timers: Solve = 0.027220296 Iter = 0.02489709 Bottom = 0.00014219 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 3.214694727e+10 MLMG: Initial residual (resid0) = 3.214694727e+10 MLMG: Final Iter. 6 resid, resid/bnorm = 1.109122314, 3.450163726e-11 MLMG: Timers: Solve = 0.217842296 Iter = 0.209813715 Bottom = 0.00337121 Done calling nodal solver Timestep 0 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.791225304 seconds MAC Proj :0.165856668 seconds Nodal Proj :0.344136534 seconds Thermal :0.183192379 seconds Reactions :1.668398982 seconds Misc :0.211675256 seconds Base State :6.000027497e-08 seconds Time to advance time step: 3.407830711 Writing plotfile 0 after all initialization ... Level 0 create thermal coeffs: Time to write plotfile: 1.115581455 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 9.529696764e-09 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 3872.434812 MLMG: Initial residual (resid0) = 3872.434812 MLMG: Final Iter. 8 resid, resid/bnorm = 1.282508606e-07, 3.311892048e-11 MLMG: Timers: Solve = 0.07553829 Iter = 0.070043415 Bottom = 0.002924369 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 7.39428878e+21 MLMG: Final Iter. 3 resid, resid/bnorm = 3.427383902e+12, 1.310739282e-13 MLMG: Timers: Solve = 0.027133029 Iter = 0.024803987 Bottom = 0.000144599 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 3872.068252 MLMG: Initial residual (resid0) = 1.244370782 MLMG: Final Iter. 6 resid, resid/bnorm = 8.884978797e-08, 2.294633828e-11 MLMG: Timers: Solve = 0.057643098 Iter = 0.052131001 Bottom = 0.002299038 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 7.380722283e+21 MLMG: Final Iter. 3 resid, resid/bnorm = 3.384434229e+12, 1.294313978e-13 MLMG: Timers: Solve = 0.028199956 Iter = 0.025646645 Bottom = 0.00013975 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 306.396391 MLMG: Initial residual (resid0) = 306.396391 MLMG: Final Iter. 6 resid, resid/bnorm = 1.055298515e-08, 3.444226323e-11 MLMG: Timers: Solve = 0.219263431 Iter = 0.210444486 Bottom = 0.003400769 Done calling nodal solver Timestep 1 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.764674983 seconds MAC Proj :0.147937993 seconds Nodal Proj :0.344424997 seconds Thermal :0.179904626 seconds Reactions :1.676536207 seconds Misc :0.193321001 seconds Base State :4.899993655e-08 seconds Time to advance time step: 3.345560972 Call to estdt for level 0 gives dt_lev = 9.530304292e-08 Minimum estdt over all levels = 9.530304292e-08 Call to estdt at beginning of step 2 gives dt =9.530304292e-08 dt_growth factor limits the new dt = 1.143563612e-08 Timestep 2 starts with TIME = 9.529696764e-09 DT = 1.143563612e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 3462.85172 MLMG: Initial residual (resid0) = 3462.85172 MLMG: Final Iter. 8 resid, resid/bnorm = 1.161296552e-07, 3.353584404e-11 MLMG: Timers: Solve = 0.073236413 Iter = 0.067726143 Bottom = 0.002861908 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 8.786804468e+21 MLMG: Final Iter. 3 resid, resid/bnorm = 4.660039516e+12, 1.782145515e-13 MLMG: Timers: Solve = 0.026646733 Iter = 0.024373768 Bottom = 0.000135243 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 3468.144518 MLMG: Initial residual (resid0) = 9.284136097 MLMG: Final Iter. 7 resid, resid/bnorm = 1.289666193e-07, 3.718605686e-11 MLMG: Timers: Solve = 0.064782447 Iter = 0.059201586 Bottom = 0.002628994 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 8.766892021e+21 MLMG: Final Iter. 3 resid, resid/bnorm = 4.62138481e+12, 1.767362741e-13 MLMG: Timers: Solve = 0.026819767 Iter = 0.024566961 Bottom = 0.000140535 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 363.9931592 MLMG: Initial residual (resid0) = 363.9931592 MLMG: Final Iter. 6 resid, resid/bnorm = 1.289093365e-08, 3.541531847e-11 MLMG: Timers: Solve = 0.21979441 Iter = 0.21170774 Bottom = 0.003327152 Done calling nodal solver Timestep 2 ends with TIME = 2.096533288e-08 DT = 1.143563612e-08 Timing summary: Advection :0.805428578 seconds MAC Proj :0.152882893 seconds Nodal Proj :0.345063841 seconds Thermal :0.175361624 seconds Reactions :1.679332277 seconds Misc :0.195509493 seconds Base State :4.999992598e-08 seconds Time to advance time step: 3.357504197 Call to estdt for level 0 gives dt_lev = 9.531080113e-08 Minimum estdt over all levels = 9.531080113e-08 Call to estdt at beginning of step 3 gives dt =9.531080113e-08 dt_growth factor limits the new dt = 1.372276334e-08 Timestep 3 starts with TIME = 2.096533288e-08 DT = 1.372276334e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 2948.940386 MLMG: Initial residual (resid0) = 2948.940386 MLMG: Final Iter. 8 resid, resid/bnorm = 1.20626261e-07, 4.090495065e-11 MLMG: Timers: Solve = 0.07340981 Iter = 0.067895979 Bottom = 0.002892583 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.042263562e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 5.772436046e+12, 2.207560896e-13 MLMG: Timers: Solve = 0.026851322 Iter = 0.024570597 Bottom = 0.000139807 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 2956.725477 MLMG: Initial residual (resid0) = 13.11341261 MLMG: Final Iter. 7 resid, resid/bnorm = 1.512362928e-07, 5.114992717e-11 MLMG: Timers: Solve = 0.065106387 Iter = 0.05959289 Bottom = 0.002560143 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.039324141e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 5.776731013e+12, 2.209203426e-13 MLMG: Timers: Solve = 0.026731711 Iter = 0.024469746 Bottom = 0.000135247 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 432.6231898 MLMG: Initial residual (resid0) = 432.6231898 MLMG: Final Iter. 6 resid, resid/bnorm = 1.510822445e-08, 3.492236386e-11 MLMG: Timers: Solve = 0.220956916 Iter = 0.212801059 Bottom = 0.003430267 Done calling nodal solver Timestep 3 ends with TIME = 3.468809622e-08 DT = 1.372276334e-08 Timing summary: Advection :0.773379235 seconds MAC Proj :0.153172081 seconds Nodal Proj :0.342999194 seconds Thermal :0.175629677 seconds Reactions :1.695322247 seconds Misc :0.197489348 seconds Base State :4.899993655e-08 seconds Time to advance time step: 3.341889834 Call to estdt for level 0 gives dt_lev = 9.532074642e-08 Minimum estdt over all levels = 9.532074642e-08 Call to estdt at beginning of step 4 gives dt =9.532074642e-08 dt_growth factor limits the new dt = 1.646731601e-08 Timestep 4 starts with TIME = 3.468809622e-08 DT = 1.646731601e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 2336.708017 MLMG: Initial residual (resid0) = 2336.708017 MLMG: Final Iter. 8 resid, resid/bnorm = 1.259787794e-07, 5.391293156e-11 MLMG: Timers: Solve = 0.074048441 Iter = 0.068570935 Bottom = 0.002892412 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.233675845e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 7.791070675e+12, 2.979550197e-13 MLMG: Timers: Solve = 0.027274958 Iter = 0.025022216 Bottom = 0.000138822 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 2348.152059 MLMG: Initial residual (resid0) = 18.2573799 MLMG: Final Iter. 7 resid, resid/bnorm = 2.131855581e-07, 9.078865112e-11 MLMG: Timers: Solve = 0.066070002 Iter = 0.060487167 Bottom = 0.002537318 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.229304008e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 7.722351198e+12, 2.95326971e-13 MLMG: Timers: Solve = 0.027266009 Iter = 0.024957318 Bottom = 0.000143604 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 512.5767955 MLMG: Initial residual (resid0) = 512.5767955 MLMG: Final Iter. 6 resid, resid/bnorm = 1.761736979e-08, 3.437020548e-11 MLMG: Timers: Solve = 0.22515568 Iter = 0.216945317 Bottom = 0.00342579 Done calling nodal solver Timestep 4 ends with TIME = 5.115541223e-08 DT = 1.646731601e-08 Timing summary: Advection :0.779560085 seconds MAC Proj :0.154838861 seconds Nodal Proj :0.349080589 seconds Thermal :0.177573687 seconds Reactions :1.764698291 seconds Misc :0.202024309 seconds Base State :5.09999154e-08 seconds Time to advance time step: 3.431680574 Call to estdt for level 0 gives dt_lev = 9.533353448e-08 Minimum estdt over all levels = 9.533353448e-08 Call to estdt at beginning of step 5 gives dt =9.533353448e-08 dt_growth factor limits the new dt = 1.976077921e-08 Timestep 5 starts with TIME = 5.115541223e-08 DT = 1.976077921e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1668.220782 MLMG: Initial residual (resid0) = 1668.220782 MLMG: Final Iter. 8 resid, resid/bnorm = 1.62495553e-07, 9.740650326e-11 MLMG: Timers: Solve = 0.075322149 Iter = 0.069762866 Bottom = 0.002978467 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.456623662e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 9.607841841e+12, 3.674340568e-13 MLMG: Timers: Solve = 0.027395502 Iter = 0.025109362 Bottom = 0.000126672 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 1684.940024 MLMG: Initial residual (resid0) = 71.11602116 MLMG: Final Iter. 8 resid, resid/bnorm = 3.206547804e-08, 1.903063468e-11 MLMG: Timers: Solve = 0.075354573 Iter = 0.069758465 Bottom = 0.003005714 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.450055918e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 9.594956939e+12, 3.669412977e-13 MLMG: Timers: Solve = 0.027280303 Iter = 0.025017792 Bottom = 0.000128533 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 604.7595071 MLMG: Initial residual (resid0) = 604.7595071 MLMG: Final Iter. 6 resid, resid/bnorm = 2.194508575e-08, 3.628729353e-11 MLMG: Timers: Solve = 0.225201352 Iter = 0.217057903 Bottom = 0.003297687 Done calling nodal solver Timestep 5 ends with TIME = 7.091619144e-08 DT = 1.976077921e-08 Timing summary: Advection :0.785190073 seconds MAC Proj :0.165439894 seconds Nodal Proj :0.348967681 seconds Thermal :0.177625482 seconds Reactions :1.818890055 seconds Misc :0.202954866 seconds Base State :4.599996828e-08 seconds Time to advance time step: 3.502999871 Call to estdt for level 0 gives dt_lev = 9.53499232e-08 Minimum estdt over all levels = 9.53499232e-08 Call to estdt at beginning of step 6 gives dt =9.53499232e-08 dt_growth factor limits the new dt = 2.371293505e-08 Timestep 6 starts with TIME = 7.091619144e-08 DT = 2.371293505e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 925.5205018 MLMG: Initial residual (resid0) = 925.5205018 MLMG: Final Iter. 9 resid, resid/bnorm = 1.65595111e-08, 1.789210619e-11 MLMG: Timers: Solve = 0.089113524 Iter = 0.083484205 Bottom = 0.003376417 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.714863853e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.064493822e+12, 2.319252965e-13 MLMG: Timers: Solve = 0.027508736 Iter = 0.025200596 Bottom = 0.000137898 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 949.6258959 MLMG: Initial residual (resid0) = 71.67495479 MLMG: Final Iter. 8 resid, resid/bnorm = 4.434289833e-08, 4.669512334e-11 MLMG: Timers: Solve = 0.074359401 Iter = 0.068782806 Bottom = 0.002972299 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.704973328e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.141803233e+12, 2.348818513e-13 MLMG: Timers: Solve = 0.027040438 Iter = 0.024778168 Bottom = 0.000147469 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 711.0822893 MLMG: Initial residual (resid0) = 711.0822893 MLMG: Final Iter. 6 resid, resid/bnorm = 2.671291766e-08, 3.756656305e-11 MLMG: Timers: Solve = 0.231991222 Iter = 0.223889049 Bottom = 0.003276933 Done calling nodal solver Timestep 6 ends with TIME = 9.462912649e-08 DT = 2.371293505e-08 Timing summary: Advection :0.794581273 seconds MAC Proj :0.179680312 seconds Nodal Proj :0.364682504 seconds Thermal :0.177519385 seconds Reactions :1.79932154 seconds Misc :0.20360426 seconds Base State :5.200035957e-08 seconds Time to advance time step: 3.523272676 Call to estdt for level 0 gives dt_lev = 9.537095711e-08 Minimum estdt over all levels = 9.537095711e-08 Call to estdt at beginning of step 7 gives dt =9.537095711e-08 dt_growth factor limits the new dt = 2.845552206e-08 Timestep 7 starts with TIME = 9.462912649e-08 DT = 2.845552206e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 676.9925932 MLMG: Initial residual (resid0) = 676.9925932 MLMG: Final Iter. 9 resid, resid/bnorm = 2.321326065e-08, 3.428879559e-11 MLMG: Timers: Solve = 0.097239624 Iter = 0.091596676 Bottom = 0.003302498 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 2.012120128e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 1.321561437e+13, 5.054066128e-13 MLMG: Timers: Solve = 0.026737377 Iter = 0.024478039 Bottom = 0.000133389 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 691.7584029 MLMG: Initial residual (resid0) = 44.05710347 MLMG: Final Iter. 8 resid, resid/bnorm = 3.633300594e-08, 5.252268102e-11 MLMG: Timers: Solve = 0.073924334 Iter = 0.068391184 Bottom = 0.003014547 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 1.997022144e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 1.209033294e+13, 4.623723156e-13 MLMG: Timers: Solve = 0.027433741 Iter = 0.025096698 Bottom = 0.000136895 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 827.0677447 MLMG: Initial residual (resid0) = 827.0677447 MLMG: Final Iter. 6 resid, resid/bnorm = 3.073860699e-08, 3.716576688e-11 MLMG: Timers: Solve = 0.232707911 Iter = 0.224503237 Bottom = 0.003289973 Done calling nodal solver Timestep 7 ends with TIME = 1.230846486e-07 DT = 2.845552206e-08 Timing summary: Advection :0.795786859 seconds MAC Proj :0.185655749 seconds Nodal Proj :0.359503041 seconds Thermal :0.17909694 seconds Reactions :1.818038184 seconds Misc :0.205670514 seconds Base State :7.900007404e-08 seconds Time to advance time step: 3.547633225 Call to estdt for level 0 gives dt_lev = 9.539807155e-08 Minimum estdt over all levels = 9.539807155e-08 Call to estdt at beginning of step 8 gives dt =9.539807155e-08 dt_growth factor limits the new dt = 3.414662647e-08 Timestep 8 starts with TIME = 1.230846486e-07 DT = 3.414662647e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 731.311611 MLMG: Initial residual (resid0) = 731.311611 MLMG: Final Iter. 9 resid, resid/bnorm = 2.958142886e-08, 4.044982798e-11 MLMG: Timers: Solve = 0.083667974 Iter = 0.078032856 Bottom = 0.003247971 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 2.35247823e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 2.052564871e+13, 7.849652917e-13 MLMG: Timers: Solve = 0.0289852 Iter = 0.026635152 Bottom = 0.000146003 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 753.4802695 MLMG: Initial residual (resid0) = 56.29847831 MLMG: Final Iter. 8 resid, resid/bnorm = 5.021422567e-08, 6.664305317e-11 MLMG: Timers: Solve = 0.075207441 Iter = 0.06952573 Bottom = 0.003064936 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 2.328684779e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 2.049128897e+13, 7.836512674e-13 MLMG: Timers: Solve = 0.028118351 Iter = 0.025703491 Bottom = 0.000141171 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 940.0866093 MLMG: Initial residual (resid0) = 940.0866093 MLMG: Final Iter. 6 resid, resid/bnorm = 3.495370282e-08, 3.718136444e-11 MLMG: Timers: Solve = 0.22542075 Iter = 0.21716694 Bottom = 0.003276609 Done calling nodal solver Timestep 8 ends with TIME = 1.57231275e-07 DT = 3.414662647e-08 Timing summary: Advection :0.809665629 seconds MAC Proj :0.173605354 seconds Nodal Proj :0.351243798 seconds Thermal :0.206379462 seconds Reactions :1.910379098 seconds Misc :0.209721612 seconds Base State :5.399988368e-08 seconds Time to advance time step: 3.664925943 Call to estdt for level 0 gives dt_lev = 9.543293932e-08 Minimum estdt over all levels = 9.543293932e-08 Call to estdt at beginning of step 9 gives dt =9.543293932e-08 dt_growth factor limits the new dt = 4.097595177e-08 Timestep 9 starts with TIME = 1.57231275e-07 DT = 4.097595177e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 2219.755231 MLMG: Initial residual (resid0) = 2219.755231 MLMG: Final Iter. 9 resid, resid/bnorm = 3.32727268e-08, 1.498936745e-11 MLMG: Timers: Solve = 0.08475181 Iter = 0.079182872 Bottom = 0.003323459 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 2.740383678e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 1.809899219e+13, 6.921623225e-13 MLMG: Timers: Solve = 0.027394118 Iter = 0.024998999 Bottom = 0.000138485 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 2253.32195 MLMG: Initial residual (resid0) = 70.26043889 MLMG: Final Iter. 8 resid, resid/bnorm = 6.901726124e-08, 3.062911682e-11 MLMG: Timers: Solve = 0.074808681 Iter = 0.069268228 Bottom = 0.003029569 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 2.730437807e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 2.02765406e+13, 7.754386152e-13 MLMG: Timers: Solve = 0.027399548 Iter = 0.025132276 Bottom = 0.000137301 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1054.509422 MLMG: Initial residual (resid0) = 1054.509422 MLMG: Final Iter. 6 resid, resid/bnorm = 3.931205583e-08, 3.727994745e-11 MLMG: Timers: Solve = 0.225155081 Iter = 0.216990468 Bottom = 0.003388153 Done calling nodal solver Timestep 9 ends with TIME = 1.982072268e-07 DT = 4.097595177e-08 Timing summary: Advection :0.807919627 seconds MAC Proj :0.174357502 seconds Nodal Proj :0.348553753 seconds Thermal :0.182027807 seconds Reactions :1.928742877 seconds Misc :0.211897391 seconds Base State :4.999992598e-08 seconds Time to advance time step: 3.657416459 Call to estdt for level 0 gives dt_lev = 9.545947407e-08 Minimum estdt over all levels = 9.545947407e-08 Call to estdt at beginning of step 10 gives dt =9.545947407e-08 dt_growth factor limits the new dt = 4.917114212e-08 Timestep 10 starts with TIME = 1.982072268e-07 DT = 4.917114212e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 2063.674887 MLMG: Initial residual (resid0) = 2063.674887 MLMG: Final Iter. 9 resid, resid/bnorm = 3.296236173e-08, 1.597265244e-11 MLMG: Timers: Solve = 0.084153955 Iter = 0.078558904 Bottom = 0.00327816 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 3.225447912e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.767150472e+13, 2.587970943e-12 MLMG: Timers: Solve = 0.027397403 Iter = 0.025122674 Bottom = 0.000138154 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 2116.623671 MLMG: Initial residual (resid0) = 84.89142245 MLMG: Final Iter. 8 resid, resid/bnorm = 9.42551992e-08, 4.453091992e-11 MLMG: Timers: Solve = 0.075436084 Iter = 0.069859832 Bottom = 0.003001263 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 3.22233896e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.628423028e+13, 2.53491721e-12 MLMG: Timers: Solve = 0.027334068 Iter = 0.025068742 Bottom = 0.000137487 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1165.461631 MLMG: Initial residual (resid0) = 1165.461631 MLMG: Final Iter. 6 resid, resid/bnorm = 4.369059114e-08, 3.748779883e-11 MLMG: Timers: Solve = 0.245325896 Iter = 0.233678942 Bottom = 0.003458524 Done calling nodal solver Timestep 10 ends with TIME = 2.473783689e-07 DT = 4.917114212e-08 Timing summary: Advection :0.808194308 seconds MAC Proj :0.174380866 seconds Nodal Proj :0.377084201 seconds Thermal :0.177692437 seconds Reactions :1.936550136 seconds Misc :0.214368016 seconds Base State :5.299989425e-08 seconds Time to advance time step: 3.692231737 Call to estdt for level 0 gives dt_lev = 9.548804265e-08 Minimum estdt over all levels = 9.548804265e-08 Call to estdt at beginning of step 11 gives dt =9.548804265e-08 dt_growth factor limits the new dt = 5.900537055e-08 Timestep 11 starts with TIME = 2.473783689e-07 DT = 5.900537055e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1009.025868 MLMG: Initial residual (resid0) = 1009.025868 MLMG: Final Iter. 9 resid, resid/bnorm = 3.11172812e-08, 3.083893305e-11 MLMG: Timers: Solve = 0.082154024 Iter = 0.076710539 Bottom = 0.00331374 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 3.802441317e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.320473873e+13, 2.417147778e-12 MLMG: Timers: Solve = 0.026820839 Iter = 0.02454691 Bottom = 0.000138781 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 970.1554543 MLMG: Initial residual (resid0) = 131.1492544 MLMG: Final Iter. 9 resid, resid/bnorm = 8.220763711e-09, 8.473656129e-12 MLMG: Timers: Solve = 0.081885501 Iter = 0.076444291 Bottom = 0.003187183 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 3.788690557e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.536081231e+13, 2.499602806e-12 MLMG: Timers: Solve = 0.026847109 Iter = 0.024582777 Bottom = 0.000133717 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1269.054352 MLMG: Initial residual (resid0) = 1269.054352 MLMG: Final Iter. 6 resid, resid/bnorm = 4.939474052e-08, 3.892247833e-11 MLMG: Timers: Solve = 0.219292044 Iter = 0.21120958 Bottom = 0.0033645 Done calling nodal solver Timestep 11 ends with TIME = 3.063837395e-07 DT = 5.900537055e-08 Timing summary: Advection :0.829086278 seconds MAC Proj :0.17873715 seconds Nodal Proj :0.34290681 seconds Thermal :0.175484551 seconds Reactions :1.874321861 seconds Misc :0.212587104 seconds Base State :5.49998731e-08 seconds Time to advance time step: 3.617051848 Call to estdt for level 0 gives dt_lev = 9.552935221e-08 Minimum estdt over all levels = 9.552935221e-08 Call to estdt at beginning of step 12 gives dt =9.552935221e-08 dt_growth factor limits the new dt = 7.080644466e-08 Timestep 12 starts with TIME = 3.063837395e-07 DT = 7.080644466e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 895.3141249 MLMG: Initial residual (resid0) = 895.3141249 MLMG: Final Iter. 9 resid, resid/bnorm = 3.053258979e-08, 3.41026562e-11 MLMG: Timers: Solve = 0.08112197 Iter = 0.075666663 Bottom = 0.003139723 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 4.484967196e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 1.487772376e+15, 5.689708978e-11 MLMG: Timers: Solve = 0.027330112 Iter = 0.02505852 Bottom = 0.000138379 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 410.7783862 MLMG: Initial residual (resid0) = 482.3381304 MLMG: Final Iter. 9 resid, resid/resid0 = 1.91893056e-08, 3.97839283e-11 MLMG: Timers: Solve = 0.082319547 Iter = 0.076825274 Bottom = 0.003347303 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 4.441461795e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 1.401147181e+15, 5.358427016e-11 MLMG: Timers: Solve = 0.027017424 Iter = 0.024750617 Bottom = 0.000142279 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1311.74867 MLMG: Initial residual (resid0) = 1311.74867 MLMG: Final Iter. 6 resid, resid/bnorm = 5.879289677e-08, 4.482024502e-11 MLMG: Timers: Solve = 0.219064957 Iter = 0.211014678 Bottom = 0.003120737 Done calling nodal solver Timestep 12 ends with TIME = 3.771901841e-07 DT = 7.080644466e-08 Timing summary: Advection :0.786896295 seconds MAC Proj :0.178111547 seconds Nodal Proj :0.340992978 seconds Thermal :0.176860057 seconds Reactions :1.898498799 seconds Misc :0.214474852 seconds Base State :4.409998837e-07 seconds Time to advance time step: 3.599734915 Call to estdt for level 0 gives dt_lev = 9.557694465e-08 Minimum estdt over all levels = 9.557694465e-08 Call to estdt at beginning of step 13 gives dt =9.557694465e-08 dt_growth factor limits the new dt = 8.496773359e-08 Timestep 13 starts with TIME = 3.771901841e-07 DT = 8.496773359e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1121.364477 MLMG: Initial residual (resid0) = 1121.364477 MLMG: Final Iter. 9 resid, resid/bnorm = 2.545914413e-08, 2.270371914e-11 MLMG: Timers: Solve = 0.082277578 Iter = 0.076875337 Bottom = 0.003335596 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.328238633e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 1.034228125e+14, 3.955213271e-12 MLMG: Timers: Solve = 0.034652192 Iter = 0.032356231 Bottom = 0.000176532 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 1215.562869 MLMG: Initial residual (resid0) = 164.8126807 MLMG: Final Iter. 9 resid, resid/bnorm = 1.433244279e-08, 1.179078693e-11 MLMG: Timers: Solve = 0.08193832 Iter = 0.076422529 Bottom = 0.003239243 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.202860905e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 9.680426788e+13, 3.702099332e-12 MLMG: Timers: Solve = 0.035212197 Iter = 0.032929134 Bottom = 0.000175713 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1498.115759 MLMG: Initial residual (resid0) = 1498.115759 MLMG: Final Iter. 6 resid, resid/bnorm = 6.250402107e-08, 4.17217566e-11 MLMG: Timers: Solve = 0.219258375 Iter = 0.210819438 Bottom = 0.003342545 Done calling nodal solver Timestep 13 ends with TIME = 4.621579177e-07 DT = 8.496773359e-08 Timing summary: Advection :0.784533223 seconds MAC Proj :0.178856374 seconds Nodal Proj :0.341218826 seconds Thermal :0.191622885 seconds Reactions :1.922011433 seconds Misc :0.217050976 seconds Base State :4.899993655e-08 seconds Time to advance time step: 3.639180305 Call to estdt for level 0 gives dt_lev = 9.561756174e-08 Minimum estdt over all levels = 9.561756174e-08 Call to estdt at beginning of step 14 gives dt =9.561756174e-08 Timestep 14 starts with TIME = 4.621579177e-07 DT = 9.561756174e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1296.296742 MLMG: Initial residual (resid0) = 1296.296742 MLMG: Final Iter. 9 resid, resid/bnorm = 2.983898639e-08, 2.301863874e-11 MLMG: Timers: Solve = 0.081992772 Iter = 0.076502628 Bottom = 0.003246237 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.954856914e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.928603821e+14, 1.884852942e-11 MLMG: Timers: Solve = 0.034693928 Iter = 0.03243999 Bottom = 0.000176305 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 872.6896169 MLMG: Initial residual (resid0) = 418.9789599 MLMG: Final Iter. 9 resid, resid/bnorm = 1.858393262e-08, 2.129500828e-11 MLMG: Timers: Solve = 0.082109404 Iter = 0.076641939 Bottom = 0.003350514 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.591455313e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.613739769e+14, 1.764439037e-11 MLMG: Timers: Solve = 0.034996834 Iter = 0.032737191 Bottom = 0.000179013 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1661.544989 MLMG: Initial residual (resid0) = 1661.544989 MLMG: Final Iter. 6 resid, resid/bnorm = 6.68269422e-08, 4.021976092e-11 MLMG: Timers: Solve = 0.219498974 Iter = 0.211384752 Bottom = 0.003362602 Done calling nodal solver Timestep 14 ends with TIME = 5.577754794e-07 DT = 9.561756174e-08 Timing summary: Advection :0.777274511 seconds MAC Proj :0.178742444 seconds Nodal Proj :0.34144864 seconds Thermal :0.191648538 seconds Reactions :1.936139653 seconds Misc :0.217850481 seconds Base State :5.09999154e-08 seconds Time to advance time step: 3.647027168 Call to estdt for level 0 gives dt_lev = 9.567387237e-08 Minimum estdt over all levels = 9.567387237e-08 Call to estdt at beginning of step 15 gives dt =9.567387237e-08 Timestep 15 starts with TIME = 5.577754794e-07 DT = 9.567387237e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1193.721011 MLMG: Initial residual (resid0) = 1193.721011 MLMG: Final Iter. 9 resid, resid/bnorm = 3.24849907e-08, 2.72132185e-11 MLMG: Timers: Solve = 0.082691439 Iter = 0.07718261 Bottom = 0.003319395 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.593320368e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.802546531e+14, 1.836644674e-11 MLMG: Timers: Solve = 0.035018374 Iter = 0.032744876 Bottom = 0.000180021 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 817.6507876 MLMG: Initial residual (resid0) = 880.7089388 MLMG: Final Iter. 9 resid, resid/resid0 = 1.902691338e-08, 2.160408796e-11 MLMG: Timers: Solve = 0.082342024 Iter = 0.076825587 Bottom = 0.003253918 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.328089736e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.475141174e+14, 1.71143458e-11 MLMG: Timers: Solve = 0.03479298 Iter = 0.032513067 Bottom = 0.000180597 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1355.222809 MLMG: Initial residual (resid0) = 1355.222809 MLMG: Final Iter. 6 resid, resid/bnorm = 6.157733878e-08, 4.543705905e-11 MLMG: Timers: Solve = 0.219102049 Iter = 0.211011313 Bottom = 0.003214281 Done calling nodal solver Timestep 15 ends with TIME = 6.534493518e-07 DT = 9.567387237e-08 Timing summary: Advection :0.788687474 seconds MAC Proj :0.179919673 seconds Nodal Proj :0.341444956 seconds Thermal :0.191580149 seconds Reactions :1.943434822 seconds Misc :0.219749715 seconds Base State :4.69999577e-08 seconds Time to advance time step: 3.668710688 Call to estdt for level 0 gives dt_lev = 9.571056553e-08 Minimum estdt over all levels = 9.571056553e-08 Call to estdt at beginning of step 16 gives dt =9.571056553e-08 Timestep 16 starts with TIME = 6.534493518e-07 DT = 9.571056553e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 960.7040972 MLMG: Initial residual (resid0) = 960.7040972 MLMG: Final Iter. 9 resid, resid/bnorm = 2.914953257e-08, 3.034184267e-11 MLMG: Timers: Solve = 0.082562878 Iter = 0.077104181 Bottom = 0.003324156 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.593814956e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.596087453e+14, 1.757688237e-11 MLMG: Timers: Solve = 0.035118711 Iter = 0.032858472 Bottom = 0.000187687 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 978.3825404 MLMG: Initial residual (resid0) = 374.842397 MLMG: Final Iter. 9 resid, resid/bnorm = 1.743438816e-08, 1.78196027e-11 MLMG: Timers: Solve = 0.081233888 Iter = 0.07572283 Bottom = 0.003270962 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.4189816e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.321252496e+14, 1.652582715e-11 MLMG: Timers: Solve = 0.034975385 Iter = 0.032716364 Bottom = 0.000179645 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1231.550857 MLMG: Initial residual (resid0) = 1231.550857 MLMG: Final Iter. 6 resid, resid/bnorm = 5.894523269e-08, 4.78626054e-11 MLMG: Timers: Solve = 0.224636234 Iter = 0.21644621 Bottom = 0.003490849 Done calling nodal solver Timestep 16 ends with TIME = 7.491599174e-07 DT = 9.571056553e-08 Timing summary: Advection :0.784203393 seconds MAC Proj :0.178429462 seconds Nodal Proj :0.346490558 seconds Thermal :0.192200395 seconds Reactions :1.932773836 seconds Misc :0.218766017 seconds Base State :5.199990483e-08 seconds Time to advance time step: 3.656805824 Call to estdt for level 0 gives dt_lev = 9.575830833e-08 Minimum estdt over all levels = 9.575830833e-08 Call to estdt at beginning of step 17 gives dt =9.575830833e-08 Timestep 17 starts with TIME = 7.491599174e-07 DT = 9.575830833e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1812.753793 MLMG: Initial residual (resid0) = 1812.753793 MLMG: Final Iter. 9 resid, resid/bnorm = 2.464366844e-08, 1.359460315e-11 MLMG: Timers: Solve = 0.084419666 Iter = 0.078795331 Bottom = 0.003424569 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.741080675e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.531491145e+14, 1.732984579e-11 MLMG: Timers: Solve = 0.035848612 Iter = 0.033582196 Bottom = 0.000179547 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 736.4775856 MLMG: Initial residual (resid0) = 1230.328214 MLMG: Final Iter. 9 resid, resid/resid0 = 2.688260769e-08, 2.184994815e-11 MLMG: Timers: Solve = 0.083438799 Iter = 0.077836497 Bottom = 0.003460762 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.190657932e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.262454394e+14, 1.630096473e-11 MLMG: Timers: Solve = 0.035544611 Iter = 0.033291676 Bottom = 0.000193731 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1403.796806 MLMG: Initial residual (resid0) = 1403.796806 MLMG: Final Iter. 6 resid, resid/bnorm = 5.714036888e-08, 4.070415934e-11 MLMG: Timers: Solve = 0.224541821 Iter = 0.216387802 Bottom = 0.003408325 Done calling nodal solver Timestep 17 ends with TIME = 8.449182257e-07 DT = 9.575830833e-08 Timing summary: Advection :0.794507384 seconds MAC Proj :0.182722699 seconds Nodal Proj :0.348590145 seconds Thermal :0.194625484 seconds Reactions :2.035861459 seconds Misc :0.224113582 seconds Base State :4.799994713e-08 seconds Time to advance time step: 3.784350761 Call to estdt for level 0 gives dt_lev = 9.578894066e-08 Minimum estdt over all levels = 9.578894066e-08 Call to estdt at beginning of step 18 gives dt =9.578894066e-08 Timestep 18 starts with TIME = 8.449182257e-07 DT = 9.578894066e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 665.9316194 MLMG: Initial residual (resid0) = 665.9316194 MLMG: Final Iter. 9 resid, resid/bnorm = 2.425420575e-08, 3.642146588e-11 MLMG: Timers: Solve = 0.081896096 Iter = 0.076418377 Bottom = 0.003042506 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.16192633e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.411060262e+14, 1.686928026e-11 MLMG: Timers: Solve = 0.034639647 Iter = 0.0323728 Bottom = 0.000170432 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 566.9129463 MLMG: Initial residual (resid0) = 567.9651174 MLMG: Final Iter. 9 resid, resid/resid0 = 1.513464554e-08, 2.664713919e-11 MLMG: Timers: Solve = 0.08342022 Iter = 0.077833207 Bottom = 0.003363256 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 4.944630524e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.131801488e+14, 1.580130697e-11 MLMG: Timers: Solve = 0.035102553 Iter = 0.032828112 Bottom = 0.000171418 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 890.9820445 MLMG: Initial residual (resid0) = 890.9820445 MLMG: Final Iter. 6 resid, resid/bnorm = 5.488209354e-08, 6.159730589e-11 MLMG: Timers: Solve = 0.219439213 Iter = 0.211317202 Bottom = 0.003285492 Done calling nodal solver Timestep 18 ends with TIME = 9.407071663e-07 DT = 9.578894066e-08 Timing summary: Advection :0.794456148 seconds MAC Proj :0.179940047 seconds Nodal Proj :0.34151649 seconds Thermal :0.191682689 seconds Reactions :1.948909058 seconds Misc :0.219570674 seconds Base State :5.500032785e-08 seconds Time to advance time step: 3.679978256 Call to estdt for level 0 gives dt_lev = 9.583007799e-08 Minimum estdt over all levels = 9.583007799e-08 Call to estdt at beginning of step 19 gives dt =9.583007799e-08 Timestep 19 starts with TIME = 9.407071663e-07 DT = 9.583007799e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1480.231162 MLMG: Initial residual (resid0) = 1480.231162 MLMG: Final Iter. 9 resid, resid/bnorm = 2.569805702e-08, 1.736084044e-11 MLMG: Timers: Solve = 0.082593502 Iter = 0.077094449 Bottom = 0.003243858 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.43512335e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.23926157e+14, 1.621226809e-11 MLMG: Timers: Solve = 0.035108251 Iter = 0.032827782 Bottom = 0.000175736 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 759.3367994 MLMG: Initial residual (resid0) = 715.3480354 MLMG: Final Iter. 9 resid, resid/bnorm = 1.394212745e-08, 1.836092688e-11 MLMG: Timers: Solve = 0.082501216 Iter = 0.076976683 Bottom = 0.003483622 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.064650233e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.010039166e+14, 1.53356496e-11 MLMG: Timers: Solve = 0.034859261 Iter = 0.03258668 Bottom = 0.000170714 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 1137.878579 MLMG: Initial residual (resid0) = 1137.878579 MLMG: Final Iter. 6 resid, resid/bnorm = 5.009496817e-08, 4.402488025e-11 MLMG: Timers: Solve = 0.218834763 Iter = 0.210721166 Bottom = 0.003415106 Done calling nodal solver Timestep 19 ends with TIME = 1.036537244e-06 DT = 9.583007799e-08 Timing summary: Advection :0.806274086 seconds MAC Proj :0.1797267 seconds Nodal Proj :0.355352652 seconds Thermal :0.191769684 seconds Reactions :1.919548406 seconds Misc :0.219224309 seconds Base State :5.199990483e-08 seconds Time to advance time step: 3.675823379 Call to estdt for level 0 gives dt_lev = 9.585646998e-08 Minimum estdt over all levels = 9.585646998e-08 Call to estdt at beginning of step 20 gives dt =9.585646998e-08 Timestep 20 starts with TIME = 1.036537244e-06 DT = 9.585646998e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1354.880928 MLMG: Initial residual (resid0) = 1354.880928 MLMG: Final Iter. 9 resid, resid/bnorm = 2.324378556e-08, 1.715559285e-11 MLMG: Timers: Solve = 0.081966742 Iter = 0.076416139 Bottom = 0.003408719 <<< STEP 4 : advect base >>> ... Level 0 create thermal coeffs: : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 5.345657735e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.181064763e+14, 1.598970521e-11 MLMG: Timers: Solve = 0.034945279 Iter = 0.032688975 Bottom = 0.000175359 <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> ... Level 0 create thermal coeffs: <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 299.7337979 MLMG: Initial residual (resid0) = 1050.094681 MLMG: Final Iter. 9 resid, resid/resid0 = 1.293341256e-08, 1.231642517e-11 MLMG: Timers: Solve = 0.083349819 Iter = 0.077861861 Bottom = 0.003492701 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> ... Level 0 create thermal coeffs: MLMG: Initial rhs = 2.614847933e+25 MLMG: Initial residual (resid0) = 4.727742204e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 3.949480127e+14, 1.510405281e-11 MLMG: Timers: Solve = 0.034886162 Iter = 0.032628984 Bottom = 0.000172251 <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> ... Level 0 create thermal coeffs: <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 933.677207 MLMG: Initial residual (resid0) = 933.677207 MLMG: Final Iter. 6 resid, resid/bnorm = 4.947625598e-08, 5.299075056e-11 MLMG: Timers: Solve = 0.22446177 Iter = 0.216318234 Bottom = 0.003468542 Done calling nodal solver Timestep 20 ends with TIME = 1.132393714e-06 DT = 9.585646998e-08 Timing summary: Advection :0.804215009 seconds MAC Proj :0.180039566 seconds Nodal Proj :0.347696687 seconds Thermal :0.191761178 seconds Reactions :1.957896728 seconds Misc :0.243060724 seconds Base State :4.499997885e-08 seconds Time to advance time step: 3.728732529 Writing plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.36599173 Writing small plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.364293293 Total Time: 87.13949696 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::amr.ref_ratio(nvals = 6) :: [2, 2, 2, 2, 2, 2] [TOP]::maestro.temperr(nvals = 3) :: [6.5e8, 6.5e8, 6.5e8] [TOP]::amr.check_file(nvals = 1) :: [flame_3d_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] AMReX (23.11-23-g9b733ec45cd9) finalized