Initializing AMReX (23.10-21-g152fdac6aa4c)... MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.10-21-g152fdac6aa4c) 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.130926728 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.252668637 Iter = 0.245364994 Bottom = 0.004172908 Done calling nodal solver Writing plotfile flame_3d_pltafter_InitProj after InitProj ... Level 0 create thermal coeffs: Time to write plotfile: 1.123549981 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.183057286 Iter = 0.175796591 Bottom = 0.003010974 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.219161662 Iter = 0.211945672 Bottom = 0.003692352 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.217859314 Iter = 0.210569248 Bottom = 0.003672687 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.26644034 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.077825648 Iter = 0.071994885 Bottom = 0.00340281 <<< 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.027855231 Iter = 0.025952757 Bottom = 0.000172515 <<< 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.075704211 Iter = 0.070302766 Bottom = 0.003616476 <<< 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.027167118 Iter = 0.025326143 Bottom = 0.000157913 <<< 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.221696828 Iter = 0.214393002 Bottom = 0.003719845 Done calling nodal solver Timestep 0 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.718470852 seconds MAC Proj :0.167971276 seconds Nodal Proj :0.340649022 seconds Thermal :0.172733477 seconds Reactions :1.668346939 seconds Misc :0.189964773 seconds Base State :4.599996828e-08 seconds Time to advance time step: 3.294010666 Writing plotfile 0 after all initialization ... Level 0 create thermal coeffs: Time to write plotfile: 1.102237894 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.0755379 Iter = 0.070210828 Bottom = 0.003500532 <<< 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.027226578 Iter = 0.025380254 Bottom = 0.000171928 <<< 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.058276144 Iter = 0.052897253 Bottom = 0.002720772 <<< 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.027073259 Iter = 0.02523491 Bottom = 0.00015914 <<< 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.219891666 Iter = 0.212500581 Bottom = 0.003644776 Done calling nodal solver Timestep 1 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.708834115 seconds MAC Proj :0.147887388 seconds Nodal Proj :0.337331188 seconds Thermal :0.167779829 seconds Reactions :1.668427609 seconds Misc :0.188880951 seconds Base State :5.400011105e-08 seconds Time to advance time step: 3.253144191 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.075442835 Iter = 0.070125179 Bottom = 0.003394403 <<< 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.027138713 Iter = 0.025308394 Bottom = 0.000160804 <<< 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.066581061 Iter = 0.061233856 Bottom = 0.003159699 <<< 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.027060776 Iter = 0.025228998 Bottom = 0.000159073 <<< 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.219518689 Iter = 0.212109132 Bottom = 0.003630296 Done calling nodal solver Timestep 2 ends with TIME = 2.096533288e-08 DT = 1.143563612e-08 Timing summary: Advection :0.71529736 seconds MAC Proj :0.155982727 seconds Nodal Proj :0.336607531 seconds Thermal :0.167659759 seconds Reactions :1.675050281 seconds Misc :0.193277725 seconds Base State :4.799994713e-08 seconds Time to advance time step: 3.248113628 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.075702718 Iter = 0.070337787 Bottom = 0.003432429 <<< 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.027177156 Iter = 0.025341302 Bottom = 0.000161427 <<< 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.066930931 Iter = 0.061585753 Bottom = 0.003008115 <<< 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.027633519 Iter = 0.025791388 Bottom = 0.000163513 <<< 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.220462424 Iter = 0.213044429 Bottom = 0.003742134 Done calling nodal solver Timestep 3 ends with TIME = 3.468809622e-08 DT = 1.372276334e-08 Timing summary: Advection :0.716698787 seconds MAC Proj :0.156612773 seconds Nodal Proj :0.337550552 seconds Thermal :0.167937714 seconds Reactions :1.683595513 seconds Misc :0.19414979 seconds Base State :4.999992598e-08 seconds Time to advance time step: 3.260758771 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.07547137 Iter = 0.070135859 Bottom = 0.003439985 <<< 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.02716359 Iter = 0.02532354 Bottom = 0.0001673 <<< 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.066877105 Iter = 0.061521689 Bottom = 0.00297505 <<< 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.027282926 Iter = 0.025451744 Bottom = 0.000164749 <<< 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.219906206 Iter = 0.212440812 Bottom = 0.003585908 Done calling nodal solver Timestep 4 ends with TIME = 5.115541223e-08 DT = 1.646731601e-08 Timing summary: Advection :0.721917812 seconds MAC Proj :0.156287246 seconds Nodal Proj :0.336747987 seconds Thermal :0.167995355 seconds Reactions :1.690781576 seconds Misc :0.194592104 seconds Base State :5.299989425e-08 seconds Time to advance time step: 3.272641093 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.076942061 Iter = 0.071618462 Bottom = 0.003616513 <<< 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.027806287 Iter = 0.025942467 Bottom = 0.000167407 <<< 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.076990028 Iter = 0.071481608 Bottom = 0.00363218 <<< 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.027491772 Iter = 0.025646667 Bottom = 0.000155642 <<< 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.220178038 Iter = 0.212804756 Bottom = 0.003395852 Done calling nodal solver Timestep 5 ends with TIME = 7.091619144e-08 DT = 1.976077921e-08 Timing summary: Advection :0.73757641 seconds MAC Proj :0.168164519 seconds Nodal Proj :0.338882103 seconds Thermal :0.170107026 seconds Reactions :1.719701863 seconds Misc :0.196185559 seconds Base State :5.400011105e-08 seconds Time to advance time step: 3.33490285 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.084745178 Iter = 0.079395552 Bottom = 0.003957576 <<< 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.027096007 Iter = 0.025251589 Bottom = 0.000160763 <<< 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.076050462 Iter = 0.07061616 Bottom = 0.003581206 <<< 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.027099675 Iter = 0.025260554 Bottom = 0.000162318 <<< 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.218226394 Iter = 0.210717512 Bottom = 0.003583582 Done calling nodal solver Timestep 6 ends with TIME = 9.462912649e-08 DT = 2.371293505e-08 Timing summary: Advection :0.73958106 seconds MAC Proj :0.174924529 seconds Nodal Proj :0.337764417 seconds Thermal :0.167835132 seconds Reactions :1.746567671 seconds Misc :0.196638954 seconds Base State :3.350000952e-07 seconds Time to advance time step: 3.36768652 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.084863651 Iter = 0.079380198 Bottom = 0.003894052 <<< 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.027684043 Iter = 0.025797623 Bottom = 0.000167188 <<< 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.075818031 Iter = 0.070449836 Bottom = 0.003557481 <<< 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.027210758 Iter = 0.025396544 Bottom = 0.000152944 <<< 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.217670622 Iter = 0.210389262 Bottom = 0.003448525 Done calling nodal solver Timestep 7 ends with TIME = 1.230846486e-07 DT = 2.845552206e-08 Timing summary: Advection :0.741395566 seconds MAC Proj :0.174754073 seconds Nodal Proj :0.334699128 seconds Thermal :0.168817372 seconds Reactions :1.759642682 seconds Misc :0.198372812 seconds Base State :4.999992598e-08 seconds Time to advance time step: 3.381961461 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.084276398 Iter = 0.078940924 Bottom = 0.003804885 <<< 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.027058824 Iter = 0.02524848 Bottom = 0.000163692 <<< 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.075720466 Iter = 0.070366398 Bottom = 0.003499682 <<< 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.027212334 Iter = 0.025356262 Bottom = 0.000160432 <<< 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.218451338 Iter = 0.211097432 Bottom = 0.003519453 Done calling nodal solver Timestep 8 ends with TIME = 1.57231275e-07 DT = 3.414662647e-08 Timing summary: Advection :0.737147315 seconds MAC Proj :0.17396189 seconds Nodal Proj :0.338615981 seconds Thermal :0.167664667 seconds Reactions :1.784770426 seconds Misc :0.199552757 seconds Base State :5.49998731e-08 seconds Time to advance time step: 3.406021687 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.085892453 Iter = 0.080428398 Bottom = 0.004022154 <<< 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.027181734 Iter = 0.025330845 Bottom = 0.000162812 <<< 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.075548278 Iter = 0.070200853 Bottom = 0.003538479 <<< 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.027228982 Iter = 0.025359847 Bottom = 0.000166344 <<< 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.218569306 Iter = 0.211184204 Bottom = 0.003491715 Done calling nodal solver Timestep 9 ends with TIME = 1.982072268e-07 DT = 4.097595177e-08 Timing summary: Advection :0.743577609 seconds MAC Proj :0.175632322 seconds Nodal Proj :0.335359584 seconds Thermal :0.167659013 seconds Reactions :1.810831974 seconds Misc :0.201281412 seconds Base State :5.299989425e-08 seconds Time to advance time step: 3.438625911 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.084374337 Iter = 0.079004206 Bottom = 0.003870128 <<< 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.027287285 Iter = 0.02540459 Bottom = 0.000161423 <<< 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.076213603 Iter = 0.070865254 Bottom = 0.003576804 <<< 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.026950135 Iter = 0.025112173 Bottom = 0.000168133 <<< 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.219237184 Iter = 0.211921681 Bottom = 0.003643322 Done calling nodal solver Timestep 10 ends with TIME = 2.473783689e-07 DT = 4.917114212e-08 Timing summary: Advection :0.73808744 seconds MAC Proj :0.174591812 seconds Nodal Proj :0.338480872 seconds Thermal :0.167775629 seconds Reactions :1.820029166 seconds Misc :0.205081396 seconds Base State :5.599986252e-08 seconds Time to advance time step: 3.448345245 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.084605248 Iter = 0.079273954 Bottom = 0.003975463 <<< 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.027092853 Iter = 0.025261986 Bottom = 0.00015965 <<< 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.083928566 Iter = 0.078583093 Bottom = 0.00375987 <<< 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.027127686 Iter = 0.025307177 Bottom = 0.000159256 <<< 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.217886767 Iter = 0.210599466 Bottom = 0.003558518 Done calling nodal solver Timestep 11 ends with TIME = 3.063837395e-07 DT = 5.900537055e-08 Timing summary: Advection :0.738449192 seconds MAC Proj :0.182468965 seconds Nodal Proj :0.334765215 seconds Thermal :0.168152879 seconds Reactions :1.846017901 seconds Misc :0.207169197 seconds Base State :4.399998943e-08 seconds Time to advance time step: 3.481327401 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.084028895 Iter = 0.07868758 Bottom = 0.003808041 <<< 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.027088924 Iter = 0.025249951 Bottom = 0.000158054 <<< 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.084493946 Iter = 0.079159379 Bottom = 0.004109455 <<< 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.027509033 Iter = 0.025676315 Bottom = 0.000159975 <<< 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.217551718 Iter = 0.210197417 Bottom = 0.003376904 Done calling nodal solver Timestep 12 ends with TIME = 3.771901841e-07 DT = 7.080644466e-08 Timing summary: Advection :0.737692918 seconds MAC Proj :0.18249759 seconds Nodal Proj :0.333729996 seconds Thermal :0.168428449 seconds Reactions :1.877734327 seconds Misc :0.211186996 seconds Base State :5.20001322e-08 seconds Time to advance time step: 3.515569794 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.0844596 Iter = 0.079121828 Bottom = 0.003969528 <<< 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.035629772 Iter = 0.033778053 Bottom = 0.000223761 <<< 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.084578538 Iter = 0.079239906 Bottom = 0.003893462 <<< 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.035475891 Iter = 0.033641099 Bottom = 0.000212014 <<< 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.218196034 Iter = 0.210780955 Bottom = 0.003598395 Done calling nodal solver Timestep 13 ends with TIME = 4.621579177e-07 DT = 8.496773359e-08 Timing summary: Advection :0.737225129 seconds MAC Proj :0.183035944 seconds Nodal Proj :0.334388189 seconds Thermal :0.184797149 seconds Reactions :1.900050013 seconds Misc :0.212349264 seconds Base State :5.60000899e-08 seconds Time to advance time step: 3.556165538 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.083906109 Iter = 0.078568433 Bottom = 0.003833579 <<< 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.035180205 Iter = 0.033342069 Bottom = 0.000212244 <<< 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.084409234 Iter = 0.079064686 Bottom = 0.003976814 <<< 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.035519692 Iter = 0.033683471 Bottom = 0.000210944 <<< 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.218248934 Iter = 0.210945427 Bottom = 0.003563441 Done calling nodal solver Timestep 14 ends with TIME = 5.577754794e-07 DT = 9.561756174e-08 Timing summary: Advection :0.739635724 seconds MAC Proj :0.182259649 seconds Nodal Proj :0.335629681 seconds Thermal :0.184628671 seconds Reactions :1.912879597 seconds Misc :0.213297833 seconds Base State :4.80001745e-08 seconds Time to advance time step: 3.57263395 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.084429724 Iter = 0.079035009 Bottom = 0.003971379 <<< 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.03568855 Iter = 0.033858204 Bottom = 0.000206194 <<< 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.084606607 Iter = 0.079251065 Bottom = 0.003936155 <<< 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.035474659 Iter = 0.033629388 Bottom = 0.000211185 <<< 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.217839099 Iter = 0.210422706 Bottom = 0.003522787 Done calling nodal solver Timestep 15 ends with TIME = 6.534493518e-07 DT = 9.567387237e-08 Timing summary: Advection :0.742034457 seconds MAC Proj :0.183056268 seconds Nodal Proj :0.334385152 seconds Thermal :0.184328341 seconds Reactions :1.915884254 seconds Misc :0.214064271 seconds Base State :4.80001745e-08 seconds Time to advance time step: 3.57804355 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.08433195 Iter = 0.078947764 Bottom = 0.003988717 <<< 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.035523262 Iter = 0.03367435 Bottom = 0.000214575 <<< 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.084377735 Iter = 0.079038144 Bottom = 0.003920052 <<< 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.035638291 Iter = 0.033819146 Bottom = 0.00021705 <<< 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.218226524 Iter = 0.210883544 Bottom = 0.00358223 Done calling nodal solver Timestep 16 ends with TIME = 7.491599174e-07 DT = 9.571056553e-08 Timing summary: Advection :0.7433453 seconds MAC Proj :0.182669125 seconds Nodal Proj :0.334642944 seconds Thermal :0.184680547 seconds Reactions :1.901547715 seconds Misc :0.21304047 seconds Base State :5.199990483e-08 seconds Time to advance time step: 3.56431425 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.084384985 Iter = 0.079060512 Bottom = 0.003989342 <<< 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.035528822 Iter = 0.033659481 Bottom = 0.000212315 <<< 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.084549575 Iter = 0.079210147 Bottom = 0.004092615 <<< 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.03556374 Iter = 0.033735098 Bottom = 0.000213609 <<< 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.21778441 Iter = 0.210385364 Bottom = 0.003641365 Done calling nodal solver Timestep 17 ends with TIME = 8.449182257e-07 DT = 9.575830833e-08 Timing summary: Advection :0.746539762 seconds MAC Proj :0.182863914 seconds Nodal Proj :0.334481003 seconds Thermal :0.184680879 seconds Reactions :1.902898393 seconds Misc :0.214467432 seconds Base State :4.900016393e-08 seconds Time to advance time step: 3.570227808 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.084102717 Iter = 0.078725396 Bottom = 0.003648466 <<< 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.035536458 Iter = 0.033732975 Bottom = 0.000213918 <<< 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.084347479 Iter = 0.078996907 Bottom = 0.003997202 <<< 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.035417784 Iter = 0.033588256 Bottom = 0.000211821 <<< 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.217856874 Iter = 0.210544808 Bottom = 0.003507218 Done calling nodal solver Timestep 18 ends with TIME = 9.407071663e-07 DT = 9.578894066e-08 Timing summary: Advection :0.749388068 seconds MAC Proj :0.182440484 seconds Nodal Proj :0.335772199 seconds Thermal :0.184832853 seconds Reactions :1.898459048 seconds Misc :0.21380626 seconds Base State :5.500010047e-08 seconds Time to advance time step: 3.568952914 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.084623169 Iter = 0.079265914 Bottom = 0.003960244 <<< 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.035659656 Iter = 0.033831332 Bottom = 0.000214026 <<< 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.086008142 Iter = 0.080555527 Bottom = 0.004137039 <<< 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.035571234 Iter = 0.033702463 Bottom = 0.000215818 <<< 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.2176533 Iter = 0.21032761 Bottom = 0.003646909 Done calling nodal solver Timestep 19 ends with TIME = 1.036537244e-06 DT = 9.583007799e-08 Timing summary: Advection :0.75297355 seconds MAC Proj :0.184702286 seconds Nodal Proj :0.334718894 seconds Thermal :0.18459749 seconds Reactions :1.901343992 seconds Misc :0.21419651 seconds Base State :5.700007932e-08 seconds Time to advance time step: 3.576819749 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.084620142 Iter = 0.079286186 Bottom = 0.004085714 <<< 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.035456595 Iter = 0.033624415 Bottom = 0.000211925 <<< 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.084596814 Iter = 0.079245434 Bottom = 0.004122474 <<< 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.035712991 Iter = 0.033863819 Bottom = 0.000211072 <<< 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.217611732 Iter = 0.210258213 Bottom = 0.00366655 Done calling nodal solver Timestep 20 ends with TIME = 1.132393714e-06 DT = 9.585646998e-08 Timing summary: Advection :0.753941335 seconds MAC Proj :0.183185161 seconds Nodal Proj :0.334181724 seconds Thermal :0.184842887 seconds Reactions :1.892108626 seconds Misc :0.214408579 seconds Base State :4.500020623e-08 seconds Time to advance time step: 3.567044597 Writing plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.289632789 Writing small plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.319935787 Total Time: 84.11535523 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.10-21-g152fdac6aa4c) finalized