Initializing AMReX (24.08-11-gb6ed5237987e)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.08-11-gb6ed5237987e) initialized Starting run at 07:05:13 UTC on 2024-08-09. Successfully read inputs file ... Castro git describe: 24.08-1-g36c0aaf32 AMReX git describe: 24.08-11-gb6ed52379 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.199137149 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001107638 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010672682 Gravity::multilevel_solve_for_new_phi() time = 0.105013568 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.010631095 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.050257143 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.006807433 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.002339371 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.066856442 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.006412203 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.067478298 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.00539801 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.062961125 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.028148797 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.081857754 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.331985256 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.017236442 Gravity::solve_for_phi() time = 0.091043413 on level 0 Castro::construct_new_gravity() time = 0.093230485 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.060858921 on level 0 Castro::construct_old_gravity() time = 0.003790925 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.066124244 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.319494845 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.015866481 Gravity::solve_for_phi() time = 0.121014329 on level 0 Castro::construct_new_gravity() time = 0.124992561 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.112990391 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.050165458 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.10803045 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.361112857 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.01669401 Gravity::solve_for_phi() time = 0.090938626 on level 0 Castro::construct_new_gravity() time = 0.09313817 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.101971297 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.004423297 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.059406808 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.352185014 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.020359221 Gravity::solve_for_phi() time = 0.096294244 on level 0 Castro::construct_new_gravity() time = 0.098407657 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.060356963 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.002714226 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.061886427 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.335585244 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.016714097 Gravity::solve_for_phi() time = 0.147106277 on level 0 Castro::construct_new_gravity() time = 0.149339608 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.078212839 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.01119419 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.063589379 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.315424697 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.016301949 Gravity::solve_for_phi() time = 0.089760879 on level 0 Castro::construct_new_gravity() time = 0.091861271 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.060268595 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.002884077 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.06256947 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.314534483 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.015997585 Gravity::solve_for_phi() time = 0.088905922 on level 0 Castro::construct_new_gravity() time = 0.091015064 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.060223374 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.002275585 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.06062693 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.316699565 on level 0  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.004318227572 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.016295093 Gravity::solve_for_phi() time = 0.08927269 on level 0 Castro::construct_new_gravity() time = 0.091400264 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.060195964 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.64852051 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.049542238 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006080519 seconds Ending run at 07:05:23 UTC on 2024-08-09. Run time = 10.33678593 Run time without initialization = 7.705088652 Average number of zones advanced per microsecond: 0.004 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 10.34 ... 10.34 ... 10.34 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7319 2.8001 2.9072 28.12% Amr::InitAmr() 1 1.708 1.708 1.708 16.52% FabArray::ParallelCopy_finish() 734 0.3668 0.6552 0.8806 8.52% Castro::computeTemp() 63 0.712 0.7673 0.8724 8.44% amrex::Copy() 523 0.5818 0.5982 0.6352 6.15% Castro::normalize_species() 62 0.4474 0.4659 0.5075 4.91% Castro::reset_internal_energy(Fab) 126 0.4287 0.4433 0.4713 4.56% MultiFab::AddProduct() 2633 0.4219 0.4337 0.4703 4.55% Castro::initMFs() 1 0.3562 0.3645 0.3866 3.74% MultiFab::Subtract() 125 0.2752 0.283 0.2998 2.90% check_for_negative_density() 11 0.1217 0.2283 0.2919 2.82% FabArray::setVal() 3453 0.2159 0.2228 0.2438 2.36% MLPoisson::Fsmooth() 3120 0.1631 0.1821 0.2175 2.10% Castro::initData() 1 0.003281 0.02836 0.1682 1.63% TagBoxArray::collate() 2 0.0001833 0.03974 0.1677 1.62% FabArray::sum() 2633 0.1384 0.1441 0.1589 1.54% Gravity::fill_multipole_BCs() 10 0.1099 0.138 0.155 1.50% FillBoundary_finish() 3240 0.08857 0.1286 0.1503 1.45% FabArray::ParallelCopy_nowait() 734 0.06445 0.07355 0.08445 0.82% FillBoundary_nowait() 3240 0.04252 0.05798 0.0741 0.72% MLPoisson::Fapply() 501 0.04835 0.05222 0.06091 0.59% MLCellLinOp::applyBC() 3630 0.02611 0.03433 0.05613 0.54% FabArray::Saxpy() 133 0.03554 0.03795 0.04693 0.45% AmrLevel::derive() 18 0.04091 0.04177 0.04259 0.41% Castro::construct_new_gravity_source() 8 0.03216 0.0337 0.03844 0.37% FabArray::Xpay() 480 0.02529 0.02732 0.03222 0.31% Castro::construct_old_gravity_source() 11 0.02445 0.02491 0.02549 0.25% MultiFab::contains_nan() 19 0.02305 0.02383 0.02507 0.24% VisMF::Write(FabArray) 4 0.01429 0.01544 0.0245 0.24% MLMG::addInterpCorrection() 390 0.0171 0.01893 0.02236 0.22% Castro::reset_internal_energy(MultiFab) 63 0.01797 0.01851 0.01913 0.19% MLMG::ResNormInf() 74 0.006248 0.01411 0.01818 0.18% Amr::writePlotFile() 2 0.007598 0.01645 0.01786 0.17% amrex::average_down 390 0.01247 0.0139 0.01687 0.16% Castro::enforce_min_density() 62 0.01394 0.01433 0.01495 0.14% Castro::wd_update() 2 0.004935 0.01035 0.01444 0.14% StateData::FillBoundary(geom) 62 0 0.007249 0.01407 0.14% Castro::evaluate_source_change() 144 0.009732 0.01013 0.01101 0.11% MLCellLinOp::defineAuxData() 10 0.007932 0.0083 0.009465 0.09% FabArray::FillBoundary() 3240 0.006974 0.007442 0.009043 0.09% FabArray::mult() 40 0.008042 0.008262 0.009039 0.09% amrex::Add() 65 0.007331 0.007617 0.00891 0.09% Amr::writeSmallPlotFile() 2 0.006617 0.008255 0.008669 0.08% MLCGSolver::bicgstab 4 0 0.0005369 0.008591 0.08% Castro::construct_new_source() 48 0.007332 0.007582 0.008554 0.08% FabArray::norminf() 143 0.005657 0.00607 0.008155 0.08% MLCellLinOp::smooth() 1560 0.005694 0.00611 0.006764 0.07% FillPatchIterator::Initialize 82 0.005677 0.005918 0.006226 0.06% Castro::construct_old_source() 66 0.005398 0.005523 0.005641 0.05% MLMG::prepareForSolve() 9 0.004994 0.005155 0.005593 0.05% FabArrayBase::getCPC() 1022 0.00388 0.004031 0.00422 0.04% Gravity::actual_solve_with_mlmg() 10 0.00347 0.003706 0.004136 0.04% Gravity::get_old_grav_vector() 11 0.003528 0.003636 0.003933 0.04% MLMG::solve() 9 0.0006799 0.001073 0.003788 0.04% Gravity::get_new_grav_vector() 9 0.002929 0.003071 0.003448 0.03% Castro::initialize_do_advance() 11 0.002754 0.003051 0.003309 0.03% MLCellLinOp::apply() 501 0.001023 0.001154 0.002472 0.02% Castro::sum_integrated_quantities() 2 0.001564 0.001998 0.002339 0.02% FabArray::ParallelCopy() 734 0.001788 0.001911 0.002087 0.02% AmrLevel::FillPatch() 82 0.00114 0.001184 0.001256 0.01% MLCellLinOp::correctionResidual() 390 0.0009439 0.0009928 0.001082 0.01% MLMG::mgVcycle() 65 0.0007813 0.0008351 0.0009091 0.01% main() 1 0.0006347 0.0007574 0.000866 0.01% Castro::save_data_for_retry() 3 0.0007703 0.0007977 0.0008222 0.01% MLMG:computeResOfCorrection() 390 0.000705 0.0007416 0.0007927 0.01% Castro::subcycle_advance_ctu() 1 8.602e-05 0.0001397 0.0006783 0.01% AmrMesh::MakeNewGrids() 2 0.000526 0.0005494 0.0005649 0.01% MLMG::actualBottomSolve() 65 6.799e-05 0.0001018 0.00056 0.01% StateDataPhysBCFunct::() 82 0.0003639 0.000438 0.0004976 0.00% Castro::Castro() 1 0.0003625 0.0004023 0.0004629 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003695 0.000391 0.0004184 0.00% Castro::clean_state() 62 0.000344 0.0003613 0.0004182 0.00% MLMG::mgVcycle_down::2 65 0.0003119 0.0003477 0.0004107 0.00% MLMG::mgVcycle_down::3 65 0.0002912 0.0003235 0.0003925 0.00% MLMG::mgVcycle_down::1 65 0.0003311 0.0003503 0.000382 0.00% MLMG::mgVcycle_down::0 65 0.0003366 0.0003515 0.000381 0.00% MLMG::mgVcycle_down::4 65 0.0002811 0.0003004 0.0003514 0.00% MLMG::mgVcycle_down::5 65 0.000294 0.0003098 0.0003493 0.00% Gravity::solve_for_phi() 8 0.0002818 0.0002969 0.0003351 0.00% FillPatchSingleLevel 82 0.0002525 0.0002661 0.0002803 0.00% MLMG::mgVcycle_up::1 65 0.0002258 0.0002448 0.0002742 0.00% Castro::do_advance_ctu() 11 0.0002234 0.0002392 0.0002735 0.00% MLMG::mgVcycle_up::0 65 0.000228 0.0002425 0.0002722 0.00% MLCellLinOp::solutionResidual() 75 0.0002239 0.0002381 0.0002634 0.00% MLMG::mgVcycle_up::2 65 0.0001922 0.0002177 0.0002561 0.00% Castro::problem_post_timestep() 2 1.158e-05 2.771e-05 0.0002406 0.00% MLMG::mgVcycle_up::3 65 0.000191 0.0002016 0.0002399 0.00% MLMG::mgVcycle_up::5 65 0.0002009 0.0002103 0.0002292 0.00% MLMG::mgVcycle_up::4 65 0.0001938 0.0002057 0.0002218 0.00% MLMG::oneIter() 65 0.0001902 0.0002027 0.0002215 0.00% Castro::do_old_sources() 11 0.0001634 0.0001725 0.0001816 0.00% Castro::initialize_advance() 1 0.0001308 0.0001451 0.0001596 0.00% Amr::coarseTimeStep() 1 0.0001309 0.0001523 0.0001591 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001179 0.0001274 0.0001469 0.00% Castro::construct_new_gravity() 8 8.011e-05 9.081e-05 0.0001408 0.00% MLMG::computeResidual() 65 0.0001213 0.0001284 0.0001394 0.00% Castro::do_new_sources() 8 0.0001004 0.0001078 0.0001307 0.00% Amr::defBaseLevel() 1 0.0001072 0.000119 0.0001305 0.00% MLMG::mgVcycle_bottom 65 0.0001111 0.000117 0.0001257 0.00% Castro::retry_advance_ctu() 11 5.599e-05 6.789e-05 0.0001233 0.00% TagBoxArray::mapPRD 2 0.0001092 0.0001142 0.0001191 0.00% Gravity::test_composite_phi() 1 9.372e-05 9.784e-05 0.0001118 0.00% Castro::check_for_nan() 19 5.425e-05 6.121e-05 9.055e-05 0.00% MLPoisson::define() 10 6.885e-05 7.454e-05 8.316e-05 0.00% Castro::construct_old_gravity() 11 5.862e-05 6.784e-05 8.167e-05 0.00% Amr::FinalizeInit() 1 8.091e-06 1.282e-05 6.787e-05 0.00% Castro::apply_source_to_state() 19 5.465e-05 5.912e-05 6.431e-05 0.00% Amr::grid_places() 2 5.468e-05 5.78e-05 6.407e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.182e-05 5.44e-05 5.895e-05 0.00% Castro::react_state() 19 4.912e-05 5.253e-05 5.678e-05 0.00% Castro::expand_state() 11 3.267e-05 3.852e-05 5.361e-05 0.00% Castro::create_source_corrector() 8 3.49e-05 3.851e-05 4.531e-05 0.00% Castro::derive() 18 4.176e-05 4.289e-05 4.468e-05 0.00% Amr::regrid() 1 1.574e-05 1.96e-05 4.038e-05 0.00% Gravity::actual_multilevel_solve() 1 2.835e-05 2.939e-05 3.097e-05 0.00% Amr::timeStep() 1 1.525e-05 1.983e-05 2.699e-05 0.00% MLMG::computeMLResidual() 9 2.13e-05 2.267e-05 2.633e-05 0.00% Castro::post_init() 1 1.178e-05 1.298e-05 1.659e-05 0.00% Amr::InitializeInit() 1 7.484e-06 9.388e-06 1.632e-05 0.00% Castro::advance() 1 7.158e-06 9.38e-06 1.373e-05 0.00% Castro::post_regrid() 2 1.049e-05 1.112e-05 1.177e-05 0.00% Amr::init() 1 6.73e-06 7.634e-06 9.314e-06 0.00% Amr::bldFineLevels() 1 6.781e-06 7.458e-06 8.629e-06 0.00% Castro::post_timestep() 1 5.691e-06 6.269e-06 6.754e-06 0.00% Amr::initialInit() 1 3.097e-06 3.458e-06 4.49e-06 0.00% Other 4669 0.03336 0.04381 0.06416 0.62% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.34 10.34 10.34 100.00% Amr::coarseTimeStep() 1 7.704 7.704 7.705 74.53% Amr::timeStep() 1 7.648 7.648 7.649 73.99% Castro::advance() 1 7.582 7.585 7.588 73.41% Castro::subcycle_advance_ctu() 1 7.532 7.535 7.537 72.92% Castro::do_advance_ctu() 11 7.507 7.51 7.513 72.68% Castro::construct_ctu_hydro_source() 11 3.385 3.402 3.433 33.21% Castro::clean_state() 62 2.968 3.093 3.376 32.66% Castro::computeTemp() 63 1.923 2.013 2.213 21.41% Amr::InitAmr() 1 1.708 1.708 1.708 16.52% Castro::reset_internal_energy(MultiFab) 63 1.206 1.246 1.34 12.97% Castro::initialize_do_advance() 11 0.8675 1.016 1.059 10.24% FabArray::ParallelCopy() 734 0.4713 0.7483 0.9693 9.38% Gravity::solve_phi_with_mlmg() 10 0.8841 0.9128 0.9297 8.99% Amr::init() 1 0.9235 0.9235 0.9235 8.93% FabArray::ParallelCopy_finish() 734 0.3668 0.6552 0.8806 8.52% Amr::initialInit() 1 0.866 0.8661 0.8661 8.38% Castro::construct_new_gravity() 8 0.7858 0.8144 0.8311 8.04% Gravity::solve_for_phi() 8 0.7669 0.7957 0.8126 7.86% Castro::evaluate_source_change() 144 0.7151 0.7362 0.804 7.78% Gravity::actual_solve_with_mlmg() 10 0.7697 0.7729 0.7751 7.50% MLMG::solve() 9 0.7418 0.7441 0.7449 7.21% Castro::do_old_sources() 11 0.6735 0.7129 0.7415 7.17% Castro::enforce_min_density() 62 0.6249 0.6417 0.6836 6.61% AmrLevel::FillPatch() 82 0.2863 0.563 0.6666 6.45% MLMG::oneIter() 65 0.6545 0.6581 0.664 6.42% MLMG::mgVcycle() 65 0.6467 0.6503 0.6549 6.34% Amr::InitializeInit() 1 0.4385 0.4733 0.636 6.15% Amr::defBaseLevel() 1 0.4385 0.4733 0.636 6.15% amrex::Copy() 523 0.5818 0.5982 0.6352 6.15% FillPatchIterator::Initialize 82 0.2036 0.4855 0.5912 5.72% FillPatchIterator::FillFromLevel0() 82 0.1947 0.4767 0.5817 5.63% FillPatchSingleLevel 82 0.1943 0.4763 0.5813 5.62% Castro::do_new_sources() 8 0.446 0.4743 0.5542 5.36% Castro::normalize_species() 62 0.4474 0.4659 0.5075 4.91% Castro::reset_internal_energy(Fab) 126 0.4287 0.4433 0.4713 4.56% MultiFab::AddProduct() 2633 0.4219 0.4337 0.4703 4.55% MLCellLinOp::smooth() 1560 0.3459 0.377 0.4481 4.34% Amr::FinalizeInit() 1 0.23 0.3928 0.4276 4.14% Castro::Castro() 1 0.3601 0.3684 0.3905 3.78% Castro::initMFs() 1 0.3563 0.3646 0.3867 3.74% Castro::expand_state() 11 0.117 0.2586 0.3315 3.21% MultiFab::Subtract() 125 0.2752 0.283 0.2998 2.90% check_for_negative_density() 11 0.1217 0.2283 0.2919 2.82% MLCellLinOp::applyBC() 3630 0.2107 0.2343 0.2785 2.69% MLMG::mgVcycle_up::1 65 0.05795 0.1517 0.2452 2.37% Castro::initData() 1 0.07681 0.1034 0.2442 2.36% FabArray::setVal() 3453 0.2159 0.2228 0.2438 2.36% Amr::bldFineLevels() 1 0.04726 0.2095 0.2434 2.35% FabArray::FillBoundary() 3240 0.1747 0.2 0.2326 2.25% MLPoisson::Fsmooth() 3120 0.1631 0.1821 0.2175 2.10% MLMG::addInterpCorrection() 390 0.03356 0.167 0.2105 2.04% Amr::grid_places() 2 0.002711 0.1654 0.2002 1.94% AmrMesh::MakeNewGrids() 2 0.002656 0.1654 0.2002 1.94% TagBoxArray::mapPRD 2 0.00113 0.1244 0.1986 1.92% Castro::create_source_corrector() 8 0.05416 0.1188 0.181 1.75% TagBoxArray::collate() 2 0.0001833 0.03974 0.1677 1.62% MLMG::mgVcycle_down::0 65 0.1537 0.1571 0.1652 1.60% FabArray::sum() 2633 0.1384 0.1441 0.1589 1.54% Gravity::fill_multipole_BCs() 10 0.1111 0.1392 0.1562 1.51% FillBoundary_finish() 3240 0.08857 0.1286 0.1503 1.45% MLMG::mgVcycle_up::2 65 0.005604 0.06273 0.1501 1.45% Castro::post_init() 1 0.1379 0.1395 0.1413 1.37% MLMG::mgVcycle_up::0 65 0.1196 0.1219 0.1253 1.21% MLCellLinOp::apply() 501 0.09197 0.0979 0.1154 1.12% Gravity::multilevel_solve_for_new_phi() 1 0.105 0.105 0.105 1.02% Gravity::actual_multilevel_solve() 1 0.1049 0.1049 0.105 1.02% FabArray::ParallelCopy_nowait() 734 0.08133 0.09115 0.1025 0.99% Amr::writePlotFile() 2 0.1001 0.1002 0.1002 0.97% MLMG::mgVcycle_up::3 65 0.005549 0.01233 0.09195 0.89% Castro::post_regrid() 2 0.08594 0.08779 0.0894 0.86% MLMG::mgVcycle_down::1 65 0.07645 0.08005 0.08626 0.83% MLMG:computeResOfCorrection() 390 0.06457 0.07126 0.08329 0.81% MLCellLinOp::correctionResidual() 390 0.06386 0.07052 0.08251 0.80% FillBoundary_nowait() 3240 0.04794 0.06392 0.08058 0.78% Castro::construct_old_gravity() 11 0.02985 0.0502 0.07642 0.74% Gravity::get_old_grav_vector() 11 0.02979 0.05013 0.07635 0.74% Castro::derive() 18 0.07106 0.07239 0.07367 0.71% AmrLevel::derive() 18 0.07102 0.07234 0.07363 0.71% Castro::construct_new_source() 48 0.05793 0.06093 0.06773 0.66% Castro::post_timestep() 1 0.06017 0.06367 0.06656 0.64% MLPoisson::Fapply() 501 0.04835 0.05222 0.06091 0.59% MLCellLinOp::solutionResidual() 75 0.05126 0.05508 0.05752 0.56% Castro::initialize_advance() 1 0.04895 0.04982 0.0508 0.49% MLMG::computeResidual() 65 0.04485 0.04569 0.04693 0.45% FabArray::Saxpy() 133 0.03554 0.03795 0.04693 0.45% MLMG::mgVcycle_down::2 65 0.006602 0.02478 0.04685 0.45% Amr::regrid() 1 0.04411 0.04519 0.04673 0.45% Castro::construct_old_source() 66 0.04241 0.04342 0.04623 0.45% Castro::construct_new_gravity_source() 8 0.03216 0.0337 0.03844 0.37% MLMG::mgVcycle_bottom 65 0.0001793 0.00255 0.03798 0.37% MLMG::actualBottomSolve() 65 6.799e-05 0.002432 0.03785 0.37% MLCGSolver::bicgstab 4 0 0.002315 0.03703 0.36% MLMG::mgVcycle_down::3 65 0.006258 0.009776 0.03253 0.31% FabArray::Xpay() 480 0.02529 0.02732 0.03222 0.31% amrex::average_down 390 0.01679 0.0217 0.03144 0.30% Castro::retry_advance_ctu() 11 0.02295 0.02358 0.02622 0.25% Castro::construct_old_gravity_source() 11 0.02445 0.02491 0.02549 0.25% Castro::check_for_nan() 19 0.02311 0.02389 0.02516 0.24% MultiFab::contains_nan() 19 0.02305 0.02383 0.02507 0.24% VisMF::Write(FabArray) 4 0.01429 0.01544 0.0245 0.24% MLPoisson::define() 10 0.01685 0.01951 0.02346 0.23% MLMG::ResNormInf() 74 0.01204 0.01921 0.02325 0.22% Castro::sum_integrated_quantities() 2 0.02148 0.02162 0.02163 0.21% Gravity::get_new_grav_vector() 9 0.01821 0.01945 0.01999 0.19% Castro::apply_source_to_state() 19 0.01474 0.01556 0.0198 0.19% MLMG::mgVcycle_down::4 65 0.006228 0.007217 0.01839 0.18% Castro::problem_post_timestep() 2 0.009668 0.01409 0.01786 0.17% Castro::wd_update() 2 0.009656 0.01407 0.01785 0.17% Castro::react_state() 19 0.01492 0.01531 0.01721 0.17% StateDataPhysBCFunct::() 82 0.0003639 0.007687 0.01454 0.14% MLMG::mgVcycle_up::4 65 0.005629 0.006375 0.01434 0.14% Gravity::test_composite_phi() 1 0.01419 0.01422 0.01425 0.14% StateData::FillBoundary(geom) 62 0 0.007249 0.01407 0.14% MLCellLinOp::defineAuxData() 10 0.01022 0.01107 0.01328 0.13% Castro::save_data_for_retry() 3 0.01175 0.012 0.01323 0.13% Amr::writeSmallPlotFile() 2 0.01311 0.01313 0.01314 0.13% MLMG::mgVcycle_down::5 65 0.006186 0.006793 0.01228 0.12% MLMG::prepareForSolve() 9 0.01013 0.01062 0.01183 0.11% MLMG::mgVcycle_up::5 65 0.005674 0.006179 0.01035 0.10% MLMG::computeMLResidual() 9 0.005477 0.008756 0.01013 0.10% FabArray::mult() 40 0.008042 0.008262 0.009039 0.09% amrex::Add() 65 0.007331 0.007617 0.00891 0.09% FabArray::norminf() 143 0.005657 0.00607 0.008155 0.08% FabArrayBase::getCPC() 1022 0.005233 0.006113 0.007922 0.08% Other 4669 0.0521 0.07146 0.1016 0.98% -------------------------------------------------------------------------------------------- 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 24 B 32 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.08-11-gb6ed5237987e) finalized