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:09:28 UTC on 2024-05-16. 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.184535112 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.00107708 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.011692667 Gravity::multilevel_solve_for_new_phi() time = 0.105624809 Gravity::fill_multipole_BCs() time = 0.010644384 ... 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.058712639 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.003888709 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.003150865 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.105215958 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -160.2179808 at index 18, 75, 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, 76, 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 = -61.58489855 at index 19, 77, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -108.2035501 at index 19, 78, 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.113047Invalid density = -15.15828486 at index 20, 80, 0 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 = Invalid density = -8.076810125 at index 20, 81, 0 -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.008833515 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.083411831 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.006532059 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.061681988 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.007440077 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.062270047 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.352640716 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01586416 Gravity::solve_for_phi() time = 0.088090824 on level 0 Castro::construct_new_gravity() time = 0.090262578 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.060851419 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.002189895 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.082363343 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.319754334 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016894162 Gravity::solve_for_phi() time = 0.122800673 on level 0 Castro::construct_new_gravity() time = 0.124988151 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.059942802 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.003791115 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.085958272 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.324122891 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015566051 Gravity::solve_for_phi() time = 0.088739602 on level 0 Castro::construct_new_gravity() time = 0.090867483 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.05983368 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.004284694 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.065303851 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.318928909 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015711251 Gravity::solve_for_phi() time = 0.089145664 on level 0 Castro::construct_new_gravity() time = 0.091274524 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.060689348 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.002247527 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.062024341 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.324037338 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016041947 Gravity::solve_for_phi() time = 0.089449684 on level 0 Castro::construct_new_gravity() time = 0.091612193 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.060635367 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.005497412 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.058548907 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.322721377 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015276352 Gravity::solve_for_phi() time = 0.087529725 on level 0 Castro::construct_new_gravity() time = 0.089662493 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.061006176 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.002527187 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.06206021 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.328012139 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015364982 Gravity::solve_for_phi() time = 0.083960231 on level 0 Castro::construct_new_gravity() time = 0.086068459 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.058930192 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.002802939 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.060005615 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.318718451 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01581185 Gravity::solve_for_phi() time = 0.08743133 on level 0 Castro::construct_new_gravity() time = 0.089566013 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.05915187 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.004491459441 [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.356635477 [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.046344414 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004006297 seconds Ending run at 07:09:38 UTC on 2024-05-16. Run time = 10.04046168 Run time without initialization = 7.407821177 Average number of zones advanced per microsecond: 0.004 Average number of zones advanced per microsecond per rank: 0.000 TinyProfiler total time across processes [min...avg...max]: 10.04 ... 10.04 ... 10.04 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8154 2.8806 2.9558 29.44% Amr::InitAmr() 1 1.715 1.715 1.715 17.08% Castro::computeTemp() 63 0.6898 0.7563 0.8426 8.39% FabArray::ParallelCopy_finish() 734 0.02746 0.4248 0.6026 6.00% amrex::Copy() 523 0.5431 0.5538 0.5813 5.79% Castro::normalize_species() 62 0.438 0.4505 0.4894 4.87% MultiFab::AddProduct() 2633 0.4253 0.4347 0.4684 4.66% Castro::reset_internal_energy(Fab) 126 0.4264 0.4339 0.4438 4.42% Castro::initMFs() 1 0.3538 0.362 0.3749 3.73% Castro::reset_internal_energy(MultiFab) 63 0.01776 0.08187 0.3049 3.04% MultiFab::Subtract() 125 0.258 0.2628 0.2696 2.68% FabArray::setVal() 3453 0.2184 0.2232 0.2345 2.34% check_for_negative_density() 11 0.0841 0.166 0.2253 2.24% MLPoisson::Fsmooth() 3120 0.1577 0.1732 0.2076 2.07% Castro::initData() 1 0.00328 0.02896 0.1705 1.70% Castro::enforce_min_density() 62 0.01379 0.04303 0.1597 1.59% FabArray::sum() 2633 0.1434 0.1483 0.1596 1.59% Gravity::fill_multipole_BCs() 10 0.1037 0.1248 0.1467 1.46% FillBoundary_finish() 3240 0.08888 0.1036 0.12 1.20% Castro::do_old_sources() 11 0.000319 0.009507 0.1084 1.08% FabArray::ParallelCopy_nowait() 734 0.06335 0.07015 0.08046 0.80% FillBoundary_nowait() 3240 0.04001 0.05255 0.06852 0.68% MLPoisson::Fapply() 501 0.04763 0.05009 0.05666 0.56% MLCellLinOp::applyBC() 3630 0.02478 0.03176 0.05135 0.51% AmrLevel::derive() 18 0.04072 0.04179 0.04525 0.45% FabArray::Saxpy() 133 0.03549 0.03649 0.04 0.40% Castro::do_new_sources() 8 0.000175 0.006254 0.03762 0.37% VisMF::Write(FabArray) 4 0.01118 0.01322 0.03426 0.34% Castro::construct_new_gravity_source() 8 0.03235 0.03264 0.0338 0.34% FabArray::Xpay() 480 0.02597 0.02729 0.03042 0.30% Amr::writePlotFile() 2 0.006529 0.02381 0.02767 0.28% Castro::construct_old_gravity_source() 11 0.02526 0.02574 0.02673 0.27% MultiFab::contains_nan() 19 0.02285 0.02318 0.0244 0.24% MLMG::addInterpCorrection() 390 0.01678 0.01823 0.02169 0.22% MLMG::ResNormInf() 74 0.006695 0.01143 0.01801 0.18% Castro::wd_update() 2 0.006713 0.01256 0.01629 0.16% amrex::average_down 390 0.01221 0.01328 0.01594 0.16% StateData::FillBoundary(geom) 62 0 0.006958 0.01362 0.14% Castro::evaluate_source_change() 144 0.009333 0.009647 0.01016 0.10% MLCellLinOp::defineAuxData() 10 0.007531 0.00798 0.009537 0.09% FabArray::mult() 40 0.007846 0.008003 0.009194 0.09% FabArray::FillBoundary() 3240 0.006833 0.007127 0.008361 0.08% amrex::Add() 65 0.00705 0.007288 0.008361 0.08% MLCGSolver::bicgstab 4 0 0.0005047 0.008074 0.08% FabArray::norminf() 143 0.005528 0.005829 0.007757 0.08% Castro::construct_new_source() 48 0.007254 0.007371 0.007721 0.08% Castro::construct_old_gravity() 11 0.0002207 0.001261 0.007527 0.07% Castro::volWgtSum() 46 0.007399 0.007441 0.007516 0.07% MLCellLinOp::smooth() 1560 0.005648 0.005901 0.006311 0.06% FillPatchIterator::Initialize 82 0.005389 0.005586 0.005879 0.06% Castro::construct_old_source() 66 0.005372 0.005504 0.005762 0.06% MLMG::prepareForSolve() 9 0.004991 0.005105 0.005621 0.06% Amr::writeSmallPlotFile() 2 0.003609 0.005198 0.005492 0.05% MLMG::solve() 9 0.0004037 0.001787 0.005303 0.05% Gravity::actual_solve_with_mlmg() 10 0.003472 0.0036 0.003984 0.04% Gravity::get_old_grav_vector() 11 0.003466 0.003564 0.003918 0.04% Castro::sum_integrated_quantities() 2 0.002333 0.003579 0.003886 0.04% Gravity::get_new_grav_vector() 9 0.003008 0.003117 0.00352 0.04% Castro::initialize_do_advance() 11 0.002481 0.002601 0.002762 0.03% MLCellLinOp::apply() 501 0.001008 0.001116 0.002312 0.02% FabArray::ParallelCopy() 734 0.001765 0.001842 0.001945 0.02% AmrLevel::FillPatch() 82 0.00108 0.001132 0.001242 0.01% main() 1 0.0008993 0.001098 0.001231 0.01% MLCellLinOp::correctionResidual() 390 0.0009253 0.0009593 0.001023 0.01% MLMG::mgVcycle() 65 0.0007903 0.0008254 0.0008764 0.01% Castro::save_data_for_retry() 3 0.000731 0.0007658 0.000803 0.01% MLMG:computeResOfCorrection() 390 0.0006901 0.0007158 0.000741 0.01% Castro::subcycle_advance_ctu() 1 8.936e-05 0.0001443 0.0007107 0.01% Castro::construct_new_gravity() 8 0.0001438 0.0002315 0.0005908 0.01% Gravity::solve_for_phi() 8 0.0003859 0.0004505 0.0005785 0.01% Castro::Castro() 1 0.0004667 0.0004925 0.000553 0.01% AmrMesh::MakeNewGrids() 2 0.0005017 0.0005144 0.0005385 0.01% MLMG::actualBottomSolve() 65 6.757e-05 9.889e-05 0.0005303 0.01% StateDataPhysBCFunct::() 82 0.0003395 0.0004166 0.0005104 0.01% FillPatchIterator::FillFromLevel0() 82 0.0003634 0.0003773 0.0003968 0.00% Castro::clean_state() 62 0.0003147 0.0003513 0.0003966 0.00% MLMG::mgVcycle_down::1 65 0.0003251 0.000339 0.0003844 0.00% MLMG::mgVcycle_down::2 65 0.0003146 0.0003408 0.0003676 0.00% MLMG::mgVcycle_down::0 65 0.0003309 0.0003385 0.0003633 0.00% MLMG::mgVcycle_down::3 65 0.0002933 0.0003168 0.0003529 0.00% MLMG::mgVcycle_down::4 65 0.0002851 0.0002986 0.0003362 0.00% MLMG::mgVcycle_down::5 65 0.0002946 0.000306 0.0003279 0.00% Castro::do_advance_ctu() 11 0.0002187 0.0002381 0.0002853 0.00% FillPatchSingleLevel 82 0.0002481 0.0002631 0.0002825 0.00% MLMG::mgVcycle_up::0 65 0.0002313 0.0002399 0.0002662 0.00% MLMG::mgVcycle_up::1 65 0.0002285 0.0002392 0.0002536 0.00% MLCellLinOp::solutionResidual() 75 0.0002206 0.0002265 0.000234 0.00% MLMG::mgVcycle_up::2 65 0.0001968 0.0002138 0.0002336 0.00% Castro::problem_post_timestep() 2 1.091e-05 2.604e-05 0.0002306 0.00% MLMG::mgVcycle_up::3 65 0.0001944 0.0002045 0.0002244 0.00% MLMG::oneIter() 65 0.0001976 0.0002035 0.0002227 0.00% MLMG::mgVcycle_up::5 65 0.000202 0.0002078 0.0002227 0.00% MLMG::mgVcycle_up::4 65 0.0001982 0.0002033 0.0002166 0.00% Gravity::update_max_rhs() 2 0.0001571 0.0001661 0.0002017 0.00% Castro::retry_advance_ctu() 11 5.03e-05 6.59e-05 0.0001796 0.00% Amr::coarseTimeStep() 1 0.0001177 0.0001402 0.0001443 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001133 0.0001207 0.0001433 0.00% Amr::defBaseLevel() 1 0.0001012 0.0001139 0.0001345 0.00% MLMG::computeResidual() 65 0.0001177 0.0001233 0.0001322 0.00% Amr::grid_places() 2 7.348e-05 9.408e-05 0.0001293 0.00% Gravity::test_composite_phi() 1 9.259e-05 0.0001032 0.0001263 0.00% MLMG::mgVcycle_bottom 65 0.000114 0.0001167 0.0001211 0.00% TagBoxArray::mapPRD 2 0.0001078 0.0001132 0.0001185 0.00% Castro::initialize_advance() 1 9.931e-05 0.0001064 0.0001175 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.034e-05 6.983e-05 9.586e-05 0.00% MLPoisson::define() 10 6.565e-05 7.051e-05 8.396e-05 0.00% Castro::check_for_nan() 19 5.036e-05 5.962e-05 8.141e-05 0.00% Castro::react_state() 19 5.016e-05 5.382e-05 6.631e-05 0.00% Amr::FinalizeInit() 1 8.448e-06 1.254e-05 6.517e-05 0.00% Castro::apply_source_to_state() 19 4.817e-05 5.208e-05 5.758e-05 0.00% Castro::derive() 18 4.048e-05 4.233e-05 4.524e-05 0.00% Castro::create_source_corrector() 8 3.467e-05 3.888e-05 4.51e-05 0.00% Castro::expand_state() 11 3.021e-05 3.322e-05 3.637e-05 0.00% Gravity::actual_multilevel_solve() 1 2.621e-05 2.735e-05 2.94e-05 0.00% Castro::post_regrid() 2 1.058e-05 1.239e-05 2.815e-05 0.00% Amr::regrid() 1 1.069e-05 1.248e-05 2.77e-05 0.00% MLMG::computeMLResidual() 9 2.014e-05 2.126e-05 2.262e-05 0.00% Amr::timeStep() 1 1.201e-05 1.54e-05 2.174e-05 0.00% Castro::post_init() 1 1.199e-05 1.28e-05 1.406e-05 0.00% Castro::advance() 1 6.927e-06 8.797e-06 1.235e-05 0.00% Amr::InitializeInit() 1 6.29e-06 7.509e-06 1.154e-05 0.00% Amr::init() 1 6.316e-06 7.21e-06 8.681e-06 0.00% Amr::bldFineLevels() 1 6.85e-06 7.41e-06 8.368e-06 0.00% Castro::post_timestep() 1 5.74e-06 6.099e-06 6.545e-06 0.00% Amr::initialInit() 1 3.008e-06 3.254e-06 4.175e-06 0.00% Other 5645 0.02875 0.0404 0.06133 0.61% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.04 10.04 10.04 100.00% Amr::coarseTimeStep() 1 7.407 7.407 7.407 73.77% Amr::timeStep() 1 7.357 7.357 7.357 73.27% Castro::advance() 1 7.29 7.293 7.296 72.67% Castro::subcycle_advance_ctu() 1 7.217 7.235 7.244 72.14% Castro::do_advance_ctu() 11 7.193 7.21 7.219 71.90% Castro::construct_ctu_hydro_source() 11 3.394 3.413 3.444 34.30% Castro::clean_state() 62 2.954 3.065 3.277 32.64% Castro::computeTemp() 63 1.929 2.023 2.253 22.44% Amr::InitAmr() 1 1.715 1.715 1.715 17.09% Castro::reset_internal_energy(MultiFab) 63 1.184 1.266 1.471 14.65% Amr::init() 1 0.9167 0.9167 0.9167 9.13% Gravity::solve_phi_with_mlmg() 10 0.8088 0.8313 0.8539 8.50% Amr::initialInit() 1 0.8533 0.8533 0.8534 8.50% Castro::initialize_do_advance() 11 0.76 0.7929 0.8185 8.15% Castro::evaluate_source_change() 144 0.724 0.7415 0.7982 7.95% Castro::do_old_sources() 11 0.7014 0.7581 0.7881 7.85% Castro::enforce_min_density() 62 0.613 0.6561 0.7584 7.55% Castro::construct_new_gravity() 8 0.7087 0.7307 0.7534 7.50% Gravity::solve_for_phi() 8 0.6912 0.7137 0.7364 7.33% Gravity::actual_solve_with_mlmg() 10 0.7028 0.7047 0.707 7.04% FabArray::ParallelCopy() 734 0.1152 0.5138 0.6978 6.95% MLMG::solve() 9 0.6756 0.6783 0.6792 6.76% Amr::InitializeInit() 1 0.4349 0.508 0.6179 6.15% Amr::defBaseLevel() 1 0.4349 0.508 0.6179 6.15% FabArray::ParallelCopy_finish() 734 0.02746 0.4248 0.6026 6.00% MLMG::oneIter() 65 0.5885 0.5936 0.5991 5.97% MLMG::mgVcycle() 65 0.5812 0.5861 0.5905 5.88% amrex::Copy() 523 0.5431 0.5538 0.5813 5.79% Castro::normalize_species() 62 0.438 0.4505 0.4894 4.87% Castro::do_new_sources() 8 0.4372 0.46 0.481 4.79% MultiFab::AddProduct() 2633 0.4253 0.4347 0.4684 4.66% AmrLevel::FillPatch() 82 0.1867 0.3522 0.4633 4.61% Castro::reset_internal_energy(Fab) 126 0.4264 0.4339 0.4438 4.42% Amr::FinalizeInit() 1 0.2354 0.3453 0.4184 4.17% MLCellLinOp::smooth() 1560 0.3108 0.3385 0.3986 3.97% FillPatchIterator::Initialize 82 0.1163 0.2805 0.393 3.91% FillPatchIterator::FillFromLevel0() 82 0.106 0.2723 0.3858 3.84% FillPatchSingleLevel 82 0.1057 0.2719 0.3854 3.84% Castro::Castro() 1 0.3578 0.3661 0.3789 3.77% Castro::initMFs() 1 0.3539 0.3622 0.375 3.74% MultiFab::Subtract() 125 0.258 0.2628 0.2696 2.68% Castro::initData() 1 0.07574 0.1407 0.2583 2.57% MLCellLinOp::applyBC() 3630 0.1798 0.2008 0.2363 2.35% FabArray::setVal() 3453 0.2184 0.2232 0.2345 2.34% Amr::bldFineLevels() 1 0.04677 0.1564 0.2281 2.27% check_for_negative_density() 11 0.0841 0.166 0.2253 2.24% MLMG::mgVcycle_up::1 65 0.05392 0.1362 0.2191 2.18% MLPoisson::Fsmooth() 3120 0.1577 0.1732 0.2076 2.07% FabArray::FillBoundary() 3240 0.1496 0.169 0.187 1.86% MLMG::addInterpCorrection() 390 0.03153 0.1506 0.1861 1.85% Amr::grid_places() 2 0.002649 0.1125 0.1856 1.85% AmrMesh::MakeNewGrids() 2 0.002574 0.1124 0.1856 1.85% TagBoxArray::mapPRD 2 0.0008198 0.1107 0.184 1.83% Castro::expand_state() 11 0.09749 0.1372 0.166 1.65% FabArray::sum() 2633 0.1434 0.1483 0.1596 1.59% Castro::post_init() 1 0.142 0.1452 0.148 1.47% Gravity::fill_multipole_BCs() 10 0.1048 0.1259 0.1478 1.47% MLMG::mgVcycle_down::0 65 0.1371 0.1397 0.1426 1.42% MLMG::mgVcycle_up::2 65 0.005412 0.05725 0.1362 1.36% MLMG::mgVcycle_up::0 65 0.1154 0.1185 0.1219 1.21% FillBoundary_finish() 3240 0.08888 0.1036 0.12 1.20% MLCellLinOp::apply() 501 0.08434 0.09168 0.1064 1.06% Amr::writePlotFile() 2 0.1056 0.1057 0.1057 1.05% Gravity::multilevel_solve_for_new_phi() 1 0.1056 0.1056 0.1057 1.05% Gravity::actual_multilevel_solve() 1 0.1055 0.1055 0.1056 1.05% FabArray::ParallelCopy_nowait() 734 0.07958 0.08715 0.1009 1.01% Castro::post_regrid() 2 0.08479 0.08761 0.09073 0.90% MLMG::mgVcycle_up::3 65 0.005379 0.01187 0.08811 0.88% Castro::derive() 18 0.06891 0.07114 0.08279 0.82% AmrLevel::derive() 18 0.06887 0.0711 0.08274 0.82% MLMG:computeResOfCorrection() 390 0.05888 0.06407 0.07569 0.75% MLCellLinOp::correctionResidual() 390 0.05817 0.06335 0.07495 0.75% FillBoundary_nowait() 3240 0.04538 0.05826 0.07456 0.74% Castro::initialize_advance() 1 0.0493 0.05831 0.07353 0.73% MLMG::mgVcycle_down::1 65 0.06037 0.06383 0.06723 0.67% Castro::post_timestep() 1 0.06033 0.06339 0.06646 0.66% Castro::construct_new_source() 48 0.05803 0.05883 0.06254 0.62% MLCellLinOp::solutionResidual() 75 0.0497 0.05601 0.05996 0.60% MLPoisson::Fapply() 501 0.04763 0.05009 0.05666 0.56% MLMG::computeResidual() 65 0.04361 0.04582 0.04745 0.47% Amr::regrid() 1 0.04356 0.04496 0.04649 0.46% Castro::construct_old_source() 66 0.04333 0.04416 0.0453 0.45% Castro::create_source_corrector() 8 0.0149 0.02932 0.04145 0.41% Castro::construct_old_gravity() 11 0.02339 0.03265 0.04033 0.40% Gravity::get_old_grav_vector() 11 0.02256 0.03139 0.04008 0.40% FabArray::Saxpy() 133 0.03549 0.03649 0.04 0.40% MLMG::mgVcycle_down::2 65 0.006312 0.02123 0.03781 0.38% MLMG::mgVcycle_bottom 65 0.0001816 0.002386 0.03537 0.35% MLMG::actualBottomSolve() 65 6.757e-05 0.002269 0.03525 0.35% MLCGSolver::bicgstab 4 0 0.002156 0.03449 0.34% VisMF::Write(FabArray) 4 0.01118 0.01322 0.03426 0.34% Castro::construct_new_gravity_source() 8 0.03235 0.03264 0.0338 0.34% FabArray::Xpay() 480 0.02597 0.02729 0.03042 0.30% MLMG::mgVcycle_down::3 65 0.006023 0.009114 0.02949 0.29% amrex::average_down 390 0.01635 0.02011 0.02799 0.28% Castro::retry_advance_ctu() 11 0.0221 0.02292 0.02703 0.27% Castro::construct_old_gravity_source() 11 0.02526 0.02574 0.02673 0.27% Castro::check_for_nan() 19 0.0229 0.02324 0.02446 0.24% MultiFab::contains_nan() 19 0.02285 0.02318 0.0244 0.24% MLMG::ResNormInf() 74 0.0119 0.01637 0.02289 0.23% Castro::sum_integrated_quantities() 2 0.02267 0.02281 0.02282 0.23% MLPoisson::define() 10 0.01497 0.01757 0.0225 0.22% Gravity::get_new_grav_vector() 9 0.01752 0.01845 0.02196 0.22% Castro::problem_post_timestep() 2 0.01113 0.01628 0.01962 0.20% Castro::wd_update() 2 0.01112 0.01626 0.0196 0.20% Castro::apply_source_to_state() 19 0.01468 0.01498 0.01758 0.18% MLMG::mgVcycle_down::4 65 0.006004 0.006843 0.01718 0.17% Castro::react_state() 19 0.01499 0.01517 0.01554 0.15% Gravity::test_composite_phi() 1 0.01405 0.01407 0.01409 0.14% StateDataPhysBCFunct::() 82 0.0003395 0.007374 0.01408 0.14% StateData::FillBoundary(geom) 62 0 0.006958 0.01362 0.14% MLMG::mgVcycle_up::4 65 0.005449 0.006053 0.01337 0.13% MLCellLinOp::defineAuxData() 10 0.009743 0.01064 0.01337 0.13% Castro::save_data_for_retry() 3 0.01101 0.01137 0.0126 0.13% MLMG::computeMLResidual() 9 0.005315 0.009628 0.01229 0.12% MLMG::prepareForSolve() 9 0.01003 0.01045 0.01222 0.12% MLMG::mgVcycle_down::5 65 0.005976 0.006448 0.01133 0.11% MLMG::mgVcycle_up::5 65 0.005488 0.005851 0.009642 0.10% FabArray::mult() 40 0.007846 0.008003 0.009194 0.09% Gravity::update_max_rhs() 2 0.002755 0.005882 0.008717 0.09% amrex::Add() 65 0.00705 0.007288 0.008361 0.08% Amr::writeSmallPlotFile() 2 0.00829 0.008293 0.008297 0.08% FabArray::norminf() 143 0.005528 0.005829 0.007757 0.08% Castro::volWgtSum() 46 0.007399 0.007441 0.007516 0.07% Other 5645 0.0467 0.06361 0.09352 0.93% -------------------------------------------------------------------------------------------- Unused ParmParse Variables: [TOP]::mg.v(nvals = 1) :: [0] Pinned Memory Usage: --------------------------------------------------------------------------------------------------------------------- Name Nalloc Nfree AvgMem min AvgMem avg AvgMem max MaxMem min MaxMem avg MaxMem max --------------------------------------------------------------------------------------------------------------------- The_Pinned_Arena::Initialize() 16 16 16 B 21 B 30 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-18-g8eff86d32d17) finalized