Initializing AMReX (24.08-9-g6f039ae18e89)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.08-9-g6f039ae18e89) initialized Starting run at 07:04:21 UTC on 2024-08-08. Successfully read inputs file ... Castro git describe: 24.08-1-g36c0aaf32 AMReX git describe: 24.08-9-g6f039ae18 Microphysics git describe: 24.08-2-g8ce3375a reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... Creating primary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Creating secondary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Generated initial model for primary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Generated initial model for secondary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Castro::numpts_1d at level 0 is 294 Setting the gravity type to PoissonGrav Initializing the data at level 0 Done initializing the level 0 data Now regridding at level lbase = 0 grid_places() time: 0.163099654 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 grid_places() time: 0.001161643 new finest: 0 Gravity::fill_multipole_BCs() time = 0.010779357 Gravity::multilevel_solve_for_new_phi() time = 0.105481971 ... norm of composite residual at level 0 1.193955541e-11 Gravity::fill_multipole_BCs() time = 0.0115031 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.04986498 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.006942375 seconds [Level 0 step 1] ADVANCE at time 0 with dt = 0.25  Beginning subcycle 1 starting at time 0 with dt = 0.25 Estimated number of subcycles remaining: 1 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Castro::construct_old_gravity() time = 0.002319734 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.064027553 on level 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Castro::construct_old_gravity() time = 0.005812125 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.063036482 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.005389497 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.061286069 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.006622049 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.062081074 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.302998926 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.015944562 Gravity::solve_for_phi() time = 0.095062429 on level 0 Castro::construct_new_gravity() time = 0.097660196 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.067429317 on level 0 Castro::construct_old_gravity() time = 0.01195337 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.077686927 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.381040716 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.016037731 Gravity::solve_for_phi() time = 0.088849558 on level 0 Castro::construct_new_gravity() time = 0.09098734 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.060979805 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.005384147 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.066920347 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.368667407 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.016821135 Gravity::solve_for_phi() time = 0.088740682 on level 0 Castro::construct_new_gravity() time = 0.090854245 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.061086225 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.002616162 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.060981264 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.318074929 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.021298629 Gravity::solve_for_phi() time = 0.101512084 on level 0 Castro::construct_new_gravity() time = 0.10371251 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.061798808 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.002344977 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.075012442 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.320170158 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.01575169 Gravity::solve_for_phi() time = 0.10442544 on level 0 Castro::construct_new_gravity() time = 0.106514753 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.061345808 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.003873365 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.087371668 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.34905871 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.01643734 Gravity::solve_for_phi() time = 0.111426384 on level 0 Castro::construct_new_gravity() time = 0.113722481 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.066708852 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.006237759 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.061780998 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.315532158 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.016830984 Gravity::solve_for_phi() time = 0.08976818 on level 0 Castro::construct_new_gravity() time = 0.091900287 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.061478267 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.003938402 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.061669273 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.313391369 on level 0  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.00449686833 [Level 0 step 1] Advanced 32768 cells Gravity::fill_multipole_BCs() time = 0.017325895 Gravity::solve_for_phi() time = 0.092655641 on level 0 Castro::construct_new_gravity() time = 0.094752538 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.061062587 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.348761759 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.049537084 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00584024 seconds Ending run at 07:04:31 UTC on 2024-08-08. Run time = 9.847169155 Run time without initialization = 7.405180383 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.847 ... 9.847 ... 9.847 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7294 2.7900 2.8465 28.91% Amr::InitAmr() 1 1.548 1.548 1.548 15.72% Castro::computeTemp() 63 0.7125 0.7633 0.8532 8.66% amrex::Copy() 523 0.5844 0.5936 0.6003 6.10% FabArray::ParallelCopy_finish() 734 0.249 0.4028 0.5905 6.00% Castro::normalize_species() 62 0.45 0.4592 0.4703 4.78% Castro::reset_internal_energy(Fab) 126 0.4305 0.4392 0.4438 4.51% MultiFab::AddProduct() 2633 0.4225 0.4315 0.4365 4.43% Castro::initMFs() 1 0.3602 0.3659 0.3926 3.99% MultiFab::Subtract() 125 0.2765 0.2816 0.2857 2.90% check_for_negative_density() 11 0.155 0.2378 0.284 2.88% FabArray::setVal() 3453 0.2161 0.2211 0.2247 2.28% MLPoisson::Fsmooth() 3120 0.1632 0.1794 0.2146 2.18% Castro::initData() 1 0.003288 0.02799 0.1632 1.66% Gravity::fill_multipole_BCs() 10 0.1145 0.1394 0.1545 1.57% FabArray::sum() 2633 0.1403 0.1429 0.1452 1.47% TagBoxArray::collate() 2 0.000217 0.02403 0.1288 1.31% FillBoundary_finish() 3240 0.08663 0.1001 0.1235 1.25% FabArray::ParallelCopy_nowait() 734 0.06621 0.07255 0.08101 0.82% FillBoundary_nowait() 3240 0.04347 0.05738 0.07278 0.74% MLPoisson::Fapply() 501 0.04854 0.05174 0.05886 0.60% MLCellLinOp::applyBC() 3630 0.02604 0.03393 0.05546 0.56% AmrLevel::derive() 18 0.04107 0.04182 0.04249 0.43% FabArray::Saxpy() 133 0.03575 0.03656 0.03981 0.40% Castro::construct_new_gravity_source() 8 0.03228 0.03284 0.03368 0.34% FabArray::Xpay() 480 0.02538 0.02695 0.03019 0.31% Castro::construct_old_gravity_source() 11 0.02458 0.02505 0.02532 0.26% VisMF::Write(FabArray) 4 0.01389 0.01516 0.02522 0.26% MultiFab::contains_nan() 19 0.02334 0.02376 0.02428 0.25% MLMG::addInterpCorrection() 390 0.01718 0.01873 0.02215 0.22% Castro::reset_internal_energy(MultiFab) 63 0.01828 0.01872 0.01931 0.20% MLMG::ResNormInf() 74 0.007164 0.01408 0.01903 0.19% Amr::writePlotFile() 2 0.007853 0.01592 0.0171 0.17% amrex::average_down 390 0.01257 0.0137 0.01674 0.17% Castro::wd_update() 2 0.004945 0.01154 0.01487 0.15% Castro::enforce_min_density() 62 0.01405 0.01434 0.01473 0.15% StateData::FillBoundary(geom) 62 0 0.0072 0.01402 0.14% Castro::evaluate_source_change() 144 0.009789 0.01014 0.01067 0.11% MLCellLinOp::defineAuxData() 10 0.007979 0.008332 0.009576 0.10% amrex::Add() 65 0.007351 0.0076 0.009099 0.09% Amr::writeSmallPlotFile() 2 0.006662 0.008498 0.008993 0.09% FabArray::FillBoundary() 3240 0.007067 0.007369 0.008831 0.09% MLCGSolver::bicgstab 4 0 0.0005237 0.00838 0.09% FabArray::mult() 40 0.008048 0.008177 0.008322 0.08% FabArray::norminf() 143 0.005646 0.006013 0.007992 0.08% Castro::construct_new_source() 48 0.007411 0.007487 0.00764 0.08% MLCellLinOp::smooth() 1560 0.005786 0.006044 0.006426 0.07% FillPatchIterator::Initialize 82 0.005764 0.005953 0.006225 0.06% Castro::construct_old_source() 66 0.005474 0.0056 0.005989 0.06% MLMG::prepareForSolve() 9 0.005014 0.005158 0.005445 0.06% Gravity::actual_solve_with_mlmg() 10 0.003574 0.003723 0.004163 0.04% FabArrayBase::getCPC() 1022 0.003899 0.004009 0.004155 0.04% Gravity::get_old_grav_vector() 11 0.003555 0.00367 0.003927 0.04% Castro::initialize_do_advance() 11 0.002697 0.002981 0.003353 0.03% Gravity::get_new_grav_vector() 9 0.002981 0.003031 0.003285 0.03% MLMG::solve() 9 0.0004488 0.001722 0.00322 0.03% MLCellLinOp::apply() 501 0.001031 0.001145 0.002378 0.02% Castro::sum_integrated_quantities() 2 0.00154 0.001787 0.002088 0.02% FabArray::ParallelCopy() 734 0.00181 0.001888 0.00204 0.02% AmrLevel::FillPatch() 82 0.001145 0.001197 0.001317 0.01% MLCellLinOp::correctionResidual() 390 0.000954 0.0009855 0.001028 0.01% main() 1 0.0006512 0.0007915 0.0009138 0.01% MLMG::mgVcycle() 65 0.0007888 0.0008248 0.0008901 0.01% Castro::save_data_for_retry() 3 0.0007586 0.0008092 0.0008686 0.01% MLMG:computeResOfCorrection() 390 0.0007156 0.0007367 0.0007784 0.01% Castro::subcycle_advance_ctu() 1 8.751e-05 0.000148 0.000683 0.01% AmrMesh::MakeNewGrids() 2 0.0005309 0.0005484 0.0005902 0.01% MLMG::actualBottomSolve() 65 6.779e-05 0.0001009 0.000569 0.01% StateDataPhysBCFunct::() 82 0.0003737 0.0004389 0.0004852 0.00% Castro::Castro() 1 0.0003612 0.00041 0.0004807 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003722 0.0003877 0.0004092 0.00% Castro::clean_state() 62 0.0003376 0.000366 0.0003954 0.00% MLMG::mgVcycle_down::1 65 0.0003375 0.0003456 0.0003726 0.00% MLMG::mgVcycle_down::0 65 0.0003367 0.0003451 0.000372 0.00% MLMG::mgVcycle_down::2 65 0.0003163 0.0003432 0.0003677 0.00% MLMG::mgVcycle_down::3 65 0.0002956 0.0003213 0.0003626 0.00% MLMG::mgVcycle_down::4 65 0.0002823 0.0002988 0.0003508 0.00% MLMG::mgVcycle_down::5 65 0.0002965 0.0003056 0.0003424 0.00% Gravity::solve_for_phi() 8 0.0002748 0.0002998 0.0003374 0.00% FillPatchSingleLevel 82 0.000258 0.0002667 0.0002897 0.00% Castro::do_advance_ctu() 11 0.0002211 0.0002444 0.0002833 0.00% MLMG::mgVcycle_up::0 65 0.0002322 0.0002408 0.0002794 0.00% Castro::problem_post_timestep() 2 1.074e-05 2.916e-05 0.000274 0.00% MLMG::mgVcycle_up::1 65 0.0002268 0.000242 0.0002622 0.00% MLCellLinOp::solutionResidual() 75 0.0002291 0.0002364 0.0002508 0.00% MLMG::mgVcycle_up::2 65 0.0001925 0.0002166 0.0002334 0.00% MLMG::mgVcycle_up::3 65 0.0001919 0.0002007 0.0002279 0.00% MLMG::mgVcycle_up::5 65 0.0002008 0.0002071 0.0002255 0.00% MLMG::mgVcycle_up::4 65 0.0001948 0.0002015 0.0002205 0.00% MLMG::oneIter() 65 0.000194 0.0002047 0.0002199 0.00% Castro::do_old_sources() 11 0.0001649 0.0001753 0.0002019 0.00% Amr::coarseTimeStep() 1 0.000154 0.0001786 0.0001945 0.00% Castro::initialize_advance() 1 0.0001216 0.0001467 0.0001837 0.00% Castro::construct_new_gravity() 8 7.674e-05 9.151e-05 0.0001712 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001239 0.0001322 0.0001529 0.00% MLMG::computeResidual() 65 0.00012 0.0001272 0.0001356 0.00% Castro::do_new_sources() 8 0.0001008 0.0001105 0.0001319 0.00% Castro::retry_advance_ctu() 11 5.335e-05 6.389e-05 0.0001295 0.00% TagBoxArray::mapPRD 2 0.0001113 0.0001164 0.0001277 0.00% MLMG::mgVcycle_bottom 65 0.0001135 0.0001161 0.000124 0.00% Amr::defBaseLevel() 1 9.892e-05 0.0001067 0.0001176 0.00% Gravity::test_composite_phi() 1 9.32e-05 9.836e-05 0.000109 0.00% Castro::construct_old_gravity() 11 6.479e-05 7.154e-05 0.0001044 0.00% Amr::grid_places() 2 5.414e-05 6.097e-05 8.464e-05 0.00% MLPoisson::define() 10 6.752e-05 7.357e-05 7.817e-05 0.00% Castro::check_for_nan() 19 5.39e-05 6.236e-05 7.416e-05 0.00% Castro::apply_source_to_state() 19 5.385e-05 5.907e-05 6.749e-05 0.00% Castro::create_source_corrector() 8 3.589e-05 4.119e-05 6.719e-05 0.00% Castro::react_state() 19 4.983e-05 5.396e-05 5.977e-05 0.00% Amr::FinalizeInit() 1 8.24e-06 1.324e-05 5.905e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.291e-05 5.428e-05 5.739e-05 0.00% Castro::derive() 18 4.169e-05 4.383e-05 5.451e-05 0.00% Amr::timeStep() 1 1.506e-05 2.538e-05 4.277e-05 0.00% Castro::expand_state() 11 3.149e-05 3.497e-05 4.037e-05 0.00% Amr::regrid() 1 1.644e-05 2.006e-05 3.947e-05 0.00% Gravity::actual_multilevel_solve() 1 2.835e-05 2.93e-05 3.173e-05 0.00% MLMG::computeMLResidual() 9 2.088e-05 2.257e-05 2.538e-05 0.00% Castro::post_regrid() 2 1.061e-05 1.186e-05 1.867e-05 0.00% Amr::InitializeInit() 1 7.649e-06 1.113e-05 1.627e-05 0.00% Castro::post_init() 1 1.23e-05 1.293e-05 1.528e-05 0.00% Castro::advance() 1 6.503e-06 8.741e-06 1.39e-05 0.00% Amr::init() 1 6.533e-06 8.413e-06 1.237e-05 0.00% Amr::bldFineLevels() 1 6.759e-06 7.587e-06 9.052e-06 0.00% Amr::initialInit() 1 3.021e-06 3.715e-06 7.728e-06 0.00% Castro::post_timestep() 1 5.58e-06 6.098e-06 6.955e-06 0.00% Other 4669 0.03362 0.04257 0.06298 0.64% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.847 9.847 9.847 100.00% Amr::coarseTimeStep() 1 7.404 7.405 7.405 75.19% Amr::timeStep() 1 7.349 7.349 7.349 74.63% Castro::advance() 1 7.281 7.284 7.289 74.02% Castro::subcycle_advance_ctu() 1 7.23 7.233 7.239 73.51% Castro::do_advance_ctu() 11 7.204 7.208 7.214 73.26% Castro::construct_ctu_hydro_source() 11 3.366 3.402 3.454 35.08% Castro::clean_state() 62 2.985 3.068 3.165 32.14% Castro::computeTemp() 63 1.935 1.998 2.093 21.25% Amr::InitAmr() 1 1.548 1.548 1.548 15.72% Castro::reset_internal_energy(MultiFab) 63 1.211 1.235 1.249 12.69% Amr::init() 1 0.8941 0.8942 0.8942 9.08% Gravity::solve_phi_with_mlmg() 10 0.8448 0.8718 0.889 9.03% Amr::initialInit() 1 0.8368 0.8369 0.8371 8.50% Castro::initialize_do_advance() 11 0.7916 0.8107 0.8342 8.47% Castro::construct_new_gravity() 8 0.7441 0.7709 0.7879 8.00% Gravity::solve_for_phi() 8 0.7263 0.7533 0.7705 7.82% Castro::evaluate_source_change() 144 0.7168 0.7314 0.7402 7.52% Castro::do_old_sources() 11 0.6414 0.7013 0.7346 7.46% Gravity::actual_solve_with_mlmg() 10 0.7284 0.7306 0.7326 7.44% MLMG::solve() 9 0.7004 0.7029 0.7036 7.15% FabArray::ParallelCopy() 734 0.3374 0.4948 0.6813 6.92% Castro::enforce_min_density() 62 0.6272 0.6384 0.6487 6.59% MLMG::oneIter() 65 0.6147 0.6179 0.622 6.32% MLMG::mgVcycle() 65 0.607 0.6101 0.6128 6.22% Amr::InitializeInit() 1 0.4446 0.4744 0.6061 6.16% Amr::defBaseLevel() 1 0.4446 0.4744 0.6061 6.16% amrex::Copy() 523 0.5844 0.5936 0.6003 6.10% FabArray::ParallelCopy_finish() 734 0.249 0.4028 0.5905 6.00% Castro::do_new_sources() 8 0.4428 0.4633 0.4878 4.95% Castro::normalize_species() 62 0.45 0.4592 0.4703 4.78% Castro::reset_internal_energy(Fab) 126 0.4305 0.4392 0.4438 4.51% MultiFab::AddProduct() 2633 0.4225 0.4315 0.4365 4.43% AmrLevel::FillPatch() 82 0.288 0.3324 0.4256 4.32% MLCellLinOp::smooth() 1560 0.3212 0.3505 0.4188 4.25% Castro::Castro() 1 0.364 0.3699 0.3968 4.03% Castro::initMFs() 1 0.3603 0.366 0.3927 3.99% Amr::FinalizeInit() 1 0.2308 0.3625 0.3923 3.98% FillPatchIterator::Initialize 82 0.2105 0.2554 0.3496 3.55% FillPatchIterator::FillFromLevel0() 82 0.2008 0.2465 0.3413 3.47% FillPatchSingleLevel 82 0.2005 0.2461 0.3409 3.46% MultiFab::Subtract() 125 0.2765 0.2816 0.2857 2.90% check_for_negative_density() 11 0.155 0.2378 0.284 2.88% MLCellLinOp::applyBC() 3630 0.1856 0.2047 0.2498 2.54% Castro::initData() 1 0.07688 0.1031 0.2367 2.40% MLMG::mgVcycle_up::1 65 0.05672 0.1444 0.2321 2.36% FabArray::setVal() 3453 0.2161 0.2211 0.2247 2.28% MLPoisson::Fsmooth() 3120 0.1632 0.1794 0.2146 2.18% Amr::bldFineLevels() 1 0.04607 0.1784 0.2088 2.12% MLMG::addInterpCorrection() 390 0.03383 0.1596 0.1981 2.01% FabArray::FillBoundary() 3240 0.156 0.1707 0.1944 1.97% Castro::expand_state() 11 0.1213 0.1377 0.1702 1.73% Amr::grid_places() 2 0.002754 0.1345 0.1643 1.67% AmrMesh::MakeNewGrids() 2 0.002699 0.1344 0.1642 1.67% TagBoxArray::mapPRD 2 0.001157 0.1091 0.1621 1.65% Gravity::fill_multipole_BCs() 10 0.1157 0.1405 0.1556 1.58% MLMG::mgVcycle_down::0 65 0.145 0.1467 0.1493 1.52% FabArray::sum() 2633 0.1403 0.1429 0.1452 1.47% MLMG::mgVcycle_up::2 65 0.005639 0.06073 0.1445 1.47% Castro::post_init() 1 0.139 0.1403 0.142 1.44% TagBoxArray::collate() 2 0.000217 0.02403 0.1288 1.31% FillBoundary_finish() 3240 0.08663 0.1001 0.1235 1.25% MLMG::mgVcycle_up::0 65 0.1124 0.1146 0.1164 1.18% MLCellLinOp::apply() 501 0.08191 0.09147 0.1123 1.14% Gravity::multilevel_solve_for_new_phi() 1 0.1054 0.1055 0.1055 1.07% Gravity::actual_multilevel_solve() 1 0.1054 0.1054 0.1054 1.07% Amr::writePlotFile() 2 0.09972 0.09981 0.09982 1.01% FabArray::ParallelCopy_nowait() 734 0.0834 0.09014 0.099 1.01% MLMG::mgVcycle_up::3 65 0.005602 0.01231 0.09221 0.94% Castro::post_regrid() 2 0.08604 0.08772 0.08903 0.90% FillBoundary_nowait() 3240 0.04893 0.06328 0.07886 0.80% MLMG:computeResOfCorrection() 390 0.06039 0.06611 0.07827 0.79% MLCellLinOp::correctionResidual() 390 0.05966 0.06537 0.07749 0.79% Castro::derive() 18 0.07141 0.0725 0.07352 0.75% AmrLevel::derive() 18 0.07136 0.07245 0.07347 0.75% MLMG::mgVcycle_down::1 65 0.06625 0.06941 0.07209 0.73% Castro::post_timestep() 1 0.06008 0.06506 0.06765 0.69% Castro::construct_new_source() 48 0.05812 0.05899 0.06046 0.61% MLPoisson::Fapply() 501 0.04854 0.05174 0.05886 0.60% MLCellLinOp::solutionResidual() 75 0.04834 0.05343 0.0566 0.57% Castro::create_source_corrector() 8 0.03101 0.04164 0.05258 0.53% Castro::initialize_advance() 1 0.04908 0.0502 0.0519 0.53% Castro::construct_old_gravity() 11 0.02675 0.03364 0.04714 0.48% Gravity::get_old_grav_vector() 11 0.02664 0.03357 0.04708 0.48% MLMG::computeResidual() 65 0.04143 0.04435 0.04678 0.48% Amr::regrid() 1 0.04421 0.04505 0.04578 0.46% Castro::construct_old_source() 66 0.04276 0.04351 0.04429 0.45% MLMG::mgVcycle_down::2 65 0.006682 0.02273 0.04041 0.41% FabArray::Saxpy() 133 0.03575 0.03656 0.03981 0.40% MLMG::mgVcycle_bottom 65 0.000182 0.002506 0.03732 0.38% MLMG::actualBottomSolve() 65 6.779e-05 0.00239 0.0372 0.38% MLCGSolver::bicgstab 4 0 0.002273 0.03637 0.37% Castro::construct_new_gravity_source() 8 0.03228 0.03284 0.03368 0.34% MLMG::mgVcycle_down::3 65 0.006298 0.00958 0.03158 0.32% FabArray::Xpay() 480 0.02538 0.02695 0.03019 0.31% amrex::average_down 390 0.01696 0.02075 0.02921 0.30% Castro::construct_old_gravity_source() 11 0.02458 0.02505 0.02532 0.26% VisMF::Write(FabArray) 4 0.01389 0.01516 0.02522 0.26% Castro::retry_advance_ctu() 11 0.02301 0.02341 0.02448 0.25% Castro::check_for_nan() 19 0.02341 0.02382 0.02435 0.25% MultiFab::contains_nan() 19 0.02334 0.02376 0.02428 0.25% MLMG::ResNormInf() 74 0.01275 0.01914 0.02397 0.24% Castro::sum_integrated_quantities() 2 0.02138 0.02155 0.02157 0.22% MLPoisson::define() 10 0.01572 0.01848 0.02123 0.22% Gravity::get_new_grav_vector() 9 0.01797 0.01835 0.01866 0.19% Castro::problem_post_timestep() 2 0.009502 0.01527 0.0183 0.19% Castro::wd_update() 2 0.009491 0.01524 0.01829 0.19% MLMG::mgVcycle_down::4 65 0.006234 0.007144 0.01812 0.18% Castro::react_state() 19 0.01496 0.01537 0.01675 0.17% Castro::apply_source_to_state() 19 0.0148 0.01506 0.01533 0.16% Gravity::test_composite_phi() 1 0.01506 0.01509 0.01511 0.15% StateDataPhysBCFunct::() 82 0.0003737 0.007639 0.01449 0.15% MLMG::mgVcycle_up::4 65 0.005641 0.006304 0.01417 0.14% StateData::FillBoundary(geom) 62 0 0.0072 0.01402 0.14% MLCellLinOp::defineAuxData() 10 0.01031 0.01108 0.01346 0.14% Amr::writeSmallPlotFile() 2 0.01302 0.01302 0.01303 0.13% Castro::save_data_for_retry() 3 0.01182 0.01202 0.01224 0.12% MLMG::mgVcycle_down::5 65 0.006249 0.006739 0.01216 0.12% MLMG::prepareForSolve() 9 0.01016 0.01064 0.01178 0.12% MLMG::computeMLResidual() 9 0.006129 0.008472 0.01093 0.11% MLMG::mgVcycle_up::5 65 0.005685 0.006111 0.01024 0.10% amrex::Add() 65 0.007351 0.0076 0.009099 0.09% FabArray::mult() 40 0.008048 0.008177 0.008322 0.08% FabArrayBase::getCPC() 1022 0.005288 0.006097 0.008038 0.08% FabArray::norminf() 143 0.005646 0.006013 0.007992 0.08% Other 4669 0.05257 0.06771 0.09704 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 19 B 27 B 33 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.08-9-g6f039ae18e89) finalized