Initializing AMReX (24.06-24-g2f93525639da)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.06-24-g2f93525639da) initialized Starting run at 07:21:09 UTC on 2024-06-28. Successfully read inputs file ... Castro git describe: 24.06-31-gc47966edd AMReX git describe: 24.06-24-g2f9352563 Microphysics git describe: 24.06-36-g01c244ad 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.167535497 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.00113003 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010569805 Gravity::multilevel_solve_for_new_phi() time = 0.098803831 Gravity::fill_multipole_BCs() time = 0.010496717 ... 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.049247054 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007152598 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.003045255 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.063630997 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 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.003888902 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.062960401 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.004130664 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.061132198 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.003635411 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.060286552 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.312472271 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015308322 Gravity::solve_for_phi() time = 0.085742798 on level 0 Castro::construct_new_gravity() time = 0.087931243 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.060211925 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.002729611 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.059462994 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.329902664 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015673646 Gravity::solve_for_phi() time = 0.086338396 on level 0 Castro::construct_new_gravity() time = 0.088524673 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.060137279 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.002195465 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.060718352 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.320291121 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015512093 Gravity::solve_for_phi() time = 0.086999653 on level 0 Castro::construct_new_gravity() time = 0.089247249 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.060286636 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.002408273 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.060086227 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.321130291 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016691483 Gravity::solve_for_phi() time = 0.087729658 on level 0 Castro::construct_new_gravity() time = 0.089938591 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.060245076 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.002247982 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.06079968 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.315776989 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015472363 Gravity::solve_for_phi() time = 0.085305159 on level 0 Castro::construct_new_gravity() time = 0.087460565 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.06001424 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.002861602 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.058617758 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.368093411 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.024964019 Gravity::solve_for_phi() time = 0.098533153 on level 0 Castro::construct_new_gravity() time = 0.100587749 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.060646802 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.003550641 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.061916348 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.322328341 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01602134 Gravity::solve_for_phi() time = 0.087732708 on level 0 Castro::construct_new_gravity() time = 0.089844043 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.060614111 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.002760271 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.061845479 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.322860688 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017551338 Gravity::solve_for_phi() time = 0.089559621 on level 0 Castro::construct_new_gravity() time = 0.091625738 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.059980283 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.004651156917 [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.105211448 [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.04935334 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.007037102 seconds Ending run at 07:21:19 UTC on 2024-06-28. Run time = 9.593288791 Run time without initialization = 7.162742035 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.593 ... 9.593 ... 9.593 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8319 2.8778 2.9433 30.68% Amr::InitAmr() 1 1.544 1.544 1.544 16.10% Castro::computeTemp() 63 0.6924 0.7528 0.8393 8.75% amrex::Copy() 523 0.536 0.5414 0.5463 5.69% MultiFab::AddProduct() 2633 0.4188 0.4278 0.4935 5.14% Castro::normalize_species() 62 0.4463 0.4513 0.4553 4.75% Castro::reset_internal_energy(Fab) 126 0.4311 0.436 0.4405 4.59% FabArray::ParallelCopy_finish() 734 0.0231 0.3136 0.4321 4.50% Castro::initMFs() 1 0.3577 0.3637 0.377 3.93% MultiFab::Subtract() 125 0.2671 0.27 0.2732 2.85% FabArray::setVal() 3453 0.2212 0.2229 0.2264 2.36% Castro::reset_internal_energy(MultiFab) 63 0.01818 0.05486 0.2172 2.26% MLPoisson::Fsmooth() 3120 0.1631 0.1756 0.2123 2.21% check_for_negative_density() 11 0.04785 0.1102 0.1762 1.84% Castro::initData() 1 0.003345 0.02784 0.1629 1.70% Gravity::fill_multipole_BCs() 10 0.1049 0.1325 0.1564 1.63% FabArray::sum() 2633 0.1416 0.1433 0.1451 1.51% FabArray::ParallelCopy_nowait() 734 0.06526 0.07248 0.08105 0.84% FillBoundary_finish() 3240 0.04707 0.05713 0.0699 0.73% FillBoundary_nowait() 3240 0.04099 0.05342 0.06721 0.70% MLPoisson::Fapply() 501 0.04808 0.05039 0.05796 0.60% MLCellLinOp::applyBC() 3630 0.02498 0.03253 0.05426 0.57% Castro::enforce_min_density() 62 0.01405 0.02182 0.05166 0.54% Castro::do_old_sources() 11 0.0003006 0.007978 0.04458 0.46% AmrLevel::derive() 18 0.04145 0.04186 0.04211 0.44% Castro::construct_new_gravity_source() 8 0.03208 0.0334 0.0419 0.44% FabArray::Saxpy() 133 0.0355 0.03617 0.03957 0.41% Castro::do_new_sources() 8 0.0001705 0.006164 0.03746 0.39% FabArray::Xpay() 480 0.02561 0.02682 0.03041 0.32% VisMF::Write(FabArray) 4 0.0152 0.01666 0.02674 0.28% Castro::construct_old_gravity_source() 11 0.02494 0.02542 0.02651 0.28% MultiFab::contains_nan() 19 0.02307 0.02345 0.02383 0.25% MLMG::addInterpCorrection() 390 0.01705 0.01832 0.0218 0.23% Amr::writePlotFile() 2 0.008737 0.01742 0.01871 0.20% amrex::average_down 390 0.01252 0.0136 0.01662 0.17% Castro::wd_update() 2 0.005059 0.01058 0.016 0.17% StateData::FillBoundary(geom) 62 0 0.007134 0.01389 0.14% Castro::evaluate_source_change() 144 0.009416 0.009613 0.0098 0.10% MLCellLinOp::defineAuxData() 10 0.007968 0.0083 0.009591 0.10% Amr::writeSmallPlotFile() 2 0.007076 0.008988 0.009467 0.10% FabArray::FillBoundary() 3240 0.006941 0.007235 0.008663 0.09% Castro::construct_old_gravity() 11 0.0002088 0.001677 0.008573 0.09% MLCGSolver::bicgstab 4 0 0.0005214 0.008342 0.09% Castro::construct_new_source() 48 0.007311 0.007474 0.008267 0.09% FabArray::mult() 40 0.008026 0.008158 0.008264 0.09% FabArray::norminf() 143 0.005516 0.005803 0.007883 0.08% Castro::volWgtSum() 46 0.007237 0.007285 0.007351 0.08% FabArrayBase::getFB() 3240 0.005227 0.005424 0.006816 0.07% MLCellLinOp::smooth() 1560 0.005699 0.005856 0.006301 0.07% FillPatchIterator::Initialize 82 0.005759 0.005995 0.006217 0.06% Castro::construct_old_source() 66 0.005448 0.005564 0.005758 0.06% MLMG::prepareForSolve() 9 0.005145 0.005271 0.005522 0.06% MLMG::ResNormInf() 74 0.002079 0.003216 0.004465 0.05% Gravity::actual_solve_with_mlmg() 10 0.003538 0.003654 0.00414 0.04% FabArrayBase::getCPC() 1022 0.003788 0.003886 0.004051 0.04% Gravity::get_old_grav_vector() 11 0.003623 0.003662 0.003719 0.04% MLMG::solve() 9 0.0004431 0.001685 0.003588 0.04% Gravity::get_new_grav_vector() 9 0.003058 0.003122 0.003326 0.03% Castro::initialize_do_advance() 11 0.00279 0.002995 0.003301 0.03% MLCellLinOp::apply() 501 0.001018 0.001122 0.00239 0.02% FabArray::ParallelCopy() 734 0.001802 0.001879 0.002068 0.02% Castro::sum_integrated_quantities() 2 0.001617 0.001841 0.002038 0.02% main() 1 0.0009283 0.001202 0.001449 0.02% AmrLevel::FillPatch() 82 0.001142 0.001173 0.00125 0.01% MLCellLinOp::correctionResidual() 390 0.0009419 0.0009589 0.0009987 0.01% MLMG::mgVcycle() 65 0.0007953 0.0008284 0.0009026 0.01% Castro::save_data_for_retry() 3 0.0007666 0.0007847 0.0008107 0.01% MLMG:computeResOfCorrection() 390 0.0006887 0.0007116 0.0007517 0.01% Castro::subcycle_advance_ctu() 1 0.0001007 0.000149 0.0006786 0.01% Gravity::solve_for_phi() 8 0.000374 0.0004499 0.0005737 0.01% MLMG::actualBottomSolve() 65 6.794e-05 0.0001028 0.0005725 0.01% Castro::Castro() 1 0.0004774 0.0005166 0.000565 0.01% AmrMesh::MakeNewGrids() 2 0.0005208 0.0005361 0.0005627 0.01% StateDataPhysBCFunct::() 82 0.000372 0.0004389 0.0004944 0.01% Castro::construct_new_gravity() 8 0.0001434 0.0002179 0.0004038 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003642 0.0003822 0.0004038 0.00% Castro::clean_state() 62 0.0003424 0.0003546 0.000377 0.00% MLMG::mgVcycle_down::3 65 0.0002936 0.0003176 0.0003675 0.00% MLMG::mgVcycle_down::2 65 0.000315 0.000341 0.0003649 0.00% MLMG::mgVcycle_down::0 65 0.0003341 0.0003424 0.0003529 0.00% MLMG::mgVcycle_down::1 65 0.00033 0.0003376 0.0003488 0.00% MLMG::mgVcycle_down::4 65 0.0002853 0.0002965 0.0003474 0.00% MLMG::mgVcycle_down::5 65 0.0002934 0.0003072 0.0003459 0.00% FillPatchSingleLevel 82 0.0002588 0.0002654 0.0002838 0.00% Castro::do_advance_ctu() 11 0.0002342 0.000249 0.0002742 0.00% Castro::problem_post_timestep() 2 1.12e-05 2.827e-05 0.0002638 0.00% MLMG::mgVcycle_up::1 65 0.0002356 0.0002431 0.0002585 0.00% MLMG::mgVcycle_up::5 65 0.0002025 0.0002103 0.0002528 0.00% MLMG::mgVcycle_up::0 65 0.0002357 0.0002418 0.0002513 0.00% MLCellLinOp::solutionResidual() 75 0.0002257 0.0002331 0.0002432 0.00% MLMG::mgVcycle_up::2 65 0.0001978 0.0002157 0.000237 0.00% MLMG::mgVcycle_up::3 65 0.0001923 0.0002006 0.0002278 0.00% MLMG::mgVcycle_up::4 65 0.0001955 0.0002017 0.0002238 0.00% MLMG::oneIter() 65 0.0001925 0.0001988 0.0002052 0.00% Amr::coarseTimeStep() 1 0.0001589 0.0001813 0.0001874 0.00% Castro::initialize_advance() 1 0.000124 0.0001431 0.00017 0.00% Amr::grid_places() 2 7.411e-05 9.444e-05 0.0001412 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001162 0.0001232 0.0001386 0.00% MLMG::computeResidual() 65 0.0001205 0.0001248 0.0001383 0.00% Castro::retry_advance_ctu() 11 5.209e-05 6.483e-05 0.000137 0.00% Gravity::test_composite_phi() 1 9.266e-05 0.0001011 0.0001294 0.00% MLMG::mgVcycle_bottom 65 0.0001135 0.0001167 0.000128 0.00% TagBoxArray::mapPRD 2 0.0001104 0.0001161 0.0001254 0.00% Amr::defBaseLevel() 1 9.835e-05 0.0001087 0.0001197 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.243e-05 6.996e-05 0.0001102 0.00% MLPoisson::define() 10 6.837e-05 7.461e-05 8.365e-05 0.00% Castro::check_for_nan() 19 5.385e-05 6.01e-05 8.257e-05 0.00% Castro::react_state() 19 5.324e-05 5.937e-05 7.359e-05 0.00% Amr::FinalizeInit() 1 8.257e-06 1.408e-05 6.495e-05 0.00% Castro::apply_source_to_state() 19 4.971e-05 5.506e-05 6.086e-05 0.00% Castro::create_source_corrector() 8 3.57e-05 3.835e-05 4.659e-05 0.00% Castro::derive() 18 4.113e-05 4.343e-05 4.569e-05 0.00% Castro::expand_state() 11 3.299e-05 3.798e-05 4.375e-05 0.00% Amr::regrid() 1 1.04e-05 1.624e-05 3.304e-05 0.00% Gravity::actual_multilevel_solve() 1 2.754e-05 2.932e-05 3.239e-05 0.00% MLMG::computeMLResidual() 9 1.969e-05 2.186e-05 2.568e-05 0.00% Amr::timeStep() 1 1.252e-05 1.637e-05 2.114e-05 0.00% Castro::advance() 1 8.018e-06 1.051e-05 2.011e-05 0.00% Castro::post_init() 1 1.159e-05 1.245e-05 1.292e-05 0.00% Amr::InitializeInit() 1 6.751e-06 8.074e-06 1.249e-05 0.00% Castro::post_regrid() 2 1.04e-05 1.1e-05 1.184e-05 0.00% Castro::post_timestep() 1 5.481e-06 6.387e-06 1.161e-05 0.00% Amr::bldFineLevels() 1 6.713e-06 7.639e-06 9.721e-06 0.00% Amr::init() 1 6.685e-06 7.609e-06 8.614e-06 0.00% Amr::initialInit() 1 2.997e-06 3.259e-06 3.557e-06 0.00% Other 1450 0.02824 0.03947 0.06215 0.65% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.593 9.593 9.593 100.00% Amr::coarseTimeStep() 1 7.162 7.162 7.162 74.66% Amr::timeStep() 1 7.105 7.105 7.105 74.06% Castro::advance() 1 7.038 7.042 7.045 73.44% Castro::subcycle_advance_ctu() 1 6.99 6.993 6.995 72.92% Castro::do_advance_ctu() 11 6.966 6.969 6.972 72.67% Castro::construct_ctu_hydro_source() 11 3.341 3.357 3.373 35.15% Castro::clean_state() 62 2.936 3.007 3.117 32.49% Castro::computeTemp() 63 1.918 1.988 2.165 22.57% Amr::InitAmr() 1 1.544 1.544 1.544 16.10% Castro::reset_internal_energy(MultiFab) 63 1.185 1.235 1.383 14.42% Amr::init() 1 0.8855 0.8856 0.8856 9.23% Amr::initialInit() 1 0.8282 0.8283 0.8285 8.64% Gravity::solve_phi_with_mlmg() 10 0.7636 0.7925 0.8187 8.53% Castro::evaluate_source_change() 144 0.7165 0.7289 0.7942 8.28% Castro::initialize_do_advance() 11 0.7456 0.7622 0.7744 8.07% Castro::construct_new_gravity() 8 0.6709 0.6991 0.7252 7.56% Gravity::solve_for_phi() 8 0.653 0.6818 0.7081 7.38% Castro::do_old_sources() 11 0.6432 0.6582 0.6688 6.97% Gravity::actual_solve_with_mlmg() 10 0.6563 0.6582 0.6606 6.89% Castro::enforce_min_density() 62 0.6112 0.627 0.6595 6.87% MLMG::solve() 9 0.6285 0.631 0.6317 6.58% Amr::InitializeInit() 1 0.4401 0.5022 0.6065 6.32% Amr::defBaseLevel() 1 0.4401 0.5022 0.6065 6.32% MLMG::oneIter() 65 0.5624 0.5636 0.5649 5.89% MLMG::mgVcycle() 65 0.5551 0.5562 0.5574 5.81% amrex::Copy() 523 0.536 0.5414 0.5463 5.69% FabArray::ParallelCopy() 734 0.1154 0.4043 0.5219 5.44% MultiFab::AddProduct() 2633 0.4188 0.4278 0.4935 5.14% Castro::do_new_sources() 8 0.4441 0.4628 0.4822 5.03% Castro::normalize_species() 62 0.4463 0.4513 0.4553 4.75% Castro::reset_internal_energy(Fab) 126 0.4311 0.436 0.4405 4.59% FabArray::ParallelCopy_finish() 734 0.0231 0.3136 0.4321 4.50% Amr::FinalizeInit() 1 0.222 0.3261 0.3882 4.05% Castro::Castro() 1 0.3619 0.3678 0.3814 3.98% Castro::initMFs() 1 0.3578 0.3638 0.3771 3.93% MLCellLinOp::smooth() 1560 0.2821 0.3072 0.3713 3.87% AmrLevel::FillPatch() 82 0.1886 0.2457 0.2749 2.87% MultiFab::Subtract() 125 0.2671 0.27 0.2732 2.85% Castro::initData() 1 0.07591 0.133 0.2427 2.53% FabArray::setVal() 3453 0.2212 0.2229 0.2264 2.36% MLMG::mgVcycle_up::1 65 0.0485 0.1316 0.215 2.24% MLPoisson::Fsmooth() 3120 0.1631 0.1756 0.2123 2.21% Amr::bldFineLevels() 1 0.0455 0.1496 0.2115 2.20% FillPatchIterator::Initialize 82 0.1184 0.1755 0.2049 2.14% FillPatchIterator::FillFromLevel0() 82 0.1068 0.1666 0.1968 2.05% FillPatchSingleLevel 82 0.1064 0.1662 0.1964 2.05% MLCellLinOp::applyBC() 3630 0.139 0.1562 0.1919 2.00% MLMG::addInterpCorrection() 390 0.03274 0.1539 0.1894 1.97% check_for_negative_density() 11 0.04785 0.1102 0.1762 1.84% Amr::grid_places() 2 0.002304 0.1066 0.1687 1.76% AmrMesh::MakeNewGrids() 2 0.002163 0.1065 0.1686 1.76% TagBoxArray::mapPRD 2 0.0007063 0.1042 0.1671 1.74% Gravity::fill_multipole_BCs() 10 0.106 0.1336 0.1575 1.64% FabArray::sum() 2633 0.1416 0.1433 0.1451 1.51% FabArray::FillBoundary() 3240 0.1084 0.1237 0.1442 1.50% MLMG::mgVcycle_up::2 65 0.005521 0.05891 0.1403 1.46% MLMG::mgVcycle_down::0 65 0.1331 0.1345 0.1361 1.42% Castro::post_init() 1 0.1307 0.134 0.1347 1.40% Castro::expand_state() 11 0.0979 0.1158 0.1316 1.37% MLMG::mgVcycle_up::0 65 0.1002 0.1004 0.1006 1.05% Amr::writePlotFile() 2 0.09914 0.09924 0.09925 1.03% MLCellLinOp::apply() 501 0.07609 0.08104 0.09884 1.03% Gravity::multilevel_solve_for_new_phi() 1 0.09875 0.09878 0.09882 1.03% Gravity::actual_multilevel_solve() 1 0.09869 0.09871 0.09875 1.03% FabArray::ParallelCopy_nowait() 734 0.08132 0.08883 0.09782 1.02% MLMG::mgVcycle_up::3 65 0.005472 0.01208 0.09137 0.95% Castro::post_regrid() 2 0.08484 0.08555 0.08885 0.93% MLMG:computeResOfCorrection() 390 0.05746 0.06202 0.07402 0.77% FillBoundary_nowait() 3240 0.04662 0.05929 0.07335 0.76% MLCellLinOp::correctionResidual() 390 0.05676 0.06131 0.07327 0.76% Castro::derive() 18 0.06953 0.07025 0.07084 0.74% AmrLevel::derive() 18 0.06948 0.07021 0.0708 0.74% FillBoundary_finish() 3240 0.04707 0.05713 0.0699 0.73% Castro::construct_new_source() 48 0.05769 0.05945 0.06757 0.70% Castro::post_timestep() 1 0.06004 0.06317 0.06726 0.70% MLMG::mgVcycle_down::1 65 0.05752 0.05958 0.06181 0.64% MLPoisson::Fapply() 501 0.04808 0.05039 0.05796 0.60% Castro::initialize_advance() 1 0.04804 0.04889 0.04996 0.52% MLCellLinOp::solutionResidual() 75 0.04413 0.04692 0.0493 0.51% Amr::regrid() 1 0.0438 0.0441 0.04569 0.48% Castro::construct_old_source() 66 0.04314 0.04385 0.04527 0.47% Castro::construct_new_gravity_source() 8 0.03208 0.0334 0.0419 0.44% FabArray::Saxpy() 133 0.0355 0.03617 0.03957 0.41% MLMG::computeResidual() 65 0.03823 0.03846 0.03914 0.41% MLMG::mgVcycle_bottom 65 0.0001814 0.002486 0.03697 0.39% MLMG::actualBottomSolve() 65 6.794e-05 0.002369 0.03684 0.38% MLMG::mgVcycle_down::2 65 0.006474 0.02085 0.03662 0.38% MLCGSolver::bicgstab 4 0 0.002251 0.03602 0.38% Castro::construct_old_gravity() 11 0.02311 0.02839 0.03463 0.36% Castro::create_source_corrector() 8 0.01526 0.02421 0.03412 0.36% Gravity::get_old_grav_vector() 11 0.02133 0.02672 0.03086 0.32% MLMG::mgVcycle_down::3 65 0.006137 0.009303 0.03071 0.32% FabArray::Xpay() 480 0.02561 0.02682 0.03041 0.32% amrex::average_down 390 0.01678 0.02005 0.02847 0.30% VisMF::Write(FabArray) 4 0.0152 0.01666 0.02674 0.28% Castro::construct_old_gravity_source() 11 0.02494 0.02542 0.02651 0.28% Castro::retry_advance_ctu() 11 0.022 0.02241 0.02406 0.25% Castro::check_for_nan() 19 0.02313 0.02351 0.0239 0.25% MultiFab::contains_nan() 19 0.02307 0.02345 0.02383 0.25% MLPoisson::define() 10 0.01589 0.01815 0.02137 0.22% Castro::sum_integrated_quantities() 2 0.02065 0.02081 0.02083 0.22% Castro::problem_post_timestep() 2 0.009184 0.01429 0.01948 0.20% Castro::wd_update() 2 0.00892 0.01427 0.01946 0.20% Gravity::get_new_grav_vector() 9 0.01763 0.01793 0.01831 0.19% MLMG::mgVcycle_down::4 65 0.006093 0.006959 0.01791 0.19% Castro::react_state() 19 0.01515 0.01551 0.01605 0.17% Castro::apply_source_to_state() 19 0.01467 0.01485 0.01511 0.16% Amr::writeSmallPlotFile() 2 0.01461 0.01461 0.01462 0.15% StateDataPhysBCFunct::() 82 0.000372 0.007573 0.01438 0.15% Gravity::test_composite_phi() 1 0.01396 0.014 0.01403 0.15% MLMG::mgVcycle_up::4 65 0.005573 0.006177 0.01399 0.15% StateData::FillBoundary(geom) 62 0 0.007134 0.01389 0.14% MLCellLinOp::defineAuxData() 10 0.01023 0.01099 0.01336 0.14% Castro::save_data_for_retry() 3 0.01082 0.01109 0.01241 0.13% MLMG::mgVcycle_down::5 65 0.006088 0.006554 0.01194 0.12% MLMG::prepareForSolve() 9 0.01028 0.01064 0.01176 0.12% MLMG::mgVcycle_up::5 65 0.005608 0.005986 0.01018 0.11% MLMG::computeMLResidual() 9 0.005484 0.007866 0.01001 0.10% MLMG::ResNormInf() 74 0.006977 0.008077 0.009264 0.10% FabArray::mult() 40 0.008026 0.008158 0.008264 0.09% FabArray::norminf() 143 0.005516 0.005803 0.007883 0.08% FabArrayBase::getCPC() 1022 0.00522 0.005927 0.007789 0.08% FabArrayBase::getFB() 3240 0.005487 0.005875 0.007693 0.08% Castro::volWgtSum() 46 0.007237 0.007285 0.007351 0.08% Other 1450 0.04664 0.06557 0.09572 1.00% -------------------------------------------------------------------------------------------- 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 22 B 31 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.06-24-g2f93525639da) finalized