Initializing AMReX (24.05-18-g8eff86d32d17)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05-18-g8eff86d32d17) initialized Starting run at 07:10:36 UTC on 2024-05-15. Successfully read inputs file ... Castro git describe: 24.05-8-g16ad4929f AMReX git describe: 24.05-18-g8eff86d32 Microphysics git describe: 24.05-5-g87839b75 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.182941991 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001072262 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010697307 Gravity::multilevel_solve_for_new_phi() time = 0.099716911 Gravity::fill_multipole_BCs() time = 0.010603854 ... 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.05022421 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007298264 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.002245799 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.067574928 on level 0 ... Entering construct_ctu_hydro_source() on level 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 = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 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.004455238 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.067293641 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.004712677 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.063633177 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.012836028 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.093733622 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.335734398 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.021333684 Gravity::solve_for_phi() time = 0.112290716 on level 0 Castro::construct_new_gravity() time = 0.1151336 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.061174043 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.004383004 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.067657381 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.342176676 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015940696 Gravity::solve_for_phi() time = 0.086178716 on level 0 Castro::construct_new_gravity() time = 0.088249989 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.060500659 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.004427352 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.082152271 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.347068433 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01598075 Gravity::solve_for_phi() time = 0.08724061 on level 0 Castro::construct_new_gravity() time = 0.089331187 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.06104745 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.007083752 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.064653588 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.33965126 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.03101454 Gravity::solve_for_phi() time = 0.111769675 on level 0 Castro::construct_new_gravity() time = 0.114645806 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.061473974 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.00364421 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.079680475 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.339597029 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.032587765 Gravity::solve_for_phi() time = 0.105132894 on level 0 Castro::construct_new_gravity() time = 0.107283066 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.061651699 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.003417421 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.068957902 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.371603727 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017454663 Gravity::solve_for_phi() time = 0.090861141 on level 0 Castro::construct_new_gravity() time = 0.094145995 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.092268798 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.003175752 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.064351975 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.357483255 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.035996535 Gravity::solve_for_phi() time = 0.110379286 on level 0 Castro::construct_new_gravity() time = 0.112418957 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.061596921 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.035143345 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.066362195 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.381277484 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.02341826 Gravity::solve_for_phi() time = 0.097238971 on level 0 Castro::construct_new_gravity() time = 0.09940976 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.061745673 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.004331112106 [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.646870762 [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.053253405 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.008486 seconds Ending run at 07:10:46 UTC on 2024-05-15. Run time = 10.01293703 Run time without initialization = 7.709704407 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.01 ... 10.01 ... 10.01 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8485 2.9100 3.0692 30.65% Amr::InitAmr() 1 1.398 1.398 1.398 13.97% Castro::computeTemp() 63 0.7095 0.7681 0.8612 8.60% FabArray::ParallelCopy_finish() 734 0.04208 0.5079 0.7442 7.43% amrex::Copy() 523 0.5479 0.5603 0.6159 6.15% Castro::normalize_species() 62 0.4426 0.457 0.4987 4.98% Castro::reset_internal_energy(Fab) 126 0.4297 0.4406 0.4728 4.72% MultiFab::AddProduct() 2633 0.4302 0.4381 0.466 4.65% Castro::initMFs() 1 0.3589 0.365 0.389 3.89% check_for_negative_density() 11 0.05907 0.2394 0.36 3.59% MultiFab::Subtract() 125 0.2611 0.266 0.2885 2.88% Castro::reset_internal_energy(MultiFab) 63 0.01788 0.07093 0.2605 2.60% FabArray::setVal() 3453 0.2189 0.2253 0.2415 2.41% Gravity::fill_multipole_BCs() 10 0.1157 0.1828 0.2132 2.13% MLPoisson::Fsmooth() 3120 0.1589 0.1752 0.2116 2.11% Castro::initData() 1 0.003282 0.029 0.1713 1.71% FabArray::sum() 2633 0.1467 0.1493 0.1595 1.59% Castro::enforce_min_density() 62 0.01389 0.03401 0.1303 1.30% FillBoundary_finish() 3240 0.07799 0.09081 0.1038 1.04% FabArray::ParallelCopy_nowait() 734 0.06525 0.0714 0.08059 0.80% FillBoundary_nowait() 3240 0.04034 0.0527 0.06748 0.67% Castro::do_old_sources() 11 0.0002986 0.01355 0.06596 0.66% MLPoisson::Fapply() 501 0.0478 0.05055 0.05845 0.58% MLCellLinOp::applyBC() 3630 0.02431 0.03198 0.05347 0.53% AmrLevel::derive() 18 0.04112 0.04235 0.04657 0.47% FabArray::Saxpy() 133 0.03588 0.0369 0.04109 0.41% Castro::construct_new_gravity_source() 8 0.03191 0.03342 0.0391 0.39% Castro::do_new_sources() 8 0.0001767 0.006466 0.03393 0.34% FabArray::Xpay() 480 0.02607 0.02758 0.03132 0.31% Castro::construct_old_gravity_source() 11 0.02552 0.02619 0.0306 0.31% VisMF::Write(FabArray) 4 0.01399 0.01589 0.02821 0.28% Amr::writePlotFile() 2 0.009752 0.02195 0.0252 0.25% MultiFab::contains_nan() 19 0.02314 0.02351 0.02466 0.25% Castro::sum_integrated_quantities() 2 0.003262 0.02221 0.02386 0.24% MLMG::addInterpCorrection() 390 0.01685 0.0183 0.02191 0.22% Castro::wd_update() 2 0.005183 0.01177 0.01881 0.19% Castro::volWgtSum() 46 0.00729 0.008038 0.01672 0.17% amrex::average_down 390 0.01227 0.01341 0.01669 0.17% StateData::FillBoundary(geom) 62 0 0.007403 0.01445 0.14% Amr::writeSmallPlotFile() 2 0.007581 0.01045 0.01092 0.11% Castro::evaluate_source_change() 144 0.009411 0.009723 0.01053 0.11% FabArray::mult() 40 0.007827 0.008086 0.009348 0.09% MLCellLinOp::defineAuxData() 10 0.007485 0.007895 0.009198 0.09% Castro::construct_new_source() 48 0.007242 0.007581 0.008863 0.09% MLMG::ResNormInf() 74 0.005416 0.007151 0.008789 0.09% FabArray::FillBoundary() 3240 0.006856 0.007171 0.008747 0.09% MLCGSolver::bicgstab 4 0 0.0005261 0.008418 0.08% FabArray::norminf() 143 0.005611 0.005945 0.00794 0.08% amrex::Add() 65 0.007106 0.007236 0.007737 0.08% Castro::construct_old_source() 66 0.005437 0.0057 0.006959 0.07% MLCellLinOp::smooth() 1560 0.005713 0.005929 0.006408 0.06% FillPatchIterator::Initialize 82 0.005541 0.005857 0.006173 0.06% Castro::construct_old_gravity() 11 0.0002229 0.002053 0.005753 0.06% MLMG::prepareForSolve() 9 0.00499 0.005177 0.005745 0.06% Gravity::actual_solve_with_mlmg() 10 0.003467 0.003674 0.004192 0.04% FabArrayBase::getCPC() 1022 0.003646 0.003829 0.004093 0.04% Gravity::get_old_grav_vector() 11 0.003507 0.003606 0.003741 0.04% MLMG::solve() 9 0.0004203 0.001964 0.003712 0.04% Gravity::get_new_grav_vector() 9 0.00294 0.003082 0.0034 0.03% Castro::initialize_do_advance() 11 0.002532 0.002856 0.003398 0.03% Castro::construct_new_gravity() 8 0.0001421 0.0006108 0.002657 0.03% MLCellLinOp::apply() 501 0.001015 0.00113 0.002439 0.02% FabArray::ParallelCopy() 734 0.001779 0.001865 0.002059 0.02% Gravity::solve_for_phi() 8 0.0003689 0.0006564 0.001904 0.02% main() 1 0.0008977 0.001145 0.001389 0.01% AmrLevel::FillPatch() 82 0.00108 0.00116 0.001285 0.01% MLCellLinOp::correctionResidual() 390 0.0009203 0.000959 0.001013 0.01% MLMG::mgVcycle() 65 0.0007864 0.0008259 0.0009175 0.01% Castro::save_data_for_retry() 3 0.0007488 0.0007881 0.0008451 0.01% MLMG:computeResOfCorrection() 390 0.0006914 0.0007191 0.0007796 0.01% Castro::subcycle_advance_ctu() 1 9.681e-05 0.0001522 0.0006484 0.01% Castro::Castro() 1 0.0004456 0.0005098 0.0005827 0.01% MLMG::actualBottomSolve() 65 6.708e-05 0.0001002 0.0005588 0.01% AmrMesh::MakeNewGrids() 2 0.0004884 0.0005156 0.000539 0.01% StateDataPhysBCFunct::() 82 0.000352 0.0004225 0.0004889 0.00% Castro::clean_state() 62 0.0003276 0.0003618 0.0004263 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003623 0.0003901 0.0004245 0.00% Castro::problem_post_timestep() 2 1.034e-05 3.623e-05 0.0004081 0.00% MLMG::mgVcycle_down::0 65 0.0003297 0.0003431 0.0003856 0.00% MLMG::mgVcycle_down::2 65 0.000313 0.0003404 0.0003717 0.00% MLMG::mgVcycle_down::3 65 0.000293 0.0003164 0.0003629 0.00% MLMG::mgVcycle_down::1 65 0.0003292 0.000338 0.0003518 0.00% MLMG::mgVcycle_down::4 65 0.0002869 0.0002985 0.0003461 0.00% MLMG::mgVcycle_down::5 65 0.0002972 0.0003066 0.0003405 0.00% Castro::do_advance_ctu() 11 0.0002211 0.0002603 0.0003235 0.00% Amr::coarseTimeStep() 1 0.0002315 0.0002847 0.0003008 0.00% FillPatchSingleLevel 82 0.0002555 0.0002694 0.0002901 0.00% MLMG::mgVcycle_up::0 65 0.0002306 0.0002403 0.0002699 0.00% MLMG::mgVcycle_up::1 65 0.0002298 0.0002399 0.0002524 0.00% MLMG::mgVcycle_up::2 65 0.0001974 0.0002169 0.0002495 0.00% MLCellLinOp::solutionResidual() 75 0.0002199 0.0002304 0.000243 0.00% MLMG::mgVcycle_up::5 65 0.0002014 0.000208 0.0002315 0.00% MLMG::mgVcycle_up::3 65 0.000196 0.000204 0.0002254 0.00% MLMG::mgVcycle_up::4 65 0.0001987 0.0002042 0.000225 0.00% MLMG::oneIter() 65 0.0001953 0.0002033 0.0002116 0.00% Castro::initialize_advance() 1 0.0001119 0.0001457 0.0001713 0.00% Amr::defBaseLevel() 1 9.057e-05 0.0001303 0.0001573 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001135 0.0001255 0.0001466 0.00% Castro::retry_advance_ctu() 11 5.043e-05 6.525e-05 0.0001398 0.00% MLMG::mgVcycle_bottom 65 0.0001171 0.0001208 0.0001313 0.00% MLMG::computeResidual() 65 0.0001191 0.0001228 0.0001271 0.00% TagBoxArray::mapPRD 2 0.0001094 0.0001145 0.0001229 0.00% Gravity::test_composite_phi() 1 9.279e-05 9.857e-05 0.0001196 0.00% Amr::grid_places() 2 7.368e-05 8.851e-05 0.0001172 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.271e-05 7.142e-05 0.0001018 0.00% Castro::check_for_nan() 19 4.94e-05 5.978e-05 0.0001015 0.00% MLPoisson::define() 10 6.53e-05 7.186e-05 9.155e-05 0.00% Castro::apply_source_to_state() 19 4.9e-05 5.693e-05 7.327e-05 0.00% Castro::react_state() 19 4.984e-05 5.653e-05 6.327e-05 0.00% Amr::FinalizeInit() 1 8.188e-06 1.35e-05 6.132e-05 0.00% Castro::create_source_corrector() 8 3.511e-05 4.22e-05 5.339e-05 0.00% Castro::derive() 18 3.974e-05 4.266e-05 5.026e-05 0.00% Castro::expand_state() 11 3.26e-05 3.765e-05 4.327e-05 0.00% Amr::timeStep() 1 1.219e-05 1.836e-05 3.287e-05 0.00% Gravity::actual_multilevel_solve() 1 2.72e-05 2.887e-05 3.065e-05 0.00% Amr::regrid() 1 1.015e-05 1.25e-05 2.982e-05 0.00% MLMG::computeMLResidual() 9 2.025e-05 2.245e-05 2.496e-05 0.00% Castro::advance() 1 6.63e-06 1.137e-05 1.834e-05 0.00% Castro::post_regrid() 2 1.062e-05 1.197e-05 1.458e-05 0.00% Amr::InitializeInit() 1 6.31e-06 7.58e-06 1.445e-05 0.00% Castro::post_init() 1 1.184e-05 1.273e-05 1.437e-05 0.00% Amr::init() 1 6.858e-06 8.005e-06 9.746e-06 0.00% Amr::bldFineLevels() 1 6.905e-06 7.878e-06 9.736e-06 0.00% Castro::post_timestep() 1 5.675e-06 6.784e-06 9e-06 0.00% Amr::initialInit() 1 2.973e-06 3.303e-06 3.655e-06 0.00% Other 4625 0.02552 0.03535 0.0586 0.59% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.01 10.01 10.01 100.00% Amr::coarseTimeStep() 1 7.709 7.709 7.709 76.99% Amr::timeStep() 1 7.647 7.647 7.647 76.37% Castro::advance() 1 7.555 7.561 7.567 75.57% Castro::subcycle_advance_ctu() 1 7.507 7.511 7.516 75.06% Castro::do_advance_ctu() 11 7.483 7.487 7.493 74.83% Castro::construct_ctu_hydro_source() 11 3.489 3.528 3.595 35.90% Castro::clean_state() 62 2.946 3.086 3.395 33.91% Castro::computeTemp() 63 1.91 2.041 2.284 22.81% Castro::reset_internal_energy(MultiFab) 63 1.193 1.272 1.47 14.68% Amr::InitAmr() 1 1.398 1.398 1.398 13.97% Castro::initialize_do_advance() 11 0.8044 0.8723 0.9182 9.17% Gravity::solve_phi_with_mlmg() 10 0.813 0.8814 0.9122 9.11% Amr::init() 1 0.9043 0.9044 0.9044 9.03% Amr::initialInit() 1 0.846 0.8461 0.8462 8.45% FabArray::ParallelCopy() 734 0.1346 0.598 0.8316 8.31% Castro::construct_new_gravity() 8 0.723 0.7895 0.8198 8.19% Gravity::solve_for_phi() 8 0.7028 0.77 0.8006 8.00% Castro::evaluate_source_change() 144 0.7337 0.747 0.7948 7.94% Castro::do_old_sources() 11 0.7036 0.7514 0.7834 7.82% Castro::enforce_min_density() 62 0.6196 0.6483 0.7525 7.52% FabArray::ParallelCopy_finish() 734 0.04208 0.5079 0.7442 7.43% Gravity::actual_solve_with_mlmg() 10 0.6947 0.6968 0.6994 6.98% MLMG::solve() 9 0.6675 0.6701 0.6709 6.70% Amr::InitializeInit() 1 0.4404 0.5065 0.6218 6.21% Amr::defBaseLevel() 1 0.4404 0.5065 0.6218 6.21% amrex::Copy() 523 0.5479 0.5603 0.6159 6.15% MLMG::oneIter() 65 0.5941 0.5953 0.5969 5.96% MLMG::mgVcycle() 65 0.5867 0.5879 0.5894 5.89% AmrLevel::FillPatch() 82 0.2111 0.4287 0.5556 5.55% Castro::do_new_sources() 8 0.4375 0.4758 0.521 5.20% Castro::normalize_species() 62 0.4426 0.457 0.4987 4.98% FillPatchIterator::Initialize 82 0.1369 0.3551 0.4845 4.84% FillPatchIterator::FillFromLevel0() 82 0.1257 0.3466 0.4769 4.76% FillPatchSingleLevel 82 0.1253 0.3462 0.4766 4.76% Castro::reset_internal_energy(Fab) 126 0.4297 0.4406 0.4728 4.72% MultiFab::AddProduct() 2633 0.4302 0.4381 0.466 4.65% Amr::FinalizeInit() 1 0.2244 0.3396 0.4057 4.05% MLCellLinOp::smooth() 1560 0.3056 0.3331 0.397 3.96% Castro::Castro() 1 0.3628 0.3693 0.3938 3.93% Castro::initMFs() 1 0.359 0.3652 0.3891 3.89% check_for_negative_density() 11 0.05907 0.2394 0.36 3.59% MultiFab::Subtract() 125 0.2611 0.266 0.2885 2.88% Castro::initData() 1 0.07627 0.1359 0.2565 2.56% FabArray::setVal() 3453 0.2189 0.2253 0.2415 2.41% Amr::bldFineLevels() 1 0.04689 0.1616 0.227 2.27% MLCellLinOp::applyBC() 3630 0.1714 0.1884 0.2249 2.25% MLMG::mgVcycle_up::1 65 0.05092 0.1358 0.2203 2.20% Castro::expand_state() 11 0.1026 0.1623 0.2152 2.15% Gravity::fill_multipole_BCs() 10 0.1168 0.1839 0.2143 2.14% MLPoisson::Fsmooth() 3120 0.1589 0.1752 0.2116 2.11% MLMG::addInterpCorrection() 390 0.03394 0.1559 0.1918 1.92% Amr::grid_places() 2 0.002634 0.1179 0.184 1.84% AmrMesh::MakeNewGrids() 2 0.002517 0.1178 0.184 1.84% TagBoxArray::mapPRD 2 0.0008113 0.1162 0.1824 1.82% FabArray::FillBoundary() 3240 0.1408 0.1564 0.1772 1.77% FabArray::sum() 2633 0.1467 0.1493 0.1595 1.59% MLMG::mgVcycle_down::0 65 0.1465 0.1491 0.1534 1.53% MLMG::mgVcycle_up::2 65 0.00545 0.05997 0.1433 1.43% Castro::post_init() 1 0.1324 0.1345 0.1362 1.36% MLMG::mgVcycle_up::0 65 0.107 0.1081 0.11 1.10% Castro::create_source_corrector() 8 0.01648 0.07043 0.1077 1.08% MLCellLinOp::apply() 501 0.08237 0.08719 0.1059 1.06% Amr::writePlotFile() 2 0.104 0.1041 0.1041 1.04% FillBoundary_finish() 3240 0.07799 0.09081 0.1038 1.04% Gravity::multilevel_solve_for_new_phi() 1 0.09969 0.09971 0.09974 1.00% Gravity::actual_multilevel_solve() 1 0.09963 0.09964 0.09966 1.00% FabArray::ParallelCopy_nowait() 734 0.08179 0.0883 0.09821 0.98% Castro::post_timestep() 1 0.08007 0.08577 0.09138 0.91% MLMG::mgVcycle_up::3 65 0.005427 0.01199 0.09057 0.90% Castro::post_regrid() 2 0.08535 0.08708 0.08925 0.89% Castro::derive() 18 0.06952 0.07187 0.08271 0.83% AmrLevel::derive() 18 0.06948 0.07183 0.08266 0.83% Castro::construct_old_gravity() 11 0.02916 0.04516 0.08044 0.80% Gravity::get_old_grav_vector() 11 0.02556 0.04311 0.08018 0.80% MLMG:computeResOfCorrection() 390 0.06014 0.0659 0.07896 0.79% MLCellLinOp::correctionResidual() 390 0.05945 0.06519 0.07818 0.78% FillBoundary_nowait() 3240 0.04574 0.05843 0.07357 0.73% Castro::construct_new_source() 48 0.05766 0.06 0.06847 0.68% MLMG::mgVcycle_down::1 65 0.06139 0.06449 0.06808 0.68% MLPoisson::Fapply() 501 0.0478 0.05055 0.05845 0.58% Castro::construct_old_source() 66 0.04388 0.04501 0.05297 0.53% MLCellLinOp::solutionResidual() 75 0.04648 0.04997 0.05221 0.52% Castro::initialize_advance() 1 0.04853 0.04987 0.05175 0.52% Amr::regrid() 1 0.04381 0.04474 0.04607 0.46% Castro::sum_integrated_quantities() 2 0.04243 0.04272 0.04275 0.43% MLMG::computeResidual() 65 0.03899 0.04027 0.0411 0.41% FabArray::Saxpy() 133 0.03588 0.0369 0.04109 0.41% Castro::construct_new_gravity_source() 8 0.03191 0.03342 0.0391 0.39% MLMG::mgVcycle_down::2 65 0.006398 0.02065 0.03684 0.37% MLMG::mgVcycle_bottom 65 0.0001846 0.002475 0.03676 0.37% MLMG::actualBottomSolve() 65 6.708e-05 0.002355 0.03663 0.37% MLCGSolver::bicgstab 4 0 0.00224 0.03583 0.36% FabArray::Xpay() 480 0.02607 0.02758 0.03132 0.31% Castro::construct_old_gravity_source() 11 0.02552 0.02619 0.0306 0.31% MLMG::mgVcycle_down::3 65 0.006067 0.009208 0.03048 0.30% VisMF::Write(FabArray) 4 0.01399 0.01589 0.02821 0.28% amrex::average_down 390 0.01642 0.01985 0.02802 0.28% Castro::check_for_nan() 19 0.02319 0.02357 0.02471 0.25% MultiFab::contains_nan() 19 0.02314 0.02351 0.02466 0.25% Castro::retry_advance_ctu() 11 0.02196 0.02261 0.02352 0.23% Castro::problem_post_timestep() 2 0.009424 0.0154 0.02212 0.22% Castro::wd_update() 2 0.009016 0.01536 0.02211 0.22% Gravity::get_new_grav_vector() 9 0.01834 0.01978 0.02078 0.21% MLPoisson::define() 10 0.01504 0.01757 0.02077 0.21% MLMG::mgVcycle_down::4 65 0.006047 0.006875 0.01772 0.18% Castro::volWgtSum() 46 0.00729 0.008038 0.01672 0.17% Castro::react_state() 19 0.0148 0.01551 0.01669 0.17% Amr::writeSmallPlotFile() 2 0.01616 0.01617 0.01617 0.16% Castro::apply_source_to_state() 19 0.01487 0.01517 0.01609 0.16% StateDataPhysBCFunct::() 82 0.000352 0.007825 0.01492 0.15% StateData::FillBoundary(geom) 62 0 0.007403 0.01445 0.14% Gravity::test_composite_phi() 1 0.01406 0.01408 0.0141 0.14% MLMG::mgVcycle_up::4 65 0.005488 0.006098 0.01398 0.14% MLMG::ResNormInf() 74 0.01066 0.01215 0.01376 0.14% MLCellLinOp::defineAuxData() 10 0.009707 0.01053 0.01293 0.13% MLMG::prepareForSolve() 9 0.01006 0.01065 0.01214 0.12% MLMG::mgVcycle_down::5 65 0.006025 0.006497 0.01187 0.12% Castro::save_data_for_retry() 3 0.01096 0.01122 0.01146 0.11% MLMG::computeMLResidual() 9 0.00657 0.009076 0.01108 0.11% MLMG::mgVcycle_up::5 65 0.005533 0.005885 0.009971 0.10% FabArray::mult() 40 0.007827 0.008086 0.009348 0.09% FabArray::norminf() 143 0.005611 0.005945 0.00794 0.08% amrex::Add() 65 0.007106 0.007236 0.007737 0.08% FabArrayBase::getCPC() 1022 0.005015 0.005758 0.007669 0.08% Other 4625 0.04395 0.06072 0.09562 0.95% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 16 B 21 B 29 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-18-g8eff86d32d17) finalized