Initializing AMReX (24.04-27-gd857968968ae)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.04-27-gd857968968ae) initialized Starting run at 21:01:10 UTC on 2024-04-26. Successfully read inputs file ... Castro git describe: 24.04-18-g6053b3e6c AMReX git describe: 24.04-27-gd85796896 Microphysics git describe: 24.04-9-g4ad225b7 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.175802035 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001075135 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.01049435 Gravity::multilevel_solve_for_new_phi() time = 0.099302272 Gravity::fill_multipole_BCs() time = 0.010469587 ... 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.073809567 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003933107 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.002201522 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.066310408 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.002355172 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.064736221 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.003289833 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.066454608 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.002570178 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.058442349 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.303974606 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.013799947 Gravity::solve_for_phi() time = 0.083426365 on level 0 Castro::construct_new_gravity() time = 0.085574974 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.059443406 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.002230591 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.058154829 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.311224677 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014211041 Gravity::solve_for_phi() time = 0.084429382 on level 0 Castro::construct_new_gravity() time = 0.086512981 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.058702493 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.002166295 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.058289576 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.314840766 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014796321 Gravity::solve_for_phi() time = 0.085473358 on level 0 Castro::construct_new_gravity() time = 0.08766143 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.060564847 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.002177365 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.063817814 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.317974937 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015821461 Gravity::solve_for_phi() time = 0.08681626 on level 0 Castro::construct_new_gravity() time = 0.089056454 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.059670838 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.002218487 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.059211774 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.317861038 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015664835 Gravity::solve_for_phi() time = 0.08566173 on level 0 Castro::construct_new_gravity() time = 0.087773232 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.059732593 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.002515196 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.059063222 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.320834357 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015399233 Gravity::solve_for_phi() time = 0.085883839 on level 0 Castro::construct_new_gravity() time = 0.088016283 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.059522472 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.003809216 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.064059516 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.320669042 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018059003 Gravity::solve_for_phi() time = 0.090015265 on level 0 Castro::construct_new_gravity() time = 0.092175913 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.058615557 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.002386747 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.06066958 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.313857067 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014906039 Gravity::solve_for_phi() time = 0.083894684 on level 0 Castro::construct_new_gravity() time = 0.086070852 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.059042346 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.004698780249 [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.031491246 [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.046679612 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.00397019 seconds Ending run at 21:01:19 UTC on 2024-04-26. Run time = 9.342666611 Run time without initialization = 7.082941289 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]: 9.343 ... 9.343 ... 9.343 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8312 2.8601 2.8917 30.95% Amr::InitAmr() 1 1.342 1.342 1.343 14.37% Castro::computeTemp() 63 0.6862 0.7387 0.819 8.77% amrex::Copy() 523 0.5821 0.5854 0.5905 6.32% FabArray::ParallelCopy_finish() 734 0.03384 0.3219 0.4691 5.02% Castro::normalize_species() 62 0.4483 0.4511 0.4566 4.89% Castro::reset_internal_energy(Fab) 126 0.4255 0.4286 0.4354 4.66% MultiFab::AddProduct() 2633 0.4227 0.4259 0.4344 4.65% Castro::initMFs() 1 0.3556 0.3634 0.3734 4.00% MultiFab::Subtract() 125 0.262 0.2638 0.2667 2.85% FabArray::setVal() 3453 0.2272 0.2293 0.2337 2.50% Castro::reset_internal_energy(MultiFab) 63 0.01813 0.05893 0.2087 2.23% MLPoisson::Fsmooth() 3120 0.1588 0.1714 0.2067 2.21% Castro::initData() 1 0.003366 0.02654 0.158 1.69% FabArray::sum() 2633 0.1464 0.1473 0.1504 1.61% Gravity::fill_multipole_BCs() 10 0.1027 0.1222 0.1403 1.50% check_for_negative_density() 11 0.02761 0.06269 0.09486 1.02% FabArray::ParallelCopy_nowait() 734 0.06346 0.07005 0.0785 0.84% FillBoundary_finish() 3240 0.05022 0.05799 0.07409 0.79% FillBoundary_nowait() 3240 0.03968 0.05301 0.06783 0.73% MLPoisson::Fapply() 501 0.04871 0.05108 0.05841 0.63% MLCellLinOp::applyBC() 3630 0.02391 0.03173 0.05252 0.56% AmrLevel::derive() 18 0.04029 0.04099 0.04594 0.49% Amr::writePlotFile() 2 0.02712 0.03939 0.04141 0.44% Castro::enforce_min_density() 62 0.01394 0.01942 0.03859 0.41% FabArray::Saxpy() 133 0.03416 0.03454 0.03768 0.40% Castro::do_old_sources() 11 0.0003138 0.0073 0.0366 0.39% Castro::do_new_sources() 8 0.0001761 0.006208 0.03584 0.38% Castro::construct_new_gravity_source() 8 0.03255 0.03308 0.03372 0.36% FabArray::Xpay() 480 0.02519 0.02628 0.02955 0.32% VisMF::Write(FabArray) 4 0.01119 0.01278 0.0276 0.30% Castro::construct_old_gravity_source() 11 0.02488 0.025 0.02512 0.27% MultiFab::contains_nan() 19 0.02282 0.02302 0.02323 0.25% MLMG::addInterpCorrection() 390 0.01728 0.01852 0.02198 0.24% amrex::average_down 390 0.01259 0.01362 0.01646 0.18% StateData::FillBoundary(geom) 62 0 0.007239 0.01396 0.15% Castro::evaluate_source_change() 144 0.009274 0.009488 0.01016 0.11% Castro::wd_update() 2 0.004823 0.007393 0.00994 0.11% MLCellLinOp::defineAuxData() 10 0.007715 0.008058 0.00921 0.10% FabArray::mult() 40 0.008481 0.008583 0.008887 0.10% FabArray::FillBoundary() 3240 0.006963 0.007263 0.00867 0.09% MLCGSolver::bicgstab 4 0 0.0005091 0.008145 0.09% FabArray::norminf() 143 0.005561 0.005758 0.007777 0.08% Castro::construct_new_source() 48 0.007379 0.007453 0.007549 0.08% Castro::volWgtSum() 46 0.007253 0.007298 0.007347 0.08% MLCellLinOp::smooth() 1560 0.005607 0.005826 0.00647 0.07% FillPatchIterator::Initialize 82 0.005509 0.005602 0.005788 0.06% Castro::construct_old_source() 66 0.005513 0.005567 0.005629 0.06% Amr::writeSmallPlotFile() 2 0.003572 0.005176 0.005628 0.06% MLMG::prepareForSolve() 9 0.005044 0.005158 0.005526 0.06% FabArrayBase::getCPC() 1022 0.00363 0.00373 0.004055 0.04% Gravity::actual_solve_with_mlmg() 10 0.003485 0.003634 0.004026 0.04% MLMG::ResNormInf() 74 0.002382 0.003155 0.003983 0.04% Gravity::get_old_grav_vector() 11 0.003633 0.003667 0.003725 0.04% MLMG::solve() 9 0.0004715 0.001325 0.00358 0.04% Castro::construct_old_gravity() 11 0.0002314 0.001138 0.003376 0.04% Gravity::get_new_grav_vector() 9 0.003079 0.003133 0.003277 0.04% Castro::initialize_do_advance() 11 0.002427 0.002574 0.002746 0.03% MLCellLinOp::apply() 501 0.001005 0.001115 0.002383 0.03% FabArray::ParallelCopy() 734 0.001819 0.001888 0.002045 0.02% Castro::sum_integrated_quantities() 2 0.001585 0.001777 0.001944 0.02% AmrLevel::FillPatch() 82 0.001104 0.001143 0.001283 0.01% main() 1 0.001009 0.001128 0.00122 0.01% MLCellLinOp::correctionResidual() 390 0.0009171 0.0009452 0.0009763 0.01% MLMG::mgVcycle() 65 0.0007958 0.0008392 0.0009431 0.01% Castro::save_data_for_retry() 3 0.0007737 0.000795 0.0008248 0.01% MLMG:computeResOfCorrection() 390 0.0007087 0.0007346 0.0007972 0.01% MLMG::mgVcycle_down::4 65 0.0002814 0.0003199 0.0007273 0.01% Castro::Castro() 1 0.0005117 0.0005568 0.0006426 0.01% Gravity::solve_for_phi() 8 0.0003914 0.0004595 0.0005956 0.01% MLMG::actualBottomSolve() 65 6.845e-05 0.0001015 0.000552 0.01% Castro::subcycle_advance_ctu() 1 9.543e-05 0.0001355 0.0005516 0.01% AmrMesh::MakeNewGrids() 2 0.0004935 0.0005111 0.0005378 0.01% StateDataPhysBCFunct::() 82 0.0003553 0.0004116 0.0004927 0.01% Castro::construct_new_gravity() 8 0.0001422 0.0002036 0.0004385 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003707 0.0003858 0.0004107 0.00% Castro::problem_post_timestep() 2 1.107e-05 3.631e-05 0.0003993 0.00% Castro::clean_state() 62 0.000324 0.0003521 0.0003942 0.00% MLMG::mgVcycle_down::2 65 0.0003098 0.0003346 0.0003799 0.00% MLMG::mgVcycle_down::0 65 0.0003324 0.0003411 0.0003614 0.00% MLMG::mgVcycle_down::1 65 0.0003275 0.00034 0.0003572 0.00% MLMG::mgVcycle_down::3 65 0.0002906 0.0003093 0.0003545 0.00% MLMG::mgVcycle_down::5 65 0.0002906 0.0003009 0.0003398 0.00% FillPatchSingleLevel 82 0.0002568 0.000267 0.0002945 0.00% MLMG::mgVcycle_up::1 65 0.0002334 0.0002461 0.0002733 0.00% Castro::do_advance_ctu() 11 0.000225 0.0002429 0.0002711 0.00% MLMG::mgVcycle_up::0 65 0.000231 0.0002379 0.0002552 0.00% MLMG::mgVcycle_up::2 65 0.0001991 0.0002195 0.0002519 0.00% MLCellLinOp::solutionResidual() 75 0.0002152 0.0002231 0.0002461 0.00% MLMG::mgVcycle_up::3 65 0.0001971 0.0002058 0.0002436 0.00% MLMG::mgVcycle_up::5 65 0.000206 0.000212 0.0002388 0.00% MLMG::mgVcycle_up::4 65 0.0002009 0.0002057 0.0002263 0.00% MLMG::oneIter() 65 0.0001957 0.0002022 0.0002246 0.00% Amr::defBaseLevel() 1 0.000106 0.0001504 0.0001766 0.00% Amr::coarseTimeStep() 1 0.0001231 0.0001526 0.000164 0.00% Amr::grid_places() 2 7.519e-05 0.0001001 0.0001575 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001138 0.0001268 0.0001557 0.00% Castro::retry_advance_ctu() 11 5.433e-05 6.656e-05 0.0001469 0.00% MLMG::computeResidual() 65 0.0001215 0.000128 0.0001442 0.00% Castro::initialize_advance() 1 0.0001126 0.0001195 0.000141 0.00% MLMG::mgVcycle_bottom 65 0.0001153 0.0001186 0.0001297 0.00% TagBoxArray::mapPRD 2 0.0001148 0.0001208 0.0001297 0.00% Gravity::test_composite_phi() 1 9.392e-05 0.0001042 0.0001282 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.669e-05 7.465e-05 0.0001002 0.00% MLPoisson::define() 10 6.635e-05 7.321e-05 8.151e-05 0.00% Amr::FinalizeInit() 1 1.147e-05 1.961e-05 7.038e-05 0.00% Castro::react_state() 19 5.239e-05 5.552e-05 6.414e-05 0.00% Castro::check_for_nan() 19 5.247e-05 5.662e-05 6.164e-05 0.00% Castro::apply_source_to_state() 19 4.961e-05 5.23e-05 5.9e-05 0.00% Castro::create_source_corrector() 8 3.634e-05 4.089e-05 5.019e-05 0.00% Castro::derive() 18 4.101e-05 4.289e-05 4.529e-05 0.00% Castro::expand_state() 11 3.238e-05 3.619e-05 4.245e-05 0.00% Amr::regrid() 1 1.077e-05 1.299e-05 3.461e-05 0.00% Gravity::actual_multilevel_solve() 1 2.736e-05 2.824e-05 3.086e-05 0.00% Amr::timeStep() 1 1.188e-05 1.411e-05 2.527e-05 0.00% Amr::init() 1 5.947e-06 9.184e-06 2.511e-05 0.00% MLMG::computeMLResidual() 9 2.003e-05 2.183e-05 2.355e-05 0.00% Amr::InitializeInit() 1 6.779e-06 9.96e-06 1.975e-05 0.00% Castro::advance() 1 6.972e-06 9.451e-06 1.614e-05 0.00% Castro::post_init() 1 1.196e-05 1.28e-05 1.409e-05 0.00% Castro::post_regrid() 2 1.106e-05 1.179e-05 1.331e-05 0.00% Amr::bldFineLevels() 1 6.84e-06 7.791e-06 1.029e-05 0.00% Castro::post_timestep() 1 5.431e-06 5.819e-06 6.804e-06 0.00% Amr::initialInit() 1 3.092e-06 3.507e-06 5.018e-06 0.00% Other 4690 0.03242 0.04103 0.06031 0.65% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.343 9.343 9.343 100.00% Amr::coarseTimeStep() 1 7.082 7.082 7.082 75.81% Amr::timeStep() 1 7.031 7.031 7.031 75.26% Castro::advance() 1 6.97 6.972 6.974 74.64% Castro::subcycle_advance_ctu() 1 6.92 6.922 6.924 74.12% Castro::do_advance_ctu() 11 6.895 6.897 6.9 73.85% Castro::construct_ctu_hydro_source() 11 3.272 3.291 3.309 35.41% Castro::clean_state() 62 2.959 3.018 3.103 33.21% Castro::computeTemp() 63 1.91 1.992 2.161 23.13% Castro::reset_internal_energy(MultiFab) 63 1.207 1.254 1.406 15.05% Amr::InitAmr() 1 1.342 1.342 1.343 14.37% Amr::init() 1 0.9163 0.9168 0.9169 9.81% Amr::initialInit() 1 0.8378 0.8383 0.8384 8.97% Gravity::solve_phi_with_mlmg() 10 0.7546 0.7755 0.7951 8.51% Castro::initialize_do_advance() 11 0.7497 0.7633 0.78 8.35% Castro::evaluate_source_change() 144 0.7309 0.7357 0.7513 8.04% Castro::construct_new_gravity() 8 0.6614 0.6818 0.7009 7.50% Gravity::solve_for_phi() 8 0.6436 0.6644 0.6838 7.32% Castro::do_old_sources() 11 0.652 0.665 0.6745 7.22% Gravity::actual_solve_with_mlmg() 10 0.6495 0.6515 0.6538 7.00% Castro::enforce_min_density() 62 0.6262 0.635 0.6503 6.96% MLMG::solve() 9 0.6228 0.6249 0.6256 6.70% Amr::InitializeInit() 1 0.4415 0.5048 0.6157 6.59% Amr::defBaseLevel() 1 0.4415 0.5048 0.6157 6.59% amrex::Copy() 523 0.5821 0.5854 0.5905 6.32% MLMG::oneIter() 65 0.5568 0.5579 0.5589 5.98% FabArray::ParallelCopy() 734 0.1236 0.411 0.5572 5.96% MLMG::mgVcycle() 65 0.5495 0.5506 0.5516 5.90% Castro::do_new_sources() 8 0.442 0.4604 0.4752 5.09% FabArray::ParallelCopy_finish() 734 0.03384 0.3219 0.4691 5.02% Castro::normalize_species() 62 0.4483 0.4511 0.4566 4.89% Castro::reset_internal_energy(Fab) 126 0.4255 0.4286 0.4354 4.66% MultiFab::AddProduct() 2633 0.4227 0.4259 0.4344 4.65% Amr::FinalizeInit() 1 0.2226 0.3335 0.3968 4.25% Castro::Castro() 1 0.3601 0.368 0.3794 4.06% Castro::initMFs() 1 0.3557 0.3635 0.3736 4.00% MLCellLinOp::smooth() 1560 0.2767 0.3023 0.3675 3.93% AmrLevel::FillPatch() 82 0.2009 0.2513 0.3054 3.27% MultiFab::Subtract() 125 0.262 0.2638 0.2667 2.85% Castro::initData() 1 0.07751 0.1351 0.253 2.71% FabArray::setVal() 3453 0.2272 0.2293 0.2337 2.50% FillPatchIterator::Initialize 82 0.1255 0.1753 0.2298 2.46% FillPatchIterator::FillFromLevel0() 82 0.1149 0.167 0.2222 2.38% FillPatchSingleLevel 82 0.1145 0.1666 0.2218 2.37% Amr::bldFineLevels() 1 0.04675 0.1572 0.221 2.37% MLMG::mgVcycle_up::1 65 0.04789 0.1309 0.2135 2.28% MLPoisson::Fsmooth() 3120 0.1588 0.1714 0.2067 2.21% MLCellLinOp::applyBC() 3630 0.1397 0.1556 0.1933 2.07% MLMG::addInterpCorrection() 390 0.03322 0.153 0.1883 2.02% Amr::grid_places() 2 0.002687 0.1135 0.1769 1.89% AmrMesh::MakeNewGrids() 2 0.002576 0.1134 0.1768 1.89% TagBoxArray::mapPRD 2 0.0008649 0.1118 0.1751 1.87% FabArray::sum() 2633 0.1464 0.1473 0.1504 1.61% FabArray::FillBoundary() 3240 0.1099 0.1239 0.143 1.53% Gravity::fill_multipole_BCs() 10 0.1039 0.1233 0.1414 1.51% MLMG::mgVcycle_up::2 65 0.005463 0.05811 0.1391 1.49% Castro::post_init() 1 0.1318 0.1328 0.1341 1.44% MLMG::mgVcycle_down::0 65 0.1316 0.1325 0.1334 1.43% Castro::expand_state() 11 0.1018 0.1148 0.1325 1.42% Amr::writePlotFile() 2 0.121 0.1211 0.1211 1.30% MLMG::mgVcycle_up::0 65 0.09863 0.09915 0.09951 1.07% Gravity::multilevel_solve_for_new_phi() 1 0.09926 0.09929 0.09932 1.06% Gravity::actual_multilevel_solve() 1 0.09919 0.09922 0.09925 1.06% MLCellLinOp::apply() 501 0.07597 0.08176 0.09878 1.06% FabArray::ParallelCopy_nowait() 734 0.08015 0.08727 0.09641 1.03% check_for_negative_density() 11 0.02761 0.06269 0.09486 1.02% MLMG::mgVcycle_up::3 65 0.005399 0.01186 0.08895 0.95% Castro::post_regrid() 2 0.08573 0.08712 0.08811 0.94% Castro::derive() 18 0.07016 0.07123 0.0764 0.82% AmrLevel::derive() 18 0.07012 0.07119 0.07636 0.82% FillBoundary_finish() 3240 0.05022 0.05799 0.07409 0.79% FillBoundary_nowait() 3240 0.04487 0.05861 0.07381 0.79% MLMG:computeResOfCorrection() 390 0.05736 0.06217 0.07351 0.79% MLCellLinOp::correctionResidual() 390 0.05665 0.06144 0.07271 0.78% Castro::post_timestep() 1 0.05772 0.05966 0.06182 0.66% MLMG::mgVcycle_down::1 65 0.05762 0.05937 0.06081 0.65% Castro::construct_new_source() 48 0.0581 0.05863 0.05943 0.64% MLPoisson::Fapply() 501 0.04871 0.05108 0.05841 0.63% Castro::initialize_advance() 1 0.04873 0.04981 0.05149 0.55% MLCellLinOp::solutionResidual() 75 0.04407 0.04697 0.04863 0.52% Amr::regrid() 1 0.04402 0.04472 0.0452 0.48% Castro::construct_old_source() 66 0.04279 0.04312 0.04348 0.47% MLMG::computeResidual() 65 0.03812 0.03854 0.03932 0.42% FabArray::Saxpy() 133 0.03416 0.03454 0.03768 0.40% MLMG::mgVcycle_down::2 65 0.00637 0.02096 0.03675 0.39% MLMG::mgVcycle_bottom 65 0.0001845 0.002424 0.03594 0.38% MLMG::actualBottomSolve() 65 6.845e-05 0.002306 0.03582 0.38% MLCGSolver::bicgstab 4 0 0.002189 0.03503 0.37% Castro::construct_new_gravity_source() 8 0.03255 0.03308 0.03372 0.36% Castro::create_source_corrector() 8 0.01668 0.02477 0.03223 0.34% MLMG::mgVcycle_down::3 65 0.006033 0.009188 0.03062 0.33% FabArray::Xpay() 480 0.02519 0.02628 0.02955 0.32% Castro::construct_old_gravity() 11 0.02316 0.02583 0.02836 0.30% amrex::average_down 390 0.01674 0.02001 0.02784 0.30% VisMF::Write(FabArray) 4 0.01119 0.01278 0.0276 0.30% Gravity::get_old_grav_vector() 11 0.02291 0.02469 0.02591 0.28% Castro::construct_old_gravity_source() 11 0.02488 0.025 0.02512 0.27% Castro::retry_advance_ctu() 11 0.02353 0.02367 0.02383 0.26% Castro::check_for_nan() 19 0.02288 0.02308 0.02329 0.25% MultiFab::contains_nan() 19 0.02282 0.02302 0.02323 0.25% Castro::sum_integrated_quantities() 2 0.02088 0.02114 0.02116 0.23% MLPoisson::define() 10 0.01523 0.01758 0.02041 0.22% Gravity::get_new_grav_vector() 9 0.01783 0.01814 0.01864 0.20% MLMG::mgVcycle_down::4 65 0.006017 0.00688 0.01747 0.19% Castro::react_state() 19 0.0153 0.01552 0.01621 0.17% Castro::apply_source_to_state() 19 0.01406 0.01425 0.01481 0.16% StateDataPhysBCFunct::() 82 0.0003553 0.00765 0.01442 0.15% Gravity::test_composite_phi() 1 0.01395 0.01398 0.01403 0.15% StateData::FillBoundary(geom) 62 0 0.007239 0.01396 0.15% MLMG::mgVcycle_up::4 65 0.005463 0.006074 0.01378 0.15% Castro::problem_post_timestep() 2 0.009091 0.01112 0.01338 0.14% Castro::wd_update() 2 0.009079 0.01109 0.01337 0.14% MLCellLinOp::defineAuxData() 10 0.009908 0.01065 0.01277 0.14% Castro::save_data_for_retry() 3 0.01198 0.01204 0.01209 0.13% MLMG::prepareForSolve() 9 0.01021 0.01061 0.01174 0.13% MLMG::mgVcycle_down::5 65 0.006 0.00646 0.01162 0.12% MLMG::mgVcycle_up::5 65 0.005515 0.005891 0.009882 0.11% MLMG::computeMLResidual() 9 0.005381 0.007833 0.009504 0.10% FabArray::mult() 40 0.008481 0.008583 0.008887 0.10% MLMG::ResNormInf() 74 0.007293 0.008034 0.00888 0.10% Amr::writeSmallPlotFile() 2 0.008276 0.008293 0.008298 0.09% FabArray::norminf() 143 0.005561 0.005758 0.007777 0.08% Castro::volWgtSum() 46 0.007253 0.007298 0.007347 0.08% FabArrayBase::getCPC() 1022 0.004896 0.00566 0.007311 0.08% Other 4690 0.05028 0.06504 0.09289 0.99% -------------------------------------------------------------------------------------------- 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 18 B 21 B 27 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.04-27-gd857968968ae) finalized