Initializing AMReX (24.05-26-ge6c93bf22695)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05-26-ge6c93bf22695) initialized Starting run at 07:15:22 UTC on 2024-05-23. Successfully read inputs file ... Castro git describe: 24.05-13-gf7309fc71 AMReX git describe: 24.05-26-ge6c93bf22 Microphysics git describe: 24.05-12-g0dd7c3c7 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.182887038 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.002993594 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.011809885 Gravity::multilevel_solve_for_new_phi() time = 0.107216428 Gravity::fill_multipole_BCs() time = 0.011268442 ... 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.052745464 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.008477482 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.00791802 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.086728601 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.036525947 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.076632059 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.006156906 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.072350753 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.006290178 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.065251554 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.323430736 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01615511 Gravity::solve_for_phi() time = 0.092157785 on level 0 Castro::construct_new_gravity() time = 0.094213653 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.061631573 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.005109157 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.092483739 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.391435732 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017890345 Gravity::solve_for_phi() time = 0.089944008 on level 0 Castro::construct_new_gravity() time = 0.092043278 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.06184027 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.002507631 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.060228292 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.330250011 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016414392 Gravity::solve_for_phi() time = 0.089843865 on level 0 Castro::construct_new_gravity() time = 0.093225704 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.061648954 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.003490288 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.064193602 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.340050475 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015642627 Gravity::solve_for_phi() time = 0.08597745 on level 0 Castro::construct_new_gravity() time = 0.088094856 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.060093184 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.002152099 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.059572388 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.324118678 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015649473 Gravity::solve_for_phi() time = 0.086066656 on level 0 Castro::construct_new_gravity() time = 0.088177797 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.059966713 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.002726328 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.060681819 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.324268307 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015525663 Gravity::solve_for_phi() time = 0.085224642 on level 0 Castro::construct_new_gravity() time = 0.087382463 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.06001393 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.002288071 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.0598615 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.325516802 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01585047 Gravity::solve_for_phi() time = 0.086026591 on level 0 Castro::construct_new_gravity() time = 0.088140597 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.060038483 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.003085221 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.060294915 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.325957578 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015233394 Gravity::solve_for_phi() time = 0.085396726 on level 0 Castro::construct_new_gravity() time = 0.087593564 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.060269353 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.004458497898 [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.408843105 [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.050324907 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.007737735 seconds Ending run at 07:15:32 UTC on 2024-05-23. Run time = 10.36497653 Run time without initialization = 7.467861103 Average number of zones advanced per microsecond: 0.004 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 10.36 ... 10.37 ... 10.37 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8566 2.9142 2.9985 28.93% Amr::InitAmr() 1 1.97 1.97 1.97 19.01% Castro::computeTemp() 63 0.7057 0.7687 0.8738 8.43% FabArray::ParallelCopy_finish() 734 0.03348 0.4261 0.6545 6.31% amrex::Copy() 523 0.5504 0.5618 0.5798 5.59% Castro::normalize_species() 62 0.4451 0.452 0.4671 4.51% MultiFab::AddProduct() 2633 0.4324 0.4407 0.4644 4.48% Castro::reset_internal_energy(Fab) 126 0.4332 0.4424 0.4589 4.43% Castro::initMFs() 1 0.3593 0.3668 0.3788 3.65% Castro::reset_internal_energy(MultiFab) 63 0.01799 0.07546 0.2869 2.77% MultiFab::Subtract() 125 0.2616 0.2667 0.2801 2.70% check_for_negative_density() 11 0.102 0.1867 0.2523 2.43% FabArray::setVal() 3453 0.221 0.2262 0.2375 2.29% MLPoisson::Fsmooth() 3120 0.1596 0.1739 0.2119 2.04% Castro::initData() 1 0.003262 0.02915 0.1708 1.65% FabArray::sum() 2633 0.1475 0.1503 0.157 1.51% Gravity::fill_multipole_BCs() 10 0.1084 0.128 0.1498 1.45% Castro::enforce_min_density() 62 0.01418 0.03207 0.0893 0.86% FabArray::ParallelCopy_nowait() 734 0.06592 0.07158 0.08068 0.78% FillBoundary_finish() 3240 0.05688 0.06748 0.07907 0.76% FillBoundary_nowait() 3240 0.0414 0.05306 0.0684 0.66% Castro::do_old_sources() 11 0.0003212 0.01331 0.06302 0.61% MLPoisson::Fapply() 501 0.04796 0.05061 0.05858 0.57% MLCellLinOp::applyBC() 3630 0.0245 0.03209 0.05282 0.51% AmrLevel::derive() 18 0.04163 0.04218 0.04304 0.42% FabArray::Saxpy() 133 0.03593 0.03684 0.04073 0.39% Castro::do_new_sources() 8 0.0001734 0.006362 0.03862 0.37% Castro::construct_new_gravity_source() 8 0.03235 0.03299 0.03405 0.33% FabArray::Xpay() 480 0.02622 0.02759 0.03153 0.30% VisMF::Write(FabArray) 4 0.01475 0.01646 0.03121 0.30% Castro::construct_old_gravity_source() 11 0.02566 0.02617 0.0272 0.26% MultiFab::contains_nan() 19 0.02306 0.02361 0.02479 0.24% Amr::writePlotFile() 2 0.008856 0.0208 0.02345 0.23% MLMG::addInterpCorrection() 390 0.01692 0.01826 0.02194 0.21% amrex::average_down 390 0.01232 0.01347 0.01646 0.16% Castro::construct_old_gravity() 11 0.0002382 0.004256 0.01479 0.14% StateData::FillBoundary(geom) 62 0 0.007105 0.01364 0.13% Castro::wd_update() 2 0.005137 0.007861 0.0126 0.12% Amr::writeSmallPlotFile() 2 0.007545 0.01077 0.01142 0.11% Castro::evaluate_source_change() 144 0.009479 0.009835 0.01055 0.10% MultiFab::max() 2 0.0001626 0.006592 0.009527 0.09% FabArray::mult() 40 0.007845 0.008093 0.00918 0.09% MLCellLinOp::defineAuxData() 10 0.007485 0.007941 0.009029 0.09% FabArray::FillBoundary() 3240 0.006926 0.007243 0.00865 0.08% MLCGSolver::bicgstab 4 0 0.000519 0.008304 0.08% FabArray::norminf() 143 0.00563 0.005895 0.008033 0.08% MLMG::ResNormInf() 74 0.003081 0.005836 0.00784 0.08% Castro::construct_new_source() 48 0.007265 0.007436 0.00764 0.07% amrex::Add() 65 0.007138 0.007279 0.007595 0.07% FillPatchIterator::Initialize 82 0.005558 0.005935 0.006362 0.06% MLCellLinOp::smooth() 1560 0.005751 0.005974 0.006336 0.06% Castro::construct_old_source() 66 0.005494 0.005587 0.005755 0.06% MLMG::prepareForSolve() 9 0.004958 0.005145 0.005554 0.05% Gravity::actual_solve_with_mlmg() 10 0.003549 0.0037 0.004163 0.04% Gravity::get_old_grav_vector() 11 0.003552 0.003652 0.003769 0.04% MLMG::solve() 9 0.000504 0.00151 0.003452 0.03% Castro::initialize_do_advance() 11 0.002666 0.002991 0.00321 0.03% Gravity::get_new_grav_vector() 9 0.003032 0.003087 0.003208 0.03% MLCellLinOp::apply() 501 0.001026 0.001138 0.002426 0.02% Castro::sum_integrated_quantities() 2 0.001627 0.0021 0.002373 0.02% Castro::save_data_for_retry() 3 0.0007637 0.001257 0.002137 0.02% Gravity::solve_for_phi() 8 0.0003978 0.0005788 0.002047 0.02% FabArray::ParallelCopy() 734 0.001796 0.001879 0.002029 0.02% Castro::construct_new_gravity() 8 0.0001456 0.0005135 0.001773 0.02% main() 1 0.001036 0.001204 0.001343 0.01% AmrMesh::MakeNewGrids() 2 0.0005035 0.0005696 0.001318 0.01% AmrLevel::FillPatch() 82 0.001125 0.001174 0.00128 0.01% MLCellLinOp::correctionResidual() 390 0.0009285 0.0009666 0.00101 0.01% MLMG::mgVcycle() 65 0.0007903 0.0008339 0.000915 0.01% MLMG:computeResOfCorrection() 390 0.0006941 0.0007246 0.0007633 0.01% Castro::subcycle_advance_ctu() 1 9.712e-05 0.0001521 0.0006185 0.01% MLMG::actualBottomSolve() 65 6.818e-05 0.000101 0.0005589 0.01% Castro::Castro() 1 0.0004463 0.0004897 0.0005564 0.01% StateDataPhysBCFunct::() 82 0.0003703 0.0004257 0.0005079 0.00% Castro::problem_post_timestep() 2 1.135e-05 3.722e-05 0.0004127 0.00% Castro::clean_state() 62 0.0003223 0.0003609 0.0004011 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003619 0.0003849 0.0004009 0.00% MLMG::mgVcycle_down::0 65 0.0003334 0.0003471 0.0003815 0.00% MLMG::mgVcycle_down::1 65 0.00033 0.0003421 0.000374 0.00% MLMG::mgVcycle_down::2 65 0.0003183 0.0003445 0.0003731 0.00% MLMG::mgVcycle_down::3 65 0.0002971 0.0003213 0.000372 0.00% MLMG::mgVcycle_down::4 65 0.0002882 0.0003018 0.0003454 0.00% MLMG::mgVcycle_down::5 65 0.0002952 0.0003111 0.0003405 0.00% FillPatchSingleLevel 82 0.000257 0.000269 0.0002962 0.00% Castro::do_advance_ctu() 11 0.0002225 0.0002543 0.0002849 0.00% MLMG::mgVcycle_up::0 65 0.0002306 0.00024 0.0002623 0.00% MLMG::mgVcycle_up::1 65 0.0002317 0.0002404 0.0002516 0.00% TagBoxArray::mapPRD 2 0.0001073 0.0001222 0.0002506 0.00% MLCellLinOp::solutionResidual() 75 0.0002196 0.000229 0.0002448 0.00% MLMG::mgVcycle_up::2 65 0.000199 0.0002168 0.0002396 0.00% MLMG::mgVcycle_up::3 65 0.0001945 0.0002051 0.0002313 0.00% MLMG::mgVcycle_up::5 65 0.0002015 0.0002093 0.0002313 0.00% MLMG::mgVcycle_up::4 65 0.0001969 0.0002068 0.000228 0.00% MLMG::oneIter() 65 0.0002 0.0002065 0.0002239 0.00% Gravity::update_max_rhs() 2 0.0001643 0.0001737 0.000187 0.00% Amr::coarseTimeStep() 1 0.0001353 0.000167 0.0001758 0.00% Amr::grid_places() 2 8.348e-05 0.0001132 0.0001684 0.00% Castro::initialize_advance() 1 0.000128 0.0001423 0.0001552 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001131 0.000124 0.000145 0.00% Castro::retry_advance_ctu() 11 5.416e-05 7.886e-05 0.0001404 0.00% MLMG::mgVcycle_bottom 65 0.0001149 0.0001209 0.000129 0.00% MLMG::computeResidual() 65 0.0001189 0.0001231 0.0001276 0.00% Amr::defBaseLevel() 1 9.806e-05 0.0001085 0.0001264 0.00% Gravity::test_composite_phi() 1 9.316e-05 9.978e-05 0.0001186 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.165e-05 7.259e-05 9.861e-05 0.00% MLPoisson::define() 10 6.371e-05 7.097e-05 8.087e-05 0.00% Castro::check_for_nan() 19 5.153e-05 5.783e-05 6.477e-05 0.00% Castro::apply_source_to_state() 19 5.023e-05 5.607e-05 6.176e-05 0.00% Castro::react_state() 19 5.299e-05 5.689e-05 6.128e-05 0.00% Castro::create_source_corrector() 8 3.688e-05 4.222e-05 5.428e-05 0.00% Castro::derive() 18 4.025e-05 4.322e-05 5.152e-05 0.00% Amr::FinalizeInit() 1 8.368e-06 1.301e-05 5.051e-05 0.00% Castro::expand_state() 11 3.351e-05 3.712e-05 4.433e-05 0.00% Amr::regrid() 1 1.083e-05 1.472e-05 4.416e-05 0.00% Gravity::actual_multilevel_solve() 1 2.662e-05 2.921e-05 3.252e-05 0.00% MLMG::computeMLResidual() 9 2.038e-05 2.262e-05 2.863e-05 0.00% Amr::timeStep() 1 1.222e-05 1.851e-05 2.798e-05 0.00% Amr::bldFineLevels() 1 7.562e-06 1.013e-05 2.51e-05 0.00% Castro::advance() 1 7.18e-06 1.048e-05 1.975e-05 0.00% Castro::post_regrid() 2 1.106e-05 1.262e-05 1.6e-05 0.00% Castro::post_init() 1 1.235e-05 1.33e-05 1.457e-05 0.00% Amr::InitializeInit() 1 6.652e-06 7.67e-06 1.276e-05 0.00% Amr::init() 1 5.922e-06 7.643e-06 9.474e-06 0.00% Castro::post_timestep() 1 5.586e-06 6.207e-06 7.026e-06 0.00% Amr::initialInit() 1 2.913e-06 3.363e-06 3.943e-06 0.00% Other 5689 0.03651 0.04585 0.06925 0.67% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.36 10.36 10.37 100.00% Amr::coarseTimeStep() 1 7.467 7.467 7.467 72.04% Amr::timeStep() 1 7.409 7.409 7.409 71.48% Castro::advance() 1 7.344 7.348 7.35 70.91% Castro::subcycle_advance_ctu() 1 7.291 7.296 7.299 70.41% Castro::do_advance_ctu() 11 7.26 7.271 7.275 70.18% Castro::construct_ctu_hydro_source() 11 3.449 3.472 3.506 33.83% Castro::clean_state() 62 2.984 3.09 3.273 31.58% Castro::computeTemp() 63 1.938 2.053 2.303 22.22% Amr::InitAmr() 1 1.97 1.97 1.97 19.01% Castro::reset_internal_energy(MultiFab) 63 1.201 1.284 1.493 14.40% Amr::init() 1 0.9266 0.9266 0.9267 8.94% Amr::initialInit() 1 0.8646 0.8646 0.8647 8.34% Castro::initialize_do_advance() 11 0.7743 0.8267 0.8607 8.30% Gravity::solve_phi_with_mlmg() 10 0.777 0.7982 0.8204 7.92% Castro::evaluate_source_change() 144 0.7374 0.7517 0.792 7.64% FabArray::ParallelCopy() 734 0.1247 0.5163 0.7419 7.16% Castro::do_old_sources() 11 0.6895 0.7158 0.7345 7.09% Castro::construct_new_gravity() 8 0.676 0.6965 0.7186 6.93% Gravity::solve_for_phi() 8 0.6585 0.6783 0.7005 6.76% Castro::enforce_min_density() 62 0.6222 0.6482 0.6976 6.73% Gravity::actual_solve_with_mlmg() 10 0.6663 0.6684 0.6708 6.47% FabArray::ParallelCopy_finish() 734 0.03348 0.4261 0.6545 6.31% MLMG::solve() 9 0.6391 0.6414 0.6422 6.20% Amr::InitializeInit() 1 0.4418 0.5121 0.6231 6.01% Amr::defBaseLevel() 1 0.4418 0.5121 0.6231 6.01% amrex::Copy() 523 0.5504 0.5618 0.5798 5.59% MLMG::oneIter() 65 0.5676 0.5693 0.5715 5.51% MLMG::mgVcycle() 65 0.56 0.5618 0.5637 5.44% Castro::do_new_sources() 8 0.4447 0.4683 0.4856 4.69% AmrLevel::FillPatch() 82 0.2019 0.3534 0.4693 4.53% Castro::normalize_species() 62 0.4451 0.452 0.4671 4.51% MultiFab::AddProduct() 2633 0.4324 0.4407 0.4644 4.48% Castro::reset_internal_energy(Fab) 126 0.4332 0.4424 0.4589 4.43% Amr::FinalizeInit() 1 0.2416 0.3525 0.4228 4.08% FillPatchIterator::Initialize 82 0.1289 0.2803 0.3977 3.84% FillPatchIterator::FillFromLevel0() 82 0.118 0.2716 0.3901 3.76% FillPatchSingleLevel 82 0.1176 0.2712 0.3897 3.76% Castro::Castro() 1 0.3633 0.3708 0.3828 3.69% Castro::initMFs() 1 0.3594 0.3669 0.379 3.66% MLCellLinOp::smooth() 1560 0.2861 0.3128 0.3773 3.64% MultiFab::Subtract() 125 0.2616 0.2667 0.2801 2.70% Castro::initData() 1 0.07729 0.14 0.254 2.45% check_for_negative_density() 11 0.102 0.1867 0.2523 2.43% FabArray::setVal() 3453 0.221 0.2262 0.2375 2.29% Amr::bldFineLevels() 1 0.0523 0.1611 0.2303 2.22% MLMG::mgVcycle_up::1 65 0.04825 0.1314 0.2151 2.07% MLPoisson::Fsmooth() 3120 0.1596 0.1739 0.2119 2.04% MLCellLinOp::applyBC() 3630 0.147 0.1657 0.1994 1.92% Castro::expand_state() 11 0.0992 0.1509 0.1913 1.85% MLMG::addInterpCorrection() 390 0.03191 0.1531 0.1882 1.82% Amr::grid_places() 2 0.004621 0.1156 0.1859 1.79% AmrMesh::MakeNewGrids() 2 0.004453 0.1155 0.1858 1.79% TagBoxArray::mapPRD 2 0.0007133 0.1119 0.1839 1.77% FabArray::sum() 2633 0.1475 0.1503 0.157 1.51% Gravity::fill_multipole_BCs() 10 0.1095 0.1291 0.1509 1.46% Castro::post_init() 1 0.1408 0.1472 0.1501 1.45% FabArray::FillBoundary() 3240 0.1159 0.1336 0.1494 1.44% MLMG::mgVcycle_up::2 65 0.005459 0.05858 0.1398 1.35% MLMG::mgVcycle_down::0 65 0.1325 0.1333 0.1354 1.31% MLMG::mgVcycle_up::0 65 0.1068 0.1077 0.1089 1.05% Gravity::multilevel_solve_for_new_phi() 1 0.1072 0.1072 0.1072 1.03% Gravity::actual_multilevel_solve() 1 0.1071 0.1071 0.1072 1.03% Amr::writePlotFile() 2 0.1036 0.1037 0.1037 1.00% MLCellLinOp::apply() 501 0.07717 0.08355 0.1006 0.97% FabArray::ParallelCopy_nowait() 734 0.08257 0.08837 0.09791 0.94% Castro::post_regrid() 2 0.08679 0.08971 0.09611 0.93% MLMG::mgVcycle_up::3 65 0.005428 0.01194 0.0895 0.86% FillBoundary_finish() 3240 0.05688 0.06748 0.07907 0.76% Castro::construct_old_gravity() 11 0.03793 0.06046 0.07526 0.73% MLMG:computeResOfCorrection() 390 0.05795 0.06291 0.07472 0.72% FillBoundary_nowait() 3240 0.04711 0.05884 0.07463 0.72% MLCellLinOp::correctionResidual() 390 0.05725 0.06218 0.07396 0.71% Castro::derive() 18 0.07051 0.07145 0.07233 0.70% AmrLevel::derive() 18 0.07046 0.07141 0.07229 0.70% Gravity::get_old_grav_vector() 11 0.02825 0.0562 0.06837 0.66% Castro::post_timestep() 1 0.05923 0.06117 0.06466 0.62% MLMG::mgVcycle_down::1 65 0.05821 0.06016 0.06235 0.60% Castro::construct_new_source() 48 0.05814 0.05922 0.06064 0.58% MLPoisson::Fapply() 501 0.04796 0.05061 0.05858 0.57% Castro::initialize_advance() 1 0.04939 0.0509 0.05519 0.53% Amr::regrid() 1 0.04695 0.04853 0.05314 0.51% MLCellLinOp::solutionResidual() 75 0.04643 0.04934 0.05109 0.49% Castro::construct_old_source() 66 0.04411 0.04487 0.04561 0.44% Castro::create_source_corrector() 8 0.0162 0.03007 0.04172 0.40% MLMG::computeResidual() 65 0.04015 0.04045 0.04113 0.40% FabArray::Saxpy() 133 0.03593 0.03684 0.04073 0.39% MLMG::mgVcycle_bottom 65 0.0001837 0.00247 0.03666 0.35% MLMG::actualBottomSolve() 65 6.818e-05 0.002349 0.03653 0.35% MLMG::mgVcycle_down::2 65 0.006426 0.0207 0.03642 0.35% MLCGSolver::bicgstab 4 0 0.002233 0.03573 0.34% Castro::construct_new_gravity_source() 8 0.03235 0.03299 0.03405 0.33% FabArray::Xpay() 480 0.02622 0.02759 0.03153 0.30% VisMF::Write(FabArray) 4 0.01475 0.01646 0.03121 0.30% MLMG::mgVcycle_down::3 65 0.006108 0.009302 0.03059 0.30% Castro::retry_advance_ctu() 11 0.02234 0.02427 0.02949 0.28% amrex::average_down 390 0.01655 0.01995 0.02747 0.26% Castro::construct_old_gravity_source() 11 0.02566 0.02617 0.0272 0.26% Castro::check_for_nan() 19 0.02311 0.02367 0.02485 0.24% MultiFab::contains_nan() 19 0.02306 0.02361 0.02479 0.24% Castro::sum_integrated_quantities() 2 0.02115 0.02146 0.02148 0.21% MLPoisson::define() 10 0.01515 0.01763 0.02107 0.20% Gravity::get_new_grav_vector() 9 0.01755 0.01859 0.01934 0.19% MLMG::mgVcycle_down::4 65 0.006069 0.006946 0.01772 0.17% Castro::save_data_for_retry() 3 0.01106 0.01248 0.01744 0.17% Amr::writeSmallPlotFile() 2 0.0166 0.01661 0.01662 0.16% Castro::apply_source_to_state() 19 0.01488 0.01519 0.01635 0.16% Castro::react_state() 19 0.01507 0.01548 0.01607 0.16% Castro::problem_post_timestep() 2 0.00915 0.01152 0.01592 0.15% Castro::wd_update() 2 0.009059 0.01149 0.0159 0.15% Gravity::test_composite_phi() 1 0.01494 0.01498 0.015 0.14% StateDataPhysBCFunct::() 82 0.0003703 0.007531 0.01411 0.14% MLMG::mgVcycle_up::4 65 0.005511 0.006144 0.01386 0.13% StateData::FillBoundary(geom) 62 0 0.007105 0.01364 0.13% MLMG::ResNormInf() 74 0.007994 0.01083 0.0129 0.12% MLCellLinOp::defineAuxData() 10 0.00968 0.01056 0.01265 0.12% MLMG::prepareForSolve() 9 0.01002 0.01054 0.01184 0.11% MLMG::mgVcycle_down::5 65 0.006037 0.006543 0.01174 0.11% Gravity::update_max_rhs() 2 0.001396 0.007842 0.01078 0.10% MLMG::computeMLResidual() 9 0.005507 0.008081 0.01029 0.10% MLMG::mgVcycle_up::5 65 0.00555 0.005935 0.009971 0.10% MultiFab::max() 2 0.0001626 0.006592 0.009527 0.09% FabArray::mult() 40 0.007845 0.008093 0.00918 0.09% FabArray::norminf() 143 0.00563 0.005895 0.008033 0.08% amrex::Add() 65 0.007138 0.007279 0.007595 0.07% Other 5689 0.05525 0.06993 0.1034 1.00% -------------------------------------------------------------------------------------------- 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 15 B 22 B 29 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-26-ge6c93bf22695) finalized