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:38:02 UTC on 2024-05-14. 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.176734166 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001038953 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010554325 Gravity::multilevel_solve_for_new_phi() time = 0.102401735 Gravity::fill_multipole_BCs() time = 0.011775762 ... 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.048116277 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003963023 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.003519184 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.066754242 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 = -9.947567811e-05 at index 0, 62, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232Invalid density = -8.076810125 at index 20, 81, 0 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.006483424 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.0668818 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.006542146 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.061736614 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.007264397 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.065566532 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.31931608 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015481981 Gravity::solve_for_phi() time = 0.090622328 on level 0 Castro::construct_new_gravity() time = 0.092690222 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.061198443 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.002370732 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.063196677 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.316206828 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014636631 Gravity::solve_for_phi() time = 0.087799801 on level 0 Castro::construct_new_gravity() time = 0.089855284 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.060957507 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.002108303 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.064728123 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.326874465 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016207508 Gravity::solve_for_phi() time = 0.089380254 on level 0 Castro::construct_new_gravity() time = 0.091422241 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.058839679 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.007124727 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.064983288 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.320634772 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015478382 Gravity::solve_for_phi() time = 0.08420989 on level 0 Castro::construct_new_gravity() time = 0.086343078 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.059576676 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.002109664 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.05842552 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.313246397 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015346223 Gravity::solve_for_phi() time = 0.084102512 on level 0 Castro::construct_new_gravity() time = 0.086202978 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.059263825 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.002040236 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.058911847 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.31594912 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01520511 Gravity::solve_for_phi() time = 0.084156923 on level 0 Castro::construct_new_gravity() time = 0.086263037 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.05918938 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.002106616 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.059331335 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.316467015 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015248198 Gravity::solve_for_phi() time = 0.083946523 on level 0 Castro::construct_new_gravity() time = 0.086062286 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.059274939 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.002179522 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.058942845 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.313693462 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015458262 Gravity::solve_for_phi() time = 0.083570867 on level 0 Castro::construct_new_gravity() time = 0.085643062 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.05925376 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.00468215758 [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.056514364 [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.046256778 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.003869137 seconds Ending run at 07:38:11 UTC on 2024-05-14. Run time = 9.331418938 Run time without initialization = 7.107605953 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.331 ... 9.331 ... 9.332 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8250 2.8604 2.9197 31.29% Amr::InitAmr() 1 1.32 1.32 1.32 14.15% Castro::computeTemp() 63 0.7105 0.7551 0.8495 9.10% amrex::Copy() 523 0.5455 0.5517 0.5737 6.15% FabArray::ParallelCopy_finish() 734 0.02469 0.3254 0.5 5.36% Castro::normalize_species() 62 0.4411 0.4468 0.4582 4.91% Castro::reset_internal_energy(Fab) 126 0.4284 0.4357 0.4449 4.77% MultiFab::AddProduct() 2633 0.4279 0.4316 0.4441 4.76% Castro::initMFs() 1 0.3602 0.3699 0.4024 4.31% MultiFab::Subtract() 125 0.2596 0.2624 0.2741 2.94% FabArray::setVal() 3453 0.2195 0.2226 0.2298 2.46% MLPoisson::Fsmooth() 3120 0.158 0.172 0.2069 2.22% Castro::reset_internal_energy(MultiFab) 63 0.01781 0.06221 0.1996 2.14% Castro::initData() 1 0.003267 0.02885 0.1706 1.83% FabArray::sum() 2633 0.1432 0.147 0.1509 1.62% Gravity::fill_multipole_BCs() 10 0.1046 0.124 0.1436 1.54% check_for_negative_density() 11 0.02975 0.08524 0.1315 1.41% FillBoundary_finish() 3240 0.05368 0.06464 0.08067 0.86% FabArray::ParallelCopy_nowait() 734 0.06206 0.06937 0.07846 0.84% FillBoundary_nowait() 3240 0.03925 0.05225 0.06675 0.72% MLPoisson::Fapply() 501 0.04754 0.04989 0.05726 0.61% MLCellLinOp::applyBC() 3630 0.02393 0.03154 0.05189 0.56% Castro::enforce_min_density() 62 0.01373 0.02022 0.05121 0.55% AmrLevel::derive() 18 0.04086 0.04149 0.04221 0.45% FabArray::Saxpy() 133 0.03571 0.0362 0.03978 0.43% Castro::construct_old_gravity_source() 11 0.02542 0.02658 0.03945 0.42% Castro::do_old_sources() 11 0.0002978 0.00747 0.03892 0.42% Castro::construct_new_gravity_source() 8 0.03216 0.03262 0.03307 0.35% FabArray::Xpay() 480 0.02612 0.0274 0.03104 0.33% Castro::do_new_sources() 8 0.0001664 0.005857 0.02963 0.32% MultiFab::contains_nan() 19 0.02291 0.02324 0.02467 0.26% VisMF::Write(FabArray) 4 0.01134 0.01254 0.02217 0.24% MLMG::addInterpCorrection() 390 0.01679 0.01801 0.02146 0.23% Amr::writePlotFile() 2 0.00659 0.01473 0.01627 0.17% amrex::average_down 390 0.01228 0.01339 0.01626 0.17% StateData::FillBoundary(geom) 62 0 0.007138 0.01498 0.16% Castro::wd_update() 2 0.005003 0.008197 0.01141 0.12% Castro::evaluate_source_change() 144 0.009365 0.009522 0.009962 0.11% Castro::construct_old_gravity() 11 0.0002085 0.001678 0.009729 0.10% MLCellLinOp::defineAuxData() 10 0.007425 0.007825 0.008896 0.10% FabArray::FillBoundary() 3240 0.006866 0.007134 0.008536 0.09% MLCGSolver::bicgstab 4 0 0.000512 0.008192 0.09% FabArray::mult() 40 0.007846 0.007947 0.00809 0.09% FabArray::norminf() 143 0.005615 0.005864 0.007875 0.08% Castro::construct_new_source() 48 0.007313 0.007403 0.007685 0.08% Castro::volWgtSum() 46 0.007318 0.007421 0.007564 0.08% MLCellLinOp::smooth() 1560 0.00568 0.005878 0.006391 0.07% Castro::construct_old_source() 66 0.005445 0.005511 0.005789 0.06% FillPatchIterator::Initialize 82 0.005439 0.005501 0.005688 0.06% Amr::writeSmallPlotFile() 2 0.003588 0.00521 0.005659 0.06% MLMG::prepareForSolve() 9 0.004925 0.005028 0.00527 0.06% MLMG::ResNormInf() 74 0.002732 0.003742 0.004657 0.05% Gravity::actual_solve_with_mlmg() 10 0.003496 0.00361 0.004038 0.04% FabArrayBase::getCPC() 1022 0.003629 0.003736 0.003894 0.04% Gravity::get_old_grav_vector() 11 0.003473 0.003576 0.003774 0.04% MLMG::solve() 9 0.0007049 0.001451 0.00341 0.04% Gravity::get_new_grav_vector() 9 0.002959 0.003074 0.003269 0.04% Castro::initialize_do_advance() 11 0.002377 0.002516 0.002762 0.03% MLCellLinOp::apply() 501 0.001018 0.001123 0.002371 0.03% Castro::sum_integrated_quantities() 2 0.001615 0.001902 0.002132 0.02% FabArray::ParallelCopy() 734 0.00176 0.001836 0.001974 0.02% main() 1 0.001007 0.001174 0.00136 0.01% AmrLevel::FillPatch() 82 0.001079 0.001102 0.001152 0.01% MLCellLinOp::correctionResidual() 390 0.0009204 0.0009515 0.001004 0.01% MLMG::mgVcycle() 65 0.0007867 0.0008258 0.0009134 0.01% Castro::save_data_for_retry() 3 0.0007549 0.0007682 0.000782 0.01% MLMG:computeResOfCorrection() 390 0.0006925 0.0007161 0.0007473 0.01% Gravity::solve_for_phi() 8 0.0003767 0.0004492 0.0006788 0.01% Castro::subcycle_advance_ctu() 1 9.13e-05 0.0001438 0.0006323 0.01% Castro::Castro() 1 0.000511 0.0005602 0.0006126 0.01% AmrMesh::MakeNewGrids() 2 0.0004961 0.0005162 0.0005368 0.01% MLMG::actualBottomSolve() 65 6.772e-05 9.788e-05 0.0005341 0.01% Castro::construct_new_gravity() 8 0.0001393 0.0002011 0.0004701 0.01% StateDataPhysBCFunct::() 82 0.00034 0.0004017 0.0004629 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003567 0.0003715 0.000385 0.00% MLMG::mgVcycle_down::2 65 0.0003127 0.0003408 0.0003822 0.00% Castro::clean_state() 62 0.0003248 0.0003407 0.0003661 0.00% MLMG::mgVcycle_down::0 65 0.0003316 0.0003398 0.000366 0.00% MLMG::mgVcycle_down::3 65 0.0002933 0.0003162 0.000363 0.00% MLMG::mgVcycle_down::1 65 0.0003277 0.0003383 0.0003573 0.00% MLMG::mgVcycle_down::4 65 0.0002864 0.0002972 0.0003422 0.00% MLMG::mgVcycle_down::5 65 0.0002974 0.000306 0.000337 0.00% FillPatchSingleLevel 82 0.0002491 0.0002572 0.0002693 0.00% Castro::problem_post_timestep() 2 1.031e-05 2.697e-05 0.0002622 0.00% MLMG::mgVcycle_up::1 65 0.0002281 0.0002365 0.0002565 0.00% MLCellLinOp::solutionResidual() 75 0.0002209 0.0002266 0.0002551 0.00% MLMG::mgVcycle_up::0 65 0.0002294 0.0002386 0.0002541 0.00% MLMG::mgVcycle_up::2 65 0.0001981 0.0002144 0.0002492 0.00% Castro::do_advance_ctu() 11 0.0002123 0.0002255 0.0002413 0.00% MLMG::mgVcycle_up::5 65 0.0002013 0.0002072 0.0002274 0.00% MLMG::mgVcycle_up::3 65 0.0001954 0.0002024 0.0002252 0.00% MLMG::mgVcycle_up::4 65 0.0001982 0.000202 0.0002242 0.00% MLMG::oneIter() 65 0.0001952 0.000201 0.0002126 0.00% Gravity::update_max_rhs() 2 0.0001584 0.0001621 0.0001679 0.00% Castro::retry_advance_ctu() 11 5.136e-05 6.905e-05 0.0001486 0.00% Amr::coarseTimeStep() 1 0.000123 0.0001438 0.0001483 0.00% Amr::defBaseLevel() 1 0.0001048 0.0001191 0.0001353 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001136 0.0001203 0.0001321 0.00% MLMG::computeResidual() 65 0.000119 0.0001218 0.0001289 0.00% TagBoxArray::mapPRD 2 0.0001095 0.0001167 0.000127 0.00% MLMG::mgVcycle_bottom 65 0.0001145 0.0001181 0.0001244 0.00% Amr::grid_places() 2 7.337e-05 9.412e-05 0.0001227 0.00% Castro::initialize_advance() 1 0.0001107 0.0001147 0.0001188 0.00% Gravity::test_composite_phi() 1 9.014e-05 9.524e-05 0.0001116 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.94e-05 6.728e-05 9.623e-05 0.00% MLPoisson::define() 10 6.65e-05 7.125e-05 8.817e-05 0.00% Amr::FinalizeInit() 1 1.085e-05 1.88e-05 7.03e-05 0.00% Castro::react_state() 19 4.923e-05 5.347e-05 5.921e-05 0.00% Castro::check_for_nan() 19 5.057e-05 5.372e-05 5.707e-05 0.00% Castro::apply_source_to_state() 19 4.874e-05 5.139e-05 5.359e-05 0.00% Castro::derive() 18 4.001e-05 4.132e-05 4.318e-05 0.00% Castro::create_source_corrector() 8 3.485e-05 3.787e-05 4.112e-05 0.00% Castro::expand_state() 11 3.087e-05 3.42e-05 3.809e-05 0.00% Amr::regrid() 1 1.061e-05 1.235e-05 2.935e-05 0.00% Gravity::actual_multilevel_solve() 1 2.583e-05 2.718e-05 2.863e-05 0.00% MLMG::computeMLResidual() 9 2.015e-05 2.11e-05 2.228e-05 0.00% Amr::timeStep() 1 1.21e-05 1.483e-05 2.064e-05 0.00% Castro::post_regrid() 2 1.085e-05 1.189e-05 1.445e-05 0.00% Castro::advance() 1 7.001e-06 9.087e-06 1.399e-05 0.00% Castro::post_init() 1 1.194e-05 1.259e-05 1.343e-05 0.00% Amr::InitializeInit() 1 6.886e-06 8.661e-06 1.263e-05 0.00% Amr::init() 1 6.296e-06 7.393e-06 1.244e-05 0.00% Amr::bldFineLevels() 1 7.067e-06 8.131e-06 9.792e-06 0.00% Castro::post_timestep() 1 5.602e-06 6.021e-06 6.425e-06 0.00% Amr::initialInit() 1 3.003e-06 3.294e-06 3.812e-06 0.00% Other 4688 0.03224 0.04405 0.06365 0.68% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.331 9.331 9.331 100.00% Amr::coarseTimeStep() 1 7.107 7.107 7.107 76.16% Amr::timeStep() 1 7.056 7.056 7.057 75.62% Castro::advance() 1 6.994 6.996 6.999 75.00% Castro::subcycle_advance_ctu() 1 6.944 6.946 6.949 74.47% Castro::do_advance_ctu() 11 6.919 6.922 6.926 74.22% Castro::construct_ctu_hydro_source() 11 3.291 3.311 3.35 35.90% Castro::clean_state() 62 2.934 3.011 3.122 33.46% Castro::computeTemp() 63 1.905 2.002 2.184 23.40% Castro::reset_internal_energy(MultiFab) 63 1.191 1.247 1.388 14.87% Amr::InitAmr() 1 1.32 1.32 1.32 14.15% Amr::init() 1 0.9031 0.9031 0.9031 9.68% Amr::initialInit() 1 0.8503 0.8503 0.8504 9.11% Gravity::solve_phi_with_mlmg() 10 0.7615 0.7823 0.8025 8.60% Castro::initialize_do_advance() 11 0.7432 0.7626 0.7804 8.36% Castro::evaluate_source_change() 144 0.7272 0.736 0.7583 8.13% Castro::construct_new_gravity() 8 0.6634 0.6837 0.7035 7.54% Castro::do_old_sources() 11 0.6474 0.6693 0.6877 7.37% Gravity::solve_for_phi() 8 0.646 0.6668 0.687 7.36% Castro::enforce_min_density() 62 0.6136 0.6265 0.6598 7.07% Gravity::actual_solve_with_mlmg() 10 0.6545 0.6564 0.6588 7.06% MLMG::solve() 9 0.628 0.6301 0.6309 6.76% Amr::InitializeInit() 1 0.4442 0.516 0.619 6.63% Amr::defBaseLevel() 1 0.4442 0.516 0.619 6.63% FabArray::ParallelCopy() 734 0.1136 0.413 0.5858 6.28% amrex::Copy() 523 0.5455 0.5517 0.5737 6.15% MLMG::oneIter() 65 0.5625 0.5631 0.5637 6.04% MLMG::mgVcycle() 65 0.5551 0.5558 0.5563 5.96% FabArray::ParallelCopy_finish() 734 0.02469 0.3254 0.5 5.36% Castro::do_new_sources() 8 0.4382 0.4593 0.4772 5.11% Castro::normalize_species() 62 0.4411 0.4468 0.4582 4.91% Castro::reset_internal_energy(Fab) 126 0.4284 0.4357 0.4449 4.77% MultiFab::AddProduct() 2633 0.4279 0.4316 0.4441 4.76% Castro::Castro() 1 0.3647 0.3744 0.4077 4.37% Amr::FinalizeInit() 1 0.2313 0.3343 0.4061 4.35% Castro::initMFs() 1 0.3604 0.37 0.4025 4.31% MLCellLinOp::smooth() 1560 0.2836 0.3087 0.3725 3.99% AmrLevel::FillPatch() 82 0.1865 0.2572 0.3197 3.43% MultiFab::Subtract() 125 0.2596 0.2624 0.2741 2.94% FillPatchIterator::Initialize 82 0.1148 0.1854 0.249 2.67% Castro::initData() 1 0.07677 0.1404 0.2473 2.65% FillPatchIterator::FillFromLevel0() 82 0.1044 0.1772 0.2416 2.59% FillPatchSingleLevel 82 0.104 0.1769 0.2412 2.58% FabArray::setVal() 3453 0.2195 0.2226 0.2298 2.46% Amr::bldFineLevels() 1 0.04704 0.1498 0.2204 2.36% MLMG::mgVcycle_up::1 65 0.04852 0.1312 0.2133 2.29% MLPoisson::Fsmooth() 3120 0.158 0.172 0.2069 2.22% MLCellLinOp::applyBC() 3630 0.1463 0.1613 0.1981 2.12% MLMG::addInterpCorrection() 390 0.03206 0.1516 0.1876 2.01% Amr::grid_places() 2 0.00292 0.1059 0.1778 1.91% AmrMesh::MakeNewGrids() 2 0.002797 0.1058 0.1777 1.90% TagBoxArray::mapPRD 2 0.0008657 0.1041 0.1762 1.89% FabArray::sum() 2633 0.1432 0.147 0.1509 1.62% FabArray::FillBoundary() 3240 0.1145 0.1297 0.15 1.61% Gravity::fill_multipole_BCs() 10 0.1057 0.1251 0.1447 1.55% Castro::post_init() 1 0.136 0.1408 0.143 1.53% Castro::expand_state() 11 0.09807 0.1164 0.1396 1.50% MLMG::mgVcycle_down::0 65 0.1325 0.1338 0.1393 1.49% MLMG::mgVcycle_up::2 65 0.005406 0.05755 0.1371 1.47% check_for_negative_density() 11 0.02975 0.08524 0.1315 1.41% Gravity::multilevel_solve_for_new_phi() 1 0.1023 0.1024 0.1024 1.10% Gravity::actual_multilevel_solve() 1 0.1023 0.1023 0.1024 1.10% MLMG::mgVcycle_up::0 65 0.09688 0.09734 0.09789 1.05% MLCellLinOp::apply() 501 0.07511 0.08052 0.09755 1.05% FabArray::ParallelCopy_nowait() 734 0.07801 0.08574 0.09568 1.03% Amr::writePlotFile() 2 0.09488 0.09496 0.09497 1.02% Castro::post_regrid() 2 0.08525 0.08752 0.09229 0.99% MLMG::mgVcycle_up::3 65 0.005418 0.01171 0.08762 0.94% FillBoundary_finish() 3240 0.05368 0.06464 0.08067 0.86% MLMG:computeResOfCorrection() 390 0.05746 0.06248 0.07389 0.79% MLCellLinOp::correctionResidual() 390 0.05676 0.06176 0.07314 0.78% FillBoundary_nowait() 3240 0.04456 0.05793 0.07273 0.78% Castro::derive() 18 0.06912 0.07006 0.0708 0.76% AmrLevel::derive() 18 0.06908 0.07001 0.07076 0.76% MLMG::mgVcycle_down::1 65 0.0615 0.06552 0.06807 0.73% Castro::post_timestep() 1 0.0578 0.06005 0.06276 0.67% Castro::construct_new_source() 48 0.05812 0.05863 0.05917 0.63% Castro::construct_old_source() 66 0.04368 0.04496 0.05756 0.62% MLPoisson::Fapply() 501 0.04754 0.04989 0.05726 0.61% Castro::initialize_advance() 1 0.0488 0.05001 0.05292 0.57% MLCellLinOp::solutionResidual() 75 0.04398 0.04655 0.04829 0.52% Amr::regrid() 1 0.04367 0.04489 0.04794 0.51% FabArray::Saxpy() 133 0.03571 0.0362 0.03978 0.43% Castro::construct_old_gravity_source() 11 0.02542 0.02658 0.03945 0.42% MLMG::computeResidual() 65 0.03758 0.0381 0.0391 0.42% MLMG::mgVcycle_down::2 65 0.006376 0.02114 0.03734 0.40% Castro::construct_old_gravity() 11 0.02495 0.03203 0.03725 0.40% Gravity::get_old_grav_vector() 11 0.02108 0.03036 0.037 0.40% Castro::create_source_corrector() 8 0.01482 0.02357 0.03696 0.40% MLMG::mgVcycle_bottom 65 0.0001824 0.002415 0.03585 0.38% MLMG::actualBottomSolve() 65 6.772e-05 0.002297 0.03572 0.38% MLCGSolver::bicgstab 4 0 0.002185 0.03495 0.37% Castro::construct_new_gravity_source() 8 0.03216 0.03262 0.03307 0.35% FabArray::Xpay() 480 0.02612 0.0274 0.03104 0.33% MLMG::mgVcycle_down::3 65 0.006087 0.009141 0.03028 0.32% amrex::average_down 390 0.01644 0.01994 0.02849 0.31% Castro::check_for_nan() 19 0.02296 0.02329 0.02472 0.26% MultiFab::contains_nan() 19 0.02291 0.02324 0.02467 0.26% Castro::retry_advance_ctu() 11 0.02199 0.02267 0.02382 0.26% VisMF::Write(FabArray) 4 0.01134 0.01254 0.02217 0.24% Castro::sum_integrated_quantities() 2 0.02083 0.02099 0.021 0.23% MLPoisson::define() 10 0.01492 0.01744 0.02062 0.22% Gravity::get_new_grav_vector() 9 0.01711 0.01757 0.01797 0.19% MLMG::mgVcycle_down::4 65 0.006083 0.006847 0.01739 0.19% Castro::react_state() 19 0.01489 0.01514 0.01558 0.17% StateDataPhysBCFunct::() 82 0.00034 0.00754 0.01544 0.17% Gravity::test_composite_phi() 1 0.01518 0.01522 0.01525 0.16% Castro::apply_source_to_state() 19 0.01474 0.01487 0.01512 0.16% StateData::FillBoundary(geom) 62 0 0.007138 0.01498 0.16% Castro::problem_post_timestep() 2 0.009391 0.01183 0.01479 0.16% Castro::wd_update() 2 0.00938 0.01181 0.01478 0.16% MLMG::mgVcycle_up::4 65 0.00547 0.006031 0.01354 0.15% Castro::save_data_for_retry() 3 0.01102 0.0113 0.01255 0.13% MLCellLinOp::defineAuxData() 10 0.009578 0.01043 0.0125 0.13% MLMG::mgVcycle_down::5 65 0.00603 0.006438 0.0115 0.12% MLMG::prepareForSolve() 9 0.009935 0.01033 0.01144 0.12% MLMG::mgVcycle_up::5 65 0.005514 0.005845 0.009819 0.11% MLMG::computeMLResidual() 9 0.005384 0.007837 0.009611 0.10% MLMG::ResNormInf() 74 0.007808 0.00867 0.009558 0.10% Gravity::update_max_rhs() 2 0.001385 0.006151 0.00843 0.09% Amr::writeSmallPlotFile() 2 0.008211 0.008218 0.008229 0.09% FabArray::mult() 40 0.007846 0.007947 0.00809 0.09% FabArray::norminf() 143 0.005615 0.005864 0.007875 0.08% Castro::volWgtSum() 46 0.007318 0.007421 0.007564 0.08% FabArrayBase::getCPC() 1022 0.004892 0.005664 0.007331 0.08% Other 4688 0.04907 0.06544 0.09169 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 19 B 21 B 25 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-18-g8eff86d32d17) finalized