Initializing AMReX (26.07-54-gaa5b8cc69640)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.07-54-gaa5b8cc69640) initialized Calling Setup() Calling ReadParameters() reading extern runtime parameters ... Calling VariableSetup() Calling BCSetup() Calling BaseStateGeometry::Init() Calling Init() Calling InitData() initdata model_File = Maximum HSE Error = 2.348271906e-14 (after putting initial model into base state arrays, and for density < base_cutoff_density) Writing plotfile double_bubble_2d_pltInitData after InitData Time to write plotfile: 0.069107629 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002857896 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile double_bubble_2d_pltafter_InitProj after InitProj Time to write plotfile: 0.065050473 Call to firstdt for level 0 gives dt_lev = 7.656411368e-08 Multiplying dt_lev by init_shrink; dt_lev = 7.656411368e-08 Minimum firstdt over all levels = 7.656411368e-08 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002776174 Iter = 0 Bottom = 0 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 2.298789937e-05 Minimum estdt over all levels = 2.298789937e-05 Call to estdt at end of istep_divu_iter = 1 gives dt = 2.298789937e-05 Multiplying dt by init_shrink; dt = 2.298789937e-05 Ignoring this new dt since it's larger than the previous dt = 7.656411368e-08 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002861415 Iter = 0 Bottom = 0 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 2.298789937e-05 Minimum estdt over all levels = 2.298789937e-05 Call to estdt at end of istep_divu_iter = 2 gives dt = 2.298789937e-05 Multiplying dt by init_shrink; dt = 2.298789937e-05 Ignoring this new dt since it's larger than the previous dt = 7.656411368e-08 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002802381 Iter = 0 Bottom = 0 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 2.298789937e-05 Minimum estdt over all levels = 2.298789937e-05 Call to estdt at end of istep_divu_iter = 3 gives dt = 2.298789937e-05 Multiplying dt by init_shrink; dt = 2.298789937e-05 Ignoring this new dt since it's larger than the previous dt = 7.656411368e-08 Writing plotfile double_bubble_2d_pltafter_DivuIter after final DivuIter Time to write plotfile: 0.065209706 Doing initial pressure iteration #1  Timestep 0 starts with TIME = 0 DT = 7.656411368e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.09220710376 MLMG: Initial residual (resid0) = 0.09220710376 MLMG: Final Iter. 9 resid, resid/bnorm = 2.455981718e-12, 2.663549355e-11 MLMG: Timers: Solve = 0.05398987 Iter = 0.050739723 Bottom = 0.000198107 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.09220710376 MLMG: Initial residual (resid0) = 0.0004260702441 MLMG: Final Iter. 7 resid, resid/bnorm = 8.282632203e-12, 8.982640019e-11 MLMG: Timers: Solve = 0.042036895 Iter = 0.039218969 Bottom = 0.000173862 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 2405121.991 MLMG: Initial residual (resid0) = 2405121.991 MLMG: Final Iter. 7 resid, resid/bnorm = 7.503622328e-07, 3.119851033e-13 MLMG: Timers: Solve = 0.097083025 Iter = 0.094265913 Bottom = 0.000123889 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.16633128 seconds MAC Proj :0.108343738 seconds Nodal Proj :0.112523277 seconds Reactions :0.013674006 seconds Misc :0.025164351 seconds Base State :5.9e-08 seconds Time to advance time step: 0.432774504 Doing initial pressure iteration #2  Timestep 0 starts with TIME = 0 DT = 7.656411368e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.04596608176 MLMG: Initial residual (resid0) = 0.04596608176 MLMG: Final Iter. 8 resid, resid/bnorm = 6.8717948e-13, 1.494970756e-11 MLMG: Timers: Solve = 0.048341767 Iter = 0.04545675 Bottom = 0.000191365 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.04596608176 MLMG: Initial residual (resid0) = 0.0002836050058 MLMG: Final Iter. 8 resid, resid/bnorm = 1.280322311e-12, 2.785363168e-11 MLMG: Timers: Solve = 0.048070527 Iter = 0.04526332 Bottom = 0.000179177 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 386248.2245 MLMG: Initial residual (resid0) = 386248.2245 MLMG: Final Iter. 6 resid, resid/bnorm = 1.884473022e-08, 4.878916983e-14 MLMG: Timers: Solve = 0.083427772 Iter = 0.080455691 Bottom = 0.000112826 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.164790225 seconds MAC Proj :0.107906443 seconds Nodal Proj :0.099044518 seconds Reactions :0.013344618 seconds Misc :0.022683382 seconds Base State :3.599999998e-08 seconds Time to advance time step: 0.408093816 Doing initial pressure iteration #3  Timestep 0 starts with TIME = 0 DT = 7.656411368e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.0433861321 MLMG: Initial residual (resid0) = 0.0433861321 MLMG: Final Iter. 8 resid, resid/bnorm = 8.726436999e-13, 2.011342468e-11 MLMG: Timers: Solve = 0.048146458 Iter = 0.045335926 Bottom = 0.000193293 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.0433861321 MLMG: Initial residual (resid0) = 0.0002834183987 MLMG: Final Iter. 8 resid, resid/bnorm = 1.283079112e-12, 2.957348465e-11 MLMG: Timers: Solve = 0.051072792 Iter = 0.047562262 Bottom = 0.000179171 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 172517.3153 MLMG: Initial residual (resid0) = 172517.3153 MLMG: Final Iter. 5 resid, resid/bnorm = 3.028708306e-08, 1.755596707e-13 MLMG: Timers: Solve = 0.069711499 Iter = 0.066872006 Bottom = 9.0658e-05 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.165764289 seconds MAC Proj :0.111061383 seconds Nodal Proj :0.085269918 seconds Reactions :0.013546938 seconds Misc :0.022339692 seconds Base State :3.100000012e-08 seconds Time to advance time step: 0.398303623 Writing plotfile 0 after all initialization Time to write plotfile: 0.066571447 Beginning main evolution  Timestep 1 starts with TIME = 0 DT = 7.656411368e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.04275829807 MLMG: Initial residual (resid0) = 0.04275829807 MLMG: Final Iter. 8 resid, resid/bnorm = 9.257219557e-13, 2.165011232e-11 MLMG: Timers: Solve = 0.048107664 Iter = 0.045272395 Bottom = 0.000193072 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.04275829807 MLMG: Initial residual (resid0) = 0.0002833664647 MLMG: Final Iter. 8 resid, resid/bnorm = 1.283870362e-12, 3.002622696e-11 MLMG: Timers: Solve = 0.048498308 Iter = 0.045684284 Bottom = 0.000174691 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 0.1841463088 MLMG: Initial residual (resid0) = 0.1841463088 MLMG: Final Iter. 7 resid, resid/bnorm = 5.730832475e-14, 3.112108253e-13 MLMG: Timers: Solve = 0.09668007 Iter = 0.093828546 Bottom = 0.000126596 Done calling nodal solver Timestep 1 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.166500473 seconds MAC Proj :0.108131132 seconds Nodal Proj :0.112642872 seconds Reactions :0.013246512 seconds Misc :0.022527115 seconds Base State :3.199999998e-08 seconds Time to advance time step: 0.429959243 Call to estdt for level 0 gives dt_lev = 3.083908518e-05 Minimum estdt over all levels = 3.083908518e-05 Call to estdt at beginning of step 2 gives dt =3.083908518e-05 dt_growth factor limits the new dt = 8.422052505e-08  Timestep 2 starts with TIME = 7.656411368e-08 DT = 8.422052505e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.1301715907 MLMG: Initial residual (resid0) = 0.1301715907 MLMG: Final Iter. 8 resid, resid/bnorm = 4.055911531e-12, 3.115819288e-11 MLMG: Timers: Solve = 0.047866569 Iter = 0.045041162 Bottom = 0.000176604 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.1301715907 MLMG: Initial residual (resid0) = 0.0008811634353 MLMG: Final Iter. 8 resid, resid/bnorm = 4.005722837e-12, 3.077263491e-11 MLMG: Timers: Solve = 0.04783159 Iter = 0.045038751 Bottom = 0.000175149 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 0.221410349 MLMG: Initial residual (resid0) = 0.221410349 MLMG: Final Iter. 7 resid, resid/bnorm = 6.311617895e-14, 2.850642674e-13 MLMG: Timers: Solve = 0.096572119 Iter = 0.093741527 Bottom = 0.000129424 Done calling nodal solver Timestep 2 ends with TIME = 1.607846387e-07 DT = 8.422052505e-08 Timing summary: Advection :0.174318469 seconds MAC Proj :0.107204911 seconds Nodal Proj :0.112692207 seconds Reactions :0.013196757 seconds Misc :0.022349344 seconds Base State :3.800000004e-08 seconds Time to advance time step: 0.430053982 Call to estdt for level 0 gives dt_lev = 3.124833366e-05 Minimum estdt over all levels = 3.124833366e-05 Call to estdt at beginning of step 3 gives dt =3.124833366e-05 dt_growth factor limits the new dt = 9.264257755e-08  Timestep 3 starts with TIME = 1.607846387e-07 DT = 9.264257755e-08 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.2067293511 MLMG: Initial residual (resid0) = 0.2067293511 MLMG: Final Iter. 8 resid, resid/bnorm = 7.954215845e-12, 3.847647081e-11 MLMG: Timers: Solve = 0.048176382 Iter = 0.045339138 Bottom = 0.000199832 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.2067293511 MLMG: Initial residual (resid0) = 0.001538514331 MLMG: Final Iter. 8 resid, resid/bnorm = 7.007030806e-12, 3.38947071e-11 MLMG: Timers: Solve = 0.048313979 Iter = 0.04527932 Bottom = 0.000170119 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 0.256751109 MLMG: Initial residual (resid0) = 0.256751109 MLMG: Final Iter. 7 resid, resid/bnorm = 6.90975055e-14, 2.691225202e-13 MLMG: Timers: Solve = 0.096281974 Iter = 0.093398169 Bottom = 0.000124411 Done calling nodal solver Timestep 3 ends with TIME = 2.534272163e-07 DT = 9.264257755e-08 Timing summary: Advection :0.175408679 seconds MAC Proj :0.107956646 seconds Nodal Proj :0.112252014 seconds Reactions :0.013256761 seconds Misc :0.022362953 seconds Base State :3.800000004e-08 seconds Time to advance time step: 0.431541338 Call to estdt for level 0 gives dt_lev = 3.131117301e-05 Minimum estdt over all levels = 3.131117301e-05 Call to estdt at beginning of step 4 gives dt =3.131117301e-05 dt_growth factor limits the new dt = 1.019068353e-07  Timestep 4 starts with TIME = 2.534272163e-07 DT = 1.019068353e-07 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.2913498304 MLMG: Initial residual (resid0) = 0.2913498304 MLMG: Final Iter. 8 resid, resid/bnorm = 1.074733965e-11, 3.688809303e-11 MLMG: Timers: Solve = 0.047880682 Iter = 0.045059695 Bottom = 0.000182725 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.2913498304 MLMG: Initial residual (resid0) = 0.002262934381 MLMG: Final Iter. 8 resid, resid/bnorm = 1.028622109e-11, 3.530539585e-11 MLMG: Timers: Solve = 0.048359748 Iter = 0.045536613 Bottom = 0.000176974 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 0.2902215403 MLMG: Initial residual (resid0) = 0.2902215403 MLMG: Final Iter. 7 resid, resid/bnorm = 7.567557692e-14, 2.607510691e-13 MLMG: Timers: Solve = 0.096558106 Iter = 0.093709566 Bottom = 0.000124619 Done calling nodal solver Timestep 4 ends with TIME = 3.553340516e-07 DT = 1.019068353e-07 Timing summary: Advection :0.174759919 seconds MAC Proj :0.107750667 seconds Nodal Proj :0.112457187 seconds Reactions :0.013502077 seconds Misc :0.022376142 seconds Base State :3.300000007e-08 seconds Time to advance time step: 0.431153902 Call to estdt for level 0 gives dt_lev = 3.131144172e-05 Minimum estdt over all levels = 3.131144172e-05 Call to estdt at beginning of step 5 gives dt =3.131144172e-05 dt_growth factor limits the new dt = 1.120975188e-07  Timestep 5 starts with TIME = 3.553340516e-07 DT = 1.120975188e-07 Cell Count: Level 0, 262144 cells <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 0.3837935459 MLMG: Initial residual (resid0) = 0.3837935459 MLMG: Final Iter. 8 resid, resid/bnorm = 1.405302524e-11, 3.661610621e-11 MLMG: Timers: Solve = 0.047744001 Iter = 0.04478695 Bottom = 0.000182936 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 0.3837935459 MLMG: Initial residual (resid0) = 0.003059812973 MLMG: Final Iter. 8 resid, resid/bnorm = 1.389730128e-11, 3.621035692e-11 MLMG: Timers: Solve = 0.048222185 Iter = 0.045391072 Bottom = 0.00018418 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 0.3244360328 MLMG: Initial residual (resid0) = 0.3244360328 MLMG: Final Iter. 7 resid, resid/bnorm = 8.271161533e-14, 2.549396706e-13 MLMG: Timers: Solve = 0.096865609 Iter = 0.09403062 Bottom = 0.000127046 Done calling nodal solver Timestep 5 ends with TIME = 4.674315704e-07 DT = 1.120975188e-07 Timing summary: Advection :0.175075395 seconds MAC Proj :0.10751296 seconds Nodal Proj :0.112873966 seconds Reactions :0.013330619 seconds Misc :0.022679794 seconds Base State :3.500000023e-08 seconds Time to advance time step: 0.431790772 Writing plotfile 5 Time to write plotfile: 0.067142446 Total Time: 3.876481054 Unused ParmParse Variables: [TOP]::amr.check_file(nvals = 1) :: [double_bubble_2d_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] [TOP]::amr.plot_file(nvals = 1) :: [double_bubble_2d_plt] [TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2] AMReX (26.07-54-gaa5b8cc69640) finalized