Initializing AMReX (24.05-32-g28b010126a1b)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05-32-g28b010126a1b) initialized Starting run at 07:40:57 UTC on 2024-05-28. Successfully read inputs file ... Castro git describe: 24.05-15-g7893b5c6d AMReX git describe: 24.05-32-g28b010126 Microphysics git describe: 24.05-22-g7b1458be 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.20433678 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001086262 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010399183 Gravity::multilevel_solve_for_new_phi() time = 0.097807925 Gravity::fill_multipole_BCs() time = 0.010612149 ... 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.049663293 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007018494 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.003534833 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.067062859 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.003403652 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.063358767 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.02612720466 at index 21, 173, 0 Invalid density = -1.592544331 at index 20, 176, 0 Invalid density = -1.592544331 at index 20, 79, 0 Invalid density = -0.02612720466 at index 21, 82, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.125 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.0625   Beginning subcycle 1 starting at time 0 with dt = 0.0625 Estimated number of subcycles remaining: 4 Castro::construct_old_gravity() time = 0.003532705 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.059531955 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -0.1060744788 at index 20, 177, 0 Invalid density = -0.1060744788 at index 20, 78, 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.002722872 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.059213369 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.308072884 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014415169 Gravity::solve_for_phi() time = 0.084515228 on level 0 Castro::construct_new_gravity() time = 0.086711589 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.060082053 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.014760925 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.095495985 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.355697558 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015023873 Gravity::solve_for_phi() time = 0.092164518 on level 0 Castro::construct_new_gravity() time = 0.094275351 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.060111538 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.002234115 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.095083706 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.321199113 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015656613 Gravity::solve_for_phi() time = 0.086500213 on level 0 Castro::construct_new_gravity() time = 0.088627298 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.060458585 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.002207838 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.060973388 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.321805024 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016364697 Gravity::solve_for_phi() time = 0.088098037 on level 0 Castro::construct_new_gravity() time = 0.090224079 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.059248685 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.002116814 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.058694365 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.315047528 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015304384 Gravity::solve_for_phi() time = 0.084914763 on level 0 Castro::construct_new_gravity() time = 0.087045537 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.0598689 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.002241846 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.061543541 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.321757318 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01641057 Gravity::solve_for_phi() time = 0.086502028 on level 0 Castro::construct_new_gravity() time = 0.088640632 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.060409299 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.002153638 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.060947591 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.323149506 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016607837 Gravity::solve_for_phi() time = 0.08654885 on level 0 Castro::construct_new_gravity() time = 0.088673034 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.060653453 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.004244786 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.06154023 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.323467754 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016385083 Gravity::solve_for_phi() time = 0.087848695 on level 0 Castro::construct_new_gravity() time = 0.089960598 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.060201177 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.004567442316 [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.233078439 [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.049196915 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006777564 seconds Ending run at 07:41:07 UTC on 2024-05-28. Run time = 9.534240169 Run time without initialization = 7.289952323 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.534 ... 9.534 ... 9.534 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8501 2.8828 2.9446 30.88% Amr::InitAmr() 1 1.318 1.318 1.318 13.82% Castro::computeTemp() 63 0.7024 0.7582 0.8458 8.87% FabArray::ParallelCopy_finish() 734 0.0584 0.4051 0.6179 6.48% amrex::Copy() 523 0.5402 0.5466 0.5601 5.87% Castro::normalize_species() 62 0.4424 0.4483 0.4645 4.87% Castro::reset_internal_energy(Fab) 126 0.4302 0.4359 0.4463 4.68% MultiFab::AddProduct() 2633 0.4246 0.4288 0.4368 4.58% Castro::initMFs() 1 0.3624 0.3684 0.3985 4.18% MultiFab::Subtract() 125 0.2721 0.2747 0.2809 2.95% Castro::reset_internal_energy(MultiFab) 63 0.01803 0.07414 0.2732 2.87% FabArray::setVal() 3453 0.2215 0.2256 0.2342 2.46% MLPoisson::Fsmooth() 3120 0.1602 0.1735 0.208 2.18% check_for_negative_density() 11 0.06307 0.1346 0.1824 1.91% Castro::initData() 1 0.00327 0.02849 0.1689 1.77% FabArray::sum() 2633 0.1417 0.1432 0.1457 1.53% Gravity::fill_multipole_BCs() 10 0.1047 0.1241 0.1443 1.51% Castro::enforce_min_density() 62 0.0138 0.03175 0.08409 0.88% FillBoundary_finish() 3240 0.05325 0.0617 0.08165 0.86% FabArray::ParallelCopy_nowait() 734 0.06557 0.07223 0.08111 0.85% FillBoundary_nowait() 3240 0.04084 0.05269 0.06766 0.71% MLPoisson::Fapply() 501 0.04804 0.05039 0.05745 0.60% MLCellLinOp::applyBC() 3630 0.02419 0.03178 0.05241 0.55% Castro::do_old_sources() 11 0.0002982 0.00985 0.04951 0.52% AmrLevel::derive() 18 0.0415 0.04189 0.0425 0.45% TagBoxArray::collate() 2 0.0001832 0.005299 0.03961 0.42% FabArray::Saxpy() 133 0.03495 0.03557 0.03884 0.41% Castro::construct_new_gravity_source() 8 0.03268 0.03308 0.03376 0.35% Castro::do_new_sources() 8 0.0001684 0.006216 0.03219 0.34% FabArray::Xpay() 480 0.02577 0.02698 0.03029 0.32% Castro::construct_old_gravity_source() 11 0.02522 0.02594 0.02919 0.31% VisMF::Write(FabArray) 4 0.01472 0.01625 0.02761 0.29% MultiFab::contains_nan() 19 0.02322 0.02368 0.02729 0.29% MLMG::addInterpCorrection() 390 0.01695 0.01821 0.02157 0.23% Amr::writePlotFile() 2 0.008257 0.01761 0.01982 0.21% Castro::construct_old_gravity() 11 0.0002254 0.005284 0.01713 0.18% amrex::average_down 390 0.01244 0.01341 0.01625 0.17% StateData::FillBoundary(geom) 62 0 0.007183 0.01406 0.15% Castro::wd_update() 2 0.005141 0.008177 0.01155 0.12% FabArray::mult() 40 0.008171 0.008435 0.01132 0.12% Castro::evaluate_source_change() 144 0.00914 0.009336 0.009713 0.10% Amr::writeSmallPlotFile() 2 0.007195 0.008814 0.009218 0.10% MLCellLinOp::defineAuxData() 10 0.00765 0.00795 0.008968 0.09% FabArray::FillBoundary() 3240 0.00684 0.007125 0.008466 0.09% MLCGSolver::bicgstab 4 0 0.0005007 0.008011 0.08% Castro::construct_new_source() 48 0.007292 0.007443 0.007947 0.08% FabArray::norminf() 143 0.005603 0.005891 0.007808 0.08% amrex::Add() 65 0.00721 0.007317 0.007728 0.08% Castro::volWgtSum() 46 0.007323 0.007374 0.007415 0.08% FillPatchIterator::Initialize 82 0.005761 0.006016 0.006278 0.07% MLCellLinOp::smooth() 1560 0.00556 0.005783 0.006225 0.07% Castro::construct_old_source() 66 0.005405 0.005591 0.006006 0.06% MLMG::prepareForSolve() 9 0.004982 0.005127 0.00536 0.06% MLMG::ResNormInf() 74 0.001864 0.003111 0.004232 0.04% Gravity::actual_solve_with_mlmg() 10 0.00348 0.003624 0.004073 0.04% Gravity::get_old_grav_vector() 11 0.00353 0.003636 0.003771 0.04% Castro::initialize_do_advance() 11 0.002666 0.003021 0.003559 0.04% Gravity::get_new_grav_vector() 9 0.003056 0.003116 0.003284 0.03% Castro::sum_integrated_quantities() 2 0.00172 0.002755 0.002963 0.03% MLCellLinOp::apply() 501 0.001004 0.001105 0.002326 0.02% MLMG::solve() 9 0.0005008 0.001327 0.002265 0.02% FabArray::ParallelCopy() 734 0.001777 0.001837 0.001968 0.02% main() 1 0.0009673 0.00119 0.001288 0.01% AmrLevel::FillPatch() 82 0.001113 0.001154 0.001227 0.01% MLCellLinOp::correctionResidual() 390 0.0009174 0.0009429 0.0009871 0.01% MLMG::mgVcycle() 65 0.0008017 0.0008361 0.0009052 0.01% Castro::save_data_for_retry() 3 0.0007956 0.0008366 0.0008706 0.01% MLMG:computeResOfCorrection() 390 0.0006987 0.0007233 0.0007615 0.01% Gravity::solve_for_phi() 8 0.0003827 0.0004687 0.0006551 0.01% Castro::Castro() 1 0.0005051 0.0005636 0.0006549 0.01% Castro::subcycle_advance_ctu() 1 0.0001094 0.0001631 0.0006483 0.01% AmrMesh::MakeNewGrids() 2 0.0005016 0.0005202 0.0005422 0.01% MLMG::actualBottomSolve() 65 6.778e-05 9.847e-05 0.0005296 0.01% StateDataPhysBCFunct::() 82 0.0003465 0.0004113 0.0004683 0.00% Castro::construct_new_gravity() 8 0.0001366 0.0002554 0.0004339 0.00% Castro::problem_post_timestep() 2 1.128e-05 3.66e-05 0.0003996 0.00% FillPatchIterator::FillFromLevel0() 82 0.000363 0.0003779 0.0003953 0.00% Castro::clean_state() 62 0.0003229 0.0003382 0.0003844 0.00% MLMG::mgVcycle_down::2 65 0.000307 0.0003376 0.0003729 0.00% MLMG::mgVcycle_down::3 65 0.0002867 0.0003108 0.0003656 0.00% MLMG::mgVcycle_down::0 65 0.0003282 0.0003404 0.0003486 0.00% MLMG::mgVcycle_down::1 65 0.0003262 0.000336 0.0003483 0.00% MLMG::mgVcycle_down::4 65 0.0002792 0.0002919 0.0003432 0.00% MLMG::mgVcycle_down::5 65 0.0002864 0.0002976 0.0003289 0.00% FillPatchSingleLevel 82 0.0002547 0.0002635 0.0002862 0.00% Castro::do_advance_ctu() 11 0.0002255 0.0002446 0.0002725 0.00% MLMG::mgVcycle_up::1 65 0.0002277 0.0002399 0.0002594 0.00% MLMG::mgVcycle_up::0 65 0.0002289 0.0002363 0.0002479 0.00% MLMG::mgVcycle_up::2 65 0.0001914 0.0002127 0.0002376 0.00% MLMG::mgVcycle_up::3 65 0.0001906 0.0002007 0.000237 0.00% MLCellLinOp::solutionResidual() 75 0.0002203 0.0002248 0.0002346 0.00% MLMG::mgVcycle_up::5 65 0.0001962 0.0002035 0.0002258 0.00% MLMG::mgVcycle_up::4 65 0.0001942 0.0002004 0.0002174 0.00% MLMG::oneIter() 65 0.0001949 0.0001991 0.0002045 0.00% Castro::retry_advance_ctu() 11 6.131e-05 7.479e-05 0.000176 0.00% Castro::initialize_advance() 1 0.0001441 0.000157 0.0001722 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001132 0.0001261 0.0001585 0.00% Amr::coarseTimeStep() 1 0.0001101 0.0001327 0.0001386 0.00% Amr::defBaseLevel() 1 0.000102 0.0001149 0.0001306 0.00% Gravity::test_composite_phi() 1 9.261e-05 9.984e-05 0.0001289 0.00% MLMG::computeResidual() 65 0.0001195 0.0001225 0.0001272 0.00% MLMG::mgVcycle_bottom 65 0.0001133 0.0001166 0.0001223 0.00% TagBoxArray::mapPRD 2 0.0001117 0.0001153 0.0001221 0.00% Amr::grid_places() 2 7.445e-05 9.159e-05 0.0001112 0.00% Castro::react_state() 19 5.003e-05 5.755e-05 8.984e-05 0.00% Amr::FinalizeInit() 1 1.031e-05 1.761e-05 8.877e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.854e-05 6.823e-05 8.72e-05 0.00% MLPoisson::define() 10 6.575e-05 7.151e-05 8.052e-05 0.00% Castro::check_for_nan() 19 4.971e-05 5.549e-05 6.637e-05 0.00% Castro::apply_source_to_state() 19 4.933e-05 5.436e-05 6.366e-05 0.00% Castro::expand_state() 11 3.406e-05 4.21e-05 5.15e-05 0.00% Castro::create_source_corrector() 8 3.517e-05 4.014e-05 4.839e-05 0.00% Castro::derive() 18 4.085e-05 4.26e-05 4.685e-05 0.00% Amr::regrid() 1 1.057e-05 1.287e-05 3.027e-05 0.00% Gravity::actual_multilevel_solve() 1 2.606e-05 2.741e-05 2.917e-05 0.00% MLMG::computeMLResidual() 9 2.01e-05 2.232e-05 2.485e-05 0.00% Castro::advance() 1 7.317e-06 1.309e-05 2.176e-05 0.00% Amr::timeStep() 1 1.157e-05 1.597e-05 2.128e-05 0.00% Castro::post_init() 1 1.166e-05 1.277e-05 1.646e-05 0.00% Castro::post_regrid() 2 1.086e-05 1.283e-05 1.519e-05 0.00% Amr::InitializeInit() 1 6.423e-06 8.893e-06 1.466e-05 0.00% Amr::bldFineLevels() 1 7.103e-06 8.274e-06 1.071e-05 0.00% Amr::init() 1 6.743e-06 7.906e-06 8.827e-06 0.00% Castro::post_timestep() 1 5.788e-06 6.262e-06 6.793e-06 0.00% Amr::initialInit() 1 2.932e-06 3.343e-06 3.909e-06 0.00% Other 5645 0.02922 0.03795 0.05646 0.59% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.534 9.534 9.534 100.00% Amr::coarseTimeStep() 1 7.289 7.289 7.289 76.45% Amr::timeStep() 1 7.233 7.233 7.233 75.86% Castro::advance() 1 7.17 7.172 7.174 75.25% Castro::subcycle_advance_ctu() 1 7.119 7.121 7.125 74.73% Castro::do_advance_ctu() 11 7.095 7.097 7.101 74.48% Castro::construct_ctu_hydro_source() 11 3.37 3.386 3.403 35.69% Castro::clean_state() 62 2.945 3.043 3.169 33.24% Castro::computeTemp() 63 1.91 2.018 2.231 23.39% Castro::reset_internal_energy(MultiFab) 63 1.192 1.259 1.456 15.27% Amr::InitAmr() 1 1.318 1.318 1.318 13.82% Amr::init() 1 0.9258 0.9259 0.9259 9.71% Amr::initialInit() 1 0.8683 0.8684 0.8685 9.11% Castro::initialize_do_advance() 11 0.7786 0.794 0.8132 8.53% Gravity::solve_phi_with_mlmg() 10 0.7628 0.784 0.8052 8.44% Castro::evaluate_source_change() 144 0.7236 0.7305 0.7439 7.80% Castro::do_old_sources() 11 0.6911 0.7117 0.7341 7.70% Castro::construct_new_gravity() 8 0.6706 0.6916 0.7126 7.47% FabArray::ParallelCopy() 734 0.1487 0.4957 0.7065 7.41% Gravity::solve_for_phi() 8 0.653 0.6743 0.6955 7.29% Castro::enforce_min_density() 62 0.6188 0.6418 0.6889 7.23% Gravity::actual_solve_with_mlmg() 10 0.6556 0.6581 0.6607 6.93% Amr::InitializeInit() 1 0.4454 0.5144 0.6481 6.80% Amr::defBaseLevel() 1 0.4454 0.5144 0.6481 6.80% MLMG::solve() 9 0.6292 0.6311 0.6318 6.63% FabArray::ParallelCopy_finish() 734 0.0584 0.4051 0.6179 6.48% MLMG::oneIter() 65 0.563 0.5639 0.5649 5.93% amrex::Copy() 523 0.5402 0.5466 0.5601 5.87% MLMG::mgVcycle() 65 0.5555 0.5564 0.5574 5.85% Castro::do_new_sources() 8 0.4408 0.463 0.4809 5.04% Castro::normalize_species() 62 0.4424 0.4483 0.4645 4.87% Castro::reset_internal_energy(Fab) 126 0.4302 0.4359 0.4463 4.68% MultiFab::AddProduct() 2633 0.4246 0.4288 0.4368 4.58% Amr::FinalizeInit() 1 0.2203 0.354 0.423 4.44% AmrLevel::FillPatch() 82 0.2182 0.3124 0.4107 4.31% Castro::Castro() 1 0.3673 0.373 0.4035 4.23% Castro::initMFs() 1 0.3625 0.3685 0.3986 4.18% MLCellLinOp::smooth() 1560 0.2818 0.3077 0.3739 3.92% FillPatchIterator::Initialize 82 0.1479 0.2411 0.3402 3.57% FillPatchIterator::FillFromLevel0() 82 0.1363 0.2322 0.3322 3.48% FillPatchSingleLevel 82 0.136 0.2319 0.3318 3.48% MultiFab::Subtract() 125 0.2721 0.2747 0.2809 2.95% Castro::initData() 1 0.07594 0.14 0.2787 2.92% Amr::bldFineLevels() 1 0.04589 0.1794 0.248 2.60% FabArray::setVal() 3453 0.2215 0.2256 0.2342 2.46% MLMG::mgVcycle_up::1 65 0.0492 0.133 0.2163 2.27% MLPoisson::Fsmooth() 3120 0.1602 0.1735 0.208 2.18% Amr::grid_places() 2 0.002698 0.1364 0.2054 2.15% AmrMesh::MakeNewGrids() 2 0.002593 0.1363 0.2054 2.15% TagBoxArray::mapPRD 2 0.0007129 0.1297 0.2035 2.13% MLCellLinOp::applyBC() 3630 0.1417 0.1589 0.1997 2.09% MLMG::addInterpCorrection() 390 0.03218 0.154 0.1898 1.99% check_for_negative_density() 11 0.06307 0.1346 0.1824 1.91% FabArray::FillBoundary() 3240 0.1122 0.1272 0.1473 1.54% FabArray::sum() 2633 0.1417 0.1432 0.1457 1.53% Castro::expand_state() 11 0.09846 0.1183 0.1454 1.53% Gravity::fill_multipole_BCs() 10 0.1058 0.1252 0.1453 1.52% MLMG::mgVcycle_up::2 65 0.005457 0.0591 0.141 1.48% MLMG::mgVcycle_down::0 65 0.1327 0.1336 0.1346 1.41% Castro::post_init() 1 0.1308 0.1317 0.1326 1.39% MLMG::mgVcycle_up::0 65 0.1011 0.1013 0.1016 1.07% Amr::writePlotFile() 2 0.09938 0.09947 0.09948 1.04% MLCellLinOp::apply() 501 0.07696 0.0811 0.09906 1.04% FabArray::ParallelCopy_nowait() 734 0.08171 0.08872 0.0981 1.03% Gravity::multilevel_solve_for_new_phi() 1 0.09777 0.0978 0.09783 1.03% Gravity::actual_multilevel_solve() 1 0.09771 0.09773 0.09776 1.03% MLMG::mgVcycle_up::3 65 0.00545 0.01179 0.08781 0.92% Castro::post_regrid() 2 0.08491 0.08578 0.08676 0.91% FillBoundary_finish() 3240 0.05325 0.0617 0.08165 0.86% FillBoundary_nowait() 3240 0.0462 0.05835 0.0737 0.77% MLMG:computeResOfCorrection() 390 0.05702 0.0619 0.07358 0.77% MLCellLinOp::correctionResidual() 390 0.05631 0.06118 0.07282 0.76% Castro::derive() 18 0.06991 0.0706 0.07247 0.76% AmrLevel::derive() 18 0.06986 0.07056 0.07243 0.76% Castro::post_timestep() 1 0.05878 0.06122 0.06345 0.67% MLMG::mgVcycle_down::1 65 0.05727 0.05881 0.06091 0.64% Castro::create_source_corrector() 8 0.01622 0.04291 0.06073 0.64% Castro::construct_new_source() 48 0.05817 0.05882 0.05955 0.62% MLPoisson::Fapply() 501 0.04804 0.05039 0.05745 0.60% Castro::initialize_advance() 1 0.04892 0.0501 0.05237 0.55% Castro::construct_old_gravity() 11 0.02698 0.03867 0.04972 0.52% MLCellLinOp::solutionResidual() 75 0.04517 0.04729 0.04924 0.52% Castro::construct_old_source() 66 0.04326 0.04434 0.04905 0.51% Amr::regrid() 1 0.04359 0.04411 0.04463 0.47% Gravity::get_old_grav_vector() 11 0.02649 0.03339 0.03982 0.42% TagBoxArray::collate() 2 0.0001832 0.005299 0.03961 0.42% MLMG::computeResidual() 65 0.03814 0.03876 0.03926 0.41% FabArray::Saxpy() 133 0.03495 0.03557 0.03884 0.41% MLMG::mgVcycle_down::2 65 0.006393 0.02093 0.03718 0.39% MLMG::mgVcycle_bottom 65 0.0001816 0.002391 0.03547 0.37% MLMG::actualBottomSolve() 65 6.778e-05 0.002275 0.03535 0.37% MLCGSolver::bicgstab 4 0 0.002162 0.03459 0.36% Castro::construct_new_gravity_source() 8 0.03268 0.03308 0.03376 0.35% MLMG::mgVcycle_down::3 65 0.006124 0.009366 0.03192 0.33% FabArray::Xpay() 480 0.02577 0.02698 0.03029 0.32% Castro::construct_old_gravity_source() 11 0.02522 0.02594 0.02919 0.31% amrex::average_down 390 0.01661 0.01983 0.02853 0.30% VisMF::Write(FabArray) 4 0.01472 0.01625 0.02761 0.29% Castro::check_for_nan() 19 0.02327 0.02373 0.02736 0.29% MultiFab::contains_nan() 19 0.02322 0.02368 0.02729 0.29% Castro::retry_advance_ctu() 11 0.0223 0.02266 0.02322 0.24% Castro::sum_integrated_quantities() 2 0.02161 0.0219 0.02193 0.23% MLPoisson::define() 10 0.01574 0.0181 0.02131 0.22% Castro::react_state() 19 0.01514 0.0158 0.02101 0.22% Gravity::get_new_grav_vector() 9 0.01757 0.0179 0.01823 0.19% MLMG::mgVcycle_down::4 65 0.006051 0.006847 0.01752 0.18% Castro::apply_source_to_state() 19 0.01446 0.01466 0.01514 0.16% Castro::problem_post_timestep() 2 0.009067 0.01187 0.01498 0.16% Castro::wd_update() 2 0.008981 0.01183 0.01497 0.16% StateDataPhysBCFunct::() 82 0.0003465 0.007594 0.0145 0.15% Amr::writeSmallPlotFile() 2 0.01418 0.01419 0.01419 0.15% StateData::FillBoundary(geom) 62 0 0.007183 0.01406 0.15% Gravity::test_composite_phi() 1 0.01402 0.01403 0.01406 0.15% MLMG::mgVcycle_up::4 65 0.005473 0.00605 0.01364 0.14% MLCellLinOp::defineAuxData() 10 0.009872 0.0106 0.01254 0.13% MLMG::mgVcycle_down::5 65 0.006016 0.006448 0.0115 0.12% MLMG::prepareForSolve() 9 0.009985 0.01049 0.01144 0.12% Castro::save_data_for_retry() 3 0.01099 0.01107 0.01133 0.12% FabArray::mult() 40 0.008171 0.008435 0.01132 0.12% MLMG::mgVcycle_up::5 65 0.00553 0.005859 0.009842 0.10% MLMG::computeMLResidual() 9 0.005786 0.007958 0.00967 0.10% MLMG::ResNormInf() 74 0.006842 0.008065 0.009211 0.10% FabArray::norminf() 143 0.005603 0.005891 0.007808 0.08% amrex::Add() 65 0.00721 0.007317 0.007728 0.08% Castro::volWgtSum() 46 0.007323 0.007374 0.007415 0.08% Other 5645 0.0483 0.06462 0.09298 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 16 B 20 B 23 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-32-g28b010126a1b) finalized