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 23:01:19 UTC on 2024-07-28. Successfully read inputs file ... Castro git describe: 24.07-39-g8949302c6 AMReX git describe: 24.07-35-g20e6f2ead Microphysics git describe: 24.07-36-gc75486b6 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.176791465 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001163054 new finest: 0 Gravity::fill_multipole_BCs() time = 0.011636254 Gravity::multilevel_solve_for_new_phi() time = 0.102842092 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.011969003 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.049177481 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003781372 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 = -8.076810125 at index 20, 174, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 20, 81, 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.002492658 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.066454434 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.007675139 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.065486259 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.003111994 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.062781502 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.004504424 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.064137551 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.303297815 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.015835142 Gravity::solve_for_phi() time = 0.088844009 on level 0 Castro::construct_new_gravity() time = 0.090982831 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.061624828 on level 0 Castro::construct_old_gravity() time = 0.002253066 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.062081184 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.307672588 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.016180866 Gravity::solve_for_phi() time = 0.087807767 on level 0 Castro::construct_new_gravity() time = 0.089983468 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.06120376 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.00232136 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.061391777 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.335002515 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.044107562 Gravity::solve_for_phi() time = 0.197884574 on level 0 Castro::construct_new_gravity() time = 0.204799399 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.091212498 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.057163887 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.062858162 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.32185396 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.018057391 Gravity::solve_for_phi() time = 0.093452967 on level 0 Castro::construct_new_gravity() time = 0.095587137 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.061604035 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.003970428 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.065055086 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.345912821 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.040476352 Gravity::solve_for_phi() time = 0.118494965 on level 0 Castro::construct_new_gravity() time = 0.120682444 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.064608272 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.005957984 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.061576332 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.330007585 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.017967656 Gravity::solve_for_phi() time = 0.091986677 on level 0 Castro::construct_new_gravity() time = 0.094132289 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.061534185 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.004262351 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.062019837 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.317761674 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.01656091 Gravity::solve_for_phi() time = 0.093851344 on level 0 Castro::construct_new_gravity() time = 0.096010717 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.061243222 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.002643736 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.062173971 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.323633453 on level 0  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.00451104808 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.019211286 Gravity::solve_for_phi() time = 0.090858888 on level 0 Castro::construct_new_gravity() time = 0.093007626 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.061303006 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.324691361 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.04782559 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003794805 seconds Ending run at 23:01:29 UTC on 2024-07-28. Run time = 9.584700054 Run time without initialization = 7.37712199 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.585 ... 9.585 ... 9.585 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7059 2.7654 2.8313 29.54% Amr::InitAmr() 1 1.3 1.3 1.3 13.57% Castro::computeTemp() 63 0.7005 0.7671 0.8628 9.00% amrex::Copy() 523 0.5819 0.5913 0.6135 6.40% FabArray::ParallelCopy_finish() 734 0.1863 0.4409 0.6077 6.34% Castro::normalize_species() 62 0.4622 0.4696 0.4949 5.16% Castro::initMFs() 1 0.3618 0.3853 0.4903 5.12% Castro::reset_internal_energy(Fab) 126 0.434 0.442 0.4541 4.74% MultiFab::AddProduct() 2633 0.4296 0.4361 0.4433 4.62% MultiFab::Subtract() 125 0.2779 0.2846 0.3242 3.38% MLPoisson::Fsmooth() 3120 0.1617 0.1794 0.2588 2.70% FabArray::setVal() 3453 0.2284 0.233 0.2426 2.53% Gravity::fill_multipole_BCs() 10 0.1522 0.1909 0.2101 2.19% check_for_negative_density() 11 0.04203 0.1164 0.1969 2.05% Castro::initData() 1 0.003361 0.02858 0.1674 1.75% FabArray::sum() 2633 0.1404 0.1426 0.1451 1.51% FillBoundary_finish() 3240 0.08057 0.1133 0.1354 1.41% FabArray::ParallelCopy_nowait() 734 0.06552 0.07193 0.08143 0.85% FillBoundary_nowait() 3240 0.04199 0.0544 0.06953 0.73% MLPoisson::Fapply() 501 0.04844 0.05209 0.06794 0.71% MLCellLinOp::applyBC() 3630 0.02526 0.0341 0.06356 0.66% TagBoxArray::collate() 2 0.0002769 0.01157 0.04545 0.47% Castro::construct_new_gravity_source() 8 0.03199 0.03353 0.04491 0.47% AmrLevel::derive() 18 0.04167 0.04208 0.04246 0.44% FabArray::Saxpy() 133 0.03419 0.03515 0.04113 0.43% FabArray::Xpay() 480 0.02626 0.02788 0.0347 0.36% Castro::construct_old_gravity_source() 11 0.02581 0.02611 0.02642 0.28% MLMG::addInterpCorrection() 390 0.01727 0.01897 0.02627 0.27% MultiFab::contains_nan() 19 0.02248 0.02302 0.02573 0.27% VisMF::Write(FabArray) 4 0.01174 0.01273 0.0217 0.23% Castro::reset_internal_energy(MultiFab) 63 0.01801 0.01852 0.01905 0.20% amrex::average_down 390 0.01259 0.01394 0.0188 0.20% MLMG::ResNormInf() 74 0.005933 0.01143 0.01832 0.19% Castro::wd_update() 2 0.00505 0.01125 0.01713 0.18% Amr::writePlotFile() 2 0.006159 0.01436 0.01557 0.16% Castro::enforce_min_density() 62 0.0138 0.01413 0.01449 0.15% StateData::FillBoundary(geom) 62 0 0.007186 0.01401 0.15% Castro::evaluate_source_change() 144 0.009841 0.01017 0.01059 0.11% MLCellLinOp::defineAuxData() 10 0.00773 0.008113 0.009379 0.10% FabArray::FillBoundary() 3240 0.006829 0.007241 0.009223 0.10% FabArray::mult() 40 0.008152 0.008328 0.009113 0.10% amrex::Add() 65 0.007371 0.007562 0.009055 0.09% MLCGSolver::bicgstab 4 0 0.0005636 0.009018 0.09% FabArray::norminf() 143 0.005593 0.005983 0.008951 0.09% Castro::construct_new_source() 48 0.007303 0.007465 0.007636 0.08% FabArrayBase::getFB() 3240 0.005243 0.005552 0.007309 0.08% MLCellLinOp::smooth() 1560 0.00558 0.005933 0.006787 0.07% FillPatchIterator::Initialize 82 0.005594 0.005787 0.006043 0.06% Castro::construct_old_source() 66 0.005475 0.0056 0.0057 0.06% MLMG::prepareForSolve() 9 0.004931 0.005147 0.005497 0.06% MLMG::solve() 9 0.0005259 0.001981 0.004553 0.05% Gravity::actual_solve_with_mlmg() 10 0.003515 0.00373 0.004503 0.05% Gravity::get_new_grav_vector() 9 0.00302 0.003152 0.004392 0.05% FabArrayBase::getCPC() 1022 0.003842 0.003981 0.004304 0.04% Gravity::get_old_grav_vector() 11 0.003572 0.003649 0.003728 0.04% Castro::initialize_do_advance() 11 0.002507 0.002621 0.002797 0.03% MLCellLinOp::apply() 501 0.001016 0.001157 0.002618 0.03% Castro::sum_integrated_quantities() 2 0.001659 0.001998 0.002211 0.02% FabArray::ParallelCopy() 734 0.001778 0.001889 0.002177 0.02% AmrLevel::FillPatch() 82 0.001144 0.001187 0.001254 0.01% MLCellLinOp::correctionResidual() 390 0.0009484 0.0009956 0.001103 0.01% MLMG::mgVcycle() 65 0.00078 0.0008232 0.000951 0.01% Gravity::solve_for_phi() 8 0.0002828 0.0003377 0.0008437 0.01% Castro::save_data_for_retry() 3 0.0007848 0.0008002 0.0008184 0.01% MLMG:computeResOfCorrection() 390 0.0007015 0.0007394 0.0008164 0.01% AmrMesh::MakeNewGrids() 2 0.0005555 0.0005742 0.0007435 0.01% main() 1 0.0006613 0.0006865 0.0007346 0.01% Castro::subcycle_advance_ctu() 1 8.181e-05 0.0001359 0.0006961 0.01% MLMG::actualBottomSolve() 65 6.875e-05 0.0001071 0.0006112 0.01% Castro::Castro() 1 0.000419 0.0004957 0.0005746 0.01% StateDataPhysBCFunct::() 82 0.0003655 0.0004317 0.00051 0.01% FillPatchIterator::FillFromLevel0() 82 0.0003689 0.0003926 0.0004338 0.00% Castro::clean_state() 62 0.0003387 0.0003688 0.0004039 0.00% MLMG::mgVcycle_down::3 65 0.0002945 0.0003211 0.0003819 0.00% MLMG::mgVcycle_down::2 65 0.0003119 0.0003448 0.000378 0.00% MLMG::mgVcycle_down::4 65 0.0002859 0.0003018 0.0003689 0.00% MLMG::mgVcycle_down::1 65 0.0003326 0.000347 0.0003675 0.00% MLMG::mgVcycle_down::0 65 0.0003338 0.0003488 0.000367 0.00% MLMG::mgVcycle_down::5 65 0.0002954 0.0003104 0.0003648 0.00% FillPatchSingleLevel 82 0.000257 0.0002715 0.0002893 0.00% Castro::do_advance_ctu() 11 0.0002189 0.0002443 0.0002733 0.00% MLMG::mgVcycle_up::1 65 0.0002318 0.0002445 0.0002652 0.00% MLMG::mgVcycle_up::0 65 0.0002332 0.0002434 0.0002617 0.00% MLCellLinOp::solutionResidual() 75 0.0002283 0.0002402 0.0002584 0.00% MLMG::mgVcycle_up::2 65 0.0001957 0.000219 0.0002473 0.00% MLMG::mgVcycle_up::4 65 0.0001977 0.0002057 0.000247 0.00% Castro::problem_post_timestep() 2 1.179e-05 2.75e-05 0.000247 0.00% MLMG::mgVcycle_up::3 65 0.0001941 0.0002049 0.0002462 0.00% MLMG::mgVcycle_up::5 65 0.0002026 0.0002107 0.0002462 0.00% MLMG::oneIter() 65 0.0001917 0.0002011 0.0002118 0.00% Castro::do_old_sources() 11 0.0001656 0.0001815 0.0001956 0.00% Amr::coarseTimeStep() 1 0.0001313 0.0001614 0.0001698 0.00% Castro::initialize_advance() 1 0.0001184 0.0001262 0.0001629 0.00% Castro::construct_new_gravity() 8 8.038e-05 9.09e-05 0.0001626 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001287 0.0001412 0.0001608 0.00% Castro::retry_advance_ctu() 11 5.197e-05 6.96e-05 0.0001558 0.00% MLMG::mgVcycle_bottom 65 0.0001139 0.0001205 0.0001344 0.00% MLMG::computeResidual() 65 0.0001208 0.0001264 0.0001323 0.00% Castro::do_new_sources() 8 0.000102 0.0001124 0.0001323 0.00% TagBoxArray::mapPRD 2 0.0001156 0.0001249 0.0001308 0.00% Amr::defBaseLevel() 1 9.219e-05 0.0001028 0.0001162 0.00% Gravity::test_composite_phi() 1 9.596e-05 9.989e-05 0.0001129 0.00% Amr::grid_places() 2 6.213e-05 6.863e-05 7.882e-05 0.00% Amr::FinalizeInit() 1 1.538e-05 3.157e-05 7.839e-05 0.00% MLPoisson::define() 10 6.673e-05 7.114e-05 7.725e-05 0.00% Castro::check_for_nan() 19 5.462e-05 6.29e-05 7.359e-05 0.00% Castro::construct_old_gravity() 11 6.163e-05 6.75e-05 7.28e-05 0.00% Castro::react_state() 19 4.751e-05 5.24e-05 6.349e-05 0.00% Castro::apply_source_to_state() 19 5.218e-05 5.572e-05 6.107e-05 0.00% Amr::regrid() 1 1.714e-05 2.227e-05 5.647e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.163e-05 5.252e-05 5.422e-05 0.00% Castro::derive() 18 3.934e-05 4.252e-05 4.798e-05 0.00% Castro::create_source_corrector() 8 3.706e-05 4.119e-05 4.509e-05 0.00% Castro::expand_state() 11 3.117e-05 3.465e-05 3.812e-05 0.00% Gravity::actual_multilevel_solve() 1 2.769e-05 2.923e-05 3.052e-05 0.00% MLMG::computeMLResidual() 9 1.981e-05 2.157e-05 2.712e-05 0.00% Amr::InitializeInit() 1 7.818e-06 1.015e-05 2.241e-05 0.00% Castro::post_regrid() 2 1.118e-05 1.307e-05 2.241e-05 0.00% Castro::post_init() 1 1.33e-05 1.482e-05 2.037e-05 0.00% Amr::timeStep() 1 1.194e-05 1.402e-05 1.822e-05 0.00% Amr::init() 1 7.293e-06 1.185e-05 1.414e-05 0.00% Amr::initialInit() 1 3.034e-06 4.861e-06 1.407e-05 0.00% Amr::bldFineLevels() 1 8.274e-06 9.823e-06 1.206e-05 0.00% Castro::advance() 1 6.755e-06 8.044e-06 1.146e-05 0.00% Castro::post_timestep() 1 5.458e-06 6.227e-06 7.652e-06 0.00% Other 1431 0.03132 0.04149 0.06 0.63% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.585 9.585 9.585 100.00% Amr::coarseTimeStep() 1 7.377 7.377 7.377 76.96% Amr::timeStep() 1 7.325 7.325 7.325 76.42% Castro::advance() 1 7.255 7.26 7.264 75.79% Castro::subcycle_advance_ctu() 1 7.204 7.208 7.213 75.25% Castro::do_advance_ctu() 11 7.175 7.182 7.187 74.99% Castro::construct_ctu_hydro_source() 11 3.236 3.262 3.293 34.36% Castro::clean_state() 62 2.995 3.098 3.234 33.74% Castro::computeTemp() 63 1.927 2.012 2.115 22.07% Amr::InitAmr() 1 1.3 1.3 1.3 13.57% Castro::reset_internal_energy(MultiFab) 63 1.223 1.245 1.274 13.29% Gravity::solve_phi_with_mlmg() 10 0.9196 0.9595 0.979 10.21% Amr::init() 1 0.907 0.9071 0.9071 9.46% Castro::construct_new_gravity() 8 0.8272 0.8628 0.8827 9.21% Castro::initialize_do_advance() 11 0.7768 0.849 0.8787 9.17% Gravity::solve_for_phi() 8 0.8034 0.8432 0.8628 9.00% Amr::initialInit() 1 0.8537 0.8537 0.8538 8.91% Gravity::actual_solve_with_mlmg() 10 0.7646 0.7667 0.7688 8.02% Castro::evaluate_source_change() 144 0.7333 0.7442 0.7564 7.89% MLMG::solve() 9 0.737 0.7393 0.7402 7.72% FabArray::ParallelCopy() 734 0.2832 0.5323 0.6971 7.27% Castro::do_old_sources() 11 0.6465 0.6757 0.6925 7.22% MLMG::oneIter() 65 0.6494 0.6563 0.6628 6.92% Castro::enforce_min_density() 62 0.6341 0.6441 0.655 6.83% MLMG::mgVcycle() 65 0.6418 0.6486 0.6535 6.82% Amr::InitializeInit() 1 0.4476 0.4965 0.6227 6.50% Amr::defBaseLevel() 1 0.4476 0.4965 0.6227 6.50% amrex::Copy() 523 0.5819 0.5913 0.6135 6.40% FabArray::ParallelCopy_finish() 734 0.1863 0.4409 0.6077 6.34% Castro::do_new_sources() 8 0.4334 0.4656 0.5168 5.39% Castro::Castro() 1 0.3663 0.39 0.4963 5.18% Castro::normalize_species() 62 0.4622 0.4696 0.4949 5.16% Castro::initMFs() 1 0.362 0.3855 0.4906 5.12% Castro::reset_internal_energy(Fab) 126 0.434 0.442 0.4541 4.74% MultiFab::AddProduct() 2633 0.4296 0.4361 0.4433 4.62% MLCellLinOp::smooth() 1560 0.3252 0.3577 0.4418 4.61% Amr::FinalizeInit() 1 0.231 0.3572 0.4061 4.24% AmrLevel::FillPatch() 82 0.2233 0.3361 0.405 4.23% FillPatchIterator::Initialize 82 0.146 0.2596 0.3294 3.44% MultiFab::Subtract() 125 0.2779 0.2846 0.3242 3.38% FillPatchIterator::FillFromLevel0() 82 0.1347 0.2509 0.3216 3.36% FillPatchSingleLevel 82 0.1343 0.2505 0.3212 3.35% MLMG::mgVcycle_up::1 65 0.05273 0.1558 0.2591 2.70% MLPoisson::Fsmooth() 3120 0.1617 0.1794 0.2588 2.70% MLCellLinOp::applyBC() 3630 0.1919 0.215 0.2465 2.57% Castro::initData() 1 0.07926 0.1051 0.244 2.55% FabArray::setVal() 3453 0.2284 0.233 0.2426 2.53% MLMG::addInterpCorrection() 390 0.03799 0.1858 0.2293 2.39% Amr::bldFineLevels() 1 0.04773 0.174 0.2229 2.33% Gravity::fill_multipole_BCs() 10 0.1534 0.1921 0.2112 2.20% FabArray::FillBoundary() 3240 0.1562 0.1809 0.2093 2.18% Castro::expand_state() 11 0.1086 0.1526 0.2019 2.11% check_for_negative_density() 11 0.04203 0.1164 0.1969 2.05% Amr::grid_places() 2 0.002881 0.1291 0.178 1.86% AmrMesh::MakeNewGrids() 2 0.002817 0.129 0.1779 1.86% TagBoxArray::mapPRD 2 0.001227 0.1161 0.1762 1.84% MLMG::mgVcycle_up::2 65 0.005459 0.07085 0.1711 1.79% MLMG::mgVcycle_down::0 65 0.1449 0.153 0.16 1.67% FabArray::sum() 2633 0.1404 0.1426 0.1451 1.51% Castro::post_init() 1 0.1371 0.1384 0.1386 1.45% FillBoundary_finish() 3240 0.08057 0.1133 0.1354 1.41% MLMG::mgVcycle_up::0 65 0.1143 0.1188 0.1219 1.27% MLCellLinOp::apply() 501 0.08825 0.09497 0.1132 1.18% MLMG::mgVcycle_up::3 65 0.005471 0.01351 0.1102 1.15% Gravity::multilevel_solve_for_new_phi() 1 0.1028 0.1028 0.1028 1.07% Gravity::actual_multilevel_solve() 1 0.1027 0.1027 0.1028 1.07% FabArray::ParallelCopy_nowait() 734 0.08292 0.08953 0.0997 1.04% Castro::create_source_corrector() 8 0.01965 0.05874 0.09851 1.03% Amr::writePlotFile() 2 0.09727 0.09737 0.09738 1.02% Castro::post_regrid() 2 0.08941 0.08964 0.09092 0.95% MLMG:computeResOfCorrection() 390 0.06378 0.0703 0.08669 0.90% MLCellLinOp::correctionResidual() 390 0.06307 0.06956 0.08587 0.90% Castro::construct_old_gravity() 11 0.02261 0.03212 0.08545 0.89% Gravity::get_old_grav_vector() 11 0.02254 0.03206 0.08538 0.89% MLMG::mgVcycle_down::1 65 0.06809 0.07257 0.08126 0.85% FillBoundary_nowait() 3240 0.04766 0.06042 0.07592 0.79% Castro::derive() 18 0.07206 0.07267 0.07312 0.76% AmrLevel::derive() 18 0.07202 0.07263 0.07308 0.76% Castro::construct_new_source() 48 0.0574 0.05963 0.07304 0.76% Castro::post_timestep() 1 0.06022 0.06475 0.06949 0.73% MLPoisson::Fapply() 501 0.04844 0.05209 0.06794 0.71% MLCellLinOp::solutionResidual() 75 0.04878 0.05361 0.05668 0.59% Castro::initialize_advance() 1 0.05097 0.05134 0.05272 0.55% Amr::regrid() 1 0.04595 0.04609 0.0468 0.49% MLMG::computeResidual() 65 0.04204 0.04368 0.04569 0.48% TagBoxArray::collate() 2 0.0002769 0.01157 0.04545 0.47% Castro::construct_new_gravity_source() 8 0.03199 0.03353 0.04491 0.47% Castro::construct_old_source() 66 0.04383 0.0444 0.04488 0.47% MLMG::mgVcycle_down::2 65 0.006437 0.02425 0.04445 0.46% FabArray::Saxpy() 133 0.03419 0.03515 0.04113 0.43% MLMG::mgVcycle_bottom 65 0.0001828 0.002732 0.04082 0.43% MLMG::actualBottomSolve() 65 6.875e-05 0.002612 0.04068 0.42% MLCGSolver::bicgstab 4 0 0.002487 0.03979 0.42% MLMG::mgVcycle_down::3 65 0.006037 0.01003 0.0367 0.38% FabArray::Xpay() 480 0.02626 0.02788 0.0347 0.36% amrex::average_down 390 0.01686 0.02111 0.03154 0.33% Castro::retry_advance_ctu() 11 0.02335 0.02478 0.03045 0.32% Castro::construct_old_gravity_source() 11 0.02581 0.02611 0.02642 0.28% Castro::check_for_nan() 19 0.02254 0.02308 0.0258 0.27% MultiFab::contains_nan() 19 0.02248 0.02302 0.02573 0.27% Gravity::get_new_grav_vector() 9 0.018 0.02042 0.02424 0.25% MLMG::ResNormInf() 74 0.01162 0.01644 0.02322 0.24% Castro::sum_integrated_quantities() 2 0.0216 0.02172 0.02174 0.23% VisMF::Write(FabArray) 4 0.01174 0.01273 0.0217 0.23% MLMG::mgVcycle_down::4 65 0.006026 0.007173 0.02155 0.22% MLPoisson::define() 10 0.01554 0.01815 0.02115 0.22% Castro::problem_post_timestep() 2 0.009643 0.01496 0.02054 0.21% Castro::wd_update() 2 0.009631 0.01493 0.02053 0.21% Castro::save_data_for_retry() 3 0.01197 0.01292 0.01825 0.19% MLMG::mgVcycle_up::4 65 0.005472 0.006331 0.01705 0.18% Castro::react_state() 19 0.01518 0.01538 0.01564 0.16% Gravity::test_composite_phi() 1 0.01551 0.01555 0.01557 0.16% Castro::apply_source_to_state() 19 0.01421 0.01443 0.01465 0.15% StateDataPhysBCFunct::() 82 0.0003655 0.007618 0.01448 0.15% StateData::FillBoundary(geom) 62 0 0.007186 0.01401 0.15% MLMG::mgVcycle_down::5 65 0.00605 0.006664 0.0136 0.14% MLCellLinOp::defineAuxData() 10 0.009999 0.01081 0.01314 0.14% MLMG::prepareForSolve() 9 0.01019 0.01095 0.01282 0.13% MLMG::computeMLResidual() 9 0.005901 0.009309 0.01214 0.13% MLMG::mgVcycle_up::5 65 0.005577 0.006042 0.01134 0.12% FabArray::mult() 40 0.008152 0.008328 0.009113 0.10% amrex::Add() 65 0.007371 0.007562 0.009055 0.09% FabArray::norminf() 143 0.005593 0.005983 0.008951 0.09% FabArrayBase::getFB() 3240 0.005517 0.00602 0.008209 0.09% FabArrayBase::getCPC() 1022 0.005251 0.00607 0.008166 0.09% Other 1431 0.05454 0.06916 0.09453 0.99% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 18 B 19 B 22 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.07-35-g20e6f2eadf0c) finalized