Initializing AMReX (24.05-10-g0d0b184f31a3)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05-10-g0d0b184f31a3) initialized Starting run at 07:10:58 UTC on 2024-05-09. Successfully read inputs file ... Castro git describe: 24.05-6-g4c1fa8f29 AMReX git describe: 24.05-10-g0d0b184f3 Microphysics git describe: 24.05-3-gca9e5320 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.168134013 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001563066 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.011649298 Gravity::multilevel_solve_for_new_phi() time = 0.100830045 Gravity::fill_multipole_BCs() time = 0.011717494 ... 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.049464939 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003959984 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.002256762 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.063378539 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 = -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 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.003742664 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.063168832 on level 0 ... 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 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.006785378 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.06110302 on level 0 ... 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 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.007167082 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.061633958 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.311189119 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015585098 Gravity::solve_for_phi() time = 0.085798331 on level 0 Castro::construct_new_gravity() time = 0.087875386 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.060526005 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.004253714 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.060773745 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.319991677 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01661542 Gravity::solve_for_phi() time = 0.087020865 on level 0 Castro::construct_new_gravity() time = 0.089066392 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.060268699 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.009060893 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.062173161 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.325727493 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015693115 Gravity::solve_for_phi() time = 0.085862626 on level 0 Castro::construct_new_gravity() time = 0.087886122 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.060190215 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.002720823 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.059885972 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.321601474 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016833623 Gravity::solve_for_phi() time = 0.087443023 on level 0 Castro::construct_new_gravity() time = 0.089465816 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.060598815 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.002709594 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.059935011 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.319805117 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016229145 Gravity::solve_for_phi() time = 0.086363182 on level 0 Castro::construct_new_gravity() time = 0.088391895 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.059993139 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.002535999 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.060181025 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.319451622 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016507041 Gravity::solve_for_phi() time = 0.086650544 on level 0 Castro::construct_new_gravity() time = 0.088690159 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.059830713 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.002188287 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.060034813 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.320642948 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017219722 Gravity::solve_for_phi() time = 0.087394614 on level 0 Castro::construct_new_gravity() time = 0.089500588 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.060107532 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.003030187 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.060224216 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.320525003 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016987814 Gravity::solve_for_phi() time = 0.087143804 on level 0 Castro::construct_new_gravity() time = 0.089258287 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.060222225 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.004660057782 [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.090355097 [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.047593776 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004022818 seconds Ending run at 07:11:07 UTC on 2024-05-09. Run time = 9.380573626 Run time without initialization = 7.142744096 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.381 ... 9.381 ... 9.381 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8313 2.8862 2.9318 31.25% Amr::InitAmr() 1 1.341 1.341 1.341 14.30% Castro::computeTemp() 63 0.6969 0.7601 0.842 8.98% amrex::Copy() 523 0.547 0.5548 0.5663 6.04% FabArray::ParallelCopy_finish() 734 0.0288 0.3172 0.5002 5.33% Castro::normalize_species() 62 0.4415 0.4497 0.4808 5.13% Castro::reset_internal_energy(Fab) 126 0.4303 0.4375 0.4485 4.78% MultiFab::AddProduct() 2633 0.4304 0.4365 0.4444 4.74% Castro::initMFs() 1 0.3634 0.3691 0.4032 4.30% MultiFab::Subtract() 125 0.2616 0.2653 0.2705 2.88% FabArray::setVal() 3453 0.2207 0.2252 0.2327 2.48% MLPoisson::Fsmooth() 3120 0.1598 0.1739 0.2126 2.27% Castro::reset_internal_energy(MultiFab) 63 0.01781 0.0597 0.1877 2.00% Castro::initData() 1 0.003333 0.02858 0.167 1.78% Gravity::fill_multipole_BCs() 10 0.1147 0.1306 0.153 1.63% FabArray::sum() 2633 0.1467 0.1488 0.1514 1.61% check_for_negative_density() 11 0.02225 0.07168 0.1265 1.35% FabArray::ParallelCopy_nowait() 734 0.06297 0.0699 0.07698 0.82% FillBoundary_finish() 3240 0.04485 0.05606 0.06716 0.72% FillBoundary_nowait() 3240 0.0408 0.05241 0.06534 0.70% Castro::enforce_min_density() 62 0.01372 0.02197 0.06181 0.66% MLPoisson::Fapply() 501 0.04782 0.05043 0.05872 0.63% MLCellLinOp::applyBC() 3630 0.02366 0.03165 0.0531 0.57% AmrLevel::derive() 18 0.04102 0.04177 0.04254 0.45% FabArray::Saxpy() 133 0.03587 0.03667 0.0409 0.44% Castro::do_new_sources() 8 0.0001689 0.006091 0.03968 0.42% Castro::do_old_sources() 11 0.0002833 0.006262 0.03559 0.38% Castro::construct_new_gravity_source() 8 0.03226 0.0331 0.03396 0.36% FabArray::Xpay() 480 0.02609 0.02747 0.03139 0.33% Castro::construct_old_gravity_source() 11 0.02539 0.02585 0.02643 0.28% MultiFab::contains_nan() 19 0.02243 0.0228 0.02325 0.25% VisMF::Write(FabArray) 4 0.01159 0.01299 0.02302 0.25% MLMG::addInterpCorrection() 390 0.01683 0.01821 0.02192 0.23% Amr::writePlotFile() 2 0.007296 0.01662 0.01908 0.20% amrex::average_down 390 0.01232 0.01343 0.01651 0.18% StateData::FillBoundary(geom) 62 0 0.006996 0.01363 0.15% Castro::wd_update() 2 0.005024 0.008458 0.01293 0.14% Castro::evaluate_source_change() 144 0.009123 0.009444 0.01001 0.11% MLCellLinOp::defineAuxData() 10 0.007527 0.007881 0.00891 0.09% FabArray::FillBoundary() 3240 0.006855 0.007147 0.008632 0.09% FabArray::mult() 40 0.008041 0.008174 0.008309 0.09% MLCGSolver::bicgstab 4 0 0.0005129 0.008206 0.09% FabArray::norminf() 143 0.005579 0.005869 0.008 0.09% Castro::construct_new_source() 48 0.007249 0.007423 0.007836 0.08% Castro::volWgtSum() 46 0.007371 0.007479 0.007592 0.08% FabArrayBase::getFB() 3240 0.00496 0.005205 0.006493 0.07% MLCellLinOp::smooth() 1560 0.005561 0.0058 0.006194 0.07% Castro::construct_old_gravity() 11 0.0002089 0.001437 0.006182 0.07% FillPatchIterator::Initialize 82 0.005423 0.005614 0.005922 0.06% Castro::construct_old_source() 66 0.005408 0.00553 0.00584 0.06% Amr::writeSmallPlotFile() 2 0.003683 0.005322 0.005683 0.06% MLMG::prepareForSolve() 9 0.004953 0.005077 0.005362 0.06% Gravity::actual_solve_with_mlmg() 10 0.003506 0.003669 0.004166 0.04% MLMG::ResNormInf() 74 0.001026 0.002736 0.004009 0.04% FabArrayBase::getCPC() 1022 0.00364 0.003747 0.003929 0.04% Gravity::get_old_grav_vector() 11 0.003547 0.003627 0.003742 0.04% MLMG::solve() 9 0.0003944 0.001818 0.003234 0.03% Gravity::get_new_grav_vector() 9 0.003014 0.003073 0.003198 0.03% Castro::initialize_do_advance() 11 0.002379 0.002577 0.002862 0.03% Castro::sum_integrated_quantities() 2 0.001647 0.00205 0.00236 0.03% MLCellLinOp::apply() 501 0.0009937 0.001108 0.002344 0.02% FabArray::ParallelCopy() 734 0.001783 0.001881 0.002049 0.02% AmrLevel::FillPatch() 82 0.001078 0.001133 0.001254 0.01% main() 1 0.0009226 0.001102 0.001245 0.01% MLCellLinOp::correctionResidual() 390 0.0009261 0.0009624 0.0009975 0.01% MLMG::mgVcycle() 65 0.0007959 0.0008413 0.0009302 0.01% Castro::save_data_for_retry() 3 0.0007514 0.0007684 0.0007918 0.01% MLMG:computeResOfCorrection() 390 0.0006848 0.0007193 0.0007738 0.01% Castro::construct_new_gravity() 8 0.0001371 0.0002661 0.0007189 0.01% Gravity::solve_for_phi() 8 0.0003717 0.0004532 0.0006978 0.01% Castro::Castro() 1 0.0004972 0.0005594 0.0006323 0.01% Castro::subcycle_advance_ctu() 1 9.052e-05 0.000139 0.0005886 0.01% MLMG::actualBottomSolve() 65 6.827e-05 9.995e-05 0.0005584 0.01% AmrMesh::MakeNewGrids() 2 0.0005062 0.0005211 0.0005394 0.01% StateDataPhysBCFunct::() 82 0.0003367 0.000408 0.0005086 0.01% FillPatchIterator::FillFromLevel0() 82 0.000362 0.0003823 0.0004081 0.00% Castro::clean_state() 62 0.0003313 0.0003541 0.0004048 0.00% MLMG::mgVcycle_down::0 65 0.0003342 0.0003449 0.0003657 0.00% MLMG::mgVcycle_down::3 65 0.0002942 0.0003145 0.0003578 0.00% MLMG::mgVcycle_down::2 65 0.0003067 0.0003338 0.0003556 0.00% MLMG::mgVcycle_down::1 65 0.000325 0.0003374 0.0003518 0.00% MLMG::mgVcycle_down::4 65 0.0002799 0.0002949 0.0003398 0.00% MLMG::mgVcycle_down::5 65 0.0002889 0.000301 0.0003391 0.00% Amr::defBaseLevel() 1 0.0001099 0.0001394 0.0003373 0.00% FillPatchSingleLevel 82 0.0002538 0.0002624 0.0002805 0.00% Castro::problem_post_timestep() 2 1.045e-05 2.759e-05 0.0002644 0.00% Castro::do_advance_ctu() 11 0.0002209 0.0002354 0.0002578 0.00% MLMG::mgVcycle_up::1 65 0.0002322 0.0002439 0.0002578 0.00% MLMG::mgVcycle_up::0 65 0.0002303 0.0002381 0.0002547 0.00% MLMG::mgVcycle_up::5 65 0.0002024 0.0002106 0.000245 0.00% MLMG::mgVcycle_up::2 65 0.0001957 0.0002173 0.0002417 0.00% MLCellLinOp::solutionResidual() 75 0.0002198 0.0002278 0.0002364 0.00% MLMG::mgVcycle_up::3 65 0.0001942 0.000205 0.0002326 0.00% MLMG::oneIter() 65 0.0001963 0.0002068 0.0002271 0.00% MLMG::mgVcycle_up::4 65 0.0001965 0.0002042 0.0002247 0.00% Amr::coarseTimeStep() 1 0.0001197 0.0001454 0.0001501 0.00% Castro::initialize_advance() 1 0.0001107 0.0001194 0.0001462 0.00% Castro::retry_advance_ctu() 11 5.213e-05 7.061e-05 0.0001461 0.00% MLMG::mgVcycle_bottom 65 0.0001172 0.0001222 0.0001387 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001133 0.0001234 0.0001375 0.00% MLMG::computeResidual() 65 0.00012 0.0001236 0.0001299 0.00% Amr::grid_places() 2 7.891e-05 9.629e-05 0.0001281 0.00% TagBoxArray::mapPRD 2 0.0001114 0.0001181 0.0001274 0.00% Gravity::test_composite_phi() 1 9.184e-05 9.74e-05 0.0001117 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.086e-05 6.99e-05 0.0001034 0.00% MLPoisson::define() 10 6.757e-05 7.317e-05 8.447e-05 0.00% Castro::check_for_nan() 19 4.945e-05 5.85e-05 8.256e-05 0.00% Castro::react_state() 19 5.148e-05 5.465e-05 6.275e-05 0.00% Amr::FinalizeInit() 1 1.063e-05 1.911e-05 6.221e-05 0.00% Castro::apply_source_to_state() 19 4.933e-05 5.282e-05 5.59e-05 0.00% Castro::derive() 18 4.106e-05 4.23e-05 4.563e-05 0.00% Castro::create_source_corrector() 8 3.42e-05 3.866e-05 4.531e-05 0.00% Castro::expand_state() 11 3.158e-05 3.491e-05 3.896e-05 0.00% Amr::regrid() 1 1.103e-05 1.347e-05 3.303e-05 0.00% Gravity::actual_multilevel_solve() 1 2.609e-05 2.767e-05 3.026e-05 0.00% MLMG::computeMLResidual() 9 2.027e-05 2.11e-05 2.203e-05 0.00% Amr::timeStep() 1 1.201e-05 1.505e-05 2.17e-05 0.00% Amr::init() 1 6.639e-06 8.454e-06 1.83e-05 0.00% Castro::post_regrid() 2 1.117e-05 1.257e-05 1.529e-05 0.00% Castro::post_init() 1 1.227e-05 1.323e-05 1.419e-05 0.00% Amr::InitializeInit() 1 7.156e-06 8.162e-06 1.324e-05 0.00% Castro::advance() 1 7.113e-06 8.454e-06 1.256e-05 0.00% Amr::bldFineLevels() 1 7.016e-06 8.769e-06 1.131e-05 0.00% Castro::post_timestep() 1 5.757e-06 6.185e-06 6.735e-06 0.00% Amr::initialInit() 1 3.034e-06 3.395e-06 6.29e-06 0.00% Other 1450 0.0275 0.03639 0.05747 0.61% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.381 9.381 9.381 100.00% Amr::coarseTimeStep() 1 7.142 7.142 7.142 76.14% Amr::timeStep() 1 7.09 7.09 7.09 75.58% Castro::advance() 1 7.026 7.029 7.032 74.96% Castro::subcycle_advance_ctu() 1 6.976 6.979 6.981 74.42% Castro::do_advance_ctu() 11 6.952 6.955 6.957 74.16% Castro::construct_ctu_hydro_source() 11 3.306 3.328 3.349 35.70% Castro::clean_state() 62 2.953 3.039 3.135 33.42% Castro::computeTemp() 63 1.919 2.013 2.194 23.39% Castro::reset_internal_energy(MultiFab) 63 1.194 1.253 1.397 14.89% Amr::InitAmr() 1 1.341 1.341 1.341 14.30% Amr::init() 1 0.8958 0.8958 0.8958 9.55% Amr::initialInit() 1 0.8416 0.8416 0.8417 8.97% Gravity::solve_phi_with_mlmg() 10 0.7658 0.783 0.8076 8.61% Castro::initialize_do_advance() 11 0.7563 0.7749 0.7992 8.52% Castro::evaluate_source_change() 144 0.734 0.7447 0.7576 8.08% Castro::construct_new_gravity() 8 0.6696 0.6861 0.71 7.57% Gravity::solve_for_phi() 8 0.6521 0.6691 0.6936 7.39% Castro::do_old_sources() 11 0.6411 0.6566 0.67 7.14% Castro::enforce_min_density() 62 0.6191 0.6358 0.669 7.13% Gravity::actual_solve_with_mlmg() 10 0.6487 0.6506 0.6528 6.96% MLMG::solve() 9 0.6221 0.6242 0.6249 6.66% Amr::InitializeInit() 1 0.4485 0.5104 0.6153 6.56% Amr::defBaseLevel() 1 0.4485 0.5104 0.6153 6.56% FabArray::ParallelCopy() 734 0.1192 0.4054 0.5863 6.25% amrex::Copy() 523 0.547 0.5548 0.5663 6.04% MLMG::oneIter() 65 0.5566 0.5577 0.5591 5.96% MLMG::mgVcycle() 65 0.5493 0.5503 0.5517 5.88% FabArray::ParallelCopy_finish() 734 0.0288 0.3172 0.5002 5.33% Castro::do_new_sources() 8 0.4429 0.4648 0.4818 5.14% Castro::normalize_species() 62 0.4415 0.4497 0.4808 5.13% Castro::reset_internal_energy(Fab) 126 0.4303 0.4375 0.4485 4.78% MultiFab::AddProduct() 2633 0.4304 0.4365 0.4444 4.74% Castro::Castro() 1 0.368 0.3737 0.4086 4.36% Castro::initMFs() 1 0.3636 0.3692 0.4034 4.30% Amr::FinalizeInit() 1 0.2263 0.3312 0.3931 4.19% MLCellLinOp::smooth() 1560 0.278 0.3032 0.3668 3.91% AmrLevel::FillPatch() 82 0.1954 0.2484 0.3286 3.50% MultiFab::Subtract() 125 0.2616 0.2653 0.2705 2.88% FillPatchIterator::Initialize 82 0.1225 0.1762 0.2574 2.74% FillPatchIterator::FillFromLevel0() 82 0.1119 0.168 0.25 2.67% FillPatchSingleLevel 82 0.1115 0.1676 0.2496 2.66% Castro::initData() 1 0.07644 0.1355 0.2459 2.62% FabArray::setVal() 3453 0.2207 0.2252 0.2327 2.48% MLMG::mgVcycle_up::1 65 0.04776 0.13 0.2132 2.27% Amr::bldFineLevels() 1 0.04739 0.1516 0.2132 2.27% MLPoisson::Fsmooth() 3120 0.1598 0.1739 0.2126 2.27% MLMG::addInterpCorrection() 390 0.03215 0.1523 0.1876 2.00% MLCellLinOp::applyBC() 3630 0.1349 0.1529 0.1862 1.99% Amr::grid_places() 2 0.002919 0.1078 0.1697 1.81% AmrMesh::MakeNewGrids() 2 0.002791 0.1077 0.1696 1.81% TagBoxArray::mapPRD 2 0.0007743 0.1055 0.1682 1.79% Gravity::fill_multipole_BCs() 10 0.1159 0.1317 0.154 1.64% FabArray::sum() 2633 0.1467 0.1488 0.1514 1.61% FabArray::FillBoundary() 3240 0.1045 0.1212 0.1424 1.52% Castro::expand_state() 11 0.09771 0.1142 0.1404 1.50% MLMG::mgVcycle_up::2 65 0.005373 0.05834 0.1393 1.48% Castro::post_init() 1 0.1343 0.1359 0.1378 1.47% MLMG::mgVcycle_down::0 65 0.1316 0.1329 0.1349 1.44% check_for_negative_density() 11 0.02225 0.07168 0.1265 1.35% Gravity::multilevel_solve_for_new_phi() 1 0.1008 0.1008 0.1009 1.08% Gravity::actual_multilevel_solve() 1 0.1007 0.1007 0.1008 1.07% MLMG::mgVcycle_up::0 65 0.09881 0.09904 0.09927 1.06% MLCellLinOp::apply() 501 0.07425 0.08006 0.09849 1.05% Amr::writePlotFile() 2 0.09756 0.09764 0.09765 1.04% FabArray::ParallelCopy_nowait() 734 0.07913 0.08635 0.09369 1.00% MLMG::mgVcycle_up::3 65 0.005351 0.01197 0.09085 0.97% Castro::post_regrid() 2 0.08561 0.08748 0.08908 0.95% MLMG:computeResOfCorrection() 390 0.05723 0.06209 0.07436 0.79% MLCellLinOp::correctionResidual() 390 0.05652 0.06137 0.07358 0.78% Castro::derive() 18 0.06938 0.07057 0.07177 0.77% AmrLevel::derive() 18 0.06933 0.07053 0.07173 0.76% FillBoundary_nowait() 3240 0.04636 0.05804 0.07116 0.76% FillBoundary_finish() 3240 0.04485 0.05606 0.06716 0.72% Castro::post_timestep() 1 0.05864 0.06128 0.06471 0.69% MLMG::mgVcycle_down::1 65 0.05817 0.05983 0.06155 0.66% Castro::construct_new_source() 48 0.05802 0.05932 0.0609 0.65% MLPoisson::Fapply() 501 0.04782 0.05043 0.05872 0.63% Castro::initialize_advance() 1 0.04886 0.05001 0.05163 0.55% MLCellLinOp::solutionResidual() 75 0.04438 0.04654 0.04885 0.52% Castro::create_source_corrector() 8 0.01507 0.03062 0.04666 0.50% Amr::regrid() 1 0.04443 0.04538 0.04638 0.49% Castro::construct_old_source() 66 0.04365 0.04435 0.04536 0.48% FabArray::Saxpy() 133 0.03587 0.03667 0.0409 0.44% Castro::construct_old_gravity() 11 0.02276 0.03182 0.04014 0.43% MLMG::computeResidual() 65 0.0382 0.03854 0.03881 0.41% Gravity::get_old_grav_vector() 11 0.02211 0.03038 0.03663 0.39% MLMG::mgVcycle_bottom 65 0.0001865 0.002456 0.03643 0.39% MLMG::actualBottomSolve() 65 6.827e-05 0.002334 0.0363 0.39% MLMG::mgVcycle_down::2 65 0.006273 0.02056 0.03627 0.39% MLCGSolver::bicgstab 4 0 0.002219 0.03551 0.38% Castro::construct_new_gravity_source() 8 0.03226 0.0331 0.03396 0.36% FabArray::Xpay() 480 0.02609 0.02747 0.03139 0.33% MLMG::mgVcycle_down::3 65 0.005997 0.009147 0.03048 0.32% amrex::average_down 390 0.01646 0.01973 0.02784 0.30% Castro::construct_old_gravity_source() 11 0.02539 0.02585 0.02643 0.28% Castro::retry_advance_ctu() 11 0.02214 0.02272 0.02391 0.25% Castro::check_for_nan() 19 0.02248 0.02286 0.02333 0.25% MultiFab::contains_nan() 19 0.02243 0.0228 0.02325 0.25% VisMF::Write(FabArray) 4 0.01159 0.01299 0.02302 0.25% Castro::sum_integrated_quantities() 2 0.02113 0.02129 0.02131 0.23% MLPoisson::define() 10 0.01516 0.01734 0.0202 0.22% Gravity::get_new_grav_vector() 9 0.01721 0.0177 0.01854 0.20% MLMG::mgVcycle_down::4 65 0.00596 0.00684 0.01769 0.19% Castro::react_state() 19 0.01495 0.01529 0.01675 0.18% Castro::problem_post_timestep() 2 0.009134 0.01209 0.01625 0.17% Castro::wd_update() 2 0.00887 0.01206 0.01624 0.17% Castro::apply_source_to_state() 19 0.01488 0.0151 0.01536 0.16% Gravity::test_composite_phi() 1 0.01514 0.01518 0.01521 0.16% StateDataPhysBCFunct::() 82 0.0003367 0.007404 0.01407 0.15% MLMG::mgVcycle_up::4 65 0.005416 0.006064 0.01387 0.15% StateData::FillBoundary(geom) 62 0 0.006996 0.01363 0.15% MLCellLinOp::defineAuxData() 10 0.009776 0.01047 0.01249 0.13% MLMG::mgVcycle_down::5 65 0.005941 0.006446 0.01174 0.13% MLMG::prepareForSolve() 9 0.009988 0.01043 0.0117 0.12% Castro::save_data_for_retry() 3 0.01105 0.01125 0.01156 0.12% MLMG::mgVcycle_up::5 65 0.005464 0.005866 0.009994 0.11% MLMG::computeMLResidual() 9 0.005464 0.007417 0.009645 0.10% MLMG::ResNormInf() 74 0.006011 0.007686 0.008898 0.09% Amr::writeSmallPlotFile() 2 0.00836 0.008363 0.008367 0.09% FabArray::mult() 40 0.008041 0.008174 0.008309 0.09% FabArray::norminf() 143 0.005579 0.005869 0.008 0.09% Castro::volWgtSum() 46 0.007371 0.007479 0.007592 0.08% FabArrayBase::getCPC() 1022 0.005049 0.005649 0.007396 0.08% FabArrayBase::getFB() 3240 0.005206 0.005628 0.007296 0.08% Other 1450 0.04541 0.0601 0.09021 0.96% -------------------------------------------------------------------------------------------- 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 20 B 23 B 30 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-10-g0d0b184f31a3) finalized