Initializing AMReX (24.07-15-gcca5aade6929)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.07-15-gcca5aade6929) initialized Starting run at 16:53:03 UTC on 2024-07-10. Successfully read inputs file ... Castro git describe: 24.07-15-g6a9b6e7c5 AMReX git describe: 24.07-15-gcca5aade6 Microphysics git describe: 24.07-10-g0b1a2310 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.193281076 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001064728 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010694598 Gravity::multilevel_solve_for_new_phi() time = 0.106113907 Gravity::fill_multipole_BCs() time = 0.010678254 ... 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.050410784 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007149305 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.003629899 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.068845073 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 = -9.947567811e-05 at index 0, 62, 0 Invalid density = Invalid density = -8.076810125 at index 20, 81, 0 -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.012287294 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.06986836 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.006512339 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.066143829 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.008240156 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.064311194 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.316183773 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017466296 Gravity::solve_for_phi() time = 0.089882767 on level 0 Castro::construct_new_gravity() time = 0.091941047 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.060909159 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.0032938 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.064934811 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.326410511 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015727235 Gravity::solve_for_phi() time = 0.086368832 on level 0 Castro::construct_new_gravity() time = 0.088424932 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.061082733 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.003133475 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.061851211 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.331978523 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018461387 Gravity::solve_for_phi() time = 0.094157821 on level 0 Castro::construct_new_gravity() time = 0.096213314 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.060141634 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.002106439 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.064650836 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.328878647 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018812765 Gravity::solve_for_phi() time = 0.093181557 on level 0 Castro::construct_new_gravity() time = 0.095242421 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.061282881 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.00311719 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.06086865 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.325562739 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015981668 Gravity::solve_for_phi() time = 0.090794745 on level 0 Castro::construct_new_gravity() time = 0.093749433 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.05985469 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.002381442 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.062702972 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.31779682 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015432807 Gravity::solve_for_phi() time = 0.085506859 on level 0 Castro::construct_new_gravity() time = 0.087558625 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.061239071 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.003370753 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.064333793 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.321583497 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016489085 Gravity::solve_for_phi() time = 0.086952728 on level 0 Castro::construct_new_gravity() time = 0.089062208 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.060466116 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.002187354 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.060849683 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.316200445 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016088264 Gravity::solve_for_phi() time = 0.086131513 on level 0 Castro::construct_new_gravity() time = 0.088185792 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.060355222 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.004568400263 [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.232275524 [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.05022821 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006840689 seconds Ending run at 16:53:13 UTC on 2024-07-10. Run time = 9.627375266 Run time without initialization = 7.290324494 Average number of zones advanced per microsecond: 0.004 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 9.627 ... 9.627 ... 9.628 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8278 2.8727 2.9677 30.83% Amr::InitAmr() 1 1.406 1.406 1.407 14.61% Castro::computeTemp() 63 0.6903 0.7523 0.8427 8.75% FabArray::ParallelCopy_finish() 734 0.05094 0.3804 0.5937 6.17% Castro::normalize_species() 62 0.4601 0.4719 0.5811 6.04% amrex::Copy() 523 0.5466 0.5532 0.5804 6.03% MultiFab::AddProduct() 2633 0.4333 0.444 0.5437 5.65% Castro::reset_internal_energy(Fab) 126 0.4298 0.4353 0.4539 4.71% Castro::initMFs() 1 0.3627 0.3668 0.3805 3.95% Castro::reset_internal_energy(MultiFab) 63 0.0179 0.08255 0.2694 2.80% MultiFab::Subtract() 125 0.2552 0.2574 0.2639 2.74% FabArray::setVal() 3453 0.224 0.2278 0.2342 2.43% MLPoisson::Fsmooth() 3120 0.1592 0.1729 0.2103 2.18% check_for_negative_density() 11 0.01891 0.1302 0.1916 1.99% Castro::initData() 1 0.003417 0.02922 0.1754 1.82% Gravity::fill_multipole_BCs() 10 0.104 0.1298 0.154 1.60% FabArray::sum() 2633 0.146 0.1474 0.1535 1.59% FillBoundary_finish() 3240 0.05615 0.06901 0.08379 0.87% Castro::enforce_min_density() 62 0.01383 0.02808 0.08138 0.85% FabArray::ParallelCopy_nowait() 734 0.06505 0.07077 0.07772 0.81% FillBoundary_nowait() 3240 0.04031 0.05295 0.06655 0.69% MLPoisson::Fapply() 501 0.04806 0.05046 0.05809 0.60% MLCellLinOp::applyBC() 3630 0.0246 0.03246 0.0543 0.56% AmrLevel::derive() 18 0.04074 0.0414 0.04254 0.44% FabArray::Saxpy() 133 0.03474 0.03539 0.03997 0.42% Castro::construct_new_gravity_source() 8 0.03205 0.0329 0.03678 0.38% Castro::do_new_sources() 8 0.0001705 0.006242 0.03638 0.38% Castro::do_old_sources() 11 0.0002971 0.007766 0.03495 0.36% FabArray::Xpay() 480 0.02672 0.02808 0.0318 0.33% VisMF::Write(FabArray) 4 0.01389 0.01582 0.02757 0.29% Castro::construct_old_gravity_source() 11 0.0251 0.02544 0.02658 0.28% MultiFab::contains_nan() 19 0.02317 0.02356 0.02646 0.27% Amr::writePlotFile() 2 0.008839 0.01976 0.02229 0.23% MLMG::addInterpCorrection() 390 0.01678 0.018 0.02156 0.22% amrex::average_down 390 0.01248 0.01355 0.01667 0.17% Castro::wd_update() 2 0.004905 0.009767 0.01444 0.15% StateData::FillBoundary(geom) 62 0 0.007181 0.01417 0.15% FabArray::mult() 40 0.008204 0.008504 0.01105 0.11% Castro::evaluate_source_change() 144 0.009407 0.009635 0.01034 0.11% Amr::writeSmallPlotFile() 2 0.006974 0.009173 0.009625 0.10% MLCellLinOp::defineAuxData() 10 0.007713 0.008162 0.009525 0.10% FabArray::FillBoundary() 3240 0.006895 0.007152 0.008607 0.09% MLCGSolver::bicgstab 4 0 0.00052 0.00832 0.09% FabArray::norminf() 143 0.005438 0.005761 0.007828 0.08% amrex::Add() 65 0.007096 0.007211 0.007766 0.08% Castro::construct_old_gravity() 11 0.0002141 0.002176 0.007713 0.08% Castro::construct_new_source() 48 0.007273 0.007323 0.007501 0.08% MLMG::ResNormInf() 74 0.002978 0.005166 0.006787 0.07% FabArrayBase::getFB() 3240 0.005108 0.005302 0.006627 0.07% MLCellLinOp::smooth() 1560 0.005539 0.005758 0.006243 0.06% FillPatchIterator::Initialize 82 0.005607 0.005941 0.006207 0.06% Castro::construct_old_source() 66 0.005367 0.005459 0.005653 0.06% MLMG::prepareForSolve() 9 0.005043 0.005194 0.005561 0.06% MLMG::solve() 9 0.0005365 0.003218 0.004669 0.05% Gravity::get_new_grav_vector() 9 0.002986 0.003163 0.004155 0.04% FabArrayBase::getCPC() 1022 0.003766 0.003904 0.004091 0.04% Gravity::actual_solve_with_mlmg() 10 0.00344 0.003584 0.004043 0.04% Gravity::get_old_grav_vector() 11 0.003537 0.003607 0.003684 0.04% Castro::initialize_do_advance() 11 0.002509 0.002865 0.003221 0.03% MLCellLinOp::apply() 501 0.0009986 0.001101 0.002361 0.02% Castro::sum_integrated_quantities() 2 0.001635 0.001979 0.002167 0.02% FabArray::ParallelCopy() 734 0.001813 0.001884 0.002042 0.02% Castro::construct_new_gravity() 8 0.000136 0.0005781 0.001807 0.02% main() 1 0.0009307 0.001224 0.001621 0.02% AmrLevel::FillPatch() 82 0.001085 0.001156 0.001243 0.01% MLCellLinOp::correctionResidual() 390 0.0009282 0.0009555 0.001005 0.01% MLMG::mgVcycle() 65 0.0007577 0.0007956 0.0008791 0.01% Castro::save_data_for_retry() 3 0.0007812 0.0008148 0.0008586 0.01% MLMG:computeResOfCorrection() 390 0.0006983 0.0007204 0.0007701 0.01% Gravity::solve_for_phi() 8 0.0003691 0.0004485 0.0007038 0.01% Castro::subcycle_advance_ctu() 1 9.419e-05 0.0001565 0.0006883 0.01% Castro::Castro() 1 0.0005054 0.000553 0.0006193 0.01% AmrMesh::MakeNewGrids() 2 0.0005223 0.0005416 0.0005612 0.01% MLMG::actualBottomSolve() 65 6.757e-05 9.939e-05 0.0005454 0.01% StateDataPhysBCFunct::() 82 0.0003456 0.0004248 0.0004763 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003707 0.0003893 0.0004092 0.00% Castro::clean_state() 62 0.0003082 0.0003375 0.0003901 0.00% MLMG::mgVcycle_down::2 65 0.0003135 0.0003393 0.0003635 0.00% MLMG::mgVcycle_down::3 65 0.0002883 0.0003118 0.0003629 0.00% MLMG::mgVcycle_down::1 65 0.0003314 0.0003412 0.0003574 0.00% MLMG::mgVcycle_down::5 65 0.0002882 0.0002991 0.0003491 0.00% MLMG::mgVcycle_down::0 65 0.000332 0.0003398 0.0003489 0.00% MLMG::mgVcycle_down::4 65 0.0002789 0.0002906 0.0003479 0.00% FillPatchSingleLevel 82 0.0002619 0.0002724 0.0002876 0.00% Castro::problem_post_timestep() 2 1.134e-05 2.869e-05 0.0002687 0.00% Castro::do_advance_ctu() 11 0.0002156 0.0002426 0.000265 0.00% MLMG::mgVcycle_up::0 65 0.0002359 0.0002429 0.0002501 0.00% MLMG::mgVcycle_up::1 65 0.0002276 0.0002398 0.0002495 0.00% MLCellLinOp::solutionResidual() 75 0.0002261 0.0002321 0.0002456 0.00% MLMG::mgVcycle_up::2 65 0.0001949 0.0002163 0.0002392 0.00% MLMG::mgVcycle_up::5 65 0.0002026 0.0002092 0.0002358 0.00% MLMG::mgVcycle_up::3 65 0.0001926 0.0002011 0.0002291 0.00% MLMG::mgVcycle_up::4 65 0.0001951 0.0002019 0.0002217 0.00% MLMG::oneIter() 65 0.0001913 0.0001958 0.0002035 0.00% Gravity::update_max_rhs() 2 0.0001582 0.0001679 0.000182 0.00% Amr::coarseTimeStep() 1 0.0001434 0.0001747 0.0001811 0.00% Castro::initialize_advance() 1 0.0001391 0.0001567 0.0001752 0.00% Amr::defBaseLevel() 1 9.594e-05 0.0001351 0.0001652 0.00% Castro::retry_advance_ctu() 11 5.342e-05 6.501e-05 0.0001541 0.00% MLMG::computeResidual() 65 0.0001187 0.0001235 0.0001457 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001109 0.0001191 0.0001319 0.00% Amr::grid_places() 2 8.278e-05 9.773e-05 0.0001294 0.00% MLMG::mgVcycle_bottom 65 0.0001162 0.0001189 0.0001279 0.00% TagBoxArray::mapPRD 2 0.0001118 0.0001173 0.0001273 0.00% Gravity::test_composite_phi() 1 9.205e-05 9.665e-05 0.000113 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.082e-05 7.012e-05 0.000112 0.00% MLPoisson::define() 10 6.423e-05 7.117e-05 9.098e-05 0.00% Castro::check_for_nan() 19 5.077e-05 5.776e-05 7.673e-05 0.00% Castro::apply_source_to_state() 19 4.973e-05 5.697e-05 7.146e-05 0.00% Castro::react_state() 19 5.102e-05 5.559e-05 6.242e-05 0.00% Amr::FinalizeInit() 1 1.2e-05 1.853e-05 6.08e-05 0.00% Castro::derive() 18 4.006e-05 4.195e-05 4.541e-05 0.00% Castro::expand_state() 11 3.425e-05 3.838e-05 4.482e-05 0.00% Castro::create_source_corrector() 8 3.268e-05 3.856e-05 4.429e-05 0.00% Amr::timeStep() 1 1.208e-05 1.6e-05 3.375e-05 0.00% Amr::regrid() 1 1.046e-05 1.256e-05 3.099e-05 0.00% Gravity::actual_multilevel_solve() 1 2.717e-05 2.843e-05 2.991e-05 0.00% MLMG::computeMLResidual() 9 1.989e-05 2.133e-05 2.276e-05 0.00% Amr::InitializeInit() 1 7.88e-06 1.159e-05 1.676e-05 0.00% Amr::init() 1 6.528e-06 8.172e-06 1.604e-05 0.00% Castro::advance() 1 7.086e-06 9.962e-06 1.536e-05 0.00% Castro::post_regrid() 2 1.07e-05 1.253e-05 1.51e-05 0.00% Castro::post_init() 1 1.115e-05 1.211e-05 1.358e-05 0.00% Amr::bldFineLevels() 1 6.68e-06 7.951e-06 9.748e-06 0.00% Amr::initialInit() 1 2.799e-06 3.586e-06 8.189e-06 0.00% Castro::post_timestep() 1 5.535e-06 6.045e-06 6.52e-06 0.00% Other 1429 0.02803 0.03949 0.0605 0.63% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.627 9.627 9.627 100.00% Amr::coarseTimeStep() 1 7.29 7.29 7.29 75.72% Amr::timeStep() 1 7.232 7.232 7.232 75.12% Castro::advance() 1 7.166 7.17 7.173 74.50% Castro::subcycle_advance_ctu() 1 7.115 7.118 7.121 73.97% Castro::do_advance_ctu() 11 7.086 7.093 7.097 73.72% Castro::construct_ctu_hydro_source() 11 3.356 3.375 3.401 35.33% Castro::clean_state() 62 2.959 3.074 3.195 33.18% Castro::computeTemp() 63 1.908 2.025 2.227 23.13% Castro::reset_internal_energy(MultiFab) 63 1.195 1.273 1.454 15.11% Amr::InitAmr() 1 1.406 1.406 1.407 14.61% Amr::init() 1 0.9295 0.9301 0.9302 9.66% Amr::initialInit() 1 0.8711 0.8717 0.8719 9.06% Castro::evaluate_source_change() 144 0.7383 0.752 0.8518 8.85% Gravity::solve_phi_with_mlmg() 10 0.7779 0.8049 0.8313 8.63% Castro::initialize_do_advance() 11 0.7834 0.7992 0.8263 8.58% Castro::construct_new_gravity() 8 0.6788 0.7049 0.7296 7.58% Gravity::solve_for_phi() 8 0.6596 0.6867 0.713 7.41% Castro::do_old_sources() 11 0.6776 0.6875 0.6982 7.25% Castro::enforce_min_density() 62 0.6173 0.6406 0.6852 7.12% FabArray::ParallelCopy() 734 0.1424 0.4697 0.6808 7.07% Gravity::actual_solve_with_mlmg() 10 0.6713 0.6733 0.6756 7.02% MLMG::solve() 9 0.6438 0.6463 0.6472 6.72% Amr::InitializeInit() 1 0.4464 0.5142 0.6381 6.63% Amr::defBaseLevel() 1 0.4464 0.5141 0.6381 6.63% FabArray::ParallelCopy_finish() 734 0.05094 0.3804 0.5937 6.17% Castro::normalize_species() 62 0.4601 0.4719 0.5811 6.04% amrex::Copy() 523 0.5466 0.5532 0.5804 6.03% MLMG::oneIter() 65 0.5738 0.5749 0.5763 5.99% MLMG::mgVcycle() 65 0.5664 0.5675 0.5689 5.91% MultiFab::AddProduct() 2633 0.4333 0.444 0.5437 5.65% Castro::do_new_sources() 8 0.4409 0.4659 0.4851 5.04% Castro::reset_internal_energy(Fab) 126 0.4298 0.4353 0.4539 4.71% Amr::FinalizeInit() 1 0.2337 0.3576 0.4254 4.42% AmrLevel::FillPatch() 82 0.2099 0.2898 0.3931 4.08% Castro::Castro() 1 0.3672 0.3715 0.3849 4.00% Castro::initMFs() 1 0.3628 0.3669 0.3806 3.95% MLCellLinOp::smooth() 1560 0.2843 0.3123 0.3778 3.92% FillPatchIterator::Initialize 82 0.1368 0.218 0.3223 3.35% FillPatchIterator::FillFromLevel0() 82 0.1266 0.2092 0.3144 3.27% FillPatchSingleLevel 82 0.1262 0.2088 0.314 3.26% MultiFab::Subtract() 125 0.2552 0.2574 0.2639 2.74% Castro::initData() 1 0.07695 0.1407 0.2632 2.73% Amr::bldFineLevels() 1 0.04862 0.1712 0.2373 2.46% FabArray::setVal() 3453 0.224 0.2278 0.2342 2.43% MLMG::mgVcycle_up::1 65 0.04817 0.1335 0.2194 2.28% MLPoisson::Fsmooth() 3120 0.1592 0.1729 0.2103 2.18% MLCellLinOp::applyBC() 3630 0.1487 0.1673 0.2059 2.14% Amr::grid_places() 2 0.002645 0.1265 0.1944 2.02% AmrMesh::MakeNewGrids() 2 0.002561 0.1265 0.1943 2.02% MLMG::addInterpCorrection() 390 0.03176 0.1559 0.193 2.00% TagBoxArray::mapPRD 2 0.0008613 0.1244 0.1928 2.00% check_for_negative_density() 11 0.01891 0.1302 0.1916 1.99% Castro::expand_state() 11 0.1028 0.1279 0.1655 1.72% FabArray::FillBoundary() 3240 0.1179 0.1349 0.1561 1.62% Gravity::fill_multipole_BCs() 10 0.1051 0.1309 0.155 1.61% FabArray::sum() 2633 0.146 0.1474 0.1535 1.59% Castro::post_init() 1 0.1386 0.1418 0.1454 1.51% MLMG::mgVcycle_down::0 65 0.1331 0.1367 0.142 1.47% MLMG::mgVcycle_up::2 65 0.005519 0.05945 0.1417 1.47% Gravity::multilevel_solve_for_new_phi() 1 0.1061 0.1061 0.1061 1.10% Gravity::actual_multilevel_solve() 1 0.106 0.106 0.1061 1.10% MLCellLinOp::apply() 501 0.07695 0.08444 0.1043 1.08% Amr::writePlotFile() 2 0.1012 0.1013 0.1013 1.05% MLMG::mgVcycle_up::0 65 0.09838 0.09876 0.09894 1.03% FabArray::ParallelCopy_nowait() 734 0.08153 0.08743 0.0947 0.98% Castro::post_regrid() 2 0.08563 0.08918 0.09236 0.96% MLMG::mgVcycle_up::3 65 0.00545 0.01203 0.0911 0.95% FillBoundary_finish() 3240 0.05615 0.06901 0.08379 0.87% MLMG:computeResOfCorrection() 390 0.05969 0.06597 0.07965 0.83% MLCellLinOp::correctionResidual() 390 0.05899 0.06525 0.07888 0.82% FillBoundary_nowait() 3240 0.0457 0.0587 0.07259 0.75% MLMG::mgVcycle_down::1 65 0.06348 0.06715 0.0714 0.74% Castro::derive() 18 0.06918 0.07028 0.07138 0.74% AmrLevel::derive() 18 0.06914 0.07023 0.07133 0.74% Castro::post_timestep() 1 0.05946 0.06261 0.06606 0.69% Castro::construct_new_source() 48 0.05764 0.05872 0.06243 0.65% MLPoisson::Fapply() 501 0.04806 0.05046 0.05809 0.60% Castro::initialize_advance() 1 0.04913 0.05114 0.05309 0.55% Castro::create_source_corrector() 8 0.01768 0.03178 0.05219 0.54% MLCellLinOp::solutionResidual() 75 0.04487 0.04763 0.05167 0.54% Amr::regrid() 1 0.04398 0.04575 0.04736 0.49% Castro::construct_old_gravity() 11 0.02238 0.03457 0.04661 0.48% Gravity::get_old_grav_vector() 11 0.02216 0.0324 0.04637 0.48% Castro::construct_old_source() 66 0.04305 0.04357 0.04458 0.46% MLMG::computeResidual() 65 0.03864 0.03922 0.04028 0.42% FabArray::Saxpy() 133 0.03474 0.03539 0.03997 0.42% MLMG::mgVcycle_down::2 65 0.006481 0.02208 0.03881 0.40% MLMG::mgVcycle_bottom 65 0.0001849 0.002487 0.03697 0.38% MLMG::actualBottomSolve() 65 6.757e-05 0.002368 0.03684 0.38% Castro::construct_new_gravity_source() 8 0.03205 0.0329 0.03678 0.38% MLCGSolver::bicgstab 4 0 0.002254 0.03606 0.37% FabArray::Xpay() 480 0.02672 0.02808 0.0318 0.33% MLMG::mgVcycle_down::3 65 0.006076 0.009248 0.03043 0.32% Castro::retry_advance_ctu() 11 0.02275 0.02367 0.03043 0.32% amrex::average_down 390 0.01672 0.02027 0.02875 0.30% VisMF::Write(FabArray) 4 0.01389 0.01582 0.02757 0.29% Castro::construct_old_gravity_source() 11 0.0251 0.02544 0.02658 0.28% Castro::check_for_nan() 19 0.02322 0.02362 0.02653 0.28% MultiFab::contains_nan() 19 0.02317 0.02356 0.02646 0.27% Castro::sum_integrated_quantities() 2 0.02082 0.02098 0.021 0.22% MLPoisson::define() 10 0.01548 0.01799 0.02079 0.22% Gravity::get_new_grav_vector() 9 0.01737 0.01858 0.01953 0.20% Castro::save_data_for_retry() 3 0.01134 0.0119 0.01893 0.20% Castro::problem_post_timestep() 2 0.009027 0.01351 0.01788 0.19% Castro::wd_update() 2 0.008909 0.01348 0.01787 0.19% MLMG::mgVcycle_down::4 65 0.006053 0.006883 0.01787 0.19% Castro::react_state() 19 0.01519 0.01572 0.01672 0.17% StateDataPhysBCFunct::() 82 0.0003456 0.007606 0.01465 0.15% Castro::apply_source_to_state() 19 0.01436 0.01449 0.01462 0.15% Amr::writeSmallPlotFile() 2 0.0144 0.0144 0.01441 0.15% Gravity::test_composite_phi() 1 0.01422 0.01425 0.01427 0.15% StateData::FillBoundary(geom) 62 0 0.007181 0.01417 0.15% MLMG::mgVcycle_up::4 65 0.005523 0.006104 0.0139 0.14% MLCellLinOp::defineAuxData() 10 0.009935 0.01081 0.0133 0.14% MLMG::mgVcycle_down::5 65 0.006038 0.006478 0.01189 0.12% MLMG::prepareForSolve() 9 0.01014 0.01057 0.01186 0.12% MLMG::ResNormInf() 74 0.008211 0.009994 0.01154 0.12% FabArray::mult() 40 0.008204 0.008504 0.01105 0.11% MLMG::computeMLResidual() 9 0.005446 0.007796 0.01069 0.11% MLMG::mgVcycle_up::5 65 0.005557 0.005907 0.01002 0.10% Gravity::update_max_rhs() 2 0.00153 0.004738 0.008298 0.09% FabArrayBase::getCPC() 1022 0.00521 0.005933 0.007844 0.08% FabArray::norminf() 143 0.005438 0.005761 0.007828 0.08% amrex::Add() 65 0.007096 0.007211 0.007766 0.08% FabArrayBase::getFB() 3240 0.005366 0.005754 0.0075 0.08% Other 1429 0.04491 0.06139 0.08935 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 17 B 21 B 31 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.07-15-gcca5aade6929) finalized