Initializing AMReX (26.09-136-g53fb957f5e13)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.09-136-g53fb957f5e13) 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.069191941 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002896694 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile double_bubble_2d_pltafter_InitProj after InitProj Time to write plotfile: 0.064788911 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.002800587 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.002796962 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.002828034 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.065985347 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.054794059 Iter = 0.051448864 Bottom = 0.000205289 <<< 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.043489498 Iter = 0.040470696 Bottom = 0.000170453 <<< 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.096606529 Iter = 0.093771442 Bottom = 0.000125818 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.166068313 seconds MAC Proj :0.110403628 seconds Nodal Proj :0.112267556 seconds Reactions :0.013667888 seconds Misc :0.025141194 seconds Base State :3.999999998e-08 seconds Time to advance time step: 0.434159742 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.048781014 Iter = 0.045961527 Bottom = 0.000199006 <<< 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.049043612 Iter = 0.046062953 Bottom = 0.000173927 <<< 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.082817998 Iter = 0.079798439 Bottom = 0.0001147 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.164017542 seconds MAC Proj :0.109463623 seconds Nodal Proj :0.098499136 seconds Reactions :0.013757954 seconds Misc :0.022407754 seconds Base State :3.100000012e-08 seconds Time to advance time step: 0.408450214 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.049498214 Iter = 0.046598717 Bottom = 0.00018766 <<< 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.04917341 Iter = 0.046301604 Bottom = 0.000167236 <<< 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.069784743 Iter = 0.066912243 Bottom = 9.3277e-05 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.16531708 seconds MAC Proj :0.110213809 seconds Nodal Proj :0.085726492 seconds Reactions :0.013194556 seconds Misc :0.022335839 seconds Base State :5.700000028e-08 seconds Time to advance time step: 0.397112219 Writing plotfile 0 after all initialization Time to write plotfile: 0.066215626 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.048931921 Iter = 0.046086324 Bottom = 0.000181235 <<< 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.049192344 Iter = 0.046362569 Bottom = 0.000173511 <<< 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.096292285 Iter = 0.093448932 Bottom = 0.000125381 Done calling nodal solver Timestep 1 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.165513778 seconds MAC Proj :0.109758504 seconds Nodal Proj :0.112211774 seconds Reactions :0.013185268 seconds Misc :0.022528107 seconds Base State :4.899999984e-08 seconds Time to advance time step: 0.430004638 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.048852763 Iter = 0.045983383 Bottom = 0.000184502 <<< 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.048937989 Iter = 0.046074523 Bottom = 0.000177612 <<< 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.096189853 Iter = 0.093326549 Bottom = 0.000130706 Done calling nodal solver Timestep 2 ends with TIME = 1.607846387e-07 DT = 8.422052505e-08 Timing summary: Advection :0.174247134 seconds MAC Proj :0.109541508 seconds Nodal Proj :0.112242206 seconds Reactions :0.013204346 seconds Misc :0.022357857 seconds Base State :5.400000003e-08 seconds Time to advance time step: 0.431877849 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.048727914 Iter = 0.045902108 Bottom = 0.00018143 <<< 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.048993267 Iter = 0.045889501 Bottom = 0.000181306 <<< 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.096172917 Iter = 0.093110933 Bottom = 0.000138198 Done calling nodal solver Timestep 3 ends with TIME = 2.534272163e-07 DT = 9.264257755e-08 Timing summary: Advection :0.17492525 seconds MAC Proj :0.109525693 seconds Nodal Proj :0.112090435 seconds Reactions :0.01324321 seconds Misc :0.022360388 seconds Base State :5.300000039e-08 seconds Time to advance time step: 0.432450102 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.048978119 Iter = 0.046124169 Bottom = 0.000177542 <<< 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.048949248 Iter = 0.046088648 Bottom = 0.000176097 <<< 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.096234321 Iter = 0.093403672 Bottom = 0.000134075 Done calling nodal solver Timestep 4 ends with TIME = 3.553340516e-07 DT = 1.019068353e-07 Timing summary: Advection :0.17502209 seconds MAC Proj :0.109561265 seconds Nodal Proj :0.112095996 seconds Reactions :0.013560088 seconds Misc :0.022383289 seconds Base State :6.100000016e-08 seconds Time to advance time step: 0.432921029 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.04880773 Iter = 0.045961327 Bottom = 0.000175071 <<< 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.048764298 Iter = 0.045906277 Bottom = 0.000172403 <<< 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.096277091 Iter = 0.093448279 Bottom = 0.000130408 Done calling nodal solver Timestep 5 ends with TIME = 4.674315704e-07 DT = 1.120975188e-07 Timing summary: Advection :0.175025167 seconds MAC Proj :0.109482002 seconds Nodal Proj :0.112300179 seconds Reactions :0.013258692 seconds Misc :0.022355256 seconds Base State :5.999999964e-08 seconds Time to advance time step: 0.432724099 Writing plotfile 5 Time to write plotfile: 0.066508766 Total Time: 3.887247459 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.09-136-g53fb957f5e13) finalized