Initializing AMReX (24.07-21-gcdacc4d908e0)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.07-21-gcdacc4d908e0) initialized Starting run at 07:33:45 UTC on 2024-07-19. Successfully read inputs file ... Castro git describe: 24.07-16-g39a7704b4 AMReX git describe: 24.07-21-gcdacc4d90 Microphysics git describe: 24.07-25-g83817c08 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 grid_places() time: 0.196586264 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001072682 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010892186 Gravity::multilevel_solve_for_new_phi() time = 0.108084502 Gravity::fill_multipole_BCs() time = 0.011544736 ... norm of composite residual at level 0 1.193955541e-11 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.049302418 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003746675 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 Castro::construct_old_gravity() time = 0.014701471 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.078468036 on level 0 ... 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 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 Castro::construct_old_gravity() time = 0.009957153 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.070309233 on level 0 ... 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 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 Castro::construct_old_gravity() time = 0.007295144 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.063257146 on level 0 ... 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 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 Castro::construct_old_gravity() time = 0.006837249 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.063674834 on level 0 ... Entering construct_ctu_hydro_source() on level 0 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.316231233 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01595332 Gravity::solve_for_phi() time = 0.089406965 on level 0 Castro::construct_new_gravity() time = 0.091500766 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.061473707 on level 0  Subcycle completed  Beginning subcycle 2 starting at time 0.03125 with dt = 0.03125 Estimated number of subcycles remaining: 7 Castro::construct_old_gravity() time = 0.00802989 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.06216464 on level 0 ... Entering construct_ctu_hydro_source() 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.985643901e+46 ener added: 2.77668231e+42 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.322423455 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016440209 Gravity::solve_for_phi() time = 0.08711723 on level 0 Castro::construct_new_gravity() time = 0.089330921 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.059666162 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  Subcycle completed  Beginning subcycle 3 starting at time 0.0625 with dt = 0.03125 Estimated number of subcycles remaining: 6 Castro::construct_old_gravity() time = 0.002383713 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.059749638 on level 0 ... Entering construct_ctu_hydro_source() 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.987377036e+46 ener added: 2.185609542e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.317593221 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014498855 Gravity::solve_for_phi() time = 0.086052352 on level 0 Castro::construct_new_gravity() time = 0.088176374 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.059742189 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  Subcycle completed  Beginning subcycle 4 starting at time 0.09375 with dt = 0.03125 Estimated number of subcycles remaining: 5 Castro::construct_old_gravity() time = 0.002080478 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.058886871 on level 0 ... Entering construct_ctu_hydro_source() 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.985020288e+46 ener added: 2.299661927e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.316702659 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01486585 Gravity::solve_for_phi() time = 0.086350236 on level 0 Castro::construct_new_gravity() time = 0.088467536 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.062702946 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  Subcycle completed  Beginning subcycle 5 starting at time 0.125 with dt = 0.03125 Estimated number of subcycles remaining: 4 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.004986335 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.062595002 on level 0 ... Entering construct_ctu_hydro_source() 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.99240328e+46 ener added: 2.331782913e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.316669224 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015292899 Gravity::solve_for_phi() time = 0.085623179 on level 0 Castro::construct_new_gravity() time = 0.088100819 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.060187255 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  Subcycle completed  Beginning subcycle 6 starting at time 0.15625 with dt = 0.03125 Estimated number of subcycles remaining: 3 Castro::construct_old_gravity() time = 0.002083944 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.059652072 on level 0 ... Entering construct_ctu_hydro_source() 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.000655828e+46 ener added: 2.414422253e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.315441919 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014888947 Gravity::solve_for_phi() time = 0.085700005 on level 0 Castro::construct_new_gravity() time = 0.087842475 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.059756836 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  Subcycle completed  Beginning subcycle 7 starting at time 0.1875 with dt = 0.03125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.003316671 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.063843576 on level 0 ... Entering construct_ctu_hydro_source() 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.993941122e+46 ener added: 2.684057583e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.318730992 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015166753 Gravity::solve_for_phi() time = 0.087098224 on level 0 Castro::construct_new_gravity() time = 0.089221968 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.059435477 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  Subcycle completed  Beginning subcycle 8 starting at time 0.21875 with dt = 0.03125 Estimated number of subcycles remaining: 1 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.002119434 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.058502711 on level 0 ... Entering construct_ctu_hydro_source() 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.993654336e+46 ener added: 2.944609967e+43 ... Leaving construct_ctu_hydro_source() on level 0 Castro::construct_ctu_hydro_source() time = 0.315266796 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015748782 Gravity::solve_for_phi() time = 0.086236169 on level 0 Castro::construct_new_gravity() time = 0.08830969 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.058949182 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  Subcycle completed Subcycling complete Zones advanced per microsecond at this level: 0.004608824638 [Level 0 step 1] Advanced 32768 cells 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] Coarse TimeStep time: 7.168018029 [STEP 1] FAB kilobyte spread across MPI nodes: [6234 ... 6234] STEP = 1 TIME = 0.25 DT = 0.25 PLOTFILE: file = wdmerger_retry_plt00001 Write plotfile time = 0.047158851 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003772788 seconds Ending run at 07:33:54 UTC on 2024-07-19. Run time = 9.475878406 Run time without initialization = 7.220032794 Average number of zones advanced per microsecond: 0.005 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 9.476 ... 9.476 ... 9.476 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8365 2.8796 2.9199 30.81% Amr::InitAmr() 1 1.324 1.324 1.324 13.97% Castro::computeTemp() 63 0.702 0.756 0.8473 8.94% FabArray::ParallelCopy_finish() 734 0.04247 0.4039 0.588 6.21% amrex::Copy() 523 0.5464 0.5543 0.5688 6.00% Castro::normalize_species() 62 0.4605 0.4674 0.4973 5.25% Castro::reset_internal_energy(Fab) 126 0.4331 0.4385 0.4486 4.73% MultiFab::AddProduct() 2633 0.4194 0.4233 0.428 4.52% Castro::initMFs() 1 0.3596 0.368 0.3849 4.06% MultiFab::Subtract() 125 0.268 0.2709 0.2743 2.90% Castro::reset_internal_energy(MultiFab) 63 0.01794 0.06571 0.2434 2.57% FabArray::setVal() 3453 0.2229 0.2272 0.2313 2.44% MLPoisson::Fsmooth() 3120 0.1627 0.1771 0.2132 2.25% Castro::initData() 1 0.003337 0.0287 0.1715 1.81% Gravity::fill_multipole_BCs() 10 0.1059 0.1257 0.1424 1.50% FabArray::sum() 2633 0.139 0.1403 0.1418 1.50% check_for_negative_density() 11 0.0514 0.08237 0.1349 1.42% FabArray::ParallelCopy_nowait() 734 0.06531 0.0707 0.08104 0.86% FillBoundary_finish() 3240 0.05595 0.06445 0.0786 0.83% FillBoundary_nowait() 3240 0.04046 0.05306 0.06749 0.71% Castro::enforce_min_density() 62 0.01373 0.02463 0.05943 0.63% MLPoisson::Fapply() 501 0.04744 0.04959 0.05722 0.60% MLCellLinOp::applyBC() 3630 0.02445 0.03251 0.05449 0.57% AmrLevel::derive() 18 0.0411 0.0417 0.04385 0.46% Castro::do_new_sources() 8 0.000175 0.007491 0.04194 0.44% FabArray::Saxpy() 133 0.03409 0.03451 0.03821 0.40% Castro::do_old_sources() 11 0.0003061 0.008156 0.03713 0.39% Castro::construct_new_gravity_source() 8 0.03234 0.03287 0.0334 0.35% FabArray::Xpay() 480 0.02636 0.02756 0.03126 0.33% Castro::construct_old_gravity_source() 11 0.02451 0.02494 0.02601 0.27% VisMF::Write(FabArray) 4 0.01117 0.01257 0.02345 0.25% MultiFab::contains_nan() 19 0.02225 0.02257 0.0234 0.25% MLMG::addInterpCorrection() 390 0.01728 0.01853 0.02203 0.23% Amr::writePlotFile() 2 0.007192 0.01648 0.01797 0.19% amrex::average_down 390 0.01258 0.01356 0.0165 0.17% StateData::FillBoundary(geom) 62 0 0.007166 0.01395 0.15% Castro::wd_update() 2 0.004852 0.007716 0.01042 0.11% Castro::evaluate_source_change() 144 0.009416 0.009667 0.01016 0.11% MLCellLinOp::defineAuxData() 10 0.007775 0.008148 0.009267 0.10% FabArray::mult() 40 0.008279 0.008478 0.008975 0.09% FabArray::FillBoundary() 3240 0.006912 0.007173 0.008673 0.09% Castro::construct_old_gravity() 11 0.0002164 0.002132 0.008371 0.09% MLCGSolver::bicgstab 4 0 0.0005072 0.008116 0.09% FabArray::norminf() 143 0.005659 0.00593 0.007972 0.08% Castro::construct_new_source() 48 0.007306 0.007393 0.007835 0.08% Castro::volWgtSum() 46 0.007295 0.007407 0.007526 0.08% FabArrayBase::getFB() 3240 0.005241 0.005428 0.006712 0.07% MLCellLinOp::smooth() 1560 0.005604 0.005838 0.006382 0.07% FillPatchIterator::Initialize 82 0.005465 0.005578 0.005769 0.06% Castro::construct_old_source() 66 0.005406 0.005498 0.005702 0.06% MLMG::prepareForSolve() 9 0.005028 0.005154 0.00534 0.06% Amr::writeSmallPlotFile() 2 0.003208 0.004843 0.005305 0.06% MLMG::ResNormInf() 74 0.003004 0.003801 0.005198 0.05% Gravity::actual_solve_with_mlmg() 10 0.003489 0.003637 0.004056 0.04% FabArrayBase::getCPC() 1022 0.003754 0.003866 0.004046 0.04% Gravity::get_old_grav_vector() 11 0.003518 0.003614 0.0039 0.04% Castro::initialize_do_advance() 11 0.002367 0.002649 0.00383 0.04% Gravity::get_new_grav_vector() 9 0.00299 0.003113 0.003537 0.04% Castro::sum_integrated_quantities() 2 0.002325 0.002676 0.002946 0.03% MLMG::solve() 9 0.0003665 0.001309 0.002634 0.03% MLCellLinOp::apply() 501 0.001015 0.001127 0.002404 0.03% FabArray::ParallelCopy() 734 0.001783 0.001849 0.001986 0.02% main() 1 0.000961 0.001195 0.001402 0.01% AmrLevel::FillPatch() 82 0.001116 0.001153 0.001207 0.01% Castro::construct_new_gravity() 8 0.0001404 0.0004316 0.001163 0.01% MLCellLinOp::correctionResidual() 390 0.0009313 0.0009572 0.001024 0.01% MLMG::mgVcycle() 65 0.0007985 0.0008278 0.0009081 0.01% Gravity::solve_for_phi() 8 0.0003832 0.0005061 0.0008047 0.01% Castro::save_data_for_retry() 3 0.0007544 0.0007719 0.0007952 0.01% MLMG:computeResOfCorrection() 390 0.0007037 0.0007309 0.0007899 0.01% Castro::Castro() 1 0.0005137 0.0005727 0.0006521 0.01% Castro::subcycle_advance_ctu() 1 9.423e-05 0.0001476 0.0006013 0.01% AmrMesh::MakeNewGrids() 2 0.0005249 0.0005463 0.0005649 0.01% MLMG::actualBottomSolve() 65 6.604e-05 9.771e-05 0.0005355 0.01% StateDataPhysBCFunct::() 82 0.0003623 0.0004245 0.0004934 0.01% FillPatchIterator::FillFromLevel0() 82 0.0003688 0.0003798 0.0003936 0.00% MLMG::mgVcycle_down::2 65 0.0003123 0.0003422 0.0003717 0.00% MLMG::mgVcycle_down::3 65 0.0002892 0.0003156 0.0003698 0.00% MLMG::mgVcycle_down::4 65 0.0002826 0.0002967 0.0003657 0.00% Castro::clean_state() 62 0.0003177 0.0003371 0.0003628 0.00% MLMG::mgVcycle_down::0 65 0.0003372 0.0003476 0.0003561 0.00% MLMG::mgVcycle_down::1 65 0.000335 0.0003441 0.0003545 0.00% MLMG::mgVcycle_down::5 65 0.000293 0.0003028 0.0003404 0.00% Castro::do_advance_ctu() 11 0.0002244 0.0002453 0.0002944 0.00% FillPatchSingleLevel 82 0.0002565 0.0002671 0.0002803 0.00% Castro::problem_post_timestep() 2 1.053e-05 2.755e-05 0.0002693 0.00% MLMG::mgVcycle_up::1 65 0.000237 0.000248 0.0002674 0.00% MLMG::mgVcycle_up::0 65 0.0002428 0.0002477 0.0002538 0.00% MLCellLinOp::solutionResidual() 75 0.0002281 0.0002361 0.0002519 0.00% MLMG::mgVcycle_up::2 65 0.0001951 0.0002169 0.0002448 0.00% MLMG::mgVcycle_up::5 65 0.0002061 0.0002113 0.000231 0.00% MLMG::mgVcycle_up::3 65 0.0001947 0.0002026 0.0002299 0.00% MLMG::mgVcycle_up::4 65 0.0001964 0.0002031 0.0002273 0.00% MLMG::oneIter() 65 0.0001956 0.0002049 0.0002141 0.00% Gravity::update_max_rhs() 2 0.0001593 0.0001638 0.0001689 0.00% Amr::coarseTimeStep() 1 0.000125 0.0001483 0.0001539 0.00% Castro::initialize_advance() 1 0.0001136 0.0001223 0.0001521 0.00% TagBoxArray::mapPRD 2 0.0001131 0.0001215 0.0001454 0.00% Castro::retry_advance_ctu() 11 5.348e-05 6.85e-05 0.0001432 0.00% MLMG::computeResidual() 65 0.0001197 0.0001243 0.0001422 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001166 0.0001254 0.0001375 0.00% Amr::grid_places() 2 8.108e-05 9.932e-05 0.000134 0.00% MLMG::mgVcycle_bottom 65 0.0001128 0.0001165 0.000127 0.00% Amr::defBaseLevel() 1 0.0001062 0.0001151 0.0001258 0.00% Gravity::test_composite_phi() 1 9.199e-05 9.65e-05 0.000112 0.00% MLPoisson::define() 10 6.759e-05 7.219e-05 7.721e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.042e-05 6.513e-05 7.68e-05 0.00% Amr::FinalizeInit() 1 1.195e-05 2.544e-05 7.383e-05 0.00% Castro::check_for_nan() 19 5.079e-05 5.669e-05 6.997e-05 0.00% Castro::react_state() 19 4.761e-05 5.315e-05 5.903e-05 0.00% Castro::apply_source_to_state() 19 4.905e-05 5.321e-05 5.84e-05 0.00% Castro::create_source_corrector() 8 3.464e-05 3.859e-05 5.083e-05 0.00% Castro::expand_state() 11 2.947e-05 3.424e-05 4.429e-05 0.00% Castro::derive() 18 3.847e-05 3.98e-05 4.159e-05 0.00% Amr::regrid() 1 1.072e-05 1.267e-05 3.182e-05 0.00% Gravity::actual_multilevel_solve() 1 2.616e-05 2.753e-05 2.859e-05 0.00% Amr::InitializeInit() 1 8.306e-06 1.421e-05 2.513e-05 0.00% Amr::timeStep() 1 1.218e-05 1.448e-05 2.255e-05 0.00% MLMG::computeMLResidual() 9 1.979e-05 2.052e-05 2.146e-05 0.00% Castro::post_regrid() 2 1.065e-05 1.161e-05 1.467e-05 0.00% Amr::bldFineLevels() 1 6.635e-06 8.356e-06 1.454e-05 0.00% Castro::post_init() 1 1.183e-05 1.256e-05 1.353e-05 0.00% Castro::advance() 1 6.728e-06 8.16e-06 1.316e-05 0.00% Amr::init() 1 6.277e-06 7.173e-06 8.165e-06 0.00% Castro::post_timestep() 1 5.486e-06 5.85e-06 6.405e-06 0.00% Amr::initialInit() 1 2.817e-06 3.131e-06 3.475e-06 0.00% Other 1448 0.02793 0.03976 0.05911 0.62% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.476 9.476 9.476 100.00% Amr::coarseTimeStep() 1 7.219 7.219 7.219 76.19% Amr::timeStep() 1 7.168 7.168 7.168 75.64% Castro::advance() 1 7.106 7.108 7.11 75.03% Castro::subcycle_advance_ctu() 1 7.054 7.057 7.06 74.51% Castro::do_advance_ctu() 11 7.027 7.032 7.034 74.23% Castro::construct_ctu_hydro_source() 11 3.309 3.329 3.363 35.49% Castro::clean_state() 62 2.97 3.046 3.175 33.51% Castro::computeTemp() 63 1.907 2.009 2.226 23.49% Castro::reset_internal_energy(MultiFab) 63 1.192 1.253 1.431 15.10% Amr::InitAmr() 1 1.324 1.324 1.324 13.98% Amr::init() 1 0.931 0.931 0.9311 9.83% Amr::initialInit() 1 0.8772 0.8773 0.8774 9.26% Castro::initialize_do_advance() 11 0.7895 0.8152 0.8413 8.88% Gravity::solve_phi_with_mlmg() 10 0.7731 0.7951 0.8133 8.58% Castro::evaluate_source_change() 144 0.7164 0.7235 0.7315 7.72% Castro::construct_new_gravity() 8 0.6698 0.692 0.7098 7.49% Castro::do_old_sources() 11 0.6449 0.6788 0.6954 7.34% Gravity::solve_for_phi() 8 0.6522 0.6741 0.6922 7.30% FabArray::ParallelCopy() 734 0.1325 0.4931 0.6752 7.13% Gravity::actual_solve_with_mlmg() 10 0.6652 0.6675 0.6699 7.07% Castro::enforce_min_density() 62 0.6145 0.6306 0.6666 7.03% MLMG::solve() 9 0.6379 0.6402 0.6409 6.76% Amr::InitializeInit() 1 0.4447 0.5114 0.64 6.75% Amr::defBaseLevel() 1 0.4447 0.5114 0.64 6.75% FabArray::ParallelCopy_finish() 734 0.04247 0.4039 0.588 6.21% MLMG::oneIter() 65 0.571 0.5725 0.5746 6.06% amrex::Copy() 523 0.5464 0.5543 0.5688 6.00% MLMG::mgVcycle() 65 0.5636 0.5651 0.5672 5.99% Castro::normalize_species() 62 0.4605 0.4674 0.4973 5.25% Castro::do_new_sources() 8 0.4432 0.462 0.4807 5.07% Castro::reset_internal_energy(Fab) 126 0.4331 0.4385 0.4486 4.73% Amr::FinalizeInit() 1 0.2373 0.3658 0.4326 4.57% MultiFab::AddProduct() 2633 0.4194 0.4233 0.428 4.52% AmrLevel::FillPatch() 82 0.2078 0.3101 0.3931 4.15% Castro::Castro() 1 0.3642 0.3726 0.3901 4.12% Castro::initMFs() 1 0.3597 0.3681 0.385 4.06% MLCellLinOp::smooth() 1560 0.2891 0.3149 0.3793 4.00% FillPatchIterator::Initialize 82 0.1359 0.2384 0.3222 3.40% FillPatchIterator::FillFromLevel0() 82 0.1253 0.2299 0.3146 3.32% FillPatchSingleLevel 82 0.1249 0.2296 0.3142 3.32% MultiFab::Subtract() 125 0.268 0.2709 0.2743 2.90% Castro::initData() 1 0.0764 0.1376 0.2735 2.89% Amr::bldFineLevels() 1 0.04846 0.175 0.2417 2.55% FabArray::setVal() 3453 0.2229 0.2272 0.2313 2.44% MLMG::mgVcycle_up::1 65 0.04912 0.1329 0.2165 2.28% MLPoisson::Fsmooth() 3120 0.1627 0.1771 0.2132 2.25% MLCellLinOp::applyBC() 3630 0.1473 0.1631 0.2004 2.12% Amr::grid_places() 2 0.002357 0.1309 0.1977 2.09% AmrMesh::MakeNewGrids() 2 0.002223 0.1308 0.1976 2.09% TagBoxArray::mapPRD 2 0.0007689 0.1292 0.1961 2.07% MLMG::addInterpCorrection() 390 0.03217 0.1547 0.1898 2.00% Castro::expand_state() 11 0.1 0.1272 0.1558 1.64% FabArray::FillBoundary() 3240 0.1158 0.1306 0.1505 1.59% Castro::post_init() 1 0.1421 0.1468 0.1493 1.58% Gravity::fill_multipole_BCs() 10 0.107 0.1268 0.1435 1.51% FabArray::sum() 2633 0.139 0.1403 0.1418 1.50% MLMG::mgVcycle_down::0 65 0.1401 0.141 0.1418 1.50% MLMG::mgVcycle_up::2 65 0.005519 0.05896 0.1404 1.48% check_for_negative_density() 11 0.0514 0.08237 0.1349 1.42% Gravity::multilevel_solve_for_new_phi() 1 0.108 0.1081 0.1081 1.14% Gravity::actual_multilevel_solve() 1 0.108 0.108 0.108 1.14% MLMG::mgVcycle_up::0 65 0.101 0.1019 0.1035 1.09% MLCellLinOp::apply() 501 0.0756 0.08088 0.09955 1.05% FabArray::ParallelCopy_nowait() 734 0.08152 0.08729 0.09832 1.04% Amr::writePlotFile() 2 0.09697 0.09705 0.09706 1.02% Castro::post_regrid() 2 0.0855 0.08806 0.09273 0.98% MLMG::mgVcycle_up::3 65 0.00547 0.01201 0.09058 0.96% FillBoundary_finish() 3240 0.05595 0.06445 0.0786 0.83% MLMG:computeResOfCorrection() 390 0.05754 0.06249 0.07515 0.79% MLCellLinOp::correctionResidual() 390 0.05682 0.06176 0.07436 0.78% FillBoundary_nowait() 3240 0.04605 0.05894 0.07366 0.78% Castro::derive() 18 0.06968 0.07048 0.0724 0.76% AmrLevel::derive() 18 0.06964 0.07044 0.07236 0.76% Castro::create_source_corrector() 8 0.01513 0.04825 0.06663 0.70% Castro::post_timestep() 1 0.05813 0.06015 0.06222 0.66% MLMG::mgVcycle_down::1 65 0.05726 0.05941 0.0612 0.65% Castro::construct_new_source() 48 0.05772 0.05839 0.05893 0.62% MLPoisson::Fapply() 501 0.04744 0.04959 0.05722 0.60% Castro::initialize_advance() 1 0.04889 0.05036 0.05303 0.56% MLCellLinOp::solutionResidual() 75 0.04464 0.04705 0.04931 0.52% Castro::construct_old_gravity() 11 0.0314 0.03913 0.04774 0.50% Amr::regrid() 1 0.04399 0.04516 0.04717 0.50% Castro::construct_old_source() 66 0.04234 0.04296 0.04436 0.47% Gravity::get_old_grav_vector() 11 0.0296 0.037 0.04265 0.45% MLMG::computeResidual() 65 0.03749 0.03879 0.04026 0.42% FabArray::Saxpy() 133 0.03409 0.03451 0.03821 0.40% MLMG::mgVcycle_down::2 65 0.006459 0.02103 0.03714 0.39% MLMG::mgVcycle_bottom 65 0.0001789 0.002458 0.03656 0.39% MLMG::actualBottomSolve() 65 6.604e-05 0.002341 0.03643 0.38% MLCGSolver::bicgstab 4 0 0.002229 0.03566 0.38% Castro::construct_new_gravity_source() 8 0.03234 0.03287 0.0334 0.35% FabArray::Xpay() 480 0.02636 0.02756 0.03126 0.33% MLMG::mgVcycle_down::3 65 0.006108 0.009228 0.0304 0.32% Castro::retry_advance_ctu() 11 0.02228 0.0234 0.0284 0.30% amrex::average_down 390 0.01678 0.01991 0.02789 0.29% Castro::construct_old_gravity_source() 11 0.02451 0.02494 0.02601 0.27% Castro::check_for_nan() 19 0.0223 0.02263 0.02346 0.25% VisMF::Write(FabArray) 4 0.01117 0.01257 0.02345 0.25% MultiFab::contains_nan() 19 0.02225 0.02257 0.0234 0.25% Castro::sum_integrated_quantities() 2 0.02168 0.02184 0.02186 0.23% MLPoisson::define() 10 0.01568 0.01822 0.02158 0.23% Gravity::get_new_grav_vector() 9 0.01778 0.01844 0.01921 0.20% MLMG::mgVcycle_down::4 65 0.006076 0.006902 0.01786 0.19% Castro::react_state() 19 0.01524 0.01577 0.01708 0.18% Castro::save_data_for_retry() 3 0.01103 0.0116 0.01691 0.18% Gravity::test_composite_phi() 1 0.01505 0.01509 0.01512 0.16% StateDataPhysBCFunct::() 82 0.0003623 0.00759 0.01444 0.15% Castro::apply_source_to_state() 19 0.01411 0.01419 0.01433 0.15% MLMG::mgVcycle_up::4 65 0.005519 0.006114 0.01398 0.15% StateData::FillBoundary(geom) 62 0 0.007166 0.01395 0.15% Castro::problem_post_timestep() 2 0.008757 0.01139 0.01382 0.15% Castro::wd_update() 2 0.00874 0.01137 0.01381 0.15% MLCellLinOp::defineAuxData() 10 0.01005 0.01081 0.01304 0.14% MLMG::mgVcycle_down::5 65 0.006077 0.006509 0.01187 0.13% MLMG::prepareForSolve() 9 0.01011 0.01053 0.01157 0.12% MLMG::ResNormInf() 74 0.00799 0.008786 0.01019 0.11% MLMG::mgVcycle_up::5 65 0.005564 0.005899 0.01005 0.11% MLMG::computeMLResidual() 9 0.005335 0.00765 0.009519 0.10% FabArray::mult() 40 0.008279 0.008478 0.008975 0.09% Gravity::update_max_rhs() 2 0.001488 0.006142 0.008693 0.09% FabArray::norminf() 143 0.005659 0.00593 0.007972 0.08% Amr::writeSmallPlotFile() 2 0.007901 0.007905 0.007912 0.08% FabArrayBase::getCPC() 1022 0.005215 0.005932 0.007807 0.08% FabArrayBase::getFB() 3240 0.005514 0.005882 0.007584 0.08% Castro::volWgtSum() 46 0.007295 0.007407 0.007526 0.08% Other 1448 0.04477 0.06165 0.08832 0.93% -------------------------------------------------------------------------------------------- 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 16 B 21 B 25 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.07-21-gcdacc4d908e0) finalized