MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.01-6-g5726371bbb1a) initialized Calling Setup() Calling ReadParameters() WARNING: burning_cutoff_density_lo not supplied in the inputs file WARNING: setting burning_cutoff_density_lo = base_cutoff_density reading extern runtime parameters ... Calling VariableSetup() Calling BCSetup() Calling BaseStateGeometry::Init() Calling Init() Calling InitData() initdata model_File = model_files/kepler_new_6.25e8.hybrid.hse.320 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.106308222 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.257027123 Iter = 0.248369366 Bottom = 0.003920813 Done calling nodal solver Writing plotfile flame_3d_pltafter_InitProj after InitProj ... Level 0 create thermal coeffs: Time to write plotfile: 1.081539631 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.991417 MLMG: Initial residual (resid0) = 1745.991417 MLMG: Final Iter. 5 resid, resid/bnorm = 4.0738445e-05, 2.333255743e-08 MLMG: Timers: Solve = 0.186001107 Iter = 0.177422143 Bottom = 0.002898683 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.638951203e-05 MLMG: Initial residual (resid0) = 8.638951203e-05 MLMG: Final Iter. 6 resid, resid/bnorm = 2.069611843e-14, 2.395674885e-10 MLMG: Timers: Solve = 0.217568989 Iter = 0.210006572 Bottom = 0.003469978 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.038642892e-06 MLMG: Initial residual (resid0) = 2.038642892e-06 MLMG: Final Iter. 6 resid, resid/bnorm = 5.303485041e-19, 2.6014782e-13 MLMG: Timers: Solve = 0.216487316 Iter = 0.208998899 Bottom = 0.003453561 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.147662878 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.663632 MLMG: Initial residual (resid0) = 3970.663632 MLMG: Final Iter. 8 resid, resid/bnorm = 3.293810309e-07, 8.295364741e-11 MLMG: Timers: Solve = 0.077280508 Iter = 0.071254643 Bottom = 0.003097514 <<< 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.43167887e+12, 1.312381812e-13 MLMG: Timers: Solve = 0.028048446 Iter = 0.025800725 Bottom = 0.000171918 <<< 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.731772 MLMG: Initial residual (resid0) = 237.2771332 MLMG: Final Iter. 8 resid, resid/bnorm = 2.836048305e-07, 7.599818203e-11 MLMG: Timers: Solve = 0.092890344 Iter = 0.089315242 Bottom = 0.003363768 <<< 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.380139262e+12, 1.292671447e-13 MLMG: Timers: Solve = 0.028529823 Iter = 0.026257311 Bottom = 0.000179125 <<< 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.214694713e+10 MLMG: Initial residual (resid0) = 3.214694713e+10 MLMG: Final Iter. 6 resid, resid/bnorm = 1.109154942, 3.450265237e-11 MLMG: Timers: Solve = 0.220569005 Iter = 0.213046964 Bottom = 0.003551189 Done calling nodal solver Timestep 0 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.7766809464 seconds MAC Proj :0.1825590134 seconds Nodal Proj :0.3470966816 seconds Thermal :0.194257021 seconds Reactions :1.58944273 seconds Misc :0.2254259586 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.394398689 Writing plotfile 0 after all initialization ... Level 0 create thermal coeffs: Time to write plotfile: 1.113813639 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.434819 MLMG: Initial residual (resid0) = 3872.434819 MLMG: Final Iter. 8 resid, resid/bnorm = 1.282558912e-07, 3.312021951e-11 MLMG: Timers: Solve = 0.075916209 Iter = 0.072377648 Bottom = 0.003165808 <<< 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.028030321 Iter = 0.025787141 Bottom = 0.000177835 <<< 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.06826 MLMG: Initial residual (resid0) = 1.244370757 MLMG: Final Iter. 6 resid, resid/bnorm = 8.885353964e-08, 2.294730714e-11 MLMG: Timers: Solve = 0.05703748 Iter = 0.053436252 Bottom = 0.002513597 <<< 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.028346906 Iter = 0.026062394 Bottom = 0.00017996 <<< 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.3963897 MLMG: Initial residual (resid0) = 306.3963897 MLMG: Final Iter. 6 resid, resid/bnorm = 1.055558485e-08, 3.445074813e-11 MLMG: Timers: Solve = 0.228212714 Iter = 0.220667922 Bottom = 0.003519144 Done calling nodal solver Timestep 1 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.7395071983 seconds MAC Proj :0.1437549591 seconds Nodal Proj :0.3533263206 seconds Thermal :0.1918227673 seconds Reactions :1.584507227 seconds Misc :0.2034182549 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.246881008 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.85173 MLMG: Initial residual (resid0) = 3462.85173 MLMG: Final Iter. 8 resid, resid/bnorm = 1.161287031e-07, 3.353556899e-11 MLMG: Timers: Solve = 0.074180575 Iter = 0.070628097 Bottom = 0.003186601 <<< 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.02785482 Iter = 0.0256041 Bottom = 0.000181208 <<< 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.14453 MLMG: Initial residual (resid0) = 9.284142909 MLMG: Final Iter. 7 resid, resid/bnorm = 1.288626009e-07, 3.715606422e-11 MLMG: Timers: Solve = 0.067147953 Iter = 0.063571767 Bottom = 0.002964224 <<< 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.028564866 Iter = 0.026269895 Bottom = 0.000187084 <<< 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.9931611 MLMG: Initial residual (resid0) = 363.9931611 MLMG: Final Iter. 6 resid, resid/bnorm = 1.288980389e-08, 3.541221449e-11 MLMG: Timers: Solve = 0.220610881 Iter = 0.213034839 Bottom = 0.003502601 Done calling nodal solver Timestep 2 ends with TIME = 2.096533288e-08 DT = 1.143563612e-08 Timing summary: Advection :0.7470943928 seconds MAC Proj :0.1519978046 seconds Nodal Proj :0.3482625484 seconds Thermal :0.1944479942 seconds Reactions :1.584976673 seconds Misc :0.2063555717 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.233490467 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.940397 MLMG: Initial residual (resid0) = 2948.940397 MLMG: Final Iter. 8 resid, resid/bnorm = 1.206280018e-07, 4.090554083e-11 MLMG: Timers: Solve = 0.07606471 Iter = 0.07253394 Bottom = 0.003190976 <<< 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.028597408 Iter = 0.026334478 Bottom = 0.000184084 <<< 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.725486 MLMG: Initial residual (resid0) = 13.11339046 MLMG: Final Iter. 7 resid, resid/bnorm = 1.512367191e-07, 5.115007121e-11 MLMG: Timers: Solve = 0.067337703 Iter = 0.063648309 Bottom = 0.002926515 <<< 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.027908654 Iter = 0.025644099 Bottom = 0.00018225 <<< 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.6231902 MLMG: Initial residual (resid0) = 432.6231902 MLMG: Final Iter. 6 resid, resid/bnorm = 1.510207692e-08, 3.490815394e-11 MLMG: Timers: Solve = 0.214987932 Iter = 0.20747685 Bottom = 0.003481895 Done calling nodal solver Timestep 3 ends with TIME = 3.468809622e-08 DT = 1.372276334e-08 Timing summary: Advection :0.7622082233 seconds MAC Proj :0.154333353 seconds Nodal Proj :0.3392760754 seconds Thermal :0.1919736862 seconds Reactions :1.590390921 seconds Misc :0.206744194 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.245291471 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.708045 MLMG: Initial residual (resid0) = 2336.708045 MLMG: Final Iter. 8 resid, resid/bnorm = 1.25987583e-07, 5.391669845e-11 MLMG: Timers: Solve = 0.074350951 Iter = 0.070836682 Bottom = 0.003205019 <<< 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.786775708e+12, 2.977907667e-13 MLMG: Timers: Solve = 0.027897986 Iter = 0.025655225 Bottom = 0.000185084 <<< 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.152088 MLMG: Initial residual (resid0) = 18.25740453 MLMG: Final Iter. 7 resid, resid/bnorm = 2.131578754e-07, 9.077686086e-11 MLMG: Timers: Solve = 0.067691914 Iter = 0.064122853 Bottom = 0.002808473 <<< 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.229304009e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 7.722351198e+12, 2.95326971e-13 MLMG: Timers: Solve = 0.028321614 Iter = 0.026073602 Bottom = 0.000193959 <<< 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.5767953 MLMG: Initial residual (resid0) = 512.5767953 MLMG: Final Iter. 6 resid, resid/bnorm = 1.761854129e-08, 3.437249102e-11 MLMG: Timers: Solve = 0.217019486 Iter = 0.209273818 Bottom = 0.00345906 Done calling nodal solver Timestep 4 ends with TIME = 5.115541223e-08 DT = 1.646731601e-08 Timing summary: Advection :0.7464287281 seconds MAC Proj :0.1527056694 seconds Nodal Proj :0.3417584896 seconds Thermal :0.1917734146 seconds Reactions :1.589448214 seconds Misc :0.2056787014 seconds Base State :1.740455627e-05 seconds Time to advance time step: 3.228127956 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.220798 MLMG: Initial residual (resid0) = 1668.220798 MLMG: Final Iter. 8 resid, resid/bnorm = 1.624680408e-07, 9.739001037e-11 MLMG: Timers: Solve = 0.075117872 Iter = 0.071559102 Bottom = 0.003153226 <<< 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.028286947 Iter = 0.026020393 Bottom = 0.00016946 <<< 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.94004 MLMG: Initial residual (resid0) = 71.11602107 MLMG: Final Iter. 8 resid, resid/bnorm = 3.196611884e-08, 1.89716655e-11 MLMG: Timers: Solve = 0.075072913 Iter = 0.071512852 Bottom = 0.003238268 <<< 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.028775575 Iter = 0.026508312 Bottom = 0.00017196 <<< 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.7595069 MLMG: Initial residual (resid0) = 604.7595069 MLMG: Final Iter. 6 resid, resid/bnorm = 2.191978909e-08, 3.624546426e-11 MLMG: Timers: Solve = 0.217008236 Iter = 0.209446194 Bottom = 0.003382811 Done calling nodal solver Timestep 5 ends with TIME = 7.091619144e-08 DT = 1.976077921e-08 Timing summary: Advection :0.7547104359 seconds MAC Proj :0.1608681679 seconds Nodal Proj :0.3418858051 seconds Thermal :0.1924395561 seconds Reactions :1.625409365 seconds Misc :0.206154108 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.281873465 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.5205218 MLMG: Initial residual (resid0) = 925.5205218 MLMG: Final Iter. 9 resid, resid/bnorm = 1.655848791e-08, 1.789100028e-11 MLMG: Timers: Solve = 0.083839461 Iter = 0.080335484 Bottom = 0.003607604 <<< 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.068788789e+12, 2.320895495e-13 MLMG: Timers: Solve = 0.028333197 Iter = 0.026099144 Bottom = 0.000178 <<< 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.6259165 MLMG: Initial residual (resid0) = 71.67495479 MLMG: Final Iter. 8 resid, resid/bnorm = 4.417317536e-08, 4.651639618e-11 MLMG: Timers: Solve = 0.07537004 Iter = 0.07179877 Bottom = 0.003346976 <<< 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.704973329e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.141803233e+12, 2.348818513e-13 MLMG: Timers: Solve = 0.028558866 Iter = 0.026249436 Bottom = 0.000176751 <<< 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.0822892 MLMG: Initial residual (resid0) = 711.0822892 MLMG: Final Iter. 6 resid, resid/bnorm = 2.67066993e-08, 3.755781814e-11 MLMG: Timers: Solve = 0.21761428 Iter = 0.210027197 Bottom = 0.003351894 Done calling nodal solver Timestep 6 ends with TIME = 9.462912649e-08 DT = 2.371293505e-08 Timing summary: Advection :0.7581651211 seconds MAC Proj :0.1698744297 seconds Nodal Proj :0.3455440998 seconds Thermal :0.1924474239 seconds Reactions :1.664015293 seconds Misc :0.2093689442 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.339787006 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.9925939 MLMG: Initial residual (resid0) = 676.9925939 MLMG: Final Iter. 9 resid, resid/bnorm = 2.321343118e-08, 3.428904745e-11 MLMG: Timers: Solve = 0.082066535 Iter = 0.078447723 Bottom = 0.003538727 <<< 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.027573985 Iter = 0.025508224 Bottom = 0.000192626 <<< 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.7584038 MLMG: Initial residual (resid0) = 44.05710321 MLMG: Final Iter. 8 resid, resid/bnorm = 3.632598577e-08, 5.251253266e-11 MLMG: Timers: Solve = 0.07351703 Iter = 0.069933552 Bottom = 0.0032116 <<< 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.027451859 Iter = 0.025374932 Bottom = 0.000164209 <<< 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.0677451 MLMG: Initial residual (resid0) = 827.0677451 MLMG: Final Iter. 6 resid, resid/bnorm = 3.071783627e-08, 3.714065318e-11 MLMG: Timers: Solve = 0.215780312 Iter = 0.208402938 Bottom = 0.003218142 Done calling nodal solver Timestep 7 ends with TIME = 1.230846486e-07 DT = 2.845552206e-08 Timing summary: Advection :0.7570135593 seconds MAC Proj :0.1662888527 seconds Nodal Proj :0.3377218246 seconds Thermal :0.1786711216 seconds Reactions :1.67113018 seconds Misc :0.2108981609 seconds Base State :2.861022949e-06 seconds Time to advance time step: 3.322097778 Call to estdt for level 0 gives dt_lev = 9.539807156e-08 Minimum estdt over all levels = 9.539807156e-08 Call to estdt at beginning of step 8 gives dt =9.539807156e-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.3116109 MLMG: Initial residual (resid0) = 731.3116109 MLMG: Final Iter. 9 resid, resid/bnorm = 2.958040568e-08, 4.044842888e-11 MLMG: Timers: Solve = 0.082599788 Iter = 0.079007518 Bottom = 0.003670271 <<< 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.027533318 Iter = 0.02546864 Bottom = 0.000172294 <<< 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.48027 MLMG: Initial residual (resid0) = 56.29847793 MLMG: Final Iter. 8 resid, resid/bnorm = 5.021658467e-08, 6.664618394e-11 MLMG: Timers: Solve = 0.073556696 Iter = 0.069972218 Bottom = 0.003222726 <<< 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.02739515 Iter = 0.025313223 Bottom = 0.000171084 <<< 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.08661 MLMG: Initial residual (resid0) = 940.08661 MLMG: Final Iter. 6 resid, resid/bnorm = 3.494798517e-08, 3.717528236e-11 MLMG: Timers: Solve = 0.216688817 Iter = 0.209272194 Bottom = 0.003429852 Done calling nodal solver Timestep 8 ends with TIME = 1.57231275e-07 DT = 3.414662647e-08 Timing summary: Advection :0.7677931786 seconds MAC Proj :0.1669569016 seconds Nodal Proj :0.3389542103 seconds Thermal :0.1805462837 seconds Reactions :1.711575508 seconds Misc :0.2191689014 seconds Base State :2.861022949e-06 seconds Time to advance time step: 3.385378838 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.755157 MLMG: Initial residual (resid0) = 2219.755157 MLMG: Final Iter. 9 resid, resid/bnorm = 3.32684067e-08, 1.498742174e-11 MLMG: Timers: Solve = 0.081980868 Iter = 0.078437348 Bottom = 0.003468396 <<< 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.027574026 Iter = 0.025470223 Bottom = 0.000184043 <<< 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.321845 MLMG: Initial residual (resid0) = 70.26043754 MLMG: Final Iter. 8 resid, resid/bnorm = 6.901615279e-08, 3.062862633e-11 MLMG: Timers: Solve = 0.073632572 Iter = 0.070066261 Bottom = 0.003162728 <<< 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.027533235 Iter = 0.025480432 Bottom = 0.000172501 <<< 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.509423 MLMG: Initial residual (resid0) = 1054.509423 MLMG: Final Iter. 6 resid, resid/bnorm = 3.930475967e-08, 3.72730284e-11 MLMG: Timers: Solve = 0.215756854 Iter = 0.208427896 Bottom = 0.003292018 Done calling nodal solver Timestep 9 ends with TIME = 1.982072268e-07 DT = 4.097595177e-08 Timing summary: Advection :0.7617650032 seconds MAC Proj :0.1663961411 seconds Nodal Proj :0.3374512196 seconds Thermal :0.1788520813 seconds Reactions :1.715215921 seconds Misc :0.2177340984 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.37780261 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.674837 MLMG: Initial residual (resid0) = 2063.674837 MLMG: Final Iter. 9 resid, resid/bnorm = 3.296190698e-08, 1.597243247e-11 MLMG: Timers: Solve = 0.082148036 Iter = 0.078559099 Bottom = 0.003555811 <<< 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.767579968e+13, 2.588135196e-12 MLMG: Timers: Solve = 0.02748936 Iter = 0.025413681 Bottom = 0.000170626 <<< 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.623674 MLMG: Initial residual (resid0) = 84.89142822 MLMG: Final Iter. 8 resid, resid/bnorm = 9.424633163e-08, 4.452673036e-11 MLMG: Timers: Solve = 0.073679405 Iter = 0.070082719 Bottom = 0.003146183 <<< 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.222338959e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.628423028e+13, 2.53491721e-12 MLMG: Timers: Solve = 0.028023029 Iter = 0.025909976 Bottom = 0.00017446 <<< 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.461632 MLMG: Initial residual (resid0) = 1165.461632 MLMG: Final Iter. 6 resid, resid/bnorm = 4.370407658e-08, 3.749936968e-11 MLMG: Timers: Solve = 0.215980355 Iter = 0.208539023 Bottom = 0.003522728 Done calling nodal solver Timestep 10 ends with TIME = 2.473783689e-07 DT = 4.917114212e-08 Timing summary: Advection :0.768153429 seconds MAC Proj :0.1665921211 seconds Nodal Proj :0.3381297588 seconds Thermal :0.180562973 seconds Reactions :1.749654055 seconds Misc :0.2235937119 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.427066326 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.111671276e-08, 3.08383697e-11 MLMG: Timers: Solve = 0.082477621 Iter = 0.078910184 Bottom = 0.003571146 <<< 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.027589818 Iter = 0.025488765 Bottom = 0.000172043 <<< 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.1554523 MLMG: Initial residual (resid0) = 131.1492543 MLMG: Final Iter. 9 resid, resid/bnorm = 8.210179203e-09, 8.46274603e-12 MLMG: Timers: Solve = 0.082836748 Iter = 0.079229602 Bottom = 0.003577184 <<< 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.788690556e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 6.536081231e+13, 2.499602806e-12 MLMG: Timers: Solve = 0.027465401 Iter = 0.025385682 Bottom = 0.000170085 <<< 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.054354 MLMG: Initial residual (resid0) = 1269.054354 MLMG: Final Iter. 6 resid, resid/bnorm = 4.938119447e-08, 3.891180414e-11 MLMG: Timers: Solve = 0.215501644 Iter = 0.208116228 Bottom = 0.003414476 Done calling nodal solver Timestep 11 ends with TIME = 3.063837395e-07 DT = 5.900537055e-08 Timing summary: Advection :0.7615644932 seconds MAC Proj :0.1760649681 seconds Nodal Proj :0.3378326893 seconds Thermal :0.1806092262 seconds Reactions :1.756725788 seconds Misc :0.2233850956 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.436552525 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.053401088e-08, 3.410424345e-11 MLMG: Timers: Solve = 0.082170036 Iter = 0.078612349 Bottom = 0.003620894 <<< 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.484967197e+22 MLMG: Final Iter. 3 resid, resid/bnorm = 1.487768081e+15, 5.689692553e-11 MLMG: Timers: Solve = 0.027440026 Iter = 0.025383681 Bottom = 0.000170918 <<< 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.7783863 MLMG: Initial residual (resid0) = 482.3381305 MLMG: Final Iter. 9 resid, resid/resid0 = 1.908621709e-08, 3.957020165e-11 MLMG: Timers: Solve = 0.082489329 Iter = 0.078914101 Bottom = 0.003570227 <<< 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.401151476e+15, 5.358443442e-11 MLMG: Timers: Solve = 0.027655278 Iter = 0.025556391 Bottom = 0.000173126 <<< 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.874686781e-08, 4.478515523e-11 MLMG: Timers: Solve = 0.213519675 Iter = 0.206253969 Bottom = 0.003244434 Done calling nodal solver Timestep 12 ends with TIME = 3.771901841e-07 DT = 7.080644466e-08 Timing summary: Advection :0.7613258362 seconds MAC Proj :0.1754560471 seconds Nodal Proj :0.3355617523 seconds Thermal :0.1790952682 seconds Reactions :1.776664972 seconds Misc :0.2266287804 seconds Base State :1.645088196e-05 seconds Time to advance time step: 3.455123663 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.364482 MLMG: Initial residual (resid0) = 1121.364482 MLMG: Final Iter. 9 resid, resid/bnorm = 2.545914413e-08, 2.270371902e-11 MLMG: Timers: Solve = 0.081845367 Iter = 0.078270597 Bottom = 0.003570729 <<< 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.034271075e+14, 3.955377524e-12 MLMG: Timers: Solve = 0.035675529 Iter = 0.033604352 Bottom = 0.000227251 <<< 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.562884 MLMG: Initial residual (resid0) = 164.8126831 MLMG: Final Iter. 9 resid, resid/bnorm = 1.432982799e-08, 1.178863568e-11 MLMG: Timers: Solve = 0.08236462 Iter = 0.078784558 Bottom = 0.003614063 <<< 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.202860904e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 9.680856285e+13, 3.702263585e-12 MLMG: Timers: Solve = 0.03579353 Iter = 0.033741852 Bottom = 0.000234668 <<< 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.115761 MLMG: Initial residual (resid0) = 1498.115761 MLMG: Final Iter. 6 resid, resid/bnorm = 6.253442431e-08, 4.174205088e-11 MLMG: Timers: Solve = 0.213086714 Iter = 0.205844966 Bottom = 0.003377602 Done calling nodal solver Timestep 13 ends with TIME = 4.621579177e-07 DT = 8.496773359e-08 Timing summary: Advection :0.7603654861 seconds MAC Proj :0.1749470234 seconds Nodal Proj :0.3350656033 seconds Thermal :0.1952881813 seconds Reactions :1.767793179 seconds Misc :0.2274036407 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.46125102 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.983864533e-08, 2.301837563e-11 MLMG: Timers: Solve = 0.082042493 Iter = 0.078460807 Bottom = 0.003576853 <<< 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.035966781 Iter = 0.033887436 Bottom = 0.000224833 <<< 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.6896171 MLMG: Initial residual (resid0) = 418.9789588 MLMG: Final Iter. 9 resid, resid/bnorm = 1.858319365e-08, 2.129416151e-11 MLMG: Timers: Solve = 0.082127494 Iter = 0.078503182 Bottom = 0.003607727 <<< 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.036014865 Iter = 0.033906062 Bottom = 0.000221126 <<< 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.544988 MLMG: Initial residual (resid0) = 1661.544988 MLMG: Final Iter. 6 resid, resid/bnorm = 6.683421816e-08, 4.022413997e-11 MLMG: Timers: Solve = 0.216097022 Iter = 0.208789399 Bottom = 0.003499477 Done calling nodal solver Timestep 14 ends with TIME = 5.577754795e-07 DT = 9.561756174e-08 Timing summary: Advection :0.7616589069 seconds MAC Proj :0.1749064922 seconds Nodal Proj :0.3380355835 seconds Thermal :0.195785284 seconds Reactions :1.819138765 seconds Misc :0.230250597 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.520155191 Call to estdt for level 0 gives dt_lev = 9.567387238e-08 Minimum estdt over all levels = 9.567387238e-08 Call to estdt at beginning of step 15 gives dt =9.567387238e-08 Timestep 15 starts with TIME = 5.577754795e-07 DT = 9.567387238e-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.248646863e-08, 2.721445659e-11 MLMG: Timers: Solve = 0.082244744 Iter = 0.078664267 Bottom = 0.003666853 <<< 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.035886406 Iter = 0.033822394 Bottom = 0.000221792 <<< 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.6508291 MLMG: Initial residual (resid0) = 880.7089387 MLMG: Final Iter. 9 resid, resid/resid0 = 1.902481017e-08, 2.160169988e-11 MLMG: Timers: Solve = 0.082726997 Iter = 0.079109393 Bottom = 0.003638062 <<< 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.328089735e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.475141174e+14, 1.71143458e-11 MLMG: Timers: Solve = 0.035833614 Iter = 0.033749394 Bottom = 0.000222085 <<< 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.22281 MLMG: Initial residual (resid0) = 1355.22281 MLMG: Final Iter. 6 resid, resid/bnorm = 6.156096788e-08, 4.542497915e-11 MLMG: Timers: Solve = 0.216878568 Iter = 0.209546319 Bottom = 0.003402352 Done calling nodal solver Timestep 15 ends with TIME = 6.534493518e-07 DT = 9.567387238e-08 Timing summary: Advection :0.7640638351 seconds MAC Proj :0.1758437157 seconds Nodal Proj :0.3393120766 seconds Thermal :0.1953532696 seconds Reactions :1.826143265 seconds Misc :0.2314968109 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.532612085 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.7040973 MLMG: Initial residual (resid0) = 960.7040973 MLMG: Final Iter. 9 resid, resid/bnorm = 2.91497031e-08, 3.034202017e-11 MLMG: Timers: Solve = 0.081985409 Iter = 0.078412473 Bottom = 0.00362802 <<< 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.593814955e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.596087453e+14, 1.757688237e-11 MLMG: Timers: Solve = 0.037099246 Iter = 0.03396452 Bottom = 0.000222043 <<< 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.3825405 MLMG: Initial residual (resid0) = 374.842397 MLMG: Final Iter. 9 resid, resid/bnorm = 1.743441658e-08, 1.781963174e-11 MLMG: Timers: Solve = 0.082268244 Iter = 0.078717517 Bottom = 0.003571937 <<< 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.418981599e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.321252496e+14, 1.652582715e-11 MLMG: Timers: Solve = 0.035577821 Iter = 0.033616143 Bottom = 0.000224462 <<< 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.550859 MLMG: Initial residual (resid0) = 1231.550859 MLMG: Final Iter. 6 resid, resid/bnorm = 5.893471666e-08, 4.785406646e-11 MLMG: Timers: Solve = 0.215033849 Iter = 0.20787731 Bottom = 0.003594353 Done calling nodal solver Timestep 16 ends with TIME = 7.491599174e-07 DT = 9.571056553e-08 Timing summary: Advection :0.8116285801 seconds MAC Proj :0.1744456291 seconds Nodal Proj :0.3364741802 seconds Thermal :0.2096042633 seconds Reactions :1.802468061 seconds Misc :0.229685545 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.564703226 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.753794 MLMG: Initial residual (resid0) = 1812.753794 MLMG: Final Iter. 9 resid, resid/bnorm = 2.464378213e-08, 1.359466587e-11 MLMG: Timers: Solve = 0.081574865 Iter = 0.078545723 Bottom = 0.003659188 <<< 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.035576154 Iter = 0.033610143 Bottom = 0.000236833 <<< 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.4776298 MLMG: Initial residual (resid0) = 1230.328214 MLMG: Final Iter. 9 resid, resid/resid0 = 2.692156431e-08, 2.188161175e-11 MLMG: Timers: Solve = 0.082062577 Iter = 0.079031476 Bottom = 0.003698231 <<< 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.190657931e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.262411444e+14, 1.630080048e-11 MLMG: Timers: Solve = 0.035794405 Iter = 0.033835395 Bottom = 0.000224125 <<< 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.796809 MLMG: Initial residual (resid0) = 1403.796809 MLMG: Final Iter. 6 resid, resid/bnorm = 5.715173756e-08, 4.071225779e-11 MLMG: Timers: Solve = 0.215271643 Iter = 0.208138562 Bottom = 0.003496854 Done calling nodal solver Timestep 17 ends with TIME = 8.449182257e-07 DT = 9.575830833e-08 Timing summary: Advection :0.7383952141 seconds MAC Proj :0.1732833385 seconds Nodal Proj :0.3357894421 seconds Thermal :0.1932866573 seconds Reactions :1.805734634 seconds Misc :0.2305841446 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.477452517 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.425488788e-08, 3.642249019e-11 MLMG: Timers: Solve = 0.081780741 Iter = 0.078764058 Bottom = 0.003523309 <<< 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.161926329e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.411060262e+14, 1.686928026e-11 MLMG: Timers: Solve = 0.035548654 Iter = 0.033560184 Bottom = 0.000219793 <<< 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.9129471 MLMG: Initial residual (resid0) = 567.9651174 MLMG: Final Iter. 9 resid, resid/resid0 = 1.513583214e-08, 2.664922841e-11 MLMG: Timers: Solve = 0.081568783 Iter = 0.078548849 Bottom = 0.003613938 <<< 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.944630522e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.131801488e+14, 1.580130697e-11 MLMG: Timers: Solve = 0.035832989 Iter = 0.033795727 Bottom = 0.000223919 <<< 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.9820512 MLMG: Initial residual (resid0) = 890.9820512 MLMG: Final Iter. 6 resid, resid/bnorm = 5.488891475e-08, 6.160496126e-11 MLMG: Timers: Solve = 0.215731978 Iter = 0.208557272 Bottom = 0.003382309 Done calling nodal solver Timestep 18 ends with TIME = 9.407071664e-07 DT = 9.578894066e-08 Timing summary: Advection :0.7354259491 seconds MAC Proj :0.1728951931 seconds Nodal Proj :0.3364503384 seconds Thermal :0.1938662529 seconds Reactions :1.793754816 seconds Misc :0.2242386341 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.457000494 Call to estdt for level 0 gives dt_lev = 9.5830078e-08 Minimum estdt over all levels = 9.5830078e-08 Call to estdt at beginning of step 19 gives dt =9.5830078e-08 Timestep 19 starts with TIME = 9.407071664e-07 DT = 9.5830078e-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.569845492e-08, 1.736110925e-11 MLMG: Timers: Solve = 0.081877284 Iter = 0.078868726 Bottom = 0.003485769 <<< 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.435123349e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.23926157e+14, 1.621226809e-11 MLMG: Timers: Solve = 0.036174824 Iter = 0.034185063 Bottom = 0.000224459 <<< 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.3367993 MLMG: Initial residual (resid0) = 715.348036 MLMG: Final Iter. 9 resid, resid/bnorm = 1.394148796e-08, 1.836008472e-11 MLMG: Timers: Solve = 0.081843992 Iter = 0.078810684 Bottom = 0.003695021 <<< 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.064650232e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.009996216e+14, 1.533548534e-11 MLMG: Timers: Solve = 0.035679822 Iter = 0.033693269 Bottom = 0.000222293 <<< 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.87858 MLMG: Initial residual (resid0) = 1137.87858 MLMG: Final Iter. 6 resid, resid/bnorm = 5.008223525e-08, 4.401369015e-11 MLMG: Timers: Solve = 0.215696021 Iter = 0.208338396 Bottom = 0.003496687 Done calling nodal solver Timestep 19 ends with TIME = 1.036537244e-06 DT = 9.5830078e-08 Timing summary: Advection :0.7394814491 seconds MAC Proj :0.173438549 seconds Nodal Proj :0.3367857933 seconds Thermal :0.1953437328 seconds Reactions :1.80825901 seconds Misc :0.224793911 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.478467226 Call to estdt for level 0 gives dt_lev = 9.585646999e-08 Minimum estdt over all levels = 9.585646999e-08 Call to estdt at beginning of step 20 gives dt =9.585646999e-08 Timestep 20 starts with TIME = 1.036537244e-06 DT = 9.585646999e-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.324404136e-08, 1.715578164e-11 MLMG: Timers: Solve = 0.081973702 Iter = 0.078954184 Bottom = 0.003668228 <<< 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.345657733e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 4.181064763e+14, 1.598970521e-11 MLMG: Timers: Solve = 0.035693072 Iter = 0.033717061 Bottom = 0.000219169 <<< 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.733798 MLMG: Initial residual (resid0) = 1050.094681 MLMG: Final Iter. 9 resid, resid/resid0 = 1.293258833e-08, 1.231564027e-11 MLMG: Timers: Solve = 0.081665116 Iter = 0.078625766 Bottom = 0.00379606 <<< 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.727742203e+22 MLMG: Final Iter. 4 resid, resid/bnorm = 3.949437177e+14, 1.510388855e-11 MLMG: Timers: Solve = 0.036383492 Iter = 0.034411272 Bottom = 0.000240125 <<< 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.6772084 MLMG: Initial residual (resid0) = 933.6772084 MLMG: Final Iter. 6 resid, resid/bnorm = 4.948662991e-08, 5.300186131e-11 MLMG: Timers: Solve = 0.213123381 Iter = 0.206063634 Bottom = 0.003494935 Done calling nodal solver Timestep 20 ends with TIME = 1.132393714e-06 DT = 9.585646999e-08 Timing summary: Advection :0.7385442257 seconds MAC Proj :0.1732878685 seconds Nodal Proj :0.333381176 seconds Thermal :0.1945786476 seconds Reactions :1.776662111 seconds Misc :0.2238726616 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.440730095 Writing plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.305018902 Writing small plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.292569876 Total Time: 82.71059442 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 (22.01-6-g5726371bbb1a) finalized