Initializing AMReX (24.07-35-g20e6f2eadf0c)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.07-35-g20e6f2eadf0c) initialized Starting run at 07:24:16 UTC on 2024-07-29. Successfully read inputs file ... Castro git describe: 24.07-34-gdbc6f6b13 AMReX git describe: 24.07-35-g20e6f2ead Microphysics git describe: 24.07-38-g83ce77ce 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.163821394 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001161171 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010372019 Gravity::multilevel_solve_for_new_phi() time = 0.099292746 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.010769841 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.047958783 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003742819 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 = -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 = -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 Castro::construct_old_gravity() time = 0.002202803 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.066620797 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.005687182 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.065455631 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, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 0 Castro::construct_old_gravity() time = 0.004215284 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.062074658 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.004367153 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.061382374 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.304800119 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.016146669 Gravity::solve_for_phi() time = 0.088493433 on level 0 Castro::construct_new_gravity() time = 0.090553555 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.060564693 on level 0 Castro::construct_old_gravity() time = 0.003828749 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.061360725 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.311801717 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.016300691 Gravity::solve_for_phi() time = 0.087302579 on level 0 Castro::construct_new_gravity() time = 0.089316078 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.06018844 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.00378353 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.062003567 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.312203043 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.016388335 Gravity::solve_for_phi() time = 0.086922607 on level 0 Castro::construct_new_gravity() time = 0.088950215 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.060210062 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.003885371 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.062179404 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.313921793 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.017912016 Gravity::solve_for_phi() time = 0.089001743 on level 0 Castro::construct_new_gravity() time = 0.091067718 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.060837305 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.004899834 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.061053629 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.312665254 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.016553215 Gravity::solve_for_phi() time = 0.087453765 on level 0 Castro::construct_new_gravity() time = 0.089481676 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.060128206 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.002202785 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.092483857 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.359246243 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.015295216 Gravity::solve_for_phi() time = 0.088366234 on level 0 Castro::construct_new_gravity() time = 0.090475024 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.061378668 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.002648012 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.060594688 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.309990718 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.015919228 Gravity::solve_for_phi() time = 0.085647695 on level 0 Castro::construct_new_gravity() time = 0.087660559 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.059186888 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.002595789 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.059198011 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.349837924 on level 0  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.004679618777 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.031725144 Gravity::solve_for_phi() time = 0.103768809 on level 0 Castro::construct_new_gravity() time = 0.105780329 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.060439573 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.063307461 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.048402589 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003810889 seconds Ending run at 07:24:26 UTC on 2024-07-29. Run time = 9.408427077 Run time without initialization = 7.116462306 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]: 9.408 ... 9.408 ... 9.408 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7267 2.7745 2.8277 30.05% Amr::InitAmr() 1 1.241 1.241 1.241 13.19% Castro::computeTemp() 63 0.7102 0.7576 0.8562 9.10% amrex::Copy() 523 0.5352 0.5447 0.5573 5.92% FabArray::ParallelCopy_finish() 734 0.1615 0.347 0.5088 5.41% Castro::normalize_species() 62 0.4475 0.4558 0.4661 4.95% Castro::reset_internal_energy(Fab) 126 0.4296 0.4377 0.4484 4.77% MultiFab::AddProduct() 2633 0.4127 0.4199 0.4269 4.54% Castro::initMFs() 1 0.3654 0.3737 0.3942 4.19% MultiFab::Subtract() 125 0.2653 0.2705 0.2865 3.05% FabArray::setVal() 3453 0.2181 0.2236 0.2308 2.45% MLPoisson::Fsmooth() 3120 0.1601 0.1751 0.2134 2.27% check_for_negative_density() 11 0.08071 0.1484 0.2094 2.23% Castro::sum_integrated_quantities() 2 0.002706 0.1594 0.1702 1.81% Castro::problem_post_timestep() 2 1.053e-05 0.01044 0.1668 1.77% Gravity::fill_multipole_BCs() 10 0.1195 0.1491 0.1653 1.76% Castro::initData() 1 0.003263 0.02812 0.1649 1.75% FabArray::sum() 2633 0.1403 0.143 0.1457 1.55% TagBoxArray::collate() 2 0.0002747 0.02544 0.1346 1.43% FabArray::ParallelCopy_nowait() 734 0.0643 0.07003 0.07735 0.82% FillBoundary_finish() 3240 0.04753 0.05763 0.0689 0.73% FillBoundary_nowait() 3240 0.04154 0.0533 0.06686 0.71% MLPoisson::Fapply() 501 0.04808 0.05067 0.05907 0.63% MLCellLinOp::applyBC() 3630 0.02404 0.03244 0.05445 0.58% AmrLevel::derive() 18 0.04074 0.0418 0.04343 0.46% FabArray::Saxpy() 133 0.03542 0.03628 0.04021 0.43% Castro::construct_new_gravity_source() 8 0.03253 0.03309 0.03361 0.36% FabArray::Xpay() 480 0.02547 0.02681 0.03078 0.33% Castro::construct_old_gravity_source() 11 0.02459 0.02506 0.02542 0.27% MultiFab::contains_nan() 19 0.02275 0.02317 0.02366 0.25% VisMF::Write(FabArray) 4 0.01137 0.01264 0.0234 0.25% MLMG::addInterpCorrection() 390 0.01698 0.01841 0.02221 0.24% Castro::reset_internal_energy(MultiFab) 63 0.01765 0.0181 0.01889 0.20% Amr::writePlotFile() 2 0.006288 0.01646 0.01888 0.20% amrex::average_down 390 0.01234 0.01348 0.01664 0.18% Castro::wd_update() 2 0.004918 0.01077 0.01544 0.16% Castro::enforce_min_density() 62 0.01357 0.01381 0.01427 0.15% StateData::FillBoundary(geom) 62 0 0.00719 0.01408 0.15% Castro::evaluate_source_change() 144 0.009371 0.009703 0.01036 0.11% MLCellLinOp::defineAuxData() 10 0.007662 0.00811 0.009248 0.10% FabArray::mult() 40 0.008251 0.008424 0.00864 0.09% FabArray::FillBoundary() 3240 0.006744 0.007102 0.008597 0.09% MLCGSolver::bicgstab 4 0 0.0005146 0.008234 0.09% FabArray::norminf() 143 0.005519 0.005842 0.007924 0.08% Castro::volWgtSum() 46 0.007284 0.00745 0.007887 0.08% Castro::construct_new_source() 48 0.007336 0.007521 0.007863 0.08% MLCellLinOp::smooth() 1560 0.00557 0.005829 0.006524 0.07% Castro::construct_old_source() 66 0.005415 0.005556 0.005856 0.06% FillPatchIterator::Initialize 82 0.005465 0.005577 0.005708 0.06% MLMG::prepareForSolve() 9 0.004833 0.004984 0.005291 0.06% Amr::writeSmallPlotFile() 2 0.003095 0.004736 0.005176 0.06% Gravity::actual_solve_with_mlmg() 10 0.003472 0.003648 0.004121 0.04% FabArrayBase::getCPC() 1022 0.003773 0.003875 0.00404 0.04% MLMG::ResNormInf() 74 0.001628 0.002418 0.003771 0.04% Gravity::get_old_grav_vector() 11 0.003516 0.003601 0.003722 0.04% MLMG::solve() 9 0.0003624 0.001785 0.003302 0.04% Gravity::get_new_grav_vector() 9 0.002924 0.003027 0.003178 0.03% Castro::initialize_do_advance() 11 0.002426 0.002567 0.002876 0.03% MLCellLinOp::apply() 501 0.0009796 0.00112 0.002419 0.03% FabArray::ParallelCopy() 734 0.001822 0.00191 0.002103 0.02% AmrLevel::FillPatch() 82 0.001084 0.001155 0.001325 0.01% MLCellLinOp::correctionResidual() 390 0.0009293 0.000968 0.001036 0.01% MLMG::mgVcycle() 65 0.000761 0.0008086 0.0008935 0.01% main() 1 0.0006261 0.0006817 0.000821 0.01% MLMG:computeResOfCorrection() 390 0.0006964 0.000736 0.0007902 0.01% Castro::save_data_for_retry() 3 0.0007345 0.000758 0.0007898 0.01% Castro::subcycle_advance_ctu() 1 8.49e-05 0.0001305 0.0006774 0.01% Castro::Castro() 1 0.0004185 0.0004973 0.0006258 0.01% AmrMesh::MakeNewGrids() 2 0.000533 0.0005556 0.0006196 0.01% MLMG::actualBottomSolve() 65 6.749e-05 0.0001004 0.0005538 0.01% StateDataPhysBCFunct::() 82 0.0003409 0.0004202 0.0004902 0.01% FillPatchIterator::FillFromLevel0() 82 0.0003657 0.0003804 0.0004055 0.00% MLMG::mgVcycle_down::1 65 0.0003326 0.0003423 0.000367 0.00% MLMG::mgVcycle_down::2 65 0.0003081 0.0003354 0.0003647 0.00% MLMG::mgVcycle_down::0 65 0.0003272 0.0003406 0.0003617 0.00% Castro::clean_state() 62 0.0003239 0.0003376 0.000359 0.00% MLMG::mgVcycle_down::3 65 0.0002875 0.0003119 0.0003586 0.00% MLMG::mgVcycle_down::5 65 0.0002903 0.000303 0.0003462 0.00% MLMG::mgVcycle_down::4 65 0.0002796 0.0002947 0.0003427 0.00% Gravity::solve_for_phi() 8 0.0002767 0.0002922 0.0003247 0.00% FillPatchSingleLevel 82 0.000253 0.0002624 0.0002741 0.00% MLMG::mgVcycle_up::0 65 0.0002335 0.0002445 0.0002738 0.00% MLMG::mgVcycle_up::1 65 0.0002304 0.0002436 0.0002627 0.00% MLCellLinOp::solutionResidual() 75 0.0002272 0.0002352 0.000244 0.00% Castro::do_advance_ctu() 11 0.0002102 0.0002244 0.00024 0.00% MLMG::mgVcycle_up::2 65 0.0001988 0.0002188 0.0002396 0.00% MLMG::mgVcycle_up::3 65 0.0001965 0.0002043 0.0002299 0.00% MLMG::mgVcycle_up::5 65 0.0002027 0.0002094 0.0002291 0.00% Castro::do_new_sources() 8 0.000101 0.000115 0.0002259 0.00% MLMG::mgVcycle_up::4 65 0.0001974 0.000207 0.0002258 0.00% MLMG::oneIter() 65 0.000186 0.000195 0.0002177 0.00% Amr::coarseTimeStep() 1 0.0001323 0.0001696 0.0001784 0.00% Castro::do_old_sources() 11 0.0001576 0.0001673 0.0001782 0.00% Castro::initialize_advance() 1 0.0001139 0.0001246 0.0001542 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001222 0.0001306 0.0001525 0.00% Castro::construct_new_gravity() 8 7.519e-05 8.46e-05 0.0001395 0.00% TagBoxArray::mapPRD 2 0.0001127 0.0001195 0.0001373 0.00% MLMG::mgVcycle_bottom 65 0.0001141 0.0001199 0.0001291 0.00% MLMG::computeResidual() 65 0.0001185 0.0001232 0.0001289 0.00% Castro::retry_advance_ctu() 11 5.874e-05 6.622e-05 0.0001238 0.00% Gravity::test_composite_phi() 1 9.045e-05 9.425e-05 0.0001081 0.00% Amr::defBaseLevel() 1 8.913e-05 9.433e-05 0.0001061 0.00% Amr::grid_places() 2 5.66e-05 6.644e-05 8.886e-05 0.00% MLPoisson::define() 10 6.783e-05 7.164e-05 7.965e-05 0.00% Castro::check_for_nan() 19 5.081e-05 5.555e-05 7.493e-05 0.00% Castro::construct_old_gravity() 11 5.866e-05 6.377e-05 7.041e-05 0.00% Amr::FinalizeInit() 1 1.565e-05 2.717e-05 6.862e-05 0.00% Castro::apply_source_to_state() 19 5.183e-05 5.432e-05 6.038e-05 0.00% Castro::react_state() 19 4.628e-05 5.031e-05 5.909e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 4.933e-05 5.049e-05 5.357e-05 0.00% Amr::regrid() 1 1.551e-05 1.989e-05 4.682e-05 0.00% Castro::derive() 18 4.037e-05 4.304e-05 4.608e-05 0.00% Castro::create_source_corrector() 8 3.481e-05 3.713e-05 4.129e-05 0.00% Castro::expand_state() 11 3.09e-05 3.289e-05 3.475e-05 0.00% Gravity::actual_multilevel_solve() 1 2.68e-05 2.766e-05 3.007e-05 0.00% Amr::init() 1 6.466e-06 1.262e-05 2.92e-05 0.00% MLMG::computeMLResidual() 9 1.949e-05 2.084e-05 2.329e-05 0.00% Amr::InitializeInit() 1 7.977e-06 9.614e-06 1.898e-05 0.00% Amr::timeStep() 1 1.173e-05 1.467e-05 1.783e-05 0.00% Castro::post_regrid() 2 1.093e-05 1.219e-05 1.615e-05 0.00% Castro::post_init() 1 1.159e-05 1.269e-05 1.567e-05 0.00% Amr::bldFineLevels() 1 7.006e-06 8.903e-06 1.448e-05 0.00% Castro::advance() 1 6.675e-06 8.141e-06 1.142e-05 0.00% Amr::initialInit() 1 3.093e-06 4.009e-06 1.06e-05 0.00% Castro::post_timestep() 1 5.533e-06 6.217e-06 7.994e-06 0.00% Other 4688 0.03289 0.0424 0.06116 0.65% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.408 9.408 9.408 100.00% Amr::coarseTimeStep() 1 7.116 7.116 7.116 75.63% Amr::timeStep() 1 7.063 7.063 7.063 75.07% Castro::advance() 1 6.995 6.999 7.003 74.43% Castro::subcycle_advance_ctu() 1 6.946 6.949 6.953 73.90% Castro::do_advance_ctu() 11 6.923 6.926 6.929 73.64% Castro::construct_ctu_hydro_source() 11 3.243 3.291 3.328 35.37% Castro::clean_state() 62 2.928 3.001 3.131 33.28% Castro::computeTemp() 63 1.899 1.956 2.07 22.00% Amr::InitAmr() 1 1.241 1.241 1.241 13.19% Castro::reset_internal_energy(MultiFab) 63 1.178 1.198 1.225 13.02% Amr::init() 1 1.05 1.05 1.05 11.16% Amr::initialInit() 1 0.9983 0.9983 0.9984 10.61% Gravity::solve_phi_with_mlmg() 10 0.7786 0.8104 0.8285 8.81% Castro::initialize_do_advance() 11 0.7446 0.7676 0.7837 8.33% Castro::evaluate_source_change() 144 0.7099 0.7223 0.7344 7.81% Castro::construct_new_gravity() 8 0.6843 0.7155 0.7333 7.79% Gravity::solve_for_phi() 8 0.6672 0.6988 0.717 7.62% Castro::do_old_sources() 11 0.6191 0.668 0.7112 7.56% Gravity::actual_solve_with_mlmg() 10 0.6573 0.6594 0.6614 7.03% Castro::enforce_min_density() 62 0.6064 0.617 0.6381 6.78% MLMG::solve() 9 0.6302 0.6325 0.6334 6.73% Amr::InitializeInit() 1 0.4476 0.4812 0.6096 6.48% Amr::defBaseLevel() 1 0.4476 0.4811 0.6096 6.48% FabArray::ParallelCopy() 734 0.2501 0.4354 0.5944 6.32% MLMG::oneIter() 65 0.5652 0.5665 0.5679 6.04% MLMG::mgVcycle() 65 0.5577 0.5589 0.5604 5.96% amrex::Copy() 523 0.5352 0.5447 0.5573 5.92% Amr::FinalizeInit() 1 0.3887 0.5172 0.5508 5.85% FabArray::ParallelCopy_finish() 734 0.1615 0.347 0.5088 5.41% Castro::do_new_sources() 8 0.4361 0.4544 0.482 5.12% Castro::normalize_species() 62 0.4475 0.4558 0.4661 4.95% Castro::reset_internal_energy(Fab) 126 0.4296 0.4377 0.4484 4.77% MultiFab::AddProduct() 2633 0.4127 0.4199 0.4269 4.54% Castro::Castro() 1 0.3697 0.3781 0.3994 4.24% Castro::initMFs() 1 0.3655 0.3738 0.3943 4.19% MLCellLinOp::smooth() 1560 0.2806 0.3058 0.371 3.94% AmrLevel::FillPatch() 82 0.1973 0.2745 0.3439 3.66% Castro::post_init() 1 0.2984 0.3005 0.3013 3.20% MultiFab::Subtract() 125 0.2653 0.2705 0.2865 3.05% FillPatchIterator::Initialize 82 0.125 0.2034 0.2743 2.92% FillPatchIterator::FillFromLevel0() 82 0.114 0.195 0.2668 2.84% FillPatchSingleLevel 82 0.1137 0.1946 0.2664 2.83% Castro::initData() 1 0.07543 0.1018 0.238 2.53% FabArray::setVal() 3453 0.2181 0.2236 0.2308 2.45% MLMG::mgVcycle_up::1 65 0.04927 0.1339 0.2191 2.33% MLPoisson::Fsmooth() 3120 0.1601 0.1751 0.2134 2.27% check_for_negative_density() 11 0.08071 0.1484 0.2094 2.23% Amr::bldFineLevels() 1 0.04605 0.1743 0.2079 2.21% MLMG::addInterpCorrection() 390 0.03318 0.1563 0.1927 2.05% MLCellLinOp::applyBC() 3630 0.1376 0.1564 0.1924 2.04% Castro::sum_integrated_quantities() 2 0.02222 0.1785 0.1889 2.01% Castro::problem_post_timestep() 2 0.009201 0.02488 0.177 1.88% Gravity::fill_multipole_BCs() 10 0.1206 0.1502 0.1663 1.77% Amr::grid_places() 2 0.002917 0.1314 0.165 1.75% AmrMesh::MakeNewGrids() 2 0.002853 0.1313 0.1649 1.75% TagBoxArray::mapPRD 2 0.001255 0.1046 0.1588 1.69% FabArray::FillBoundary() 3240 0.1067 0.124 0.1481 1.57% Castro::expand_state() 11 0.1 0.1202 0.1464 1.56% FabArray::sum() 2633 0.1403 0.143 0.1457 1.55% MLMG::mgVcycle_up::2 65 0.005477 0.05996 0.1428 1.52% MLMG::mgVcycle_down::0 65 0.1323 0.1337 0.1347 1.43% TagBoxArray::collate() 2 0.0002747 0.02544 0.1346 1.43% MLCellLinOp::apply() 501 0.07647 0.08236 0.1018 1.08% MLMG::mgVcycle_up::0 65 0.09924 0.09995 0.101 1.07% Gravity::multilevel_solve_for_new_phi() 1 0.09924 0.09926 0.0993 1.06% Gravity::actual_multilevel_solve() 1 0.09919 0.09921 0.09925 1.05% Amr::writePlotFile() 2 0.09664 0.09672 0.09674 1.03% FabArray::ParallelCopy_nowait() 734 0.08059 0.08645 0.09398 1.00% MLMG::mgVcycle_up::3 65 0.005416 0.01218 0.0927 0.99% Castro::post_regrid() 2 0.08439 0.08523 0.08736 0.93% MLMG:computeResOfCorrection() 390 0.05743 0.06375 0.07691 0.82% MLCellLinOp::correctionResidual() 390 0.05673 0.06301 0.07612 0.81% FillBoundary_nowait() 3240 0.04711 0.05923 0.07306 0.78% Castro::derive() 18 0.06859 0.07015 0.07302 0.78% AmrLevel::derive() 18 0.06855 0.07011 0.07297 0.78% FillBoundary_finish() 3240 0.04753 0.05763 0.0689 0.73% Castro::post_timestep() 1 0.06018 0.06458 0.06784 0.72% MLMG::mgVcycle_down::1 65 0.05841 0.06042 0.06213 0.66% Castro::construct_new_source() 48 0.05827 0.05921 0.06028 0.64% MLPoisson::Fapply() 501 0.04808 0.05067 0.05907 0.63% Castro::initialize_advance() 1 0.04813 0.04936 0.05116 0.54% MLCellLinOp::solutionResidual() 75 0.04434 0.04653 0.04921 0.52% Amr::regrid() 1 0.04345 0.04404 0.04612 0.49% Castro::construct_old_source() 66 0.04259 0.04348 0.04463 0.47% Castro::create_source_corrector() 8 0.01565 0.02859 0.04409 0.47% FabArray::Saxpy() 133 0.03542 0.03628 0.04021 0.43% MLMG::computeResidual() 65 0.03755 0.03844 0.03901 0.41% MLMG::mgVcycle_bottom 65 0.0001816 0.002496 0.03709 0.39% MLMG::actualBottomSolve() 65 6.749e-05 0.002376 0.03696 0.39% MLMG::mgVcycle_down::2 65 0.006366 0.02068 0.03622 0.39% MLCGSolver::bicgstab 4 0 0.002261 0.03617 0.38% Castro::construct_new_gravity_source() 8 0.03253 0.03309 0.03361 0.36% Castro::construct_old_gravity() 11 0.02237 0.02766 0.03198 0.34% Gravity::get_old_grav_vector() 11 0.02231 0.0276 0.03192 0.34% FabArray::Xpay() 480 0.02547 0.02681 0.03078 0.33% MLMG::mgVcycle_down::3 65 0.006023 0.009226 0.0305 0.32% amrex::average_down 390 0.01667 0.02003 0.02811 0.30% Castro::construct_old_gravity_source() 11 0.02459 0.02506 0.02542 0.27% Castro::retry_advance_ctu() 11 0.02166 0.02219 0.02427 0.26% Castro::check_for_nan() 19 0.0228 0.02322 0.02372 0.25% MultiFab::contains_nan() 19 0.02275 0.02317 0.02366 0.25% VisMF::Write(FabArray) 4 0.01137 0.01264 0.0234 0.25% MLPoisson::define() 10 0.01543 0.018 0.02063 0.22% Castro::wd_update() 2 0.009188 0.01444 0.01882 0.20% MLMG::mgVcycle_down::4 65 0.006006 0.00694 0.0181 0.19% Gravity::get_new_grav_vector() 9 0.01712 0.01759 0.01797 0.19% Castro::react_state() 19 0.01451 0.01501 0.01741 0.19% Castro::apply_source_to_state() 19 0.01466 0.015 0.01604 0.17% StateDataPhysBCFunct::() 82 0.0003409 0.007611 0.01457 0.15% Gravity::test_composite_phi() 1 0.01423 0.01425 0.01427 0.15% MLMG::mgVcycle_up::4 65 0.005475 0.006143 0.01408 0.15% StateData::FillBoundary(geom) 62 0 0.00719 0.01408 0.15% MLCellLinOp::defineAuxData() 10 0.00988 0.01077 0.013 0.14% Castro::save_data_for_retry() 3 0.0108 0.01112 0.01272 0.14% MLMG::mgVcycle_down::5 65 0.006008 0.006522 0.01201 0.13% MLMG::prepareForSolve() 9 0.009878 0.01036 0.01153 0.12% MLMG::mgVcycle_up::5 65 0.005552 0.005975 0.01025 0.11% MLMG::computeMLResidual() 9 0.005479 0.007514 0.009981 0.11% FabArray::mult() 40 0.008251 0.008424 0.00864 0.09% MLMG::ResNormInf() 74 0.006577 0.007303 0.00857 0.09% FabArray::norminf() 143 0.005519 0.005842 0.007924 0.08% Castro::volWgtSum() 46 0.007284 0.00745 0.007887 0.08% FabArrayBase::getCPC() 1022 0.005157 0.005925 0.007804 0.08% Amr::writeSmallPlotFile() 2 0.007777 0.007787 0.007792 0.08% Other 4688 0.05109 0.06753 0.09374 1.00% -------------------------------------------------------------------------------------------- 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 20 B 28 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.07-35-g20e6f2eadf0c) finalized