Initializing AMReX (24.05-32-g28b010126a1b)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05-32-g28b010126a1b) initialized Starting run at 07:09:09 UTC on 2024-05-29. Successfully read inputs file ... Castro git describe: 24.05-15-g7893b5c6d AMReX git describe: 24.05-32-g28b010126 Microphysics git describe: 24.05-22-g7b1458be reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... Creating primary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Creating secondary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Generated initial model for primary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Generated initial model for secondary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Castro::numpts_1d at level 0 is 294 Setting the gravity type to PoissonGrav Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.179037758 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001051804 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010613314 Gravity::multilevel_solve_for_new_phi() time = 0.100113761 Gravity::fill_multipole_BCs() time = 0.010582224 ... norm of composite residual at level 0 1.193955541e-11 TIME= 0 MASS = 2.549996568e+33 TIME= 0 XMOM = 0 TIME= 0 YMOM = 0 TIME= 0 ZMOM = 0 TIME= 0 ANG MOM X = 0 TIME= 0 ANG MOM Y = 0 TIME= 0 ANG MOM Z = -3.115378115e+34 TIME= 0 RHO*e = 1.325915365e+50 TIME= 0 RHO*K = 3.09960185e+49 TIME= 0 RHO*E = 1.63587555e+50 TIME= 0 RHO*PHI = -5.434813569e+50 TIME= 0 TOTAL ENERGY= -1.081531234e+50 TIME= 0 CENTER OF MASS X-LOC = 316439254.4 TIME= 0 CENTER OF MASS X-VEL = 0 TIME= 0 CENTER OF MASS Y-LOC = -4.854673355e-07 TIME= 0 CENTER OF MASS Y-VEL = 0 TIME= 0 CENTER OF MASS Z-LOC = 0 TIME= 0 CENTER OF MASS Z-VEL = 0 TIME= 0 MAXIMUM TEMPERATURE = 10000000 TIME= 0 MAXIMUM DENSITY = 4476474 TIME= 0 MAXIMUM T_S / T_E = 0 INITIAL GRIDS Level 0 32 grids 32768 cells 100 % of domain smallest grid: 32 x 32 biggest grid: 32 x 32 PLOTFILE: file = wdmerger_retry_plt00000 Write plotfile time = 0.102371904 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.004309155 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 Castro::construct_old_gravity() time = 0.002231546 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -5.111789233e+47 ener added: 3.959219145e+48 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.174137392e+41 ymom added: 6.491992098e+28 zmom added: 0 eint added: 0 ener added: 6.931398274e+47 Castro::do_old_sources() time = 0.066001654 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.003637875 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -5.024842118e+47 ener added: 6.151153052e+47 Contributions to the state from old-time sources: mass added: 0 xmom added: -5.870686959e+40 ymom added: 3.245996049e+28 zmom added: 0 eint added: 0 ener added: 3.465699137e+47 Castro::do_old_sources() time = 0.065006839 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Castro::construct_old_gravity() time = 0.003895985 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -2.93534348e+40 ymom added: 1.622998024e+28 zmom added: 0 eint added: 0 ener added: 1.732849568e+47 Castro::do_old_sources() time = 0.062243576 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 78, 0 Invalid density = -0.1060744788 at index 20, 177, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.0625 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.03125   Beginning subcycle 1 starting at time 0 with dt = 0.03125 Estimated number of subcycles remaining: 8 Castro::construct_old_gravity() time = 0.004341409 on level 0 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46767174e+40 ymom added: 8.114990122e+27 zmom added: 0 eint added: 0 ener added: 8.664247842e+46 Castro::do_old_sources() time = 0.059732631 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.314382477 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014514157 Gravity::solve_for_phi() time = 0.090336399 on level 0 Castro::construct_new_gravity() time = 0.092394812 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.674192048e+44 ener added: 3.847280946e+42 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.752957884e+35 ymom added: -3.585455572e+27 zmom added: 0 eint added: 0 ener added: 6.209484931e+43 Castro::do_new_sources() time = 0.061203153 on level 0  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 Castro::construct_old_gravity() time = 0.002893828 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.985023871e+46 ener added: 5.45426343e+40 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467706799e+40 ymom added: 9.393353618e+26 zmom added: 0 eint added: 0 ener added: 8.722516755e+46 Castro::do_old_sources() time = 0.070913536 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.985643901e+46 ener added: 2.77668231e+42 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.331805191 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017877365 Gravity::solve_for_phi() time = 0.088013616 on level 0 Castro::construct_new_gravity() time = 0.090793376 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.270322473e+44 ener added: 4.896470657e+43 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.769111125e+35 ymom added: 3.317462454e+23 zmom added: 0 eint added: 0 ener added: 6.005402878e+43 Castro::do_new_sources() time = 0.060176878 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.334150457e+29 ener added: 1.334150457e+29  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 Castro::construct_old_gravity() time = 0.002247396 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.984358138e+46 ener added: 6.992681121e+42 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467742181e+40 ymom added: 9.469667061e+26 zmom added: 0 eint added: 0 ener added: 8.781433933e+46 Castro::do_old_sources() time = 0.064331589 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.987377036e+46 ener added: 2.185609542e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.337595695 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015209192 Gravity::solve_for_phi() time = 0.091373223 on level 0 Castro::construct_new_gravity() time = 0.093489847 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.025861471e+43 ener added: 1.182920356e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.788465098e+35 ymom added: -1.334068531e+24 zmom added: 0 eint added: 0 ener added: 6.020728567e+43 Castro::do_new_sources() time = 0.060301457 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 8.567892228e+28 ener added: 8.567892228e+28  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 Castro::construct_old_gravity() time = 0.005778798 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983924894e+46 ener added: 1.155580676e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.46777795e+40 ymom added: 9.403931719e+26 zmom added: 0 eint added: 0 ener added: 8.841709465e+46 Castro::do_old_sources() time = 0.058605257 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.985020288e+46 ener added: 2.299661927e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.326136737 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014989104 Gravity::solve_for_phi() time = 0.095975022 on level 0 Castro::construct_new_gravity() time = 0.098018239 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -4.577058949e+43 ener added: 1.487868819e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.809569034e+35 ymom added: 1.794499263e+23 zmom added: 0 eint added: 0 ener added: 4.688674119e+43 Castro::do_new_sources() time = 0.05849146 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.793059048e+29 ener added: 4.793059048e+29  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.50359098e+29 ener added: 2.50359098e+29 Castro::construct_old_gravity() time = 0.002121745 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.98357515e+46 ener added: 1.522753583e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467814142e+40 ymom added: 9.411487506e+26 zmom added: 0 eint added: 0 ener added: 8.90325471e+46 Castro::do_old_sources() time = 0.058330021 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.99240328e+46 ener added: 2.331782913e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.324486578 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014681774 Gravity::solve_for_phi() time = 0.090228375 on level 0 Castro::construct_new_gravity() time = 0.09247893 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.032940322e+43 ener added: 1.61912462e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.835708706e+35 ymom added: 3.356421978e+24 zmom added: 0 eint added: 0 ener added: 2.183890438e+43 Castro::do_new_sources() time = 0.058498748 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.823721814e+28 ener added: 7.823721814e+28  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 Castro::construct_old_gravity() time = 0.003081449 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.983089144e+46 ener added: 2.018395698e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467850856e+40 ymom added: 9.456822224e+26 zmom added: 0 eint added: 0 ener added: 8.965576069e+46 Castro::do_old_sources() time = 0.182690074 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.000655828e+46 ener added: 2.414422253e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.437121514 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015136427 Gravity::solve_for_phi() time = 0.103301163 on level 0 Castro::construct_new_gravity() time = 0.105376349 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -8.829626388e+43 ener added: 1.712554021e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.854387385e+35 ymom added: 3.424896292e+24 zmom added: 0 eint added: 0 ener added: -1.318403855e+43 Castro::do_new_sources() time = 0.058791334 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 1.416634103e+29 ener added: 1.416634103e+29  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.002134909 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982961815e+46 ener added: 2.174330548e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467887944e+40 ymom added: 9.567136705e+26 zmom added: 0 eint added: 0 ener added: 9.028440984e+46 Castro::do_old_sources() time = 0.091208499 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993941122e+46 ener added: 2.684057583e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.326492059 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01500515 Gravity::solve_for_phi() time = 0.084377705 on level 0 Castro::construct_new_gravity() time = 0.086503278 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -1.15291724e+44 ener added: 1.741677821e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.882731577e+35 ymom added: -2.475700629e+24 zmom added: 0 eint added: 0 ener added: -4.024106959e+43 Castro::do_new_sources() time = 0.058971678 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 4.554229645e+29 ener added: 4.554229645e+29  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 2.882423826e+28 ener added: 2.882423826e+28 Castro::construct_old_gravity() time = 0.002136966 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -6.982570067e+46 ener added: 2.527968049e+43 Contributions to the state from old-time sources: mass added: 0 xmom added: -1.467925598e+40 ymom added: 9.539935874e+26 zmom added: 0 eint added: 0 ener added: 9.091984554e+46 Castro::do_old_sources() time = 0.059204661 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 6.993654336e+46 ener added: 2.944609967e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.326793025 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015067701 Gravity::solve_for_phi() time = 0.085067207 on level 0 Castro::construct_new_gravity() time = 0.087217065 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: -8.770686754e+43 ener added: 1.764875055e+44 Contributions to the state from new-time sources: mass added: 0 xmom added: -1.944888844e+35 ymom added: -9.954748546e+24 zmom added: 0 eint added: 0 ener added: -1.240335825e+43 Castro::do_new_sources() time = 0.058988608 on level 0 Contributions to the state from negative energy resets: mass added: 0 xmom added: 0 ymom added: 0 zmom added: 0 eint added: 7.411946981e+29 ener added: 7.411946981e+29  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.004478744977 [Level 0 step 1] Advanced 32768 cells TIME= 0.25 MASS = 2.549996568e+33 TIME= 0.25 XMOM = -4.340682746e+38 TIME= 0.25 YMOM = 1.137357411e+28 TIME= 0.25 ZMOM = 0 TIME= 0.25 ANG MOM X = 0 TIME= 0.25 ANG MOM Y = 0 TIME= 0.25 ANG MOM Z = -1.467966168e+38 TIME= 0.25 RHO*e = 1.325895685e+50 TIME= 0.25 RHO*K = 3.170941956e+49 TIME= 0.25 RHO*E = 1.642989881e+50 TIME= 0.25 RHO*PHI = -5.449013736e+50 TIME= 0.25 TOTAL ENERGY= -1.081516988e+50 TIME= 0.25 CENTER OF MASS X-LOC = 316424366.8 TIME= 0.25 CENTER OF MASS X-VEL = -170223.0819 TIME= 0.25 CENTER OF MASS Y-LOC = 1.213668339e-06 TIME= 0.25 CENTER OF MASS Y-VEL = 4.460231145e-06 TIME= 0.25 CENTER OF MASS Z-LOC = 0 TIME= 0.25 CENTER OF MASS Z-VEL = 0 TIME= 0.25 MAXIMUM TEMPERATURE = 58192817.11 TIME= 0.25 MAXIMUM DENSITY = 4458827.313 TIME= 0.25 MAXIMUM T_S / T_E = 0 [STEP 1] Coarse TimeStep time: 7.376360418 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.046862357 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00395084 seconds Ending run at 07:09:19 UTC on 2024-05-29. Run time = 10.23311916 Run time without initialization = 7.427963009 Average number of zones advanced per microsecond: 0.004 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 10.23 ... 10.23 ... 10.23 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8263 2.8857 2.9971 29.29% Amr::InitAmr() 1 1.857 1.857 1.857 18.15% Castro::computeTemp() 63 0.6954 0.7543 0.8299 8.11% FabArray::ParallelCopy_finish() 734 0.04522 0.4079 0.5949 5.81% amrex::Copy() 523 0.5374 0.5476 0.5742 5.61% Castro::normalize_species() 62 0.44 0.4489 0.4876 4.76% Castro::reset_internal_energy(Fab) 126 0.4292 0.4371 0.4605 4.50% MultiFab::AddProduct() 2633 0.422 0.4288 0.4444 4.34% Castro::initMFs() 1 0.3553 0.361 0.3672 3.59% check_for_negative_density() 11 0.1128 0.2111 0.3524 3.44% MultiFab::Subtract() 125 0.2703 0.2746 0.2889 2.82% Castro::reset_internal_energy(MultiFab) 63 0.01778 0.06705 0.2358 2.30% FabArray::setVal() 3453 0.219 0.2237 0.2326 2.27% MLPoisson::Fsmooth() 3120 0.1595 0.175 0.2113 2.06% Castro::initData() 1 0.003261 0.02874 0.1701 1.66% FabArray::sum() 2633 0.1406 0.143 0.1481 1.45% Castro::do_old_sources() 11 0.0003016 0.03679 0.1445 1.41% Gravity::fill_multipole_BCs() 10 0.1036 0.1216 0.1415 1.38% FillBoundary_finish() 3240 0.07212 0.08802 0.1037 1.01% FabArray::ParallelCopy_nowait() 734 0.06308 0.07169 0.08218 0.80% Castro::enforce_min_density() 62 0.01347 0.0339 0.08168 0.80% Amr::writePlotFile() 2 0.06153 0.06977 0.0711 0.69% FillBoundary_nowait() 3240 0.03857 0.0527 0.06794 0.66% MLPoisson::Fapply() 501 0.04796 0.05096 0.05833 0.57% MLCellLinOp::applyBC() 3630 0.02421 0.03183 0.05293 0.52% AmrLevel::derive() 18 0.0412 0.0417 0.04243 0.41% FabArray::Saxpy() 133 0.03467 0.0355 0.03901 0.38% Castro::construct_old_gravity_source() 11 0.0251 0.02612 0.03682 0.36% Castro::construct_new_gravity_source() 8 0.03267 0.03323 0.03582 0.35% Castro::do_new_sources() 8 0.0001775 0.0055 0.03363 0.33% FabArray::Xpay() 480 0.02568 0.02711 0.03056 0.30% MultiFab::contains_nan() 19 0.02297 0.02338 0.02391 0.23% VisMF::Write(FabArray) 4 0.01162 0.01261 0.02195 0.21% MLMG::addInterpCorrection() 390 0.01685 0.01818 0.02163 0.21% amrex::average_down 390 0.01226 0.01341 0.01628 0.16% StateData::FillBoundary(geom) 62 0 0.007068 0.01367 0.13% Castro::wd_update() 2 0.004403 0.008751 0.01364 0.13% Castro::evaluate_source_change() 144 0.009093 0.009324 0.009858 0.10% MLMG::solve() 9 0.001528 0.006092 0.009692 0.09% MLCellLinOp::defineAuxData() 10 0.007438 0.007794 0.008924 0.09% FabArray::FillBoundary() 3240 0.00682 0.007139 0.008605 0.08% MLMG::ResNormInf() 74 0.002444 0.006077 0.008398 0.08% FabArray::mult() 40 0.007937 0.008137 0.008368 0.08% MLCGSolver::bicgstab 4 0 0.0005036 0.008058 0.08% amrex::Add() 65 0.007192 0.007298 0.008023 0.08% FabArray::norminf() 143 0.005577 0.005891 0.007855 0.08% Castro::construct_new_source() 48 0.007235 0.007376 0.007659 0.07% Castro::construct_old_source() 66 0.005383 0.005594 0.007246 0.07% MLCellLinOp::smooth() 1560 0.00558 0.005794 0.006251 0.06% Amr::writeSmallPlotFile() 2 0.004095 0.005725 0.006022 0.06% MLMG::prepareForSolve() 9 0.004959 0.00519 0.005855 0.06% FillPatchIterator::Initialize 82 0.005472 0.005555 0.005684 0.06% Gravity::actual_solve_with_mlmg() 10 0.003445 0.003646 0.004112 0.04% Gravity::get_old_grav_vector() 11 0.003501 0.003566 0.003639 0.04% Gravity::get_new_grav_vector() 9 0.002986 0.003089 0.00331 0.03% Castro::initialize_do_advance() 11 0.002409 0.002552 0.002729 0.03% Castro::sum_integrated_quantities() 2 0.001752 0.002329 0.002596 0.03% MLCellLinOp::apply() 501 0.0009868 0.001106 0.00235 0.02% Castro::construct_old_gravity() 11 0.0002113 0.001004 0.002322 0.02% FabArray::ParallelCopy() 734 0.001757 0.001836 0.00198 0.02% Gravity::solve_for_phi() 8 0.0003858 0.0005591 0.001499 0.01% Castro::construct_new_gravity() 8 0.0001383 0.0004575 0.001493 0.01% main() 1 0.0009053 0.001054 0.001186 0.01% AmrLevel::FillPatch() 82 0.001097 0.001128 0.001173 0.01% MLCellLinOp::correctionResidual() 390 0.0009063 0.0009369 0.0009974 0.01% MLMG::mgVcycle() 65 0.0007982 0.0008346 0.0008983 0.01% Castro::save_data_for_retry() 3 0.0007367 0.0007533 0.0007668 0.01% MLMG:computeResOfCorrection() 390 0.000696 0.0007217 0.0007635 0.01% Castro::subcycle_advance_ctu() 1 8.79e-05 0.0001326 0.000605 0.01% MLMG::actualBottomSolve() 65 6.693e-05 0.0001001 0.0005571 0.01% AmrMesh::MakeNewGrids() 2 0.0004947 0.0005103 0.0005407 0.01% Castro::Castro() 1 0.0004536 0.0004798 0.0005129 0.01% StateDataPhysBCFunct::() 82 0.0003391 0.0004098 0.0004666 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003591 0.0003736 0.0003989 0.00% Castro::problem_post_timestep() 2 1.02e-05 3.503e-05 0.0003842 0.00% Castro::clean_state() 62 0.0003187 0.0003445 0.0003827 0.00% MLMG::mgVcycle_down::2 65 0.0003043 0.0003345 0.0003746 0.00% MLMG::mgVcycle_down::3 65 0.0002865 0.0003098 0.0003597 0.00% MLMG::mgVcycle_down::4 65 0.00028 0.0002933 0.000355 0.00% MLMG::mgVcycle_down::1 65 0.0003223 0.0003333 0.0003549 0.00% MLMG::mgVcycle_down::0 65 0.0003287 0.0003386 0.0003523 0.00% MLMG::mgVcycle_down::5 65 0.0002856 0.000297 0.0003405 0.00% FillPatchSingleLevel 82 0.0002474 0.0002605 0.0002801 0.00% MLMG::mgVcycle_up::1 65 0.0002262 0.000238 0.0002517 0.00% Castro::do_advance_ctu() 11 0.0002134 0.0002315 0.0002494 0.00% MLMG::mgVcycle_up::0 65 0.0002287 0.0002341 0.0002441 0.00% MLMG::mgVcycle_up::5 65 0.0001976 0.000204 0.0002435 0.00% MLMG::mgVcycle_up::2 65 0.0001914 0.0002135 0.0002365 0.00% MLCellLinOp::solutionResidual() 75 0.0002142 0.0002213 0.0002356 0.00% MLMG::mgVcycle_up::3 65 0.0001883 0.0001987 0.0002237 0.00% MLMG::mgVcycle_up::4 65 0.0001915 0.0001987 0.0002223 0.00% MLMG::oneIter() 65 0.0001944 0.0002006 0.0002094 0.00% Amr::defBaseLevel() 1 9.59e-05 0.0001256 0.0001538 0.00% Castro::retry_advance_ctu() 11 4.958e-05 5.964e-05 0.0001524 0.00% Amr::coarseTimeStep() 1 0.0001121 0.0001365 0.0001415 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001139 0.0001214 0.0001319 0.00% MLMG::computeResidual() 65 0.0001188 0.0001221 0.0001258 0.00% Amr::grid_places() 2 7.444e-05 9.073e-05 0.0001236 0.00% MLMG::mgVcycle_bottom 65 0.0001135 0.0001171 0.0001227 0.00% TagBoxArray::mapPRD 2 0.0001074 0.0001117 0.0001209 0.00% Gravity::test_composite_phi() 1 9.356e-05 9.684e-05 0.0001109 0.00% Castro::initialize_advance() 1 9.642e-05 0.000102 0.0001067 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.173e-05 6.67e-05 8.174e-05 0.00% Castro::check_for_nan() 19 5.106e-05 5.915e-05 8.009e-05 0.00% MLPoisson::define() 10 6.563e-05 7.01e-05 7.862e-05 0.00% Amr::FinalizeInit() 1 8.2e-06 1.246e-05 6.402e-05 0.00% Castro::react_state() 19 4.883e-05 5.252e-05 5.628e-05 0.00% Castro::apply_source_to_state() 19 4.858e-05 5.143e-05 5.466e-05 0.00% Castro::derive() 18 3.989e-05 4.223e-05 4.918e-05 0.00% Castro::create_source_corrector() 8 3.47e-05 3.912e-05 4.433e-05 0.00% Castro::expand_state() 11 2.915e-05 3.285e-05 4.34e-05 0.00% Amr::regrid() 1 1.019e-05 1.251e-05 2.942e-05 0.00% Gravity::actual_multilevel_solve() 1 2.683e-05 2.743e-05 2.79e-05 0.00% MLMG::computeMLResidual() 9 2.05e-05 2.219e-05 2.685e-05 0.00% Amr::timeStep() 1 1.218e-05 1.432e-05 2.026e-05 0.00% Castro::post_init() 1 1.202e-05 1.283e-05 1.738e-05 0.00% Castro::advance() 1 7.009e-06 9.284e-06 1.412e-05 0.00% Amr::InitializeInit() 1 6.252e-06 7.273e-06 1.259e-05 0.00% Castro::post_regrid() 2 1.035e-05 1.097e-05 1.178e-05 0.00% Amr::init() 1 6.326e-06 7.065e-06 9.603e-06 0.00% Amr::bldFineLevels() 1 6.794e-06 7.324e-06 7.927e-06 0.00% Castro::post_timestep() 1 5.73e-06 6.305e-06 7.012e-06 0.00% Amr::initialInit() 1 2.835e-06 3.064e-06 3.447e-06 0.00% Other 5693 0.03612 0.04446 0.0645 0.63% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.23 10.23 10.23 100.00% Amr::coarseTimeStep() 1 7.428 7.428 7.428 72.58% Amr::timeStep() 1 7.376 7.376 7.376 72.08% Castro::advance() 1 7.311 7.314 7.316 71.50% Castro::subcycle_advance_ctu() 1 7.262 7.265 7.266 71.01% Castro::do_advance_ctu() 11 7.239 7.241 7.243 70.78% Castro::construct_ctu_hydro_source() 11 3.413 3.464 3.574 34.92% Castro::clean_state() 62 2.933 3.04 3.266 31.91% Castro::computeTemp() 63 1.906 2.008 2.215 21.64% Amr::InitAmr() 1 1.857 1.857 1.857 18.15% Castro::reset_internal_energy(MultiFab) 63 1.187 1.253 1.417 13.85% Amr::init() 1 0.9474 0.9474 0.9474 9.26% Gravity::solve_phi_with_mlmg() 10 0.8002 0.8195 0.84 8.21% Amr::initialInit() 1 0.8397 0.8397 0.8398 8.21% Castro::do_old_sources() 11 0.68 0.7706 0.8376 8.19% Castro::initialize_do_advance() 11 0.757 0.7745 0.7959 7.78% Castro::evaluate_source_change() 144 0.7193 0.7303 0.7567 7.39% Castro::construct_new_gravity() 8 0.7066 0.7252 0.7452 7.28% Gravity::solve_for_phi() 8 0.6891 0.7075 0.7279 7.11% Castro::enforce_min_density() 62 0.6141 0.6455 0.7214 7.05% Gravity::actual_solve_with_mlmg() 10 0.6941 0.6961 0.6984 6.82% FabArray::ParallelCopy() 734 0.1363 0.4978 0.6867 6.71% MLMG::solve() 9 0.6678 0.6699 0.6707 6.55% Amr::InitializeInit() 1 0.4391 0.5037 0.6166 6.03% Amr::defBaseLevel() 1 0.4391 0.5037 0.6166 6.03% FabArray::ParallelCopy_finish() 734 0.04522 0.4079 0.5949 5.81% MLMG::oneIter() 65 0.5906 0.5916 0.5936 5.80% MLMG::mgVcycle() 65 0.5832 0.5841 0.5861 5.73% amrex::Copy() 523 0.5374 0.5476 0.5742 5.61% Castro::normalize_species() 62 0.44 0.4489 0.4876 4.76% Castro::do_new_sources() 8 0.4391 0.4593 0.4755 4.65% Castro::reset_internal_energy(Fab) 126 0.4292 0.4371 0.4605 4.50% AmrLevel::FillPatch() 82 0.2091 0.3274 0.4588 4.48% MultiFab::AddProduct() 2633 0.422 0.4288 0.4444 4.34% Amr::FinalizeInit() 1 0.2231 0.3361 0.4006 3.92% MLCellLinOp::smooth() 1560 0.3022 0.3291 0.3961 3.87% FillPatchIterator::Initialize 82 0.1379 0.2566 0.3886 3.80% FillPatchIterator::FillFromLevel0() 82 0.1271 0.2482 0.3811 3.72% FillPatchSingleLevel 82 0.1268 0.2478 0.3807 3.72% Castro::Castro() 1 0.3593 0.365 0.3712 3.63% Castro::initMFs() 1 0.3554 0.3611 0.3673 3.59% check_for_negative_density() 11 0.1128 0.2111 0.3524 3.44% MultiFab::Subtract() 125 0.2703 0.2746 0.2889 2.82% Castro::initData() 1 0.07733 0.1375 0.2558 2.50% FabArray::setVal() 3453 0.219 0.2237 0.2326 2.27% MLMG::mgVcycle_up::1 65 0.05223 0.1379 0.2234 2.18% Amr::bldFineLevels() 1 0.04657 0.1592 0.223 2.18% MLCellLinOp::applyBC() 3630 0.1644 0.1854 0.2223 2.17% MLPoisson::Fsmooth() 3120 0.1595 0.175 0.2113 2.06% MLMG::addInterpCorrection() 390 0.03223 0.1557 0.1937 1.89% Amr::grid_places() 2 0.002593 0.1155 0.1801 1.76% AmrMesh::MakeNewGrids() 2 0.002469 0.1154 0.18 1.76% TagBoxArray::mapPRD 2 0.0007739 0.1137 0.1785 1.74% FabArray::FillBoundary() 3240 0.1363 0.1535 0.1722 1.68% Amr::writePlotFile() 2 0.1498 0.1499 0.1499 1.47% FabArray::sum() 2633 0.1406 0.143 0.1481 1.45% MLMG::mgVcycle_down::0 65 0.1404 0.1438 0.1461 1.43% Gravity::fill_multipole_BCs() 10 0.1047 0.1227 0.1426 1.39% MLMG::mgVcycle_up::2 65 0.005424 0.05886 0.1409 1.38% Castro::expand_state() 11 0.1034 0.1166 0.1388 1.36% Castro::post_init() 1 0.1326 0.1334 0.1352 1.32% MLMG::mgVcycle_up::0 65 0.1021 0.1031 0.1055 1.03% FillBoundary_finish() 3240 0.07212 0.08802 0.1037 1.01% MLCellLinOp::apply() 501 0.08043 0.08823 0.1036 1.01% Gravity::multilevel_solve_for_new_phi() 1 0.1001 0.1001 0.1001 0.98% Gravity::actual_multilevel_solve() 1 0.1 0.1 0.1001 0.98% FabArray::ParallelCopy_nowait() 734 0.07893 0.08807 0.09897 0.97% MLMG::mgVcycle_up::3 65 0.005404 0.01181 0.08863 0.87% Castro::post_regrid() 2 0.08529 0.0871 0.08792 0.86% MLMG:computeResOfCorrection() 390 0.06085 0.06603 0.07801 0.76% MLCellLinOp::correctionResidual() 390 0.06014 0.06531 0.07725 0.75% FillBoundary_nowait() 3240 0.0438 0.05837 0.07395 0.72% Castro::derive() 18 0.06936 0.07016 0.07157 0.70% AmrLevel::derive() 18 0.06932 0.07012 0.07153 0.70% MLMG::mgVcycle_down::1 65 0.06596 0.06858 0.07092 0.69% Castro::post_timestep() 1 0.05999 0.06223 0.06542 0.64% Castro::construct_new_source() 48 0.05799 0.05892 0.06182 0.60% Castro::construct_old_source() 66 0.04297 0.04461 0.05923 0.58% MLPoisson::Fapply() 501 0.04796 0.05096 0.05833 0.57% MLCellLinOp::solutionResidual() 75 0.04538 0.05039 0.05376 0.53% Castro::initialize_advance() 1 0.04821 0.04918 0.04977 0.49% Castro::create_source_corrector() 8 0.0156 0.03042 0.04905 0.48% Amr::regrid() 1 0.0439 0.04468 0.04503 0.44% MLMG::mgVcycle_down::2 65 0.006297 0.02246 0.04108 0.40% MLMG::computeResidual() 65 0.03789 0.03939 0.04103 0.40% FabArray::Saxpy() 133 0.03467 0.0355 0.03901 0.38% Castro::construct_old_gravity_source() 11 0.0251 0.02612 0.03682 0.36% Castro::construct_new_gravity_source() 8 0.03267 0.03323 0.03582 0.35% MLMG::mgVcycle_bottom 65 0.0001804 0.002406 0.0357 0.35% MLMG::actualBottomSolve() 65 6.693e-05 0.002289 0.03558 0.35% MLCGSolver::bicgstab 4 0 0.002174 0.03478 0.34% Castro::construct_old_gravity() 11 0.02262 0.02667 0.03414 0.33% Gravity::get_old_grav_vector() 11 0.02241 0.02567 0.03387 0.33% FabArray::Xpay() 480 0.02568 0.02711 0.03056 0.30% MLMG::mgVcycle_down::3 65 0.006024 0.009205 0.03028 0.30% amrex::average_down 390 0.0164 0.02026 0.02884 0.28% Castro::check_for_nan() 19 0.02303 0.02344 0.02399 0.23% MultiFab::contains_nan() 19 0.02297 0.02338 0.02391 0.23% Castro::retry_advance_ctu() 11 0.02188 0.02235 0.02299 0.22% VisMF::Write(FabArray) 4 0.01162 0.01261 0.02195 0.21% Castro::sum_integrated_quantities() 2 0.02108 0.02137 0.02139 0.21% MLPoisson::define() 10 0.01476 0.01715 0.0204 0.20% Gravity::get_new_grav_vector() 9 0.0173 0.01812 0.01878 0.18% MLMG::mgVcycle_down::4 65 0.006006 0.006836 0.01761 0.17% Castro::problem_post_timestep() 2 0.009886 0.01247 0.01713 0.17% Castro::wd_update() 2 0.009873 0.01243 0.01711 0.17% Castro::react_state() 19 0.01486 0.01519 0.01614 0.16% Castro::apply_source_to_state() 19 0.01438 0.01462 0.01551 0.15% StateDataPhysBCFunct::() 82 0.0003391 0.007478 0.01411 0.14% Gravity::test_composite_phi() 1 0.01403 0.01405 0.01406 0.14% MLMG::mgVcycle_up::4 65 0.005434 0.006063 0.01381 0.13% StateData::FillBoundary(geom) 62 0 0.007068 0.01367 0.13% MLMG::prepareForSolve() 9 0.01001 0.01079 0.01349 0.13% MLMG::ResNormInf() 74 0.007674 0.01102 0.01334 0.13% MLMG::computeMLResidual() 9 0.006792 0.0104 0.01327 0.13% MLCellLinOp::defineAuxData() 10 0.009612 0.01039 0.01248 0.12% MLMG::mgVcycle_down::5 65 0.005999 0.006434 0.01163 0.11% Castro::save_data_for_retry() 3 0.0108 0.01105 0.0113 0.11% MLMG::mgVcycle_up::5 65 0.005484 0.005847 0.00988 0.10% Amr::writeSmallPlotFile() 2 0.008812 0.008819 0.008826 0.09% FabArray::mult() 40 0.007937 0.008137 0.008368 0.08% amrex::Add() 65 0.007192 0.007298 0.008023 0.08% FabArray::norminf() 143 0.005577 0.005891 0.007855 0.08% Other 5693 0.05517 0.06952 0.1014 0.99% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 16 B 20 B 26 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-32-g28b010126a1b) finalized