MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.01-22-ge34f5c67fec9) 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.160024166 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.255925159 Iter = 0.247702114 Bottom = 0.003901565 Done calling nodal solver Writing plotfile flame_3d_pltafter_InitProj after InitProj ... Level 0 create thermal coeffs: Time to write plotfile: 1.112236023 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.184936393 Iter = 0.177036724 Bottom = 0.002951224 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.250091877 Iter = 0.237214305 Bottom = 0.004124356 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.212540962 Iter = 0.205575798 Bottom = 0.003411853 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.125135183 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.073683364 Iter = 0.068342251 Bottom = 0.003138102 <<< 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.027601943 Iter = 0.025604891 Bottom = 0.000172958 <<< 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.073828365 Iter = 0.069136755 Bottom = 0.003287852 <<< 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.026841648 Iter = 0.024893429 Bottom = 0.000185793 <<< 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.213113007 Iter = 0.206150052 Bottom = 0.003427227 Done calling nodal solver Timestep 0 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.7608020306 seconds MAC Proj :0.1604700089 seconds Nodal Proj :0.333170414 seconds Thermal :0.1769940853 seconds Reactions :1.559479237 seconds Misc :0.2212486267 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.288615227 Writing plotfile 0 after all initialization ... Level 0 create thermal coeffs: Time to write plotfile: 1.083586693 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.074098325 Iter = 0.069437591 Bottom = 0.003168768 <<< 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.027063316 Iter = 0.025134847 Bottom = 0.000177542 <<< 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.059028533 Iter = 0.05426859 Bottom = 0.002612639 <<< 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.028167863 Iter = 0.026150477 Bottom = 0.000177876 <<< 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.214477472 Iter = 0.207416183 Bottom = 0.003517477 Done calling nodal solver Timestep 1 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.7451739311 seconds MAC Proj :0.1458978653 seconds Nodal Proj :0.3353559971 seconds Thermal :0.1798746586 seconds Reactions :1.544775248 seconds Misc :0.2039833069 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.190225601 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.074445993 Iter = 0.069786259 Bottom = 0.003162559 <<< 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.027178441 Iter = 0.025239722 Bottom = 0.000172542 <<< 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.065800112 Iter = 0.061086003 Bottom = 0.002973477 <<< 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.027120857 Iter = 0.025138097 Bottom = 0.00017225 <<< 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.212893672 Iter = 0.205865508 Bottom = 0.003418186 Done calling nodal solver Timestep 2 ends with TIME = 2.096533288e-08 DT = 1.143563612e-08 Timing summary: Advection :0.7433214188 seconds MAC Proj :0.1528806686 seconds Nodal Proj :0.3320214748 seconds Thermal :0.176271677 seconds Reactions :1.546830654 seconds Misc :0.2068994045 seconds Base State :2.861022949e-06 seconds Time to advance time step: 3.16263485 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.073518613 Iter = 0.068829337 Bottom = 0.003156892 <<< 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.026808523 Iter = 0.024881012 Bottom = 0.000170586 <<< 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.065335734 Iter = 0.060632333 Bottom = 0.002809183 <<< 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.026976899 Iter = 0.025057096 Bottom = 0.000170585 <<< 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.212446627 Iter = 0.20541613 Bottom = 0.00344702 Done calling nodal solver Timestep 3 ends with TIME = 3.468809622e-08 DT = 1.372276334e-08 Timing summary: Advection :0.7444589138 seconds MAC Proj :0.1514430046 seconds Nodal Proj :0.3320724964 seconds Thermal :0.1753163338 seconds Reactions :1.562522411 seconds Misc :0.208013773 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.178247213 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.074244867 Iter = 0.069619258 Bottom = 0.003270059 <<< 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.027030024 Iter = 0.025086429 Bottom = 0.000172459 <<< 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.066129989 Iter = 0.060440083 Bottom = 0.002778641 <<< 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.027073399 Iter = 0.025126638 Bottom = 0.000170127 <<< 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.216122857 Iter = 0.209120151 Bottom = 0.003478727 Done calling nodal solver Timestep 4 ends with TIME = 5.115541223e-08 DT = 1.646731601e-08 Timing summary: Advection :0.7473516464 seconds MAC Proj :0.152971983 seconds Nodal Proj :0.3354296684 seconds Thermal :0.1751914024 seconds Reactions :1.564337254 seconds Misc :0.2082846165 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.187983036 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.073454902 Iter = 0.068751127 Bottom = 0.003179141 <<< 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.026900606 Iter = 0.024965053 Bottom = 0.00018146 <<< 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.073360569 Iter = 0.068719668 Bottom = 0.003239143 <<< 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.026993814 Iter = 0.025055721 Bottom = 0.000166585 <<< 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.215069469 Iter = 0.208058598 Bottom = 0.003274309 Done calling nodal solver Timestep 5 ends with TIME = 7.091619144e-08 DT = 1.976077921e-08 Timing summary: Advection :0.7543013096 seconds MAC Proj :0.1594843864 seconds Nodal Proj :0.3342256546 seconds Thermal :0.1752068996 seconds Reactions :1.58994031 seconds Misc :0.208997488 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.226555824 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.082150367 Iter = 0.0775173 Bottom = 0.003582436 <<< 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.027008148 Iter = 0.025085095 Bottom = 0.000169126 <<< 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.076533711 Iter = 0.07170356 Bottom = 0.003350769 <<< 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.029065785 Iter = 0.02709294 Bottom = 0.000184209 <<< 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.21238908 Iter = 0.20536 Bottom = 0.003372436 Done calling nodal solver Timestep 6 ends with TIME = 9.462912649e-08 DT = 2.371293505e-08 Timing summary: Advection :0.7634267807 seconds MAC Proj :0.1714806557 seconds Nodal Proj :0.3325893879 seconds Thermal :0.1810457706 seconds Reactions :1.61317873 seconds Misc :0.2133293152 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.279504299 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.085473551 Iter = 0.080718275 Bottom = 0.0043234 <<< 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.027691235 Iter = 0.025737391 Bottom = 0.000173583 <<< 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.077062922 Iter = 0.072300354 Bottom = 0.003957439 <<< 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.027671985 Iter = 0.025719849 Bottom = 0.000169666 <<< 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.218258695 Iter = 0.210977781 Bottom = 0.003476186 Done calling nodal solver Timestep 7 ends with TIME = 1.230846486e-07 DT = 2.845552206e-08 Timing summary: Advection :0.7982754707 seconds MAC Proj :0.1757206917 seconds Nodal Proj :0.3422210217 seconds Thermal :0.1791372299 seconds Reactions :1.630722761 seconds Misc :0.2142920494 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.344830275 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.087082268 Iter = 0.082324075 Bottom = 0.004425315 <<< 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.027917861 Iter = 0.02593535 Bottom = 0.000179792 <<< 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.074292449 Iter = 0.069539465 Bottom = 0.003749812 <<< 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.027538818 Iter = 0.025584682 Bottom = 0.000201085 <<< 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.219220242 Iter = 0.211940286 Bottom = 0.003564268 Done calling nodal solver Timestep 8 ends with TIME = 1.57231275e-07 DT = 3.414662647e-08 Timing summary: Advection :0.8018352985 seconds MAC Proj :0.1746134758 seconds Nodal Proj :0.3432450294 seconds Thermal :0.1800091267 seconds Reactions :1.672626257 seconds Misc :0.2189440727 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.395724773 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.084305587 Iter = 0.079569811 Bottom = 0.004050187 <<< 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.033884144 Iter = 0.031425464 Bottom = 0.00026721 <<< 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.07634571 Iter = 0.071570142 Bottom = 0.003819727 <<< 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.028646782 Iter = 0.02668498 Bottom = 0.000181293 <<< 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.216201018 Iter = 0.208952228 Bottom = 0.003603769 Done calling nodal solver Timestep 9 ends with TIME = 1.982072268e-07 DT = 4.097595177e-08 Timing summary: Advection :0.8101613522 seconds MAC Proj :0.1741809845 seconds Nodal Proj :0.3401134014 seconds Thermal :0.1906638145 seconds Reactions :1.677231312 seconds Misc :0.227189064 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.424020052 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.081746531 Iter = 0.076974755 Bottom = 0.004123231 <<< 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.027594734 Iter = 0.025628474 Bottom = 0.000184334 <<< 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.075046995 Iter = 0.070233885 Bottom = 0.00367402 <<< 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.027564484 Iter = 0.025579224 Bottom = 0.000179251 <<< 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.219734662 Iter = 0.211544742 Bottom = 0.003668437 Done calling nodal solver Timestep 10 ends with TIME = 2.473783689e-07 DT = 4.917114212e-08 Timing summary: Advection :0.8048856258 seconds MAC Proj :0.1700179577 seconds Nodal Proj :0.3434228897 seconds Thermal :0.1779170036 seconds Reactions :1.70854187 seconds Misc :0.2214651108 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.430763245 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.082659327 Iter = 0.077898135 Bottom = 0.004125316 <<< 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.027688027 Iter = 0.025655598 Bottom = 0.0001745 <<< 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.084638839 Iter = 0.079879854 Bottom = 0.004232106 <<< 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.028601573 Iter = 0.02656452 Bottom = 0.000190085 <<< 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.215878558 Iter = 0.208632893 Bottom = 0.003996063 Done calling nodal solver Timestep 11 ends with TIME = 3.063837395e-07 DT = 5.900537055e-08 Timing summary: Advection :0.8019382954 seconds MAC Proj :0.1805052757 seconds Nodal Proj :0.3391320705 seconds Thermal :0.1814582348 seconds Reactions :1.730305195 seconds Misc :0.2275028229 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.465299606 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.084508797 Iter = 0.079757062 Bottom = 0.004281899 <<< 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.027603817 Iter = 0.025623099 Bottom = 0.000177375 <<< 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.084417671 Iter = 0.07964952 Bottom = 0.004432898 <<< 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.026548604 Iter = 0.024563176 Bottom = 0.000175084 <<< 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.212437997 Iter = 0.205200915 Bottom = 0.003421269 Done calling nodal solver Timestep 12 ends with TIME = 3.771901841e-07 DT = 7.080644466e-08 Timing summary: Advection :0.8019940853 seconds MAC Proj :0.1822304726 seconds Nodal Proj :0.3349387646 seconds Thermal :0.178504467 seconds Reactions :1.74542141 seconds Misc :0.2280173302 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.475589037 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.085603261 Iter = 0.08086836 Bottom = 0.004229313 <<< 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.034643982 Iter = 0.032643554 Bottom = 0.000231085 <<< 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.085544635 Iter = 0.080795359 Bottom = 0.004503483 <<< 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.036615242 Iter = 0.03467044 Bottom = 0.000262417 <<< 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.219018616 Iter = 0.211072656 Bottom = 0.003620146 Done calling nodal solver Timestep 13 ends with TIME = 4.621579177e-07 DT = 8.496773359e-08 Timing summary: Advection :0.800246954 seconds MAC Proj :0.1847879887 seconds Nodal Proj :0.3430473804 seconds Thermal :0.1951408386 seconds Reactions :1.769273281 seconds Misc :0.2297084332 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.526675701 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.084963215 Iter = 0.080056647 Bottom = 0.004145107 <<< 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.035757905 Iter = 0.033769102 Bottom = 0.000263043 <<< 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.085712761 Iter = 0.080886985 Bottom = 0.00442269 <<< 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.034804191 Iter = 0.032825804 Bottom = 0.000249043 <<< 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.215734806 Iter = 0.208488725 Bottom = 0.00363727 Done calling nodal solver Timestep 14 ends with TIME = 5.577754795e-07 DT = 9.561756174e-08 Timing summary: Advection :0.8025708199 seconds MAC Proj :0.1839387417 seconds Nodal Proj :0.3383970261 seconds Thermal :0.1938409805 seconds Reactions :1.78020215 seconds Misc :0.231659174 seconds Base State :2.861022949e-06 seconds Time to advance time step: 3.53507781 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.085525469 Iter = 0.080785984 Bottom = 0.004416609 <<< 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.034756482 Iter = 0.032783555 Bottom = 0.000232628 <<< 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.085867637 Iter = 0.080921652 Bottom = 0.004398815 <<< 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.034680648 Iter = 0.032687054 Bottom = 0.000230542 <<< 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.219087741 Iter = 0.211228616 Bottom = 0.00387448 Done calling nodal solver Timestep 15 ends with TIME = 6.534493518e-07 DT = 9.567387238e-08 Timing summary: Advection :0.8025403023 seconds MAC Proj :0.18471241 seconds Nodal Proj :0.3434481621 seconds Thermal :0.1927378178 seconds Reactions :1.784221649 seconds Misc :0.2322664261 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.54442811 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.084652047 Iter = 0.079915146 Bottom = 0.004286317 <<< 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.036912077 Iter = 0.034936608 Bottom = 0.000249878 <<< 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.085177342 Iter = 0.08042744 Bottom = 0.004403066 <<< 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.034358646 Iter = 0.03263822 Bottom = 0.000225586 <<< 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.217536858 Iter = 0.210761988 Bottom = 0.003712228 Done calling nodal solver Timestep 16 ends with TIME = 7.491599174e-07 DT = 9.571056553e-08 Timing summary: Advection :0.8549349308 seconds MAC Proj :0.1830508709 seconds Nodal Proj :0.3393788338 seconds Thermal :0.1976544857 seconds Reactions :1.777191877 seconds Misc :0.2315924168 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.588476658 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.085418467 Iter = 0.080768942 Bottom = 0.004439068 <<< 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.035320693 Iter = 0.033581476 Bottom = 0.000225835 <<< 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.086956143 Iter = 0.082236076 Bottom = 0.00431794 <<< 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.035954155 Iter = 0.033727768 Bottom = 0.000228875 <<< 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.216868438 Iter = 0.21014886 Bottom = 0.003971982 Done calling nodal solver Timestep 17 ends with TIME = 8.449182257e-07 DT = 9.575830833e-08 Timing summary: Advection :0.7622251511 seconds MAC Proj :0.1850328445 seconds Nodal Proj :0.3364500999 seconds Thermal :0.1859130859 seconds Reactions :1.775635004 seconds Misc :0.2308802605 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.480565786 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.085342093 Iter = 0.0806439 Bottom = 0.004074688 <<< 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.036382991 Iter = 0.034676399 Bottom = 0.000258501 <<< 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.083490124 Iter = 0.078849849 Bottom = 0.00422169 <<< 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.036437491 Iter = 0.034742149 Bottom = 0.000247503 <<< 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.218337529 Iter = 0.211521242 Bottom = 0.003609436 Done calling nodal solver Timestep 18 ends with TIME = 9.407071664e-07 DT = 9.578894066e-08 Timing summary: Advection :0.7666275501 seconds MAC Proj :0.1812393665 seconds Nodal Proj :0.3370821476 seconds Thermal :0.1882209778 seconds Reactions :1.767836332 seconds Misc :0.2297592163 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.475148678 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.084687089 Iter = 0.080001063 Bottom = 0.004048316 <<< 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.03632049 Iter = 0.034622023 Bottom = 0.000267126 <<< 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.082054242 Iter = 0.077384799 Bottom = 0.004292106 <<< 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.034688565 Iter = 0.032990639 Bottom = 0.000227209 <<< 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.212333205 Iter = 0.205622751 Bottom = 0.003684355 Done calling nodal solver Timestep 19 ends with TIME = 1.036537244e-06 DT = 9.5830078e-08 Timing summary: Advection :0.7662293911 seconds MAC Proj :0.1789529324 seconds Nodal Proj :0.3299484253 seconds Thermal :0.1840672493 seconds Reactions :1.766027689 seconds Misc :0.229695797 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.459273338 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.082578078 Iter = 0.077849968 Bottom = 0.004328232 <<< 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.034271771 Iter = 0.032591304 Bottom = 0.000238377 <<< 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.082256784 Iter = 0.077605883 Bottom = 0.004403148 <<< 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.034704773 Iter = 0.033017847 Bottom = 0.000227001 <<< 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.212681873 Iter = 0.205998379 Bottom = 0.003740311 Done calling nodal solver Timestep 20 ends with TIME = 1.132393714e-06 DT = 9.585646999e-08 Timing summary: Advection :0.7630152702 seconds MAC Proj :0.1771130562 seconds Nodal Proj :0.3293297291 seconds Thermal :0.1818537712 seconds Reactions :1.765399456 seconds Misc :0.2300188541 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.451102495 Writing plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.292636871 Writing small plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.29177928 Total Time: 82.5226295 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-22-ge34f5c67fec9) finalized