Initializing AMReX (24.07-37-g83ecf62e0b43)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.07-37-g83ecf62e0b43) initialized Starting run at 07:34:48 UTC on 2024-07-31. Successfully read inputs file ... Castro git describe: 24.07-35-gc7c6acad2 AMReX git describe: 24.07-37-g83ecf62e0 Microphysics git describe: 24.07-39-g727afa98 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 Now regridding at level lbase = 0 grid_places() time: 0.190207266 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001129034 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010772996 Gravity::multilevel_solve_for_new_phi() time = 0.098958783 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.01060499 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.049829522 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.005528717 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 ... 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 = -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 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 Castro::construct_old_gravity() time = 0.053811058 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.069493634 on level 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 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.007146713 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.072223449 on level 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 ... 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 Castro::construct_old_gravity() time = 0.004884846 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.05897364 on level 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 ... Entering construct_ctu_hydro_source() on level 0 Castro::construct_old_gravity() time = 0.003262791 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.058604637 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at level 0 Castro::construct_ctu_hydro_source() time = 0.300188225 on level 0  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014434855 Gravity::solve_for_phi() time = 0.089345889 on level 0 Castro::construct_new_gravity() time = 0.091428816 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.059381202 on level 0 Castro::construct_old_gravity() time = 0.026068357 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.059316208 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.343461518 on level 0  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.044359586 Gravity::solve_for_phi() time = 0.116825387 on level 0 Castro::construct_new_gravity() time = 0.119007727 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.059572229 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 Castro::construct_old_gravity() time = 0.007901961 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.076871175 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.363064317 on level 0  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.014652751 Gravity::solve_for_phi() time = 0.085411496 on level 0 Castro::construct_new_gravity() time = 0.087476667 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.058302506 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 Castro::construct_old_gravity() time = 0.003325384 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.058804238 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.30906295 on level 0  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.01506762 Gravity::solve_for_phi() time = 0.085113199 on level 0 Castro::construct_new_gravity() time = 0.087130961 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.059201914 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 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.003500192 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.058267891 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.305448692 on level 0  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.01463293 Gravity::solve_for_phi() time = 0.085252765 on level 0 Castro::construct_new_gravity() time = 0.087282581 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.058701443 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 Castro::construct_old_gravity() time = 0.007388466 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.062848001 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.315405783 on level 0  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016843326 Gravity::solve_for_phi() time = 0.087472121 on level 0 Castro::construct_new_gravity() time = 0.08951144 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.060546177 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 Castro::construct_old_gravity() time = 0.004169017 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.060907695 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.312095661 on level 0  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 ... Entering construct_ctu_hydro_source() on level 0 Gravity::fill_multipole_BCs() time = 0.016484484 Gravity::solve_for_phi() time = 0.087491886 on level 0 Castro::construct_new_gravity() time = 0.089524973 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.06042831 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 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.00311398 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.060612959 on level 0 ... Leaving construct_ctu_hydro_source() on level 0 ... new-time level Poisson gravity solve at 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 Castro::construct_ctu_hydro_source() time = 0.313766855 on level 0  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.004581460012 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.017098741 Gravity::solve_for_phi() time = 0.089198297 on level 0 Castro::construct_new_gravity() time = 0.091227206 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.060286982 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 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 TIME = 0.25 DT = 0.25 [STEP 1] Coarse TimeStep time: 7.21304939 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.048708879 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00376709 seconds Ending run at 07:34:58 UTC on 2024-07-31. Run time = 10.02487307 Run time without initialization = 7.266407259 Average number of zones advanced per microsecond: 0.005 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 10.02 ... 10.02 ... 10.02 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7230 2.7666 2.8200 28.13% Amr::InitAmr() 1 1.817 1.817 1.817 18.12% Castro::computeTemp() 63 0.6919 0.7583 0.8468 8.45% FabArray::ParallelCopy_finish() 734 0.3067 0.5058 0.6975 6.96% amrex::Copy() 523 0.5365 0.5467 0.5599 5.58% Castro::normalize_species() 62 0.4482 0.4602 0.5028 5.02% Castro::reset_internal_energy(Fab) 126 0.43 0.4381 0.4536 4.52% MultiFab::AddProduct() 2633 0.4141 0.4223 0.4499 4.49% Castro::initMFs() 1 0.3567 0.3619 0.3825 3.82% MultiFab::Subtract() 125 0.2657 0.2695 0.2758 2.75% FabArray::setVal() 3453 0.2182 0.2233 0.2342 2.34% MLPoisson::Fsmooth() 3120 0.1596 0.1732 0.2097 2.09% check_for_negative_density() 11 0.06873 0.1385 0.1869 1.86% TagBoxArray::collate() 2 0.0002464 0.03755 0.1737 1.73% Gravity::fill_multipole_BCs() 10 0.1311 0.158 0.1729 1.72% Castro::initData() 1 0.003283 0.02811 0.165 1.65% FabArray::sum() 2633 0.1409 0.1438 0.1524 1.52% FabArray::ParallelCopy_nowait() 734 0.0632 0.07224 0.09849 0.98% FillBoundary_finish() 3240 0.04902 0.06402 0.078 0.78% FillBoundary_nowait() 3240 0.04063 0.05291 0.06817 0.68% MLPoisson::Fapply() 501 0.04808 0.05028 0.05798 0.58% MLCellLinOp::applyBC() 3630 0.02397 0.03197 0.05354 0.53% AmrLevel::derive() 18 0.04113 0.04175 0.04317 0.43% Castro::sum_integrated_quantities() 2 0.002607 0.03866 0.04141 0.41% FabArray::Saxpy() 133 0.03547 0.03605 0.03967 0.40% Castro::problem_post_timestep() 2 1.115e-05 0.002413 0.03841 0.38% Castro::construct_new_gravity_source() 8 0.03249 0.03303 0.03434 0.34% FabArray::Xpay() 480 0.02543 0.02678 0.03038 0.30% VisMF::Write(FabArray) 4 0.01163 0.01381 0.02683 0.27% MultiFab::contains_nan() 19 0.02275 0.02342 0.02598 0.26% Castro::construct_old_gravity_source() 11 0.02456 0.02488 0.02535 0.25% MLMG::addInterpCorrection() 390 0.017 0.01825 0.02189 0.22% Amr::writePlotFile() 2 0.007175 0.01777 0.01999 0.20% Castro::reset_internal_energy(MultiFab) 63 0.01768 0.01821 0.01909 0.19% amrex::average_down 390 0.01238 0.01348 0.01646 0.16% Castro::wd_update() 2 0.004939 0.0106 0.01615 0.16% Castro::enforce_min_density() 62 0.01355 0.01387 0.01465 0.15% StateData::FillBoundary(geom) 62 0 0.007286 0.01395 0.14% Castro::evaluate_source_change() 144 0.009329 0.00969 0.01019 0.10% MLCellLinOp::defineAuxData() 10 0.007646 0.007968 0.00913 0.09% FabArray::mult() 40 0.008244 0.00848 0.009026 0.09% Castro::volWgtSum() 46 0.007364 0.007605 0.008831 0.09% FabArray::FillBoundary() 3240 0.006769 0.007032 0.00843 0.08% MLCGSolver::bicgstab 4 0 0.0005094 0.008151 0.08% Castro::construct_new_source() 48 0.007352 0.007491 0.007894 0.08% FabArray::norminf() 143 0.005482 0.005763 0.007776 0.08% Amr::writeSmallPlotFile() 2 0.003609 0.006173 0.006943 0.07% MLCellLinOp::smooth() 1560 0.005493 0.005723 0.006148 0.06% Castro::construct_old_source() 66 0.005476 0.005567 0.005882 0.06% FillPatchIterator::Initialize 82 0.005426 0.005573 0.005772 0.06% MLMG::prepareForSolve() 9 0.00484 0.004934 0.005234 0.05% MLMG::ResNormInf() 74 0.002143 0.003375 0.004425 0.04% Gravity::get_old_grav_vector() 11 0.003506 0.003681 0.004224 0.04% Gravity::actual_solve_with_mlmg() 10 0.003447 0.003596 0.004049 0.04% FabArrayBase::getCPC() 1022 0.003721 0.003845 0.004019 0.04% Castro::initialize_do_advance() 11 0.002413 0.002627 0.00332 0.03% Gravity::get_new_grav_vector() 9 0.002928 0.002983 0.003089 0.03% MLMG::solve() 9 0.0003636 0.001369 0.002995 0.03% MLCellLinOp::apply() 501 0.0009825 0.001099 0.002332 0.02% FabArray::ParallelCopy() 734 0.001836 0.001901 0.002047 0.02% AmrLevel::FillPatch() 82 0.001098 0.001148 0.00121 0.01% MLCellLinOp::correctionResidual() 390 0.0009284 0.0009568 0.001027 0.01% MLMG::mgVcycle() 65 0.000763 0.0007999 0.0008926 0.01% Castro::save_data_for_retry() 3 0.0007433 0.0007766 0.0008208 0.01% main() 1 0.0005948 0.0006685 0.0007861 0.01% MLMG:computeResOfCorrection() 390 0.0006979 0.0007242 0.000776 0.01% Castro::subcycle_advance_ctu() 1 8.074e-05 0.0001874 0.0006118 0.01% AmrMesh::MakeNewGrids() 2 0.0005279 0.0005431 0.0005674 0.01% MLMG::actualBottomSolve() 65 6.696e-05 9.868e-05 0.000548 0.01% StateDataPhysBCFunct::() 82 0.0003427 0.0004287 0.0004756 0.00% Castro::Castro() 1 0.0003529 0.0003797 0.0004591 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003644 0.0003844 0.0004069 0.00% Castro::clean_state() 62 0.0003115 0.000341 0.0003805 0.00% MLMG::mgVcycle_down::3 65 0.0002889 0.0003099 0.0003727 0.00% MLMG::mgVcycle_down::2 65 0.0003083 0.000333 0.0003706 0.00% MLMG::mgVcycle_down::4 65 0.0002826 0.0002944 0.0003589 0.00% MLMG::mgVcycle_down::1 65 0.0003321 0.0003388 0.0003492 0.00% MLMG::mgVcycle_down::0 65 0.0003278 0.0003339 0.0003451 0.00% MLMG::mgVcycle_down::5 65 0.0002911 0.0002998 0.0003348 0.00% Gravity::solve_for_phi() 8 0.0002716 0.0002894 0.000312 0.00% Castro::do_advance_ctu() 11 0.0002216 0.0002449 0.000299 0.00% FillPatchSingleLevel 82 0.0002492 0.0002625 0.0002845 0.00% MLMG::mgVcycle_up::0 65 0.0002345 0.0002417 0.0002572 0.00% MLMG::mgVcycle_up::1 65 0.0002277 0.0002417 0.0002544 0.00% MLMG::mgVcycle_up::2 65 0.0001978 0.0002169 0.0002543 0.00% Castro::retry_advance_ctu() 11 5.005e-05 7.58e-05 0.0002375 0.00% MLCellLinOp::solutionResidual() 75 0.0002203 0.0002272 0.0002367 0.00% MLMG::mgVcycle_up::3 65 0.0001948 0.000203 0.000231 0.00% MLMG::mgVcycle_up::5 65 0.0002026 0.000208 0.0002245 0.00% MLMG::mgVcycle_up::4 65 0.0001976 0.000204 0.0002181 0.00% Castro::do_old_sources() 11 0.0001597 0.0001767 0.0002119 0.00% MLMG::oneIter() 65 0.0001867 0.0001915 0.0001976 0.00% Amr::coarseTimeStep() 1 0.0001173 0.0001502 0.0001566 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001211 0.0001303 0.0001456 0.00% Castro::construct_new_gravity() 8 7.71e-05 8.661e-05 0.0001425 0.00% Castro::initialize_advance() 1 0.0001001 0.0001116 0.0001327 0.00% Castro::apply_source_to_state() 19 5.152e-05 5.988e-05 0.0001316 0.00% MLMG::mgVcycle_bottom 65 0.0001139 0.0001171 0.0001266 0.00% MLMG::computeResidual() 65 0.0001183 0.0001203 0.0001253 0.00% Castro::do_new_sources() 8 9.951e-05 0.0001096 0.0001214 0.00% TagBoxArray::mapPRD 2 0.0001102 0.0001132 0.0001183 0.00% Amr::defBaseLevel() 1 9.424e-05 0.0001056 0.0001179 0.00% Gravity::test_composite_phi() 1 9.017e-05 9.336e-05 0.0001066 0.00% MLPoisson::define() 10 6.724e-05 7.289e-05 7.951e-05 0.00% Castro::construct_old_gravity() 11 6.027e-05 6.465e-05 7.521e-05 0.00% Amr::grid_places() 2 5.358e-05 5.659e-05 6.476e-05 0.00% Amr::FinalizeInit() 1 8.058e-06 1.209e-05 6.45e-05 0.00% Castro::check_for_nan() 19 4.923e-05 5.42e-05 5.795e-05 0.00% Castro::create_source_corrector() 8 3.271e-05 3.834e-05 5.37e-05 0.00% Castro::react_state() 19 4.635e-05 4.863e-05 5.353e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.017e-05 5.151e-05 5.255e-05 0.00% Castro::derive() 18 4.016e-05 4.175e-05 4.306e-05 0.00% Amr::regrid() 1 1.512e-05 1.884e-05 4.04e-05 0.00% Castro::expand_state() 11 3.041e-05 3.427e-05 3.727e-05 0.00% Castro::post_regrid() 2 1.061e-05 1.218e-05 2.845e-05 0.00% Gravity::actual_multilevel_solve() 1 2.711e-05 2.77e-05 2.844e-05 0.00% Amr::timeStep() 1 1.244e-05 1.678e-05 2.569e-05 0.00% MLMG::computeMLResidual() 9 1.886e-05 2.082e-05 2.253e-05 0.00% Castro::advance() 1 6.928e-06 8.992e-06 1.953e-05 0.00% Castro::post_init() 1 1.176e-05 1.279e-05 1.792e-05 0.00% Amr::InitializeInit() 1 6.814e-06 7.948e-06 1.481e-05 0.00% Amr::init() 1 7.727e-06 1.036e-05 1.287e-05 0.00% Castro::post_timestep() 1 5.698e-06 6.193e-06 9.8e-06 0.00% Amr::bldFineLevels() 1 6.688e-06 7.577e-06 8.79e-06 0.00% Amr::initialInit() 1 2.866e-06 3.612e-06 6.723e-06 0.00% Other 4688 0.0339 0.04344 0.06228 0.62% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.02 10.02 10.02 100.00% Amr::coarseTimeStep() 1 7.266 7.266 7.266 72.48% Amr::timeStep() 1 7.213 7.213 7.213 71.95% Castro::advance() 1 7.144 7.149 7.153 71.35% Castro::subcycle_advance_ctu() 1 7.055 7.091 7.103 70.85% Castro::do_advance_ctu() 11 7.031 7.067 7.079 70.62% Castro::construct_ctu_hydro_source() 11 3.244 3.271 3.313 33.05% Castro::clean_state() 62 2.912 3.012 3.098 30.90% Castro::computeTemp() 63 1.872 1.958 2.049 20.44% Amr::InitAmr() 1 1.817 1.817 1.817 18.12% Castro::reset_internal_energy(MultiFab) 63 1.18 1.2 1.231 12.28% Castro::initialize_do_advance() 11 0.8568 0.903 0.9478 9.45% Amr::init() 1 0.9415 0.9415 0.9415 9.39% Amr::initialInit() 1 0.8857 0.8857 0.8858 8.84% Gravity::solve_phi_with_mlmg() 10 0.7922 0.8211 0.8363 8.34% FabArray::ParallelCopy() 734 0.4076 0.5963 0.7836 7.82% Castro::evaluate_source_change() 144 0.7118 0.726 0.7729 7.71% Castro::construct_new_gravity() 8 0.698 0.7266 0.7418 7.40% Gravity::solve_for_phi() 8 0.6812 0.71 0.7253 7.23% FabArray::ParallelCopy_finish() 734 0.3067 0.5058 0.6975 6.96% Castro::do_old_sources() 11 0.6442 0.6694 0.6907 6.89% Castro::enforce_min_density() 62 0.6077 0.6205 0.6685 6.67% Gravity::actual_solve_with_mlmg() 10 0.6592 0.6613 0.6632 6.62% MLMG::solve() 9 0.6313 0.6336 0.6343 6.33% Amr::InitializeInit() 1 0.4375 0.4681 0.6261 6.25% Amr::defBaseLevel() 1 0.4375 0.4681 0.6261 6.25% MLMG::oneIter() 65 0.5659 0.5666 0.5674 5.66% MLMG::mgVcycle() 65 0.5584 0.5591 0.5599 5.58% amrex::Copy() 523 0.5365 0.5467 0.5599 5.58% Castro::normalize_species() 62 0.4482 0.4602 0.5028 5.02% AmrLevel::FillPatch() 82 0.3574 0.4204 0.5016 5.00% Castro::do_new_sources() 8 0.4265 0.4518 0.4763 4.75% Castro::reset_internal_energy(Fab) 126 0.43 0.4381 0.4536 4.52% MultiFab::AddProduct() 2633 0.4141 0.4223 0.4499 4.49% Amr::FinalizeInit() 1 0.2595 0.4176 0.4482 4.47% FillPatchIterator::Initialize 82 0.2858 0.349 0.4318 4.31% FillPatchIterator::FillFromLevel0() 82 0.28 0.3405 0.4243 4.23% FillPatchSingleLevel 82 0.2796 0.3401 0.424 4.23% Castro::Castro() 1 0.3605 0.3657 0.3864 3.85% Castro::initMFs() 1 0.3569 0.362 0.3826 3.82% MLCellLinOp::smooth() 1560 0.2832 0.3081 0.373 3.72% MultiFab::Subtract() 125 0.2657 0.2695 0.2758 2.75% Castro::initData() 1 0.07572 0.1011 0.2383 2.38% FabArray::setVal() 3453 0.2182 0.2233 0.2342 2.34% Amr::bldFineLevels() 1 0.04571 0.2035 0.2339 2.33% Castro::expand_state() 11 0.1308 0.1806 0.2295 2.29% MLMG::mgVcycle_up::1 65 0.04887 0.1322 0.2155 2.15% MLPoisson::Fsmooth() 3120 0.1596 0.1732 0.2097 2.09% MLCellLinOp::applyBC() 3630 0.14 0.1618 0.1986 1.98% Amr::grid_places() 2 0.002701 0.1607 0.1913 1.91% AmrMesh::MakeNewGrids() 2 0.002648 0.1607 0.1913 1.91% MLMG::addInterpCorrection() 390 0.0322 0.1542 0.1892 1.89% TagBoxArray::mapPRD 2 0.001122 0.1218 0.1892 1.89% check_for_negative_density() 11 0.06873 0.1385 0.1869 1.86% Gravity::fill_multipole_BCs() 10 0.1322 0.1591 0.174 1.74% TagBoxArray::collate() 2 0.0002464 0.03755 0.1737 1.73% Castro::post_init() 1 0.1696 0.1714 0.172 1.72% FabArray::sum() 2633 0.1409 0.1438 0.1524 1.52% FabArray::FillBoundary() 3240 0.109 0.1299 0.1467 1.46% MLMG::mgVcycle_up::2 65 0.005462 0.05891 0.1408 1.40% MLMG::mgVcycle_down::0 65 0.1331 0.1349 0.1374 1.37% Castro::create_source_corrector() 8 0.04513 0.09789 0.1282 1.28% FabArray::ParallelCopy_nowait() 734 0.07945 0.08858 0.1151 1.15% MLMG::mgVcycle_up::0 65 0.1001 0.1004 0.1008 1.01% MLCellLinOp::apply() 501 0.07705 0.08303 0.1008 1.01% Gravity::multilevel_solve_for_new_phi() 1 0.09891 0.09893 0.09896 0.99% Amr::writePlotFile() 2 0.09882 0.09893 0.09894 0.99% Gravity::actual_multilevel_solve() 1 0.09886 0.09888 0.09891 0.99% MLMG::mgVcycle_up::3 65 0.005444 0.01194 0.0903 0.90% Castro::initialize_advance() 1 0.04826 0.05813 0.08983 0.90% Castro::post_regrid() 2 0.08482 0.08537 0.08721 0.87% Castro::construct_old_gravity() 11 0.02753 0.04083 0.08651 0.86% Gravity::get_old_grav_vector() 11 0.02746 0.04077 0.08645 0.86% FillBoundary_finish() 3240 0.04902 0.06402 0.078 0.78% MLMG:computeResOfCorrection() 390 0.05898 0.06376 0.07581 0.76% MLCellLinOp::correctionResidual() 390 0.05826 0.06304 0.07504 0.75% FillBoundary_nowait() 3240 0.04613 0.05881 0.075 0.75% Castro::derive() 18 0.06935 0.07003 0.07147 0.71% AmrLevel::derive() 18 0.06931 0.06999 0.07143 0.71% Castro::post_timestep() 1 0.05986 0.06403 0.06876 0.69% MLMG::mgVcycle_down::1 65 0.05885 0.06247 0.06473 0.65% Castro::construct_new_source() 48 0.05822 0.05902 0.06043 0.60% Castro::sum_integrated_quantities() 2 0.022 0.05791 0.06031 0.60% MLPoisson::Fapply() 501 0.04808 0.05028 0.05798 0.58% MLCellLinOp::solutionResidual() 75 0.04426 0.04714 0.04944 0.49% Castro::problem_post_timestep() 2 0.009573 0.01672 0.04882 0.49% Amr::regrid() 1 0.04352 0.04392 0.04504 0.45% Castro::construct_old_source() 66 0.04265 0.04315 0.04366 0.44% FabArray::Saxpy() 133 0.03547 0.03605 0.03967 0.40% MLMG::computeResidual() 65 0.03789 0.03874 0.03917 0.39% MLMG::mgVcycle_down::2 65 0.006358 0.02068 0.03674 0.37% MLMG::mgVcycle_bottom 65 0.0001815 0.002451 0.03644 0.36% MLMG::actualBottomSolve() 65 6.696e-05 0.002334 0.03631 0.36% MLCGSolver::bicgstab 4 0 0.00222 0.03552 0.35% Castro::construct_new_gravity_source() 8 0.03249 0.03303 0.03434 0.34% FabArray::Xpay() 480 0.02543 0.02678 0.03038 0.30% MLMG::mgVcycle_down::3 65 0.006038 0.009208 0.03032 0.30% amrex::average_down 390 0.01672 0.02009 0.0287 0.29% VisMF::Write(FabArray) 4 0.01163 0.01381 0.02683 0.27% Castro::check_for_nan() 19 0.0228 0.02348 0.02604 0.26% MultiFab::contains_nan() 19 0.02275 0.02342 0.02598 0.26% Castro::construct_old_gravity_source() 11 0.02456 0.02488 0.02535 0.25% Castro::retry_advance_ctu() 11 0.0217 0.02238 0.02408 0.24% MLPoisson::define() 10 0.01624 0.0186 0.02142 0.21% Castro::wd_update() 2 0.009561 0.01431 0.01953 0.19% MLMG::mgVcycle_down::4 65 0.00601 0.006852 0.01784 0.18% Gravity::get_new_grav_vector() 9 0.01716 0.01745 0.01768 0.18% Castro::apply_source_to_state() 19 0.01468 0.01489 0.01608 0.16% Castro::react_state() 19 0.01444 0.01481 0.01544 0.15% StateDataPhysBCFunct::() 82 0.0003427 0.007714 0.01442 0.14% Gravity::test_composite_phi() 1 0.0141 0.01413 0.01414 0.14% StateData::FillBoundary(geom) 62 0 0.007286 0.01395 0.14% MLMG::mgVcycle_up::4 65 0.005493 0.006064 0.01394 0.14% MLCellLinOp::defineAuxData() 10 0.009865 0.01059 0.01284 0.13% Castro::save_data_for_retry() 3 0.01079 0.01122 0.0127 0.13% MLMG::mgVcycle_down::5 65 0.006004 0.006438 0.01185 0.12% MLMG::prepareForSolve() 9 0.009906 0.01027 0.01141 0.11% MLMG::mgVcycle_up::5 65 0.005512 0.005871 0.0101 0.10% MLMG::computeMLResidual() 9 0.005357 0.007752 0.009569 0.10% Amr::writeSmallPlotFile() 2 0.009534 0.009542 0.009548 0.10% MLMG::ResNormInf() 74 0.006989 0.008205 0.009211 0.09% FabArray::mult() 40 0.008244 0.00848 0.009026 0.09% Castro::volWgtSum() 46 0.007364 0.007605 0.008831 0.09% FabArray::norminf() 143 0.005482 0.005763 0.007776 0.08% FabArrayBase::getCPC() 1022 0.005165 0.005876 0.007722 0.08% Other 4688 0.05454 0.07086 0.0976 0.97% -------------------------------------------------------------------------------------------- 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 15 B 18 B 28 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.07-37-g83ecf62e0b43) finalized