Initializing AMReX (26.05-63-gc2eb6db7ee79)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.05-63-gc2eb6db7ee79) 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.068741462 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.003104807 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile double_bubble_2d_pltafter_InitProj after InitProj Time to write plotfile: 0.065010675 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.00280057 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.002802456 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.002940659 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.067540354 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.054583527 Iter = 0.051189578 Bottom = 0.000195749 <<< 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.042917297 Iter = 0.040058903 Bottom = 0.000168452 <<< 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.096604714 Iter = 0.093782182 Bottom = 0.000123996 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.16449033 seconds MAC Proj :0.109557047 seconds Nodal Proj :0.112104678 seconds Reactions :0.013296877 seconds Misc :0.025017357 seconds Base State :6.900000005e-08 seconds Time to advance time step: 0.431394441 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.048961805 Iter = 0.046126091 Bottom = 0.000192684 <<< 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.048532312 Iter = 0.045695827 Bottom = 0.000168691 <<< 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.083165947 Iter = 0.080309205 Bottom = 0.000106658 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.163493555 seconds MAC Proj :0.109003846 seconds Nodal Proj :0.098752505 seconds Reactions :0.013211737 seconds Misc :0.022310408 seconds Base State :3.100000001e-08 seconds Time to advance time step: 0.407084347 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.048386761 Iter = 0.045568374 Bottom = 0.00018085 <<< 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.048217274 Iter = 0.045363588 Bottom = 0.000175547 <<< 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.069888714 Iter = 0.067014462 Bottom = 8.983e-05 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.163669562 seconds MAC Proj :0.108104392 seconds Nodal Proj :0.08543909 seconds Reactions :0.013486645 seconds Misc :0.02230228 seconds Base State :3.200000021e-08 seconds Time to advance time step: 0.393313236 Writing plotfile 0 after all initialization Time to write plotfile: 0.066354121 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.04800742 Iter = 0.045191091 Bottom = 0.000182121 <<< 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.048310395 Iter = 0.045460674 Bottom = 0.000175785 <<< 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.097046822 Iter = 0.094194839 Bottom = 0.00012947 Done calling nodal solver Timestep 1 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.164252317 seconds MAC Proj :0.107797683 seconds Nodal Proj :0.113001399 seconds Reactions :0.01321263 seconds Misc :0.022534079 seconds Base State :5.500000011e-08 seconds Time to advance time step: 0.427886092 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.04824289 Iter = 0.045429957 Bottom = 0.00017171 <<< 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.048596649 Iter = 0.045738718 Bottom = 0.000164712 <<< 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.096742589 Iter = 0.09388679 Bottom = 0.000123919 Done calling nodal solver Timestep 2 ends with TIME = 1.607846387e-07 DT = 8.422052505e-08 Timing summary: Advection :0.173316232 seconds MAC Proj :0.10841292 seconds Nodal Proj :0.11263866 seconds Reactions :0.013266517 seconds Misc :0.022351169 seconds Base State :3.39999997e-08 seconds Time to advance time step: 0.430281373 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.048414904 Iter = 0.045611208 Bottom = 0.000176238 <<< 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.048160971 Iter = 0.045275276 Bottom = 0.000181878 <<< 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.096444892 Iter = 0.093478076 Bottom = 0.00012882 Done calling nodal solver Timestep 3 ends with TIME = 2.534272163e-07 DT = 9.264257755e-08 Timing summary: Advection :0.173982617 seconds MAC Proj :0.108054465 seconds Nodal Proj :0.112337613 seconds Reactions :0.013519827 seconds Misc :0.022326992 seconds Base State :5.100000022e-08 seconds Time to advance time step: 0.430496121 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.048423905 Iter = 0.045580803 Bottom = 0.000174896 <<< 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.047783402 Iter = 0.044868698 Bottom = 0.000173246 <<< 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.096824749 Iter = 0.093841521 Bottom = 0.000129694 Done calling nodal solver Timestep 4 ends with TIME = 3.553340516e-07 DT = 1.019068353e-07 Timing summary: Advection :0.173361653 seconds MAC Proj :0.107858204 seconds Nodal Proj :0.112780303 seconds Reactions :0.013238228 seconds Misc :0.022343675 seconds Base State :3.300000007e-08 seconds Time to advance time step: 0.429882912 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.048228134 Iter = 0.045412269 Bottom = 0.000171754 <<< 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.048374878 Iter = 0.045556834 Bottom = 0.000169578 <<< 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.096647724 Iter = 0.093810062 Bottom = 0.000128362 Done calling nodal solver Timestep 5 ends with TIME = 4.674315704e-07 DT = 1.120975188e-07 Timing summary: Advection :0.174817182 seconds MAC Proj :0.108114063 seconds Nodal Proj :0.11256861 seconds Reactions :0.013518457 seconds Misc :0.022343172 seconds Base State :4.00000002e-08 seconds Time to advance time step: 0.431653871 Writing plotfile 5 Time to write plotfile: 0.066923841 Total Time: 3.86535522 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.05-63-gc2eb6db7ee79) finalized