Initializing AMReX (24.08-9-g6f039ae18e89)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.08-9-g6f039ae18e89) initialized Starting run at 07:58:11 UTC on 2024-08-07. Successfully read inputs file ... Castro git describe: 24.08-1-g36c0aaf32 AMReX git describe: 24.08-9-g6f039ae18 Microphysics git describe: 24.08-2-g8ce3375a 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.177534445 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.00125717 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010694266 Gravity::multilevel_solve_for_new_phi() time = 0.102857548 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.010585757 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.048812046 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003734918 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 = -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 Castro::construct_old_gravity() time = 0.002722891 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.068547771 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 = -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 Castro::construct_old_gravity() time = 0.00555469 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.067264966 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.006098512 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.063042603 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.015862991 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.068490907 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.305216119 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.015963406 Gravity::solve_for_phi() time = 0.08874518 on level 0 Castro::construct_new_gravity() time = 0.09087544 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.060437002 on level 0 Castro::construct_old_gravity() time = 0.003213776 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.0613639 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.308452031 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.016312532 Gravity::solve_for_phi() time = 0.090809164 on level 0 Castro::construct_new_gravity() time = 0.092886916 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.059555273 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.003974468 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.062531181 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.311976951 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.014987188 Gravity::solve_for_phi() time = 0.093259239 on level 0 Castro::construct_new_gravity() time = 0.09538389 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.060421695 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.00239319 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.08987958 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.347414605 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.015409137 Gravity::solve_for_phi() time = 0.091197438 on level 0 Castro::construct_new_gravity() time = 0.093401005 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.05976576 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.012393105 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.086769847 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.352947145 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.017346479 Gravity::solve_for_phi() time = 0.120361766 on level 0 Castro::construct_new_gravity() time = 0.122441104 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.067506497 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.013053741 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.102783834 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.413518185 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.016080496 Gravity::solve_for_phi() time = 0.090048355 on level 0 Castro::construct_new_gravity() time = 0.092132049 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.059676239 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.047287563 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.104949461 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.328518183 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.016243007 Gravity::solve_for_phi() time = 0.089422361 on level 0 Castro::construct_new_gravity() time = 0.095149769 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.059675001 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.00260136 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.060074391 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.310717631 on level 0  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.004494326612 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.015931144 Gravity::solve_for_phi() time = 0.087504718 on level 0 Castro::construct_new_gravity() time = 0.089592249 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.05941432 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.353716208 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.049396627 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003753004 seconds Ending run at 07:58:21 UTC on 2024-08-07. Run time = 9.70876276 Run time without initialization = 7.407811657 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]: 9.709 ... 9.709 ... 9.709 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7212 2.7664 2.8802 29.67% Amr::InitAmr() 1 1.397 1.397 1.397 14.39% Castro::computeTemp() 63 0.7115 0.7565 0.8323 8.57% FabArray::ParallelCopy_finish() 734 0.2183 0.4697 0.6862 7.07% amrex::Copy() 523 0.5818 0.594 0.6547 6.74% Castro::normalize_species() 62 0.4487 0.459 0.5173 5.33% MultiFab::AddProduct() 2633 0.4226 0.4298 0.4735 4.88% Castro::reset_internal_energy(Fab) 126 0.4303 0.4368 0.45 4.63% Castro::initMFs() 1 0.3611 0.369 0.3819 3.93% check_for_negative_density() 11 0.08677 0.2415 0.3331 3.43% MultiFab::Subtract() 125 0.2742 0.2795 0.2957 3.05% FabArray::setVal() 3453 0.2164 0.2202 0.2376 2.45% MLPoisson::Fsmooth() 3120 0.1632 0.1793 0.2145 2.21% Castro::initData() 1 0.003312 0.02842 0.1684 1.73% FabArray::sum() 2633 0.14 0.1424 0.1582 1.63% Gravity::fill_multipole_BCs() 10 0.1023 0.1343 0.147 1.51% TagBoxArray::collate() 2 0.0003637 0.03054 0.1467 1.51% FillBoundary_finish() 3240 0.07225 0.09241 0.1116 1.15% FabArray::ParallelCopy_nowait() 734 0.06375 0.07153 0.07932 0.82% FillBoundary_nowait() 3240 0.04349 0.05682 0.07184 0.74% MLPoisson::Fapply() 501 0.04832 0.05102 0.05862 0.60% MLCellLinOp::applyBC() 3630 0.02534 0.03369 0.05518 0.57% AmrLevel::derive() 18 0.04094 0.04153 0.04266 0.44% FabArray::Saxpy() 133 0.03537 0.03617 0.03972 0.41% Castro::construct_new_gravity_source() 8 0.03203 0.03266 0.03487 0.36% FabArray::Xpay() 480 0.02525 0.02671 0.03092 0.32% Castro::construct_old_gravity_source() 11 0.02438 0.02488 0.02696 0.28% VisMF::Write(FabArray) 4 0.01143 0.0125 0.02406 0.25% MultiFab::contains_nan() 19 0.0229 0.02342 0.02396 0.25% MLMG::addInterpCorrection() 390 0.01706 0.01849 0.02207 0.23% Castro::reset_internal_energy(MultiFab) 63 0.01776 0.01813 0.01887 0.19% Amr::writePlotFile() 2 0.006075 0.01633 0.01793 0.18% amrex::average_down 390 0.01245 0.0136 0.01666 0.17% StateData::FillBoundary(geom) 62 0 0.007261 0.01458 0.15% Castro::enforce_min_density() 62 0.01362 0.01385 0.01412 0.15% Castro::wd_update() 2 0.004571 0.01064 0.01341 0.14% Castro::evaluate_source_change() 144 0.009584 0.01 0.0107 0.11% MLCellLinOp::defineAuxData() 10 0.00774 0.00824 0.009657 0.10% FabArray::FillBoundary() 3240 0.006971 0.007335 0.008972 0.09% MLMG::ResNormInf() 74 0.003447 0.00737 0.008853 0.09% FabArray::mult() 40 0.007975 0.008114 0.008458 0.09% MLCGSolver::bicgstab 4 0 0.0005249 0.008398 0.09% FabArray::norminf() 143 0.005614 0.006017 0.007908 0.08% amrex::Add() 65 0.007257 0.007385 0.007494 0.08% Castro::volWgtSum() 46 0.007272 0.007378 0.007494 0.08% Castro::construct_new_source() 48 0.007256 0.007333 0.007452 0.08% Castro::construct_old_source() 66 0.0054 0.005599 0.006783 0.07% MLCellLinOp::smooth() 1560 0.005745 0.006014 0.006456 0.07% FillPatchIterator::Initialize 82 0.00548 0.005605 0.005765 0.06% MLMG::prepareForSolve() 9 0.004946 0.005074 0.005454 0.06% Amr::writeSmallPlotFile() 2 0.003018 0.004603 0.005032 0.05% MLMG::solve() 9 0.0004623 0.002057 0.004208 0.04% FabArrayBase::getCPC() 1022 0.003818 0.003937 0.004112 0.04% Gravity::actual_solve_with_mlmg() 10 0.003459 0.003606 0.004049 0.04% Gravity::get_old_grav_vector() 11 0.003474 0.003538 0.003603 0.04% Castro::sum_integrated_quantities() 2 0.001486 0.003276 0.003586 0.04% Gravity::get_new_grav_vector() 9 0.00293 0.002982 0.003086 0.03% Castro::initialize_do_advance() 11 0.002374 0.002533 0.00273 0.03% MLCellLinOp::apply() 501 0.001013 0.001137 0.002419 0.02% FabArray::ParallelCopy() 734 0.001785 0.001869 0.002052 0.02% AmrLevel::FillPatch() 82 0.001105 0.001154 0.001262 0.01% MLCellLinOp::correctionResidual() 390 0.0009465 0.0009783 0.001054 0.01% MLMG::mgVcycle() 65 0.0007818 0.0008205 0.0009161 0.01% main() 1 0.0006363 0.0006905 0.0008701 0.01% Castro::save_data_for_retry() 3 0.0007432 0.0007691 0.0008009 0.01% MLMG:computeResOfCorrection() 390 0.0007029 0.0007276 0.0007614 0.01% Castro::subcycle_advance_ctu() 1 7.801e-05 0.0001324 0.0007359 0.01% AmrMesh::MakeNewGrids() 2 0.0005329 0.0005507 0.0006253 0.01% MLMG::actualBottomSolve() 65 6.78e-05 0.0001009 0.0005531 0.01% Castro::Castro() 1 0.0004081 0.0004685 0.000514 0.01% StateDataPhysBCFunct::() 82 0.0003717 0.0004341 0.0005046 0.01% Castro::clean_state() 62 0.0003121 0.000345 0.0003996 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003619 0.0003817 0.0003975 0.00% MLMG::mgVcycle_down::1 65 0.0003337 0.0003465 0.0003718 0.00% MLMG::mgVcycle_down::3 65 0.000291 0.0003169 0.0003702 0.00% MLMG::mgVcycle_down::2 65 0.0003128 0.0003412 0.0003671 0.00% MLMG::mgVcycle_down::0 65 0.0003316 0.0003431 0.0003657 0.00% MLMG::mgVcycle_down::4 65 0.0002814 0.0002962 0.0003557 0.00% MLMG::mgVcycle_down::5 65 0.0002906 0.000304 0.0003444 0.00% Gravity::solve_for_phi() 8 0.0002721 0.000285 0.0003098 0.00% Castro::do_advance_ctu() 11 0.0002043 0.0002248 0.0002856 0.00% FillPatchSingleLevel 82 0.0002506 0.0002622 0.0002738 0.00% MLMG::mgVcycle_up::1 65 0.0002225 0.0002382 0.0002664 0.00% MLMG::mgVcycle_up::0 65 0.0002286 0.0002365 0.0002621 0.00% MLMG::mgVcycle_up::2 65 0.0001906 0.0002138 0.0002605 0.00% Castro::problem_post_timestep() 2 1.057e-05 2.703e-05 0.0002563 0.00% MLCellLinOp::solutionResidual() 75 0.0002243 0.0002336 0.0002458 0.00% MLMG::mgVcycle_up::5 65 0.0001992 0.0002071 0.0002266 0.00% MLMG::mgVcycle_up::3 65 0.0001898 0.0001991 0.0002203 0.00% MLMG::mgVcycle_up::4 65 0.0001919 0.0002003 0.0002195 0.00% MLMG::oneIter() 65 0.0001911 0.0001987 0.0002166 0.00% Castro::do_old_sources() 11 0.000158 0.00017 0.0001976 0.00% Amr::FinalizeInit() 1 1.305e-05 3.47e-05 0.0001782 0.00% Amr::coarseTimeStep() 1 0.0001452 0.0001687 0.0001746 0.00% Castro::construct_new_gravity() 8 7.956e-05 9.062e-05 0.0001596 0.00% Amr::defBaseLevel() 1 0.0001181 0.0001405 0.0001567 0.00% Castro::initialize_advance() 1 0.000112 0.0001213 0.0001566 0.00% TagBoxArray::mapPRD 2 0.0001145 0.000122 0.0001409 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001154 0.0001241 0.0001325 0.00% MLMG::computeResidual() 65 0.0001211 0.0001244 0.0001321 0.00% MLMG::mgVcycle_bottom 65 0.0001105 0.0001154 0.0001222 0.00% Castro::retry_advance_ctu() 11 5.06e-05 6.44e-05 0.000121 0.00% Castro::do_new_sources() 8 9.949e-05 0.0001041 0.0001106 0.00% Gravity::test_composite_phi() 1 9.246e-05 9.608e-05 0.0001102 0.00% Amr::grid_places() 2 5.825e-05 6.571e-05 7.992e-05 0.00% MLPoisson::define() 10 6.73e-05 7.159e-05 7.817e-05 0.00% Castro::check_for_nan() 19 5.221e-05 5.788e-05 7.476e-05 0.00% Castro::construct_old_gravity() 11 5.925e-05 6.368e-05 7.22e-05 0.00% Castro::apply_source_to_state() 19 5.23e-05 5.532e-05 5.987e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.04e-05 5.227e-05 5.46e-05 0.00% Castro::react_state() 19 4.801e-05 4.995e-05 5.338e-05 0.00% Castro::derive() 18 4.079e-05 4.396e-05 4.893e-05 0.00% Amr::regrid() 1 1.56e-05 1.99e-05 4.428e-05 0.00% Castro::create_source_corrector() 8 3.469e-05 3.772e-05 4.18e-05 0.00% Castro::expand_state() 11 3.065e-05 3.347e-05 3.639e-05 0.00% Gravity::actual_multilevel_solve() 1 2.739e-05 2.885e-05 3.136e-05 0.00% Amr::timeStep() 1 1.392e-05 1.623e-05 2.579e-05 0.00% Amr::InitializeInit() 1 7.881e-06 1.034e-05 2.45e-05 0.00% MLMG::computeMLResidual() 9 1.986e-05 2.158e-05 2.303e-05 0.00% Amr::bldFineLevels() 1 7.025e-06 9.196e-06 1.67e-05 0.00% Castro::post_regrid() 2 1.045e-05 1.163e-05 1.387e-05 0.00% Castro::post_init() 1 1.169e-05 1.24e-05 1.339e-05 0.00% Castro::advance() 1 6.644e-06 8.082e-06 1.248e-05 0.00% Amr::init() 1 6.213e-06 7.103e-06 8.096e-06 0.00% Castro::post_timestep() 1 5.389e-06 5.911e-06 6.634e-06 0.00% Amr::initialInit() 1 2.928e-06 3.566e-06 6.322e-06 0.00% Other 4623 0.02598 0.03676 0.05668 0.58% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.709 9.709 9.709 100.00% Amr::coarseTimeStep() 1 7.407 7.407 7.407 76.29% Amr::timeStep() 1 7.354 7.354 7.354 75.74% Castro::advance() 1 7.287 7.289 7.292 75.11% Castro::subcycle_advance_ctu() 1 7.236 7.239 7.243 74.60% Castro::do_advance_ctu() 11 7.211 7.215 7.219 74.35% Castro::construct_ctu_hydro_source() 11 3.299 3.378 3.458 35.61% Castro::clean_state() 62 2.974 3.054 3.217 33.14% Castro::computeTemp() 63 1.924 1.986 2.076 21.38% Amr::InitAmr() 1 1.397 1.397 1.397 14.39% Castro::reset_internal_energy(MultiFab) 63 1.211 1.23 1.301 13.40% Amr::init() 1 0.9036 0.9036 0.9036 9.31% Castro::initialize_do_advance() 11 0.7656 0.8246 0.8821 9.09% Gravity::solve_phi_with_mlmg() 10 0.8188 0.8511 0.8658 8.92% Amr::initialInit() 1 0.8507 0.8507 0.8508 8.76% Castro::evaluate_source_change() 144 0.717 0.7289 0.8047 8.29% Castro::do_old_sources() 11 0.6508 0.7117 0.7919 8.16% FabArray::ParallelCopy() 734 0.3073 0.5606 0.775 7.98% Castro::construct_new_gravity() 8 0.7218 0.7548 0.7698 7.93% Gravity::solve_for_phi() 8 0.7039 0.7361 0.7509 7.73% Gravity::actual_solve_with_mlmg() 10 0.7126 0.715 0.717 7.39% Castro::enforce_min_density() 62 0.6252 0.6365 0.6914 7.12% MLMG::solve() 9 0.6846 0.6873 0.6881 7.09% FabArray::ParallelCopy_finish() 734 0.2183 0.4697 0.6862 7.07% amrex::Copy() 523 0.5818 0.594 0.6547 6.74% Amr::InitializeInit() 1 0.4439 0.4785 0.6198 6.38% Amr::defBaseLevel() 1 0.4439 0.4785 0.6198 6.38% MLMG::oneIter() 65 0.6107 0.6134 0.6148 6.33% MLMG::mgVcycle() 65 0.6031 0.6058 0.6073 6.25% AmrLevel::FillPatch() 82 0.251 0.3933 0.5277 5.44% Castro::normalize_species() 62 0.4487 0.459 0.5173 5.33% Castro::do_new_sources() 8 0.4389 0.4582 0.476 4.90% MultiFab::AddProduct() 2633 0.4226 0.4298 0.4735 4.88% FillPatchIterator::Initialize 82 0.1737 0.3162 0.4522 4.66% Castro::reset_internal_energy(Fab) 126 0.4303 0.4368 0.45 4.63% FillPatchIterator::FillFromLevel0() 82 0.1628 0.3078 0.4446 4.58% FillPatchSingleLevel 82 0.1624 0.3074 0.4442 4.58% MLCellLinOp::smooth() 1560 0.3155 0.3451 0.4108 4.23% Amr::FinalizeInit() 1 0.2309 0.3722 0.4068 4.19% Castro::Castro() 1 0.3654 0.3734 0.3862 3.98% Castro::initMFs() 1 0.3612 0.3691 0.382 3.93% check_for_negative_density() 11 0.08677 0.2415 0.3331 3.43% MultiFab::Subtract() 125 0.2742 0.2795 0.2957 3.05% Castro::initData() 1 0.07698 0.1038 0.2424 2.50% MLCellLinOp::applyBC() 3630 0.1768 0.1961 0.2377 2.45% FabArray::setVal() 3453 0.2164 0.2202 0.2376 2.45% MLMG::mgVcycle_up::1 65 0.05555 0.1435 0.2325 2.39% Amr::bldFineLevels() 1 0.04621 0.1885 0.2236 2.30% MLPoisson::Fsmooth() 3120 0.1632 0.1793 0.2145 2.21% MLMG::addInterpCorrection() 390 0.03411 0.1607 0.2005 2.06% FabArray::FillBoundary() 3240 0.1377 0.1624 0.1873 1.93% Amr::grid_places() 2 0.002882 0.1441 0.1788 1.84% AmrMesh::MakeNewGrids() 2 0.002823 0.1441 0.1787 1.84% TagBoxArray::mapPRD 2 0.001141 0.1123 0.1724 1.78% Castro::expand_state() 11 0.1043 0.1324 0.1716 1.77% FabArray::sum() 2633 0.14 0.1424 0.1582 1.63% MLMG::mgVcycle_down::0 65 0.1417 0.1443 0.1489 1.53% Gravity::fill_multipole_BCs() 10 0.1036 0.1354 0.1482 1.53% MLMG::mgVcycle_up::2 65 0.005581 0.06141 0.1468 1.51% TagBoxArray::collate() 2 0.0003637 0.03054 0.1467 1.51% Castro::post_init() 1 0.1364 0.1398 0.142 1.46% MLMG::mgVcycle_up::0 65 0.1098 0.1108 0.1122 1.16% FillBoundary_finish() 3240 0.07225 0.09241 0.1116 1.15% MLCellLinOp::apply() 501 0.08084 0.08748 0.108 1.11% Gravity::multilevel_solve_for_new_phi() 1 0.1028 0.1028 0.1029 1.06% Gravity::actual_multilevel_solve() 1 0.1028 0.1028 0.1028 1.06% Amr::writePlotFile() 2 0.0985 0.09858 0.09859 1.02% FabArray::ParallelCopy_nowait() 734 0.0809 0.08904 0.09713 1.00% MLMG::mgVcycle_up::3 65 0.00556 0.01231 0.09271 0.95% Castro::post_regrid() 2 0.08591 0.08817 0.0915 0.94% Castro::construct_old_gravity() 11 0.02667 0.0609 0.09007 0.93% Gravity::get_old_grav_vector() 11 0.02661 0.06084 0.09001 0.93% Castro::create_source_corrector() 8 0.02634 0.06218 0.08582 0.88% MLMG:computeResOfCorrection() 390 0.06007 0.0664 0.07879 0.81% MLCellLinOp::correctionResidual() 390 0.05936 0.06567 0.07803 0.80% FillBoundary_nowait() 3240 0.04942 0.06266 0.0779 0.80% Castro::derive() 18 0.07119 0.07217 0.07611 0.78% AmrLevel::derive() 18 0.07115 0.07213 0.07607 0.78% MLMG::mgVcycle_down::1 65 0.06904 0.07184 0.07528 0.78% Castro::post_timestep() 1 0.06131 0.06428 0.06687 0.69% Castro::construct_new_source() 48 0.05753 0.0583 0.0607 0.63% MLPoisson::Fapply() 501 0.04832 0.05102 0.05862 0.60% MLCellLinOp::solutionResidual() 75 0.04555 0.04889 0.05259 0.54% Castro::initialize_advance() 1 0.04886 0.04997 0.05186 0.53% Amr::regrid() 1 0.04432 0.04563 0.04815 0.50% Castro::construct_old_source() 66 0.04238 0.04321 0.04539 0.47% MLMG::computeResidual() 65 0.03924 0.04052 0.04446 0.46% MLMG::mgVcycle_down::2 65 0.00659 0.02254 0.04112 0.42% FabArray::Saxpy() 133 0.03537 0.03617 0.03972 0.41% MLMG::mgVcycle_bottom 65 0.0001794 0.002507 0.03732 0.38% MLMG::actualBottomSolve() 65 6.78e-05 0.002391 0.0372 0.38% MLCGSolver::bicgstab 4 0 0.002275 0.0364 0.37% Castro::construct_new_gravity_source() 8 0.03203 0.03266 0.03487 0.36% MLMG::mgVcycle_down::3 65 0.006216 0.009499 0.03159 0.33% FabArray::Xpay() 480 0.02525 0.02671 0.03092 0.32% amrex::average_down 390 0.01675 0.02058 0.02985 0.31% Castro::construct_old_gravity_source() 11 0.02438 0.02488 0.02696 0.28% Castro::retry_advance_ctu() 11 0.02253 0.0232 0.0242 0.25% VisMF::Write(FabArray) 4 0.01143 0.0125 0.02406 0.25% Castro::check_for_nan() 19 0.02296 0.02348 0.02403 0.25% MultiFab::contains_nan() 19 0.0229 0.02342 0.02396 0.25% Castro::sum_integrated_quantities() 2 0.02289 0.02306 0.02308 0.24% Gravity::get_new_grav_vector() 9 0.01774 0.01949 0.02177 0.22% MLPoisson::define() 10 0.01694 0.01867 0.02094 0.22% MLMG::mgVcycle_down::4 65 0.006217 0.007104 0.01817 0.19% Castro::apply_source_to_state() 19 0.01466 0.01503 0.01789 0.18% Castro::problem_post_timestep() 2 0.01008 0.01433 0.01689 0.17% Castro::wd_update() 2 0.01007 0.0143 0.01688 0.17% Castro::react_state() 19 0.01463 0.01487 0.0155 0.16% StateDataPhysBCFunct::() 82 0.0003717 0.007695 0.01506 0.16% StateData::FillBoundary(geom) 62 0 0.007261 0.01458 0.15% MLMG::mgVcycle_up::4 65 0.005633 0.006275 0.0142 0.15% Gravity::test_composite_phi() 1 0.01416 0.01418 0.01419 0.15% MLMG::ResNormInf() 74 0.009135 0.0124 0.01375 0.14% MLCellLinOp::defineAuxData() 10 0.01005 0.01096 0.01351 0.14% Castro::save_data_for_retry() 3 0.01165 0.01185 0.01211 0.12% MLMG::mgVcycle_down::5 65 0.006186 0.006693 0.0121 0.12% MLMG::prepareForSolve() 9 0.01002 0.01048 0.01175 0.12% MLMG::computeMLResidual() 9 0.005454 0.007779 0.01052 0.11% MLMG::mgVcycle_up::5 65 0.005675 0.006093 0.0103 0.11% FabArray::mult() 40 0.007975 0.008114 0.008458 0.09% FabArrayBase::getCPC() 1022 0.005301 0.006001 0.007924 0.08% FabArray::norminf() 143 0.005614 0.006017 0.007908 0.08% Amr::writeSmallPlotFile() 2 0.007705 0.007707 0.00771 0.08% amrex::Add() 65 0.007257 0.007385 0.007494 0.08% Castro::volWgtSum() 46 0.007272 0.007378 0.007494 0.08% Other 4623 0.04472 0.0643 0.09313 0.96% -------------------------------------------------------------------------------------------- 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 18 B 25 B 38 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.08-9-g6f039ae18e89) finalized