Initializing AMReX (26.04-109-g60df5ae33414)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.04-109-g60df5ae33414) 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.072967499 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.002831384 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile double_bubble_2d_pltafter_InitProj after InitProj Time to write plotfile: 0.065028014 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.002775809 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.002907431 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.002802709 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.066291721 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.053963538 Iter = 0.050588176 Bottom = 0.000204712 <<< 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.042427724 Iter = 0.039608409 Bottom = 0.000168136 <<< 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.096710711 Iter = 0.093814005 Bottom = 0.000123233 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.165515488 seconds MAC Proj :0.108647263 seconds Nodal Proj :0.112226044 seconds Reactions :0.013390386 seconds Misc :0.025250788 seconds Base State :3.100000001e-08 seconds Time to advance time step: 0.431605797 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.04788545 Iter = 0.044921217 Bottom = 0.000191746 <<< 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.048047855 Iter = 0.045217034 Bottom = 0.000179058 <<< 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.082985881 Iter = 0.0801535 Bottom = 0.000113229 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.163668085 seconds MAC Proj :0.107648743 seconds Nodal Proj :0.098861155 seconds Reactions :0.013233285 seconds Misc :0.022575992 seconds Base State :3.199999998e-08 seconds Time to advance time step: 0.406322262 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.047857048 Iter = 0.045005627 Bottom = 0.000196644 <<< 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.048233297 Iter = 0.045413823 Bottom = 0.000182723 <<< 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.069834047 Iter = 0.066960706 Bottom = 9.5094e-05 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.165027251 seconds MAC Proj :0.107699532 seconds Nodal Proj :0.085571266 seconds Reactions :0.01328156 seconds Misc :0.022319038 seconds Base State :3.500000001e-08 seconds Time to advance time step: 0.394279246 Writing plotfile 0 after all initialization Time to write plotfile: 0.066563783 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.047590828 Iter = 0.044770833 Bottom = 0.000189985 <<< 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.047725313 Iter = 0.044914564 Bottom = 0.000181227 <<< 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.096302599 Iter = 0.093455896 Bottom = 0.000127916 Done calling nodal solver Timestep 1 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.165083589 seconds MAC Proj :0.106987033 seconds Nodal Proj :0.112363412 seconds Reactions :0.013227696 seconds Misc :0.022844442 seconds Base State :3.199999998e-08 seconds Time to advance time step: 0.427399096 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.047682231 Iter = 0.044691651 Bottom = 0.000180134 <<< 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.048027168 Iter = 0.04520559 Bottom = 0.00018163 <<< 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.096471274 Iter = 0.093614329 Bottom = 0.000131547 Done calling nodal solver Timestep 2 ends with TIME = 1.607846387e-07 DT = 8.422052505e-08 Timing summary: Advection :0.173535262 seconds MAC Proj :0.107373341 seconds Nodal Proj :0.112376665 seconds Reactions :0.013553801 seconds Misc :0.022329464 seconds Base State :3.400000015e-08 seconds Time to advance time step: 0.429480485 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.048031647 Iter = 0.045193103 Bottom = 0.000181579 <<< 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.048149048 Iter = 0.045301999 Bottom = 0.000182693 <<< 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.096711787 Iter = 0.093849132 Bottom = 0.000133266 Done calling nodal solver Timestep 3 ends with TIME = 2.534272163e-07 DT = 9.264257755e-08 Timing summary: Advection :0.174320598 seconds MAC Proj :0.107897367 seconds Nodal Proj :0.112753317 seconds Reactions :0.013228089 seconds Misc :0.022489266 seconds Base State :3.199999998e-08 seconds Time to advance time step: 0.430984264 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.073147889 Iter = 0.070323034 Bottom = 0.000184247 <<< 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.051502246 Iter = 0.047431267 Bottom = 0.000187783 <<< 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.097795657 Iter = 0.0949842 Bottom = 0.000155699 Done calling nodal solver Timestep 4 ends with TIME = 3.553340516e-07 DT = 1.019068353e-07 Timing summary: Advection :0.208312225 seconds MAC Proj :0.138178182 seconds Nodal Proj :0.114126383 seconds Reactions :0.014461701 seconds Misc :0.023026329 seconds Base State :3.09999999e-08 seconds Time to advance time step: 0.498468194 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.048237522 Iter = 0.045343289 Bottom = 0.000221041 <<< 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.047848085 Iter = 0.044951787 Bottom = 0.000182711 <<< 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.096581333 Iter = 0.093748708 Bottom = 0.000129705 Done calling nodal solver Timestep 5 ends with TIME = 4.674315704e-07 DT = 1.120975188e-07 Timing summary: Advection :0.174332247 seconds MAC Proj :0.107900619 seconds Nodal Proj :0.112572665 seconds Reactions :0.013281257 seconds Misc :0.022501768 seconds Base State :3.899999967e-08 seconds Time to advance time step: 0.430893545 Writing plotfile 5 Time to write plotfile: 0.066662715 Total Time: 3.940222511 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.04-109-g60df5ae33414) finalized