Initializing AMReX (24.07-33-g80b0c7ea3dbc)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.07-33-g80b0c7ea3dbc) initialized Starting run at 21:30:58 UTC on 2024-07-26. Successfully read inputs file ... Castro git describe: 24.07-27-g48044aec5 AMReX git describe: 24.07-33-g80b0c7ea3 Microphysics git describe: 24.07-34-gecbc595e 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.178424997 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001072568 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010524964 Gravity::multilevel_solve_for_new_phi() time = 0.102084511 Gravity::fill_multipole_BCs() time = 0.010510245 ... 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.047234154 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.00376025 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.002141692 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.065723238 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.006379659 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.065499151 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.003880713 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.06209868 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.007886463 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.063298562 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.314627775 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015431994 Gravity::solve_for_phi() time = 0.087603311 on level 0 Castro::construct_new_gravity() time = 0.089741734 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.060673563 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.002353289 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.061238759 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.316763081 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015532047 Gravity::solve_for_phi() time = 0.090308827 on level 0 Castro::construct_new_gravity() time = 0.0924604 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.060586121 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.002337039 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.059839462 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.318330837 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015686035 Gravity::solve_for_phi() time = 0.093987233 on level 0 Castro::construct_new_gravity() time = 0.096066257 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.06192083 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.003502044 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.060116505 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.317889842 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016028222 Gravity::solve_for_phi() time = 0.089122371 on level 0 Castro::construct_new_gravity() time = 0.091200564 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.059836082 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.002205651 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.060871802 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.317811419 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016181904 Gravity::solve_for_phi() time = 0.087736347 on level 0 Castro::construct_new_gravity() time = 0.089795641 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.05999821 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.005375227 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.059923088 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.319362559 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017710499 Gravity::solve_for_phi() time = 0.094509606 on level 0 Castro::construct_new_gravity() time = 0.096577987 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.059855337 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.002511542 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.060676095 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.319951894 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017119364 Gravity::solve_for_phi() time = 0.092584331 on level 0 Castro::construct_new_gravity() time = 0.094648854 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.059645333 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.005129215 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.060285923 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.319269003 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01607149 Gravity::solve_for_phi() time = 0.08797293 on level 0 Castro::construct_new_gravity() time = 0.090040681 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.059422734 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.004654061325 [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.099569602 [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.048482183 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003747457 seconds Ending run at 21:31:08 UTC on 2024-07-26. Run time = 10.00834429 Run time without initialization = 7.152677391 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]: 10.01 ... 10.01 ... 10.01 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.7811 2.8509 2.9018 28.99% Amr::InitAmr() 1 1.971 1.971 1.971 19.69% Castro::computeTemp() 63 0.6981 0.7447 0.8163 8.16% amrex::Copy() 523 0.5327 0.544 0.5549 5.54% FabArray::ParallelCopy_finish() 734 0.03478 0.3451 0.5461 5.46% Castro::normalize_species() 62 0.449 0.4587 0.47 4.70% Castro::reset_internal_energy(Fab) 126 0.4217 0.4315 0.4385 4.38% MultiFab::AddProduct() 2633 0.4219 0.4312 0.4379 4.38% Castro::initMFs() 1 0.3503 0.3563 0.3624 3.62% MultiFab::Subtract() 125 0.2663 0.2713 0.2763 2.76% Castro::reset_internal_energy(MultiFab) 63 0.01779 0.06069 0.2481 2.48% FabArray::setVal() 3453 0.2193 0.2248 0.2309 2.31% MLPoisson::Fsmooth() 3120 0.157 0.1728 0.2114 2.11% check_for_negative_density() 11 0.0354 0.09922 0.1678 1.68% Castro::initData() 1 0.003296 0.02812 0.1664 1.66% Gravity::fill_multipole_BCs() 10 0.1025 0.1233 0.1484 1.48% FabArray::sum() 2633 0.1388 0.1415 0.1436 1.43% FillBoundary_finish() 3240 0.06641 0.07672 0.0887 0.89% FabArray::ParallelCopy_nowait() 734 0.06353 0.07103 0.08144 0.81% FillBoundary_nowait() 3240 0.04086 0.05291 0.06848 0.68% MLCellLinOp::applyBC() 3630 0.02756 0.0358 0.05985 0.60% MLPoisson::Fapply() 501 0.04806 0.05137 0.05965 0.60% Castro::enforce_min_density() 62 0.01369 0.02264 0.05707 0.57% AmrLevel::derive() 18 0.03984 0.04078 0.04191 0.42% FabArray::Saxpy() 133 0.03587 0.03691 0.04117 0.41% Castro::do_new_sources() 8 0.00017 0.006218 0.03802 0.38% Castro::construct_new_gravity_source() 8 0.03199 0.03283 0.03371 0.34% Castro::do_old_sources() 11 0.0003042 0.006893 0.03268 0.33% FabArray::Xpay() 480 0.02555 0.02707 0.03105 0.31% Castro::construct_old_gravity_source() 11 0.02491 0.02543 0.02597 0.26% MultiFab::contains_nan() 19 0.0226 0.02322 0.02521 0.25% VisMF::Write(FabArray) 4 0.01151 0.01264 0.02278 0.23% MLMG::addInterpCorrection() 390 0.01707 0.01859 0.02227 0.22% Amr::writePlotFile() 2 0.007076 0.01672 0.01898 0.19% amrex::average_down 390 0.01232 0.01353 0.01663 0.17% Castro::wd_update() 2 0.004943 0.008744 0.01484 0.15% StateData::FillBoundary(geom) 62 0 0.007062 0.01382 0.14% Castro::evaluate_source_change() 144 0.009433 0.009813 0.01045 0.10% MLCellLinOp::defineAuxData() 10 0.007674 0.008103 0.009334 0.09% FabArray::mult() 40 0.008209 0.008422 0.008847 0.09% FabArray::FillBoundary() 3240 0.006846 0.007179 0.008736 0.09% MLCGSolver::bicgstab 4 0 0.0005197 0.008314 0.08% Castro::construct_new_source() 48 0.007226 0.007525 0.008275 0.08% FabArray::norminf() 143 0.005624 0.005951 0.008091 0.08% FabArrayBase::getFB() 3240 0.005105 0.005385 0.00675 0.07% MLCellLinOp::smooth() 1560 0.005676 0.005957 0.006507 0.07% Castro::construct_old_source() 66 0.005354 0.005527 0.00611 0.06% MLMG::ResNormInf() 74 0.001791 0.003659 0.005849 0.06% FillPatchIterator::Initialize 82 0.005409 0.005549 0.0057 0.06% MLMG::prepareForSolve() 9 0.004801 0.004992 0.005248 0.05% Amr::writeSmallPlotFile() 2 0.003045 0.004679 0.005106 0.05% MLMG::solve() 9 0.0003935 0.001755 0.004216 0.04% Gravity::actual_solve_with_mlmg() 10 0.003422 0.003608 0.004143 0.04% FabArrayBase::getCPC() 1022 0.003712 0.003873 0.004076 0.04% Gravity::get_old_grav_vector() 11 0.003497 0.003611 0.003684 0.04% Gravity::get_new_grav_vector() 9 0.002991 0.003082 0.003234 0.03% Castro::initialize_do_advance() 11 0.002393 0.002625 0.003006 0.03% Castro::construct_old_gravity() 11 0.0002165 0.0008948 0.002887 0.03% MLCellLinOp::apply() 501 0.001016 0.001138 0.002448 0.02% Castro::sum_integrated_quantities() 2 0.001561 0.002013 0.002441 0.02% FabArray::ParallelCopy() 734 0.001808 0.001902 0.002093 0.02% AmrLevel::FillPatch() 82 0.001095 0.001135 0.001244 0.01% Castro::construct_new_gravity() 8 0.0001482 0.0003338 0.00115 0.01% MLCellLinOp::correctionResidual() 390 0.0009389 0.0009731 0.001016 0.01% main() 1 0.0006168 0.0006935 0.0009095 0.01% MLMG::mgVcycle() 65 0.0007718 0.0008102 0.0008942 0.01% MLMG:computeResOfCorrection() 390 0.0007034 0.000737 0.000796 0.01% Castro::save_data_for_retry() 3 0.000733 0.0007626 0.0007917 0.01% Gravity::solve_for_phi() 8 0.0003803 0.0004706 0.0006754 0.01% MLMG::actualBottomSolve() 65 6.741e-05 0.0001024 0.0005788 0.01% AmrMesh::MakeNewGrids() 2 0.000522 0.0005366 0.0005554 0.01% Castro::subcycle_advance_ctu() 1 9.217e-05 0.0001355 0.0005452 0.01% StateDataPhysBCFunct::() 82 0.0003498 0.0004159 0.0005032 0.01% Castro::Castro() 1 0.0003495 0.0003736 0.0004094 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003559 0.0003759 0.0003987 0.00% Castro::clean_state() 62 0.0003182 0.0003415 0.0003878 0.00% MLMG::mgVcycle_down::3 65 0.0002906 0.0003172 0.0003754 0.00% MLMG::mgVcycle_down::1 65 0.0003311 0.0003413 0.0003676 0.00% MLMG::mgVcycle_down::0 65 0.0003262 0.000337 0.0003623 0.00% MLMG::mgVcycle_down::2 65 0.0003075 0.0003387 0.0003622 0.00% MLMG::mgVcycle_down::4 65 0.0002793 0.000295 0.0003472 0.00% MLMG::mgVcycle_down::5 65 0.000288 0.0003016 0.0003418 0.00% FillPatchSingleLevel 82 0.0002518 0.0002612 0.0002709 0.00% Castro::do_advance_ctu() 11 0.0002205 0.0002407 0.0002699 0.00% MLMG::mgVcycle_up::0 65 0.000234 0.0002435 0.0002606 0.00% Castro::problem_post_timestep() 2 1.032e-05 2.755e-05 0.0002581 0.00% MLMG::mgVcycle_up::1 65 0.0002302 0.0002419 0.0002561 0.00% MLCellLinOp::solutionResidual() 75 0.0002252 0.0002322 0.0002479 0.00% MLMG::mgVcycle_up::2 65 0.0001937 0.0002164 0.0002429 0.00% MLMG::mgVcycle_up::3 65 0.0001905 0.0002029 0.0002412 0.00% MLMG::mgVcycle_up::5 65 0.0001987 0.0002079 0.0002315 0.00% MLMG::mgVcycle_up::4 65 0.0001933 0.0002045 0.0002282 0.00% MLMG::oneIter() 65 0.0001888 0.0001999 0.0002152 0.00% Amr::coarseTimeStep() 1 0.0001157 0.0001414 0.0001478 0.00% Amr::defBaseLevel() 1 0.0001087 0.0001228 0.0001424 0.00% Castro::retry_advance_ctu() 11 4.99e-05 6.165e-05 0.0001401 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001179 0.0001257 0.0001361 0.00% MLMG::computeResidual() 65 0.0001195 0.000124 0.0001348 0.00% MLMG::mgVcycle_bottom 65 0.0001113 0.0001162 0.0001241 0.00% TagBoxArray::mapPRD 2 0.0001083 0.0001154 0.0001233 0.00% Amr::grid_places() 2 7.857e-05 8.94e-05 0.0001177 0.00% Gravity::test_composite_phi() 1 9.087e-05 9.619e-05 0.0001134 0.00% Castro::initialize_advance() 1 9.893e-05 0.0001033 0.0001111 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.011e-05 6.804e-05 9.329e-05 0.00% MLPoisson::define() 10 6.593e-05 7.129e-05 7.694e-05 0.00% Castro::check_for_nan() 19 5.039e-05 5.563e-05 6.327e-05 0.00% Castro::apply_source_to_state() 19 5.106e-05 5.437e-05 6.049e-05 0.00% Castro::react_state() 19 4.971e-05 5.401e-05 5.867e-05 0.00% Amr::FinalizeInit() 1 8.47e-06 1.175e-05 5.523e-05 0.00% Castro::derive() 18 3.998e-05 4.22e-05 4.594e-05 0.00% Castro::create_source_corrector() 8 3.461e-05 3.832e-05 4.45e-05 0.00% Amr::timeStep() 1 1.221e-05 1.515e-05 4.024e-05 0.00% Castro::expand_state() 11 3.014e-05 3.285e-05 3.766e-05 0.00% Gravity::actual_multilevel_solve() 1 2.699e-05 2.795e-05 3.022e-05 0.00% Amr::regrid() 1 1.055e-05 1.262e-05 2.946e-05 0.00% MLMG::computeMLResidual() 9 1.929e-05 2.085e-05 2.225e-05 0.00% Amr::InitializeInit() 1 6.945e-06 8.547e-06 1.541e-05 0.00% Amr::init() 1 6.413e-06 1.009e-05 1.494e-05 0.00% Castro::post_init() 1 1.178e-05 1.254e-05 1.359e-05 0.00% Castro::post_regrid() 2 1.052e-05 1.127e-05 1.258e-05 0.00% Castro::advance() 1 6.773e-06 8.584e-06 1.208e-05 0.00% Amr::bldFineLevels() 1 6.912e-06 7.85e-06 1.188e-05 0.00% Castro::post_timestep() 1 5.447e-06 6.018e-06 6.823e-06 0.00% Amr::initialInit() 1 2.886e-06 3.987e-06 6.53e-06 0.00% Other 1496 0.03493 0.04452 0.06411 0.64% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.01 10.01 10.01 100.00% Amr::coarseTimeStep() 1 7.152 7.152 7.152 71.46% Amr::timeStep() 1 7.1 7.1 7.1 70.94% Castro::advance() 1 7.034 7.038 7.041 70.35% Castro::subcycle_advance_ctu() 1 6.985 6.989 6.992 69.86% Castro::do_advance_ctu() 11 6.962 6.965 6.969 69.63% Castro::construct_ctu_hydro_source() 11 3.298 3.319 3.34 33.37% Castro::clean_state() 62 2.908 3.016 3.117 31.15% Castro::computeTemp() 63 1.875 1.985 2.157 21.55% Amr::InitAmr() 1 1.971 1.971 1.971 19.69% Castro::reset_internal_energy(MultiFab) 63 1.175 1.24 1.407 14.06% Amr::init() 1 0.8848 0.8848 0.8848 8.84% Gravity::solve_phi_with_mlmg() 10 0.7901 0.8112 0.8371 8.36% Amr::initialInit() 1 0.8334 0.8334 0.8335 8.33% Castro::initialize_do_advance() 11 0.754 0.7742 0.8035 8.03% Castro::evaluate_source_change() 144 0.7178 0.7328 0.744 7.43% Castro::construct_new_gravity() 8 0.6933 0.7144 0.7397 7.39% Gravity::solve_for_phi() 8 0.6761 0.6972 0.723 7.22% Gravity::actual_solve_with_mlmg() 10 0.6838 0.6861 0.6886 6.88% Castro::do_old_sources() 11 0.6378 0.66 0.6746 6.74% MLMG::solve() 9 0.6563 0.6587 0.6596 6.59% Castro::enforce_min_density() 62 0.6107 0.6311 0.6539 6.53% FabArray::ParallelCopy() 734 0.1267 0.4345 0.633 6.32% Amr::InitializeInit() 1 0.4311 0.4949 0.608 6.07% Amr::defBaseLevel() 1 0.4311 0.4949 0.608 6.07% MLMG::oneIter() 65 0.5859 0.5883 0.5911 5.91% MLMG::mgVcycle() 65 0.5787 0.581 0.5837 5.83% amrex::Copy() 523 0.5327 0.544 0.5549 5.54% FabArray::ParallelCopy_finish() 734 0.03478 0.3451 0.5461 5.46% Castro::do_new_sources() 8 0.4341 0.4626 0.4819 4.82% Castro::normalize_species() 62 0.449 0.4587 0.47 4.70% Castro::reset_internal_energy(Fab) 126 0.4217 0.4315 0.4385 4.38% MultiFab::AddProduct() 2633 0.4219 0.4312 0.4379 4.38% Amr::FinalizeInit() 1 0.2254 0.3385 0.4023 4.02% MLCellLinOp::smooth() 1560 0.2947 0.3226 0.3906 3.90% Castro::Castro() 1 0.3541 0.3602 0.3663 3.66% Castro::initMFs() 1 0.3504 0.3564 0.3625 3.62% AmrLevel::FillPatch() 82 0.2005 0.2599 0.3541 3.54% FillPatchIterator::Initialize 82 0.1289 0.1892 0.2849 2.85% FillPatchIterator::FillFromLevel0() 82 0.118 0.1809 0.2774 2.77% FillPatchSingleLevel 82 0.1177 0.1805 0.277 2.77% MultiFab::Subtract() 125 0.2663 0.2713 0.2763 2.76% Castro::initData() 1 0.0757 0.1335 0.2521 2.52% FabArray::setVal() 3453 0.2193 0.2248 0.2309 2.31% MLMG::mgVcycle_up::1 65 0.05122 0.1386 0.2266 2.26% Amr::bldFineLevels() 1 0.04658 0.1591 0.2217 2.22% MLCellLinOp::applyBC() 3630 0.1597 0.1785 0.2164 2.16% MLPoisson::Fsmooth() 3120 0.157 0.1728 0.2114 2.11% MLMG::addInterpCorrection() 390 0.03308 0.1602 0.1974 1.97% Amr::grid_places() 2 0.002677 0.1157 0.1795 1.79% AmrMesh::MakeNewGrids() 2 0.002596 0.1157 0.1794 1.79% TagBoxArray::mapPRD 2 0.0008112 0.1138 0.1779 1.78% check_for_negative_density() 11 0.0354 0.09922 0.1678 1.68% FabArray::FillBoundary() 3240 0.1276 0.1427 0.1591 1.59% Castro::expand_state() 11 0.1 0.1214 0.1553 1.55% Gravity::fill_multipole_BCs() 10 0.1037 0.1244 0.1495 1.49% MLMG::mgVcycle_up::2 65 0.005443 0.06145 0.1471 1.47% FabArray::sum() 2633 0.1388 0.1415 0.1436 1.43% MLMG::mgVcycle_down::0 65 0.1354 0.1374 0.1413 1.41% Castro::post_init() 1 0.1347 0.1361 0.138 1.38% MLMG::mgVcycle_up::0 65 0.1031 0.1051 0.1064 1.06% MLCellLinOp::apply() 501 0.08027 0.0862 0.1051 1.05% Gravity::multilevel_solve_for_new_phi() 1 0.102 0.1021 0.1021 1.02% Gravity::actual_multilevel_solve() 1 0.102 0.102 0.102 1.02% FabArray::ParallelCopy_nowait() 734 0.07969 0.0875 0.0987 0.99% Amr::writePlotFile() 2 0.09601 0.09612 0.09613 0.96% MLMG::mgVcycle_up::3 65 0.005407 0.01226 0.09396 0.94% FillBoundary_finish() 3240 0.06641 0.07672 0.0887 0.89% Castro::post_regrid() 2 0.0847 0.08664 0.08807 0.88% MLMG:computeResOfCorrection() 390 0.05915 0.0648 0.07717 0.77% MLCellLinOp::correctionResidual() 390 0.05844 0.06407 0.07637 0.76% FillBoundary_nowait() 3240 0.04663 0.05876 0.07465 0.75% Castro::derive() 18 0.06774 0.06926 0.07055 0.70% AmrLevel::derive() 18 0.0677 0.06921 0.07051 0.70% MLMG::mgVcycle_down::1 65 0.06426 0.06637 0.06972 0.70% Castro::post_timestep() 1 0.0588 0.06153 0.06598 0.66% Castro::construct_new_source() 48 0.05787 0.05931 0.06094 0.61% MLPoisson::Fapply() 501 0.04806 0.05137 0.05965 0.60% MLCellLinOp::solutionResidual() 75 0.04674 0.04956 0.05193 0.52% Castro::create_source_corrector() 8 0.01557 0.02899 0.05099 0.51% Castro::initialize_advance() 1 0.04779 0.04888 0.04971 0.50% Amr::regrid() 1 0.04325 0.04444 0.04549 0.45% Castro::construct_old_source() 66 0.04308 0.04403 0.04541 0.45% FabArray::Saxpy() 133 0.03587 0.03691 0.04117 0.41% MLMG::computeResidual() 65 0.03937 0.04038 0.04103 0.41% Castro::construct_old_gravity() 11 0.0231 0.02995 0.03915 0.39% Gravity::get_old_grav_vector() 11 0.02184 0.02905 0.03883 0.39% MLMG::mgVcycle_down::2 65 0.006389 0.0213 0.0379 0.38% MLMG::mgVcycle_bottom 65 0.0001799 0.002528 0.03766 0.38% MLMG::actualBottomSolve() 65 6.741e-05 0.002412 0.03754 0.38% MLCGSolver::bicgstab 4 0 0.002294 0.03671 0.37% Castro::construct_new_gravity_source() 8 0.03199 0.03283 0.03371 0.34% MLMG::mgVcycle_down::3 65 0.00609 0.009434 0.03151 0.31% FabArray::Xpay() 480 0.02555 0.02707 0.03105 0.31% amrex::average_down 390 0.01654 0.02015 0.02857 0.29% Castro::construct_old_gravity_source() 11 0.02491 0.02543 0.02597 0.26% Castro::retry_advance_ctu() 11 0.02183 0.02264 0.02575 0.26% Castro::check_for_nan() 19 0.02265 0.02327 0.02526 0.25% MultiFab::contains_nan() 19 0.0226 0.02322 0.02521 0.25% VisMF::Write(FabArray) 4 0.01151 0.01264 0.02278 0.23% MLPoisson::define() 10 0.01592 0.01836 0.02172 0.22% Castro::sum_integrated_quantities() 2 0.02101 0.02115 0.02116 0.21% MLMG::mgVcycle_down::4 65 0.00604 0.006983 0.01848 0.18% Gravity::get_new_grav_vector() 9 0.01741 0.01782 0.01846 0.18% Castro::problem_post_timestep() 2 0.009196 0.01266 0.01835 0.18% Castro::wd_update() 2 0.009172 0.01263 0.01834 0.18% Castro::apply_source_to_state() 19 0.01482 0.01518 0.01603 0.16% Castro::react_state() 19 0.01493 0.01534 0.01599 0.16% MLMG::mgVcycle_up::4 65 0.005484 0.006183 0.01447 0.14% StateDataPhysBCFunct::() 82 0.0003498 0.007478 0.01426 0.14% Gravity::test_composite_phi() 1 0.01406 0.01409 0.01411 0.14% Castro::save_data_for_retry() 3 0.01084 0.0113 0.01402 0.14% StateData::FillBoundary(geom) 62 0 0.007062 0.01382 0.14% MLCellLinOp::defineAuxData() 10 0.009881 0.01075 0.01312 0.13% MLMG::mgVcycle_down::5 65 0.006007 0.006563 0.01235 0.12% MLMG::prepareForSolve() 9 0.009835 0.01045 0.0116 0.12% MLMG::computeMLResidual() 9 0.006145 0.008559 0.0111 0.11% MLMG::ResNormInf() 74 0.00685 0.008636 0.01074 0.11% MLMG::mgVcycle_up::5 65 0.00555 0.005991 0.01038 0.10% FabArray::mult() 40 0.008209 0.008422 0.008847 0.09% FabArray::norminf() 143 0.005624 0.005951 0.008091 0.08% FabArrayBase::getCPC() 1022 0.005156 0.005914 0.007868 0.08% Amr::writeSmallPlotFile() 2 0.00774 0.007744 0.007751 0.08% FabArrayBase::getFB() 3240 0.005362 0.005848 0.007644 0.08% Other 1496 0.05301 0.06941 0.09784 0.98% -------------------------------------------------------------------------------------------- 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 13 B 18 B 21 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.07-33-g80b0c7ea3dbc) finalized