Initializing AMReX (26.09-13-gb105892be1b2)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.09-13-gb105892be1b2) 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.069055667 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.003016417 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile double_bubble_2d_pltafter_InitProj after InitProj Time to write plotfile: 0.065126094 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.002849176 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.002804708 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.002846476 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.06686585 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.054035379 Iter = 0.050462912 Bottom = 0.000200461 <<< 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.04231943 Iter = 0.039443146 Bottom = 0.000173788 <<< 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.194466889e-07, 2.991310593e-13 MLMG: Timers: Solve = 0.0954119 Iter = 0.092504001 Bottom = 0.000127142 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.16388759 seconds MAC Proj :0.108210329 seconds Nodal Proj :0.110996168 seconds Reactions :0.013372198 seconds Misc :0.025196046 seconds Base State :5.5e-08 seconds Time to advance time step: 0.428386933 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.871759563e-13, 1.49496309e-11 MLMG: Timers: Solve = 0.047659897 Iter = 0.044851542 Bottom = 0.000201051 <<< 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.0002836044351 MLMG: Final Iter. 8 resid, resid/bnorm = 1.280322636e-12, 2.785363875e-11 MLMG: Timers: Solve = 0.048265623 Iter = 0.045399034 Bottom = 0.000179587 <<< 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.909938874e-08, 4.944848294e-14 MLMG: Timers: Solve = 0.081982009 Iter = 0.07914807 Bottom = 0.000112722 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.164301565 seconds MAC Proj :0.107521452 seconds Nodal Proj :0.097725116 seconds Reactions :0.013202866 seconds Misc :0.022526255 seconds Base State :3.199999998e-08 seconds Time to advance time step: 0.405555333 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.726384415e-13, 2.011330348e-11 MLMG: Timers: Solve = 0.047242833 Iter = 0.044439804 Bottom = 0.00018861 <<< 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.0002834175439 MLMG: Final Iter. 8 resid, resid/bnorm = 1.283079816e-12, 2.95735009e-11 MLMG: Timers: Solve = 0.047519103 Iter = 0.044683124 Bottom = 0.000186461 <<< 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 = 2.992373993e-08, 1.734535451e-13 MLMG: Timers: Solve = 0.068980572 Iter = 0.066108482 Bottom = 9.628e-05 Done calling nodal solver Timestep 0 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.16414096 seconds MAC Proj :0.10635916 seconds Nodal Proj :0.08471393 seconds Reactions :0.013194613 seconds Misc :0.022339906 seconds Base State :5.800000014e-08 seconds Time to advance time step: 0.391033927 Writing plotfile 0 after all initialization Time to write plotfile: 0.066138295 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.257211425e-13, 2.16500933e-11 MLMG: Timers: Solve = 0.047452554 Iter = 0.044457924 Bottom = 0.000194953 <<< 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.0002833663214 MLMG: Final Iter. 8 resid, resid/bnorm = 1.283870579e-12, 3.002623204e-11 MLMG: Timers: Solve = 0.046904174 Iter = 0.044058545 Bottom = 0.000181313 <<< 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.567421524e-14, 3.023368515e-13 MLMG: Timers: Solve = 0.09513286 Iter = 0.092207251 Bottom = 0.000133423 Done calling nodal solver Timestep 1 ends with TIME = 7.656411368e-08 DT = 7.656411368e-08 Timing summary: Advection :0.164534022 seconds MAC Proj :0.106158864 seconds Nodal Proj :0.110907026 seconds Reactions :0.013235848 seconds Misc :0.022525335 seconds Base State :6.099999994e-08 seconds Time to advance time step: 0.424762561 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.048280976 Iter = 0.045434277 Bottom = 0.000180125 <<< 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.047871432 Iter = 0.045001383 Bottom = 0.000171771 <<< 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.050715484e-14, 2.732806082e-13 MLMG: Timers: Solve = 0.096478829 Iter = 0.093490443 Bottom = 0.000129491 Done calling nodal solver Timestep 2 ends with TIME = 1.607846387e-07 DT = 8.422052505e-08 Timing summary: Advection :0.173588257 seconds MAC Proj :0.108110534 seconds Nodal Proj :0.112294838 seconds Reactions :0.013225076 seconds Misc :0.022353742 seconds Base State :5.9e-08 seconds Time to advance time step: 0.429863478 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.954215628e-12, 3.847646977e-11 MLMG: Timers: Solve = 0.047372754 Iter = 0.044557063 Bottom = 0.000171874 <<< 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.007030373e-12, 3.389470501e-11 MLMG: Timers: Solve = 0.047809406 Iter = 0.044962649 Bottom = 0.000181211 <<< 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.722400414e-14, 2.618255649e-13 MLMG: Timers: Solve = 0.094963391 Iter = 0.092127803 Bottom = 0.000134245 Done calling nodal solver Timestep 3 ends with TIME = 2.534272163e-07 DT = 9.264257755e-08 Timing summary: Advection :0.173630707 seconds MAC Proj :0.106766906 seconds Nodal Proj :0.110846616 seconds Reactions :0.013205651 seconds Misc :0.022613369 seconds Base State :4.399999964e-08 seconds Time to advance time step: 0.42734917 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.074734052e-11, 3.688809601e-11 MLMG: Timers: Solve = 0.047289058 Iter = 0.044394254 Bottom = 0.000181055 <<< 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.028621936e-11, 3.530538989e-11 MLMG: Timers: Solve = 0.047549524 Iter = 0.044717705 Bottom = 0.000179948 <<< 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.438494265e-14, 2.563040034e-13 MLMG: Timers: Solve = 0.09518381 Iter = 0.092347175 Bottom = 0.000133352 Done calling nodal solver Timestep 4 ends with TIME = 3.553340516e-07 DT = 1.019068353e-07 Timing summary: Advection :0.175677611 seconds MAC Proj :0.106377857 seconds Nodal Proj :0.11110187 seconds Reactions :0.013507728 seconds Misc :0.022334071 seconds Base State :5.600000019e-08 seconds Time to advance time step: 0.429306091 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.40530248e-11, 3.661610508e-11 MLMG: Timers: Solve = 0.04770067 Iter = 0.044782361 Bottom = 0.000179137 <<< 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.389729998e-11, 3.621035353e-11 MLMG: Timers: Solve = 0.047981938 Iter = 0.045154024 Bottom = 0.000175105 <<< 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.146261443e-14, 2.510899105e-13 MLMG: Timers: Solve = 0.09523408 Iter = 0.092381616 Bottom = 0.000129082 Done calling nodal solver Timestep 5 ends with TIME = 4.674315704e-07 DT = 1.120975188e-07 Timing summary: Advection :0.174776488 seconds MAC Proj :0.107238294 seconds Nodal Proj :0.111147233 seconds Reactions :0.013457883 seconds Misc :0.022356119 seconds Base State :5.800000036e-08 seconds Time to advance time step: 0.429274019 Writing plotfile 5 Time to write plotfile: 0.065905972 Total Time: 3.8463629 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-13-gb105892be1b2) finalized