MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.01-9-gaa079d4dc93b) initialized Calling Setup() Calling ReadParameters() reading extern runtime parameters ... Calling VariableSetup() Calling BCSetup() Calling BaseStateGeometry::Init() Calling Init() Calling InitData() initdata model_File = model_files/kepler_new_6.25e8.hybrid.hse.320 cutoff densities: low density cutoff (for mapping the model) = 3000000 buoyancy cutoff density (for zeroing rho - rho_0, centrifugal term) = 15000000 anelastic cutoff = 3000000 Writing plotfile flame_3d_pltInitData after InitData ... Level 0 create thermal coeffs: Time to write plotfile: 1.717870712 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.52299787 Iter = 0.507195454 Bottom = 0.007196412 Done calling nodal solver Writing plotfile flame_3d_pltafter_InitProj after InitProj ... Level 0 create thermal coeffs: Time to write plotfile: 1.79661274 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.379516576 Iter = 0.363913369 Bottom = 0.005311028 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.467920956 Iter = 0.453644923 Bottom = 0.006687369 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.459082576 Iter = 0.444306416 Bottom = 0.006450326 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.801627636 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.150831 Iter = 0.142509789 Bottom = 0.005492027 <<< 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.055951919 Iter = 0.052640693 Bottom = 0.000283959 <<< 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.147152772 Iter = 0.142349289 Bottom = 0.005679363 <<< 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.054886997 Iter = 0.051468604 Bottom = 0.000276251 <<< 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.426889282 Iter = 0.414267841 Bottom = 0.005465446 Done calling nodal solver Timestep 0 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :1.310234785 seconds MAC Proj :0.3161220551 seconds Nodal Proj :0.6338288784 seconds Thermal :0.3280050755 seconds Reactions :2.621032953 seconds Misc :0.325933218 seconds Base State :3.576278687e-06 seconds Time to advance time step: 5.5770998 Writing plotfile 0 after all initialization ... Level 0 create thermal coeffs: Time to write plotfile: 1.320343494 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.077601449 Iter = 0.073920971 Bottom = 0.003328185 <<< 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.029209029 Iter = 0.026811516 Bottom = 0.000197209 <<< 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.058009555 Iter = 0.054277618 Bottom = 0.002575223 <<< 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.028432024 Iter = 0.026195512 Bottom = 0.000173375 <<< 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.222983088 Iter = 0.215137713 Bottom = 0.003552977 Done calling nodal solver Timestep 1 ends with TIME = 9.529696764e-09 DT = 9.529696764e-09 Timing summary: Advection :0.7914729118 seconds MAC Proj :0.1467313766 seconds Nodal Proj :0.3507814407 seconds Thermal :0.1951975822 seconds Reactions :1.604038954 seconds Misc :0.2093498707 seconds Base State :3.814697266e-06 seconds Time to advance time step: 3.328635693 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.129129928 Iter = 0.125070782 Bottom = 0.005560361 <<< 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.054783288 Iter = 0.051797397 Bottom = 0.000321377 <<< 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.068581068 Iter = 0.064906925 Bottom = 0.003647312 <<< 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.02785023 Iter = 0.025830428 Bottom = 0.000180375 <<< 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.221913582 Iter = 0.214381667 Bottom = 0.00388152 Done calling nodal solver Timestep 2 ends with TIME = 2.096533288e-08 DT = 1.143563612e-08 Timing summary: Advection :0.9344067574 seconds MAC Proj :0.2099168301 seconds Nodal Proj :0.3468141556 seconds Thermal :0.2191309929 seconds Reactions :1.60546422 seconds Misc :0.2251782417 seconds Base State :3.099441528e-06 seconds Time to advance time step: 3.541331291 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.077867325 Iter = 0.074752309 Bottom = 0.003996645 <<< 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.029875615 Iter = 0.027607061 Bottom = 0.000204835 <<< 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.067378937 Iter = 0.064237921 Bottom = 0.003483687 <<< 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.028694109 Iter = 0.026674474 Bottom = 0.000209001 <<< 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.225627518 Iter = 0.216977514 Bottom = 0.003901729 Done calling nodal solver Timestep 3 ends with TIME = 3.468809622e-08 DT = 1.372276334e-08 Timing summary: Advection :0.7425508499 seconds MAC Proj :0.1553161144 seconds Nodal Proj :0.3520736694 seconds Thermal :0.18754673 seconds Reactions :1.61301446 seconds Misc :0.2106897831 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.261583805 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.078085285 Iter = 0.074740018 Bottom = 0.004077062 <<< 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.028382565 Iter = 0.026291888 Bottom = 0.000181835 <<< 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.069173738 Iter = 0.066030597 Bottom = 0.003545602 <<< 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.028481108 Iter = 0.026447597 Bottom = 0.000174835 <<< 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.223656299 Iter = 0.216142592 Bottom = 0.003812145 Done calling nodal solver Timestep 4 ends with TIME = 5.115541223e-08 DT = 1.646731601e-08 Timing summary: Advection :0.7426800728 seconds MAC Proj :0.1574885845 seconds Nodal Proj :0.3486461639 seconds Thermal :0.18440938 seconds Reactions :1.623684645 seconds Misc :0.2071447372 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.26443553 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.077550782 Iter = 0.074388891 Bottom = 0.004021897 <<< 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.028303524 Iter = 0.026282721 Bottom = 0.000172293 <<< 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.081206218 Iter = 0.078041243 Bottom = 0.004170313 <<< 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.028036564 Iter = 0.025983762 Bottom = 0.000172917 <<< 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.224159885 Iter = 0.215618507 Bottom = 0.003807771 Done calling nodal solver Timestep 5 ends with TIME = 7.091619144e-08 DT = 1.976077921e-08 Timing summary: Advection :0.7534031868 seconds MAC Proj :0.1689071655 seconds Nodal Proj :0.3497064114 seconds Thermal :0.1866602898 seconds Reactions :1.651454687 seconds Misc :0.2091977596 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.319723368 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.086534038 Iter = 0.083396812 Bottom = 0.004503149 <<< 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.029471446 Iter = 0.027422019 Bottom = 0.000199959 <<< 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.07697032 Iter = 0.073824304 Bottom = 0.004003396 <<< 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.030115283 Iter = 0.027478478 Bottom = 0.000202543 <<< 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.224982265 Iter = 0.217031014 Bottom = 0.003897145 Done calling nodal solver Timestep 6 ends with TIME = 9.462912649e-08 DT = 2.371293505e-08 Timing summary: Advection :0.7545611858 seconds MAC Proj :0.1737346649 seconds Nodal Proj :0.3510582447 seconds Thermal :0.1882214546 seconds Reactions :1.669339418 seconds Misc :0.2116436958 seconds Base State :4.053115845e-06 seconds Time to advance time step: 3.34893465 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.092980572 Iter = 0.089829888 Bottom = 0.004577065 <<< 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.028579275 Iter = 0.026474222 Bottom = 0.000191709 <<< 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.07729828 Iter = 0.074191347 Bottom = 0.004054562 <<< 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.028729859 Iter = 0.026717682 Bottom = 0.000181709 <<< 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.227723987 Iter = 0.218917774 Bottom = 0.003703602 Done calling nodal solver Timestep 7 ends with TIME = 1.230846486e-07 DT = 2.845552206e-08 Timing summary: Advection :0.7624540329 seconds MAC Proj :0.1806726456 seconds Nodal Proj :0.3522479534 seconds Thermal :0.1858315468 seconds Reactions :1.692599773 seconds Misc :0.2152552605 seconds Base State :3.814697266e-06 seconds Time to advance time step: 3.389446974 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.085904159 Iter = 0.082752976 Bottom = 0.004463233 <<< 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.02885836 Iter = 0.026807058 Bottom = 0.000186877 <<< 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.076985863 Iter = 0.073882763 Bottom = 0.004028894 <<< 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.029660989 Iter = 0.027402019 Bottom = 0.000200126 <<< 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.225313099 Iter = 0.215873967 Bottom = 0.003707895 Done calling nodal solver Timestep 8 ends with TIME = 1.57231275e-07 DT = 3.414662647e-08 Timing summary: Advection :0.7621042728 seconds MAC Proj :0.1729960442 seconds Nodal Proj :0.350902319 seconds Thermal :0.1881546974 seconds Reactions :1.72401762 seconds Misc :0.2153103352 seconds Base State :3.337860107e-06 seconds Time to advance time step: 3.413893461 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.087640877 Iter = 0.084568444 Bottom = 0.004443024 <<< 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.029098361 Iter = 0.027041142 Bottom = 0.000190709 <<< 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.076630569 Iter = 0.073493385 Bottom = 0.003957104 <<< 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.028816109 Iter = 0.026583556 Bottom = 0.000181043 <<< 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.223344672 Iter = 0.215789299 Bottom = 0.003736688 Done calling nodal solver Timestep 9 ends with TIME = 1.982072268e-07 DT = 4.097595177e-08 Timing summary: Advection :0.7601108551 seconds MAC Proj :0.174352169 seconds Nodal Proj :0.3472864628 seconds Thermal :0.1860382557 seconds Reactions :1.742019892 seconds Misc :0.222826004 seconds Base State :3.576278687e-06 seconds Time to advance time step: 3.43300724 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.088098629 Iter = 0.085012446 Bottom = 0.004484816 <<< 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.028935277 Iter = 0.026933349 Bottom = 0.000190127 <<< 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.077371656 Iter = 0.074072514 Bottom = 0.003952853 <<< 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.028139523 Iter = 0.026070304 Bottom = 0.000185959 <<< 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.226209812 Iter = 0.218001977 Bottom = 0.00392427 Done calling nodal solver Timestep 10 ends with TIME = 2.473783689e-07 DT = 4.917114212e-08 Timing summary: Advection :0.7641949654 seconds MAC Proj :0.1756689548 seconds Nodal Proj :0.3524496555 seconds Thermal :0.1836385727 seconds Reactions :1.755758524 seconds Misc :0.2215001583 seconds Base State :3.814697266e-06 seconds Time to advance time step: 3.453623056 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.083925024 Iter = 0.080934341 Bottom = 0.004342522 <<< 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.0286529 Iter = 0.026486847 Bottom = 0.000189125 <<< 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.083407438 Iter = 0.080399964 Bottom = 0.004306313 <<< 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.027748771 Iter = 0.02576951 Bottom = 0.000176794 <<< 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.479428725 Iter = 0.464913066 Bottom = 0.007691624 Done calling nodal solver Timestep 11 ends with TIME = 3.063837395e-07 DT = 5.900537055e-08 Timing summary: Advection :0.8006560802 seconds MAC Proj :0.1770009995 seconds Nodal Proj :0.6887471676 seconds Thermal :0.1850781441 seconds Reactions :1.784157753 seconds Misc :0.225194931 seconds Base State :3.814697266e-06 seconds Time to advance time step: 3.861491919 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.17493546 Iter = 0.17024006 Bottom = 0.007965082 <<< 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.056717256 Iter = 0.053091362 Bottom = 0.000386293 <<< 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.171564276 Iter = 0.166739 Bottom = 0.007643416 <<< 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.056988841 Iter = 0.053386114 Bottom = 0.000388752 <<< 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.460991003 Iter = 0.446687345 Bottom = 0.006744284 Done calling nodal solver Timestep 12 ends with TIME = 3.771901841e-07 DT = 7.080644466e-08 Timing summary: Advection :1.295273542 seconds MAC Proj :0.3628418446 seconds Nodal Proj :0.6859264374 seconds Thermal :0.3286674023 seconds Reactions :3.152982235 seconds Misc :0.3479578495 seconds Base State :4.053115845e-06 seconds Time to advance time step: 6.174438 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.168258842 Iter = 0.163554734 Bottom = 0.007784081 <<< 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.074739351 Iter = 0.071359625 Bottom = 0.00050692 <<< 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.168469217 Iter = 0.163119565 Bottom = 0.007927207 <<< 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.075786356 Iter = 0.072326755 Bottom = 0.000532836 <<< 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.472807523 Iter = 0.457989904 Bottom = 0.007514581 Done calling nodal solver Timestep 13 ends with TIME = 4.621579177e-07 DT = 8.496773359e-08 Timing summary: Advection :1.288070917 seconds MAC Proj :0.3532485962 seconds Nodal Proj :0.6945073605 seconds Thermal :0.3681623936 seconds Reactions :3.167862654 seconds Misc :0.3816711903 seconds Base State :4.053115845e-06 seconds Time to advance time step: 6.254281998 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.164935533 Iter = 0.160489967 Bottom = 0.007778832 <<< 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.074446599 Iter = 0.071011665 Bottom = 0.000560087 <<< 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.165372784 Iter = 0.160543009 Bottom = 0.007532664 <<< 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.07327351 Iter = 0.06983645 Bottom = 0.000497043 <<< 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.45469397 Iter = 0.44075023 Bottom = 0.007038495 Done calling nodal solver Timestep 14 ends with TIME = 5.577754795e-07 DT = 9.561756174e-08 Timing summary: Advection :1.283554792 seconds MAC Proj :0.3461792469 seconds Nodal Proj :0.6748688221 seconds Thermal :0.362844944 seconds Reactions :2.896687984 seconds Misc :0.3293771744 seconds Base State :3.576278687e-06 seconds Time to advance time step: 5.894285202 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.170422102 Iter = 0.165679328 Bottom = 0.007738582 <<< 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.072874633 Iter = 0.069490906 Bottom = 0.000497462 <<< 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.167342462 Iter = 0.162695937 Bottom = 0.007546831 <<< 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.077810867 Iter = 0.074394557 Bottom = 0.000550045 <<< 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.465486652 Iter = 0.451392786 Bottom = 0.006997036 Done calling nodal solver Timestep 15 ends with TIME = 6.534493518e-07 DT = 9.567387238e-08 Timing summary: Advection :1.286169291 seconds MAC Proj :0.3537392616 seconds Nodal Proj :0.6866848469 seconds Thermal :0.3687601089 seconds Reactions :2.833623171 seconds Misc :0.3335988522 seconds Base State :3.576278687e-06 seconds Time to advance time step: 5.863332748 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.167577338 Iter = 0.162834896 Bottom = 0.007691583 <<< 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.075787522 Iter = 0.072533214 Bottom = 0.000518128 <<< 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.166306623 Iter = 0.16137493 Bottom = 0.00723537 <<< 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.075534188 Iter = 0.072219255 Bottom = 0.000471962 <<< 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.4396781 Iter = 0.425603901 Bottom = 0.006804162 Done calling nodal solver Timestep 16 ends with TIME = 7.491599174e-07 DT = 9.571056553e-08 Timing summary: Advection :1.297283411 seconds MAC Proj :0.3521404266 seconds Nodal Proj :0.6479489803 seconds Thermal :0.3694243431 seconds Reactions :2.804925442 seconds Misc :0.3326721191 seconds Base State :4.291534424e-06 seconds Time to advance time step: 5.805166483 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.16475624 Iter = 0.159662713 Bottom = 0.00731308 <<< 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.073060175 Iter = 0.069709241 Bottom = 0.000465544 <<< 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.165072867 Iter = 0.160531051 Bottom = 0.007300874 <<< 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.071849127 Iter = 0.068501401 Bottom = 0.000481461 <<< 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.439471015 Iter = 0.425644984 Bottom = 0.006878287 Done calling nodal solver Timestep 17 ends with TIME = 8.449182257e-07 DT = 9.575830833e-08 Timing summary: Advection :1.281953096 seconds MAC Proj :0.3457894325 seconds Nodal Proj :0.6422810555 seconds Thermal :0.3543417454 seconds Reactions :2.806772947 seconds Misc :0.3336398602 seconds Base State :3.814697266e-06 seconds Time to advance time step: 5.765524626 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.160176299 Iter = 0.155778068 Bottom = 0.006845202 <<< 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.070079576 Iter = 0.067007769 Bottom = 0.000427127 <<< 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.157743119 Iter = 0.153481055 Bottom = 0.006812828 <<< 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.070633496 Iter = 0.067490562 Bottom = 0.000435295 <<< 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.399679895 Iter = 0.387350955 Bottom = 0.005773613 Done calling nodal solver Timestep 18 ends with TIME = 9.407071664e-07 DT = 9.578894066e-08 Timing summary: Advection :1.222531796 seconds MAC Proj :0.3328726292 seconds Nodal Proj :0.5822906494 seconds Thermal :0.3371317387 seconds Reactions :2.452083826 seconds Misc :0.2974100113 seconds Base State :3.814697266e-06 seconds Time to advance time step: 5.224946022 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.161371018 Iter = 0.157036829 Bottom = 0.006919871 <<< 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.0726813 Iter = 0.069539617 Bottom = 0.000441793 <<< 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.152258512 Iter = 0.147969907 Bottom = 0.006603118 <<< 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.063729837 Iter = 0.061085531 Bottom = 0.00039421 <<< 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.44646144 Iter = 0.433365496 Bottom = 0.006731577 Done calling nodal solver Timestep 19 ends with TIME = 1.036537244e-06 DT = 9.5830078e-08 Timing summary: Advection :1.183233023 seconds MAC Proj :0.3278591633 seconds Nodal Proj :0.6553001404 seconds Thermal :0.3135666847 seconds Reactions :2.427036285 seconds Misc :0.2411322594 seconds Base State :3.814697266e-06 seconds Time to advance time step: 5.148620129 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.149211079 Iter = 0.144870598 Bottom = 0.005363573 <<< 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.037273071 Iter = 0.035018893 Bottom = 0.00026325 <<< 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.088647093 Iter = 0.084998407 Bottom = 0.004833484 <<< 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.038148326 Iter = 0.035996565 Bottom = 0.000267251 <<< 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.229865214 Iter = 0.222018673 Bottom = 0.004015229 Done calling nodal solver Timestep 20 ends with TIME = 1.132393714e-06 DT = 9.585646999e-08 Timing summary: Advection :0.9482433796 seconds MAC Proj :0.2508742809 seconds Nodal Proj :0.3572545052 seconds Thermal :0.2054429054 seconds Reactions :2.114015102 seconds Misc :0.235316515 seconds Base State :4.053115845e-06 seconds Time to advance time step: 4.111582994 Writing plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.345907211 Writing small plotfile 20 ... Level 0 create thermal coeffs: Time to write plotfile: 1.341055393 Total Time: 110.8869588 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-9-gaa079d4dc93b) finalized