MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.01-13-g0278677bbc8c) 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.107061625 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.25605866 Iter = 0.247804031 Bottom = 0.004202773 Done calling nodal solver Writing plotfile flame_3d_pltafter_InitProj after InitProj ... Level 0 create thermal coeffs: Time to write plotfile: 1.095941544 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.17993699 Iter = 0.172096739 Bottom = 0.0030846 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.219373291 Iter = 0.212420503 Bottom = 0.003797482 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.213961928 Iter = 0.207003306 Bottom = 0.003719562 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.084888697 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.074364868 Iter = 0.069235881 Bottom = 0.003727853 <<< 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.026497187 Iter = 0.024801136 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.071233642 Iter = 0.068617586 Bottom = 0.003851354 <<< 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.026304395 Iter = 0.024598761 Bottom = 0.000165417 <<< 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.21312684 Iter = 0.206217218 Bottom = 0.003711728 Done calling nodal solver Timestep 0 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.7163758278 seconds MAC Proj :0.156167984 seconds Nodal Proj :0.3300659657 seconds Thermal :0.1658394337 seconds Reactions :1.549788952 seconds Misc :0.2242267132 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.179795742 Writing plotfile 0 after all initialization ... Level 0 create thermal coeffs: Time to write plotfile: 1.088063955 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.071157308 Iter = 0.068540085 Bottom = 0.003616812 <<< 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.026313395 Iter = 0.024622761 Bottom = 0.000168668 <<< 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.054611759 Iter = 0.051980411 Bottom = 0.002945184 <<< 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.026602896 Iter = 0.02485397 Bottom = 0.000169543 <<< 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.213325757 Iter = 0.206219844 Bottom = 0.003715731 Done calling nodal solver Timestep 1 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.6773819923 seconds MAC Proj :0.1345148087 seconds Nodal Proj :0.3306117058 seconds Thermal :0.1666812897 seconds Reactions :1.543816566 seconds Misc :0.1934328079 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.081871033 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.071814187 Iter = 0.069204256 Bottom = 0.003661646 <<< 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.026501771 Iter = 0.024785595 Bottom = 0.000168917 <<< 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.063810059 Iter = 0.061164711 Bottom = 0.003383395 <<< 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.026365728 Iter = 0.024654803 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 = 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.213318549 Iter = 0.206376719 Bottom = 0.003708105 Done calling nodal solver Timestep 2 ends with TIME = 2.096533288e-08 DT = 1.143563612e-08 Timing summary: Advection :0.6809980869 seconds MAC Proj :0.1444804668 seconds Nodal Proj :0.3304507732 seconds Thermal :0.1661198139 seconds Reactions :1.553560019 seconds Misc :0.1941268444 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.070109606 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.071355643 Iter = 0.06871667 Bottom = 0.00368448 <<< 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.026566938 Iter = 0.024861428 Bottom = 0.000177834 <<< 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.063144014 Iter = 0.06055375 Bottom = 0.003371852 <<< 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.026548688 Iter = 0.024837387 Bottom = 0.000186335 <<< 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.213090465 Iter = 0.205984134 Bottom = 0.003761522 Done calling nodal solver Timestep 3 ends with TIME = 3.468809622e-08 DT = 1.372276334e-08 Timing summary: Advection :0.6838884354 seconds MAC Proj :0.1432425976 seconds Nodal Proj :0.3308656216 seconds Thermal :0.1665296555 seconds Reactions :1.557762384 seconds Misc :0.1949930191 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.077648401 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.071962812 Iter = 0.069207006 Bottom = 0.003930396 <<< 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.026307145 Iter = 0.024578968 Bottom = 0.000174292 <<< 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.063447266 Iter = 0.060781959 Bottom = 0.00327506 <<< 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.026533063 Iter = 0.024821679 Bottom = 0.000172793 <<< 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.212961881 Iter = 0.205924842 Bottom = 0.003606603 Done calling nodal solver Timestep 4 ends with TIME = 5.115541223e-08 DT = 1.646731601e-08 Timing summary: Advection :0.6890413761 seconds MAC Proj :0.1440997124 seconds Nodal Proj :0.3301017284 seconds Thermal :0.1660766602 seconds Reactions :1.571270704 seconds Misc :0.1959123611 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.096881151 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.071958396 Iter = 0.06932359 Bottom = 0.003742604 <<< 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.026666897 Iter = 0.024912637 Bottom = 0.000164167 <<< 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.072013813 Iter = 0.069287464 Bottom = 0.003821354 <<< 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.026650438 Iter = 0.024915012 Bottom = 0.000167377 <<< 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.214026553 Iter = 0.207012222 Bottom = 0.003515353 Done calling nodal solver Timestep 5 ends with TIME = 7.091619144e-08 DT = 1.976077921e-08 Timing summary: Advection :0.6993534565 seconds MAC Proj :0.1528606415 seconds Nodal Proj :0.3310964108 seconds Thermal :0.1667714119 seconds Reactions :1.594016314 seconds Misc :0.1964857578 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.140958309 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.080370193 Iter = 0.077787428 Bottom = 0.00420748 <<< 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.026717939 Iter = 0.02492697 Bottom = 0.000189626 <<< 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.071507894 Iter = 0.068829629 Bottom = 0.003850439 <<< 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.026366317 Iter = 0.024653683 Bottom = 0.000177668 <<< 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.213076967 Iter = 0.20612372 Bottom = 0.003671564 Done calling nodal solver Timestep 6 ends with TIME = 9.462912649e-08 DT = 2.371293505e-08 Timing summary: Advection :0.6996078491 seconds MAC Proj :0.1606543064 seconds Nodal Proj :0.3301625252 seconds Thermal :0.1661307812 seconds Reactions :1.611179829 seconds Misc :0.197527647 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.165629864 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.081250339 Iter = 0.07862449 Bottom = 0.004133484 <<< 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.026402609 Iter = 0.024680849 Bottom = 0.000171168 <<< 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.072524956 Iter = 0.069904816 Bottom = 0.003805939 <<< 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.026576777 Iter = 0.024879684 Bottom = 0.000159626 <<< 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.213033342 Iter = 0.206079886 Bottom = 0.003500771 Done calling nodal solver Timestep 7 ends with TIME = 1.230846486e-07 DT = 2.845552206e-08 Timing summary: Advection :0.7077946663 seconds MAC Proj :0.1625165939 seconds Nodal Proj :0.3296208382 seconds Thermal :0.1662330627 seconds Reactions :1.645870447 seconds Misc :0.2018361092 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.214240074 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.080049833 Iter = 0.077450651 Bottom = 0.0042919 <<< 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.026747987 Iter = 0.024993643 Bottom = 0.000175376 <<< 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.071379783 Iter = 0.068759643 Bottom = 0.003841147 <<< 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.026474568 Iter = 0.024765433 Bottom = 0.000168625 <<< 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.214340058 Iter = 0.207279601 Bottom = 0.003559354 Done calling nodal solver Timestep 8 ends with TIME = 1.57231275e-07 DT = 3.414662647e-08 Timing summary: Advection :0.7094581127 seconds MAC Proj :0.1602227688 seconds Nodal Proj :0.3309831619 seconds Thermal :0.1667823792 seconds Reactions :1.660230637 seconds Misc :0.2034132481 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.231459856 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.081708842 Iter = 0.079004701 Bottom = 0.004139442 <<< 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.026755444 Iter = 0.025047851 Bottom = 0.000173669 <<< 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.07116849 Iter = 0.068540266 Bottom = 0.00369248 <<< 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.02653161 Iter = 0.024775224 Bottom = 0.000171001 <<< 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.212735466 Iter = 0.205753009 Bottom = 0.003695313 Done calling nodal solver Timestep 9 ends with TIME = 1.982072268e-07 DT = 4.097595177e-08 Timing summary: Advection :0.7153139114 seconds MAC Proj :0.1617555618 seconds Nodal Proj :0.3299438953 seconds Thermal :0.1668922901 seconds Reactions :1.698982954 seconds Misc :0.207034111 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.280299902 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.08058971 Iter = 0.077998862 Bottom = 0.004211233 <<< 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.026442359 Iter = 0.024757433 Bottom = 0.000174417 <<< 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.072118037 Iter = 0.069508606 Bottom = 0.003739563 <<< 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.026595069 Iter = 0.024852017 Bottom = 0.000175376 <<< 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.21700024 Iter = 0.209896866 Bottom = 0.003733979 Done calling nodal solver Timestep 10 ends with TIME = 2.473783689e-07 DT = 4.917114212e-08 Timing summary: Advection :0.7104008198 seconds MAC Proj :0.1614353657 seconds Nodal Proj :0.334246397 seconds Thermal :0.1659011841 seconds Reactions :1.692052364 seconds Misc :0.208565712 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.27298069 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.082021094 Iter = 0.079296328 Bottom = 0.004223192 <<< 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.026275817 Iter = 0.024566307 Bottom = 0.000172085 <<< 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.080049 Iter = 0.077424193 Bottom = 0.004219357 <<< 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.026358734 Iter = 0.024656682 Bottom = 0.000167667 <<< 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.213311635 Iter = 0.206319012 Bottom = 0.003718563 Done calling nodal solver Timestep 11 ends with TIME = 3.063837395e-07 DT = 5.900537055e-08 Timing summary: Advection :0.7111592293 seconds MAC Proj :0.1708936691 seconds Nodal Proj :0.329605341 seconds Thermal :0.1662597656 seconds Reactions :1.726373672 seconds Misc :0.2124712467 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.317129374 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.080301043 Iter = 0.077685069 Bottom = 0.004225236 <<< 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.026767778 Iter = 0.02502856 Bottom = 0.000174208 <<< 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.080614044 Iter = 0.078008362 Bottom = 0.004170858 <<< 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.026550651 Iter = 0.024847517 Bottom = 0.000170418 <<< 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.213577095 Iter = 0.206545847 Bottom = 0.003451312 Done calling nodal solver Timestep 12 ends with TIME = 3.771901841e-07 DT = 7.080644466e-08 Timing summary: Advection :0.7099220753 seconds MAC Proj :0.1696152687 seconds Nodal Proj :0.3312602043 seconds Thermal :0.1671881676 seconds Reactions :1.746636152 seconds Misc :0.2156496048 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.340646267 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.08076667 Iter = 0.078164905 Bottom = 0.004279942 <<< 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.035171117 Iter = 0.033480941 Bottom = 0.000228293 <<< 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.080717044 Iter = 0.078063779 Bottom = 0.004307108 <<< 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.034582365 Iter = 0.032874355 Bottom = 0.000228293 <<< 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.213222426 Iter = 0.206278804 Bottom = 0.003745813 Done calling nodal solver Timestep 13 ends with TIME = 4.621579177e-07 DT = 8.496773359e-08 Timing summary: Advection :0.7069494724 seconds MAC Proj :0.1701929569 seconds Nodal Proj :0.3298244476 seconds Thermal :0.1824786663 seconds Reactions :1.769278765 seconds Misc :0.2166447639 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.375750542 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.080081665 Iter = 0.077466734 Bottom = 0.004249024 <<< 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.034703073 Iter = 0.03298573 Bottom = 0.000222877 <<< 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.084061564 Iter = 0.081428674 Bottom = 0.004285026 <<< 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.03485699 Iter = 0.03314469 Bottom = 0.000225251 <<< 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.213830139 Iter = 0.206838515 Bottom = 0.003699938 Done calling nodal solver Timestep 14 ends with TIME = 5.577754795e-07 DT = 9.561756174e-08 Timing summary: Advection :0.7086145878 seconds MAC Proj :0.1728856564 seconds Nodal Proj :0.3311359882 seconds Thermal :0.1830284595 seconds Reactions :1.781971693 seconds Misc :0.2190141678 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.397019863 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.080567794 Iter = 0.077934404 Bottom = 0.004295816 <<< 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.034857241 Iter = 0.033152106 Bottom = 0.000243376 <<< 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.080191083 Iter = 0.077571027 Bottom = 0.004255775 <<< 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.03443253 Iter = 0.032709604 Bottom = 0.000238335 <<< 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.21263334 Iter = 0.205679384 Bottom = 0.003632854 Done calling nodal solver Timestep 15 ends with TIME = 6.534493518e-07 DT = 9.567387238e-08 Timing summary: Advection :0.7102925777 seconds MAC Proj :0.169479847 seconds Nodal Proj :0.3298053741 seconds Thermal :0.1832511425 seconds Reactions :1.785492659 seconds Misc :0.2192845345 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.397992134 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.079920456 Iter = 0.077314192 Bottom = 0.004289649 <<< 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.034311613 Iter = 0.032613061 Bottom = 0.000221167 <<< 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.080414002 Iter = 0.077787945 Bottom = 0.004287483 <<< 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.035065993 Iter = 0.033320399 Bottom = 0.000239252 <<< 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.216505112 Iter = 0.209256529 Bottom = 0.003846689 Done calling nodal solver Timestep 16 ends with TIME = 7.491599174e-07 DT = 9.571056553e-08 Timing summary: Advection :0.7138454914 seconds MAC Proj :0.1691300869 seconds Nodal Proj :0.3385374546 seconds Thermal :0.1822328568 seconds Reactions :1.776566505 seconds Misc :0.2190709114 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.399823904 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.082276804 Iter = 0.079078868 Bottom = 0.004385982 <<< 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.03593483 Iter = 0.033891236 Bottom = 0.000248501 <<< 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.0831426 Iter = 0.079905623 Bottom = 0.004480941 <<< 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.035810621 Iter = 0.033770193 Bottom = 0.000256793 <<< 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.212732049 Iter = 0.205473591 Bottom = 0.003757521 Done calling nodal solver Timestep 17 ends with TIME = 8.449182257e-07 DT = 9.575830833e-08 Timing summary: Advection :0.7434825897 seconds MAC Proj :0.1756622791 seconds Nodal Proj :0.3346393108 seconds Thermal :0.1989870071 seconds Reactions :1.773914337 seconds Misc :0.2237384319 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.450787783 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.082724181 Iter = 0.079468204 Bottom = 0.004147066 <<< 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.035864122 Iter = 0.033835694 Bottom = 0.000245086 <<< 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.082475388 Iter = 0.079265287 Bottom = 0.0043974 <<< 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.035679746 Iter = 0.033610358 Bottom = 0.00025067 <<< 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.213403511 Iter = 0.206147053 Bottom = 0.003734564 Done calling nodal solver Timestep 18 ends with TIME = 9.407071664e-07 DT = 9.578894066e-08 Timing summary: Advection :0.7443621159 seconds MAC Proj :0.1753623486 seconds Nodal Proj :0.3345906734 seconds Thermal :0.198335886 seconds Reactions :1.766933441 seconds Misc :0.2237842083 seconds Base State :7.390975952e-06 seconds Time to advance time step: 3.443743229 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.082486846 Iter = 0.079236453 Bottom = 0.004130772 <<< 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.035424953 Iter = 0.033395649 Bottom = 0.000244753 <<< 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.08282314 Iter = 0.079613621 Bottom = 0.004352732 <<< 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.035732579 Iter = 0.033711942 Bottom = 0.000248251 <<< 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.212848549 Iter = 0.205580716 Bottom = 0.003807606 Done calling nodal solver Timestep 19 ends with TIME = 1.036537244e-06 DT = 9.5830078e-08 Timing summary: Advection :0.7457613945 seconds MAC Proj :0.1755020618 seconds Nodal Proj :0.3344187737 seconds Thermal :0.198839426 seconds Reactions :1.768179178 seconds Misc :0.2236001492 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.44666791 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.082085678 Iter = 0.078903034 Bottom = 0.004391024 <<< 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.035485661 Iter = 0.033468941 Bottom = 0.00024321 <<< 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.082338345 Iter = 0.079148244 Bottom = 0.004456315 <<< 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.035445744 Iter = 0.033434024 Bottom = 0.000242168 <<< 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.212490756 Iter = 0.205217797 Bottom = 0.003796938 Done calling nodal solver Timestep 20 ends with TIME = 1.132393714e-06 DT = 9.585646999e-08 Timing summary: Advection :0.7444000244 seconds MAC Proj :0.1746275425 seconds Nodal Proj :0.3339738846 seconds Thermal :0.198176384 seconds Reactions :1.763564587 seconds Misc :0.2247238159 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.439849854 Writing plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.294322014 Writing small plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.293727636 Total Time: 80.51255369 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-13-g0278677bbc8c) finalized