Initializing AMReX (24.04-22-g64d2360b209c)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.04-22-g64d2360b209c) initialized Starting run at 01:18:13 UTC on 2024-04-20. Successfully read inputs file ... Castro git describe: 24.04-23-gc71c4dfd3 AMReX git describe: 24.04-22-g64d2360b2 Microphysics git describe: 24.04-9-g4ad225b7 reading extern runtime parameters ... 13 Species: He4 C12 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Ni56 Successfully read inputs file ... Creating primary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Creating secondary with CO core with mass fractions C = 0.5and O = 0.5 and a He shell of solar mass =0. Generated initial model for primary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Generated initial model for secondary WD of mass 0.64 solar masses, central density 4.500e+06 g cm**-3, and radius 8.734e+08 cm. Castro::numpts_1d at level 0 is 294 Setting the gravity type to PoissonGrav Initializing the data at level 0 Done initializing the level 0 data grid_places() time: 0.167551192 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001050365 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010549278 Gravity::multilevel_solve_for_new_phi() time = 0.101225094 Gravity::fill_multipole_BCs() time = 0.01057836 ... 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.054046179 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.006658677 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.002231589 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.069005009 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.004157375 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.071726527 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.003896423 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.06651485 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.003249836 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.059118746 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.311964596 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014187845 Gravity::solve_for_phi() time = 0.085010637 on level 0 Castro::construct_new_gravity() time = 0.087170568 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.060319774 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.002286308 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.060190604 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.319127742 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.014424187 Gravity::solve_for_phi() time = 0.085434782 on level 0 Castro::construct_new_gravity() time = 0.087559512 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.060169954 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.002216671 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.060463091 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.337028214 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016766908 Gravity::solve_for_phi() time = 0.091765399 on level 0 Castro::construct_new_gravity() time = 0.093893653 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.062880176 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.006023382 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.066740969 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.339080715 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016582803 Gravity::solve_for_phi() time = 0.088285307 on level 0 Castro::construct_new_gravity() time = 0.090374109 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.063732079 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.006239076 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.062594278 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.338176845 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018942563 Gravity::solve_for_phi() time = 0.089604192 on level 0 Castro::construct_new_gravity() time = 0.091766224 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.059760126 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.002291542 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.060133347 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.314962839 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015004883 Gravity::solve_for_phi() time = 0.085036872 on level 0 Castro::construct_new_gravity() time = 0.087127737 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.060028271 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.002318016 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.059540213 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.316615435 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01517175 Gravity::solve_for_phi() time = 0.084305576 on level 0 Castro::construct_new_gravity() time = 0.086499049 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.059985545 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.002455701 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.063639435 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.353011808 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017021942 Gravity::solve_for_phi() time = 0.087139894 on level 0 Castro::construct_new_gravity() time = 0.089331247 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.059962141 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.004586047608 [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.203707189 [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.050483911 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.006862529 seconds Ending run at 01:18:23 UTC on 2024-04-20. Run time = 9.538246948 Run time without initialization = 7.262033744 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.538 ... 9.538 ... 9.538 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8301 2.8722 2.9508 30.94% Amr::InitAmr() 1 1.376 1.377 1.377 14.43% Castro::computeTemp() 63 0.6852 0.7393 0.8232 8.63% amrex::Copy() 523 0.581 0.5867 0.6088 6.38% Castro::normalize_species() 62 0.4466 0.4509 0.4724 4.95% FabArray::ParallelCopy_finish() 734 0.02112 0.3273 0.4641 4.87% MultiFab::AddProduct() 2633 0.4372 0.4414 0.4551 4.77% Castro::reset_internal_energy(Fab) 126 0.4246 0.4284 0.4346 4.56% Castro::initMFs() 1 0.354 0.3634 0.3765 3.95% MultiFab::Subtract() 125 0.2708 0.2733 0.2787 2.92% FabArray::setVal() 3453 0.2271 0.2297 0.2407 2.52% check_for_negative_density() 11 0.06015 0.1625 0.238 2.49% MLPoisson::Fsmooth() 3120 0.1605 0.1739 0.2163 2.27% Castro::reset_internal_energy(MultiFab) 63 0.01802 0.05371 0.2127 2.23% Castro::initData() 1 0.003365 0.02716 0.152 1.59% FabArray::sum() 2633 0.1454 0.1469 0.1516 1.59% Gravity::fill_multipole_BCs() 10 0.1043 0.1267 0.1453 1.52% FabArray::ParallelCopy_nowait() 734 0.06387 0.07036 0.07914 0.83% FillBoundary_finish() 3240 0.04759 0.06093 0.06875 0.72% FillBoundary_nowait() 3240 0.03984 0.05262 0.0667 0.70% MLPoisson::Fapply() 501 0.04892 0.05114 0.0596 0.62% Castro::enforce_min_density() 62 0.01394 0.02391 0.0555 0.58% Castro::do_old_sources() 11 0.0002855 0.01256 0.05485 0.58% MLCellLinOp::applyBC() 3630 0.02433 0.03222 0.05454 0.57% AmrLevel::derive() 18 0.04067 0.04098 0.04145 0.43% FabArray::Saxpy() 133 0.03538 0.03596 0.04096 0.43% Castro::construct_new_gravity_source() 8 0.03265 0.03293 0.03334 0.35% Castro::do_new_sources() 8 0.00017 0.00584 0.03198 0.34% FabArray::Xpay() 480 0.0259 0.02706 0.03092 0.32% VisMF::Write(FabArray) 4 0.01399 0.01562 0.02787 0.29% Castro::construct_old_gravity_source() 11 0.02546 0.02566 0.02582 0.27% Amr::writePlotFile() 2 0.008653 0.02199 0.02389 0.25% MultiFab::contains_nan() 19 0.02318 0.02333 0.0237 0.25% MLMG::addInterpCorrection() 390 0.01701 0.01824 0.02196 0.23% TagBoxArray::collate() 2 9.587e-05 0.003441 0.02081 0.22% amrex::average_down 390 0.01253 0.01354 0.01666 0.17% StateData::FillBoundary(geom) 62 0 0.007357 0.01413 0.15% Castro::construct_old_gravity() 11 0.0002295 0.001547 0.01072 0.11% Castro::wd_update() 2 0.004627 0.008151 0.01031 0.11% Castro::evaluate_source_change() 144 0.009337 0.009579 0.01014 0.11% MLCellLinOp::defineAuxData() 10 0.007815 0.008162 0.009347 0.10% Amr::writeSmallPlotFile() 2 0.007169 0.009056 0.009312 0.10% FabArray::mult() 40 0.008399 0.008494 0.008727 0.09% FabArray::FillBoundary() 3240 0.006897 0.007147 0.008622 0.09% MLCGSolver::bicgstab 4 0 0.0005224 0.008358 0.09% FabArray::norminf() 143 0.005633 0.005884 0.007978 0.08% Castro::construct_new_source() 48 0.007291 0.007465 0.007945 0.08% Castro::volWgtSum() 46 0.00738 0.007459 0.007922 0.08% Castro::construct_old_source() 66 0.005488 0.0057 0.007226 0.08% MLCellLinOp::smooth() 1560 0.005864 0.006134 0.006616 0.07% FabArrayBase::getFB() 3240 0.004976 0.005164 0.006431 0.07% Castro::sum_integrated_quantities() 2 0.001654 0.005524 0.006145 0.06% FillPatchIterator::Initialize 82 0.005621 0.005862 0.006095 0.06% MLMG::prepareForSolve() 9 0.005081 0.005258 0.005616 0.06% Gravity::actual_solve_with_mlmg() 10 0.003484 0.003624 0.004076 0.04% FabArrayBase::getCPC() 1022 0.003678 0.003774 0.003893 0.04% Gravity::get_old_grav_vector() 11 0.003598 0.003645 0.003729 0.04% MLMG::ResNormInf() 74 0.001872 0.0026 0.003346 0.04% Gravity::get_new_grav_vector() 9 0.003031 0.003108 0.0033 0.03% Castro::initialize_do_advance() 11 0.002893 0.003022 0.003177 0.03% MLMG::solve() 9 0.0003855 0.001562 0.003121 0.03% MLCellLinOp::apply() 501 0.001051 0.001157 0.002433 0.03% FabArray::ParallelCopy() 734 0.001808 0.001879 0.002062 0.02% main() 1 0.000938 0.001212 0.001316 0.01% AmrLevel::FillPatch() 82 0.001104 0.001148 0.001259 0.01% MLCellLinOp::correctionResidual() 390 0.000931 0.0009527 0.0009836 0.01% MLMG::mgVcycle() 65 0.0007962 0.0008304 0.0009438 0.01% Castro::save_data_for_retry() 3 0.0007955 0.0008313 0.0008742 0.01% MLMG:computeResOfCorrection() 390 0.0007139 0.0007403 0.0007826 0.01% Gravity::solve_for_phi() 8 0.0003815 0.0004842 0.0006814 0.01% Castro::subcycle_advance_ctu() 1 0.0001064 0.0001544 0.0006588 0.01% Castro::Castro() 1 0.0004725 0.0005265 0.0005851 0.01% Castro::construct_new_gravity() 8 0.0001414 0.0002882 0.0005531 0.01% MLMG::actualBottomSolve() 65 6.825e-05 0.0001006 0.0005525 0.01% StateDataPhysBCFunct::() 82 0.0003459 0.0004224 0.0005442 0.01% AmrMesh::MakeNewGrids() 2 0.0004987 0.0005199 0.0005379 0.01% FillPatchIterator::FillFromLevel0() 82 0.0003753 0.0003907 0.0004286 0.00% Castro::clean_state() 62 0.0003324 0.00035 0.0003839 0.00% MLMG::mgVcycle_down::3 65 0.0002951 0.0003156 0.0003642 0.00% MLMG::mgVcycle_down::2 65 0.0003093 0.000335 0.0003631 0.00% MLMG::mgVcycle_down::1 65 0.000332 0.0003434 0.000356 0.00% MLMG::mgVcycle_down::0 65 0.0003336 0.0003398 0.0003502 0.00% MLMG::mgVcycle_down::4 65 0.0002859 0.0003003 0.0003478 0.00% MLMG::mgVcycle_down::5 65 0.0002979 0.0003048 0.000338 0.00% Castro::problem_post_timestep() 2 1.11e-05 2.872e-05 0.0002829 0.00% FillPatchSingleLevel 82 0.0002574 0.000266 0.0002824 0.00% Castro::do_advance_ctu() 11 0.0002189 0.000247 0.0002749 0.00% MLMG::mgVcycle_up::1 65 0.0002335 0.0002429 0.0002596 0.00% MLMG::mgVcycle_up::2 65 0.0001972 0.0002195 0.0002532 0.00% MLMG::mgVcycle_up::0 65 0.0002324 0.0002413 0.0002526 0.00% MLCellLinOp::solutionResidual() 75 0.0002207 0.0002253 0.0002334 0.00% MLMG::mgVcycle_up::5 65 0.0002024 0.0002085 0.0002323 0.00% MLMG::mgVcycle_up::3 65 0.0001956 0.0002038 0.0002269 0.00% MLMG::oneIter() 65 0.0001974 0.0002032 0.0002213 0.00% MLMG::mgVcycle_up::4 65 0.0001988 0.0002032 0.0002212 0.00% Amr::coarseTimeStep() 1 0.000143 0.0001755 0.0001804 0.00% Castro::initialize_advance() 1 0.00013 0.0001516 0.0001748 0.00% Amr::defBaseLevel() 1 9.182e-05 0.0001252 0.0001423 0.00% Castro::retry_advance_ctu() 11 5.725e-05 6.617e-05 0.0001403 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001136 0.0001211 0.0001374 0.00% MLMG::computeResidual() 65 0.0001226 0.000125 0.0001335 0.00% MLMG::mgVcycle_bottom 65 0.0001161 0.0001186 0.0001254 0.00% TagBoxArray::mapPRD 2 0.0001078 0.0001158 0.0001245 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.387e-05 7.514e-05 0.0001238 0.00% Gravity::test_composite_phi() 1 9.261e-05 9.943e-05 0.0001135 0.00% Amr::grid_places() 2 7.639e-05 8.962e-05 0.0001094 0.00% MLPoisson::define() 10 6.599e-05 7.066e-05 7.666e-05 0.00% Castro::check_for_nan() 19 4.994e-05 5.653e-05 6.769e-05 0.00% Castro::react_state() 19 5.094e-05 5.469e-05 6.269e-05 0.00% Castro::apply_source_to_state() 19 5.038e-05 5.476e-05 5.937e-05 0.00% Amr::FinalizeInit() 1 9.123e-06 1.517e-05 5.378e-05 0.00% Castro::create_source_corrector() 8 3.445e-05 3.93e-05 5.353e-05 0.00% Castro::derive() 18 4.092e-05 4.312e-05 4.872e-05 0.00% Castro::expand_state() 11 3.269e-05 3.667e-05 4.234e-05 0.00% Gravity::actual_multilevel_solve() 1 2.918e-05 3.073e-05 3.342e-05 0.00% Amr::regrid() 1 1.095e-05 1.297e-05 3.068e-05 0.00% Amr::timeStep() 1 1.196e-05 1.758e-05 2.975e-05 0.00% MLMG::computeMLResidual() 9 2.172e-05 2.333e-05 2.607e-05 0.00% Castro::advance() 1 7.13e-06 1.093e-05 2.434e-05 0.00% Amr::InitializeInit() 1 7.108e-06 8.688e-06 1.505e-05 0.00% Castro::post_regrid() 2 1.107e-05 1.2e-05 1.458e-05 0.00% Castro::post_init() 1 1.185e-05 1.274e-05 1.412e-05 0.00% Amr::init() 1 6.781e-06 7.772e-06 9.831e-06 0.00% Amr::bldFineLevels() 1 6.85e-06 7.609e-06 9.298e-06 0.00% Castro::post_timestep() 1 5.097e-06 5.837e-06 6.625e-06 0.00% Amr::initialInit() 1 3.173e-06 3.676e-06 6.133e-06 0.00% Other 1448 0.02753 0.03603 0.05707 0.60% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.538 9.538 9.538 100.00% Amr::coarseTimeStep() 1 7.261 7.261 7.261 76.13% Amr::timeStep() 1 7.204 7.204 7.204 75.52% Castro::advance() 1 7.14 7.142 7.145 74.91% Castro::subcycle_advance_ctu() 1 7.089 7.091 7.095 74.38% Castro::do_advance_ctu() 11 7.062 7.066 7.07 74.12% Castro::construct_ctu_hydro_source() 11 3.376 3.406 3.437 36.04% Castro::clean_state() 62 2.976 3.045 3.165 33.18% Castro::computeTemp() 63 1.918 1.999 2.15 22.54% Castro::reset_internal_energy(MultiFab) 63 1.213 1.26 1.41 14.78% Amr::InitAmr() 1 1.377 1.377 1.377 14.43% Amr::init() 1 0.8989 0.899 0.899 9.42% Amr::initialInit() 1 0.8375 0.8375 0.8376 8.78% Gravity::solve_phi_with_mlmg() 10 0.7635 0.7872 0.8081 8.47% Castro::initialize_do_advance() 11 0.7545 0.7811 0.7946 8.33% Castro::evaluate_source_change() 144 0.7435 0.7509 0.7744 8.12% Castro::construct_new_gravity() 8 0.6684 0.6917 0.7125 7.47% Gravity::solve_for_phi() 8 0.6503 0.6741 0.695 7.29% Castro::do_old_sources() 11 0.6461 0.6778 0.6916 7.25% Castro::enforce_min_density() 62 0.6367 0.6522 0.6788 7.12% Gravity::actual_solve_with_mlmg() 10 0.6567 0.6586 0.6609 6.93% MLMG::solve() 9 0.6296 0.6319 0.6326 6.63% amrex::Copy() 523 0.581 0.5867 0.6088 6.38% Amr::InitializeInit() 1 0.4415 0.5012 0.6079 6.37% Amr::defBaseLevel() 1 0.4415 0.5012 0.6079 6.37% MLMG::oneIter() 65 0.5642 0.5649 0.5657 5.93% MLMG::mgVcycle() 65 0.5569 0.5575 0.5583 5.85% FabArray::ParallelCopy() 734 0.1153 0.4172 0.5535 5.80% Castro::do_new_sources() 8 0.4436 0.4639 0.4869 5.10% Castro::normalize_species() 62 0.4466 0.4509 0.4724 4.95% FabArray::ParallelCopy_finish() 734 0.02112 0.3273 0.4641 4.87% MultiFab::AddProduct() 2633 0.4372 0.4414 0.4551 4.77% Castro::reset_internal_energy(Fab) 126 0.4246 0.4284 0.4346 4.56% Amr::FinalizeInit() 1 0.2297 0.3363 0.396 4.15% Castro::Castro() 1 0.3582 0.368 0.381 3.99% Castro::initMFs() 1 0.3541 0.3635 0.3766 3.95% MLCellLinOp::smooth() 1560 0.2824 0.3073 0.3728 3.91% AmrLevel::FillPatch() 82 0.1988 0.2651 0.3033 3.18% MultiFab::Subtract() 125 0.2708 0.2733 0.2787 2.92% Castro::initData() 1 0.07807 0.1319 0.2453 2.57% FabArray::setVal() 3453 0.2271 0.2297 0.2407 2.52% check_for_negative_density() 11 0.06015 0.1625 0.238 2.49% FillPatchIterator::Initialize 82 0.1197 0.1888 0.2275 2.39% FillPatchIterator::FillFromLevel0() 82 0.1086 0.1802 0.2194 2.30% FillPatchSingleLevel 82 0.1082 0.1798 0.219 2.30% MLPoisson::Fsmooth() 3120 0.1605 0.1739 0.2163 2.27% MLMG::mgVcycle_up::1 65 0.04833 0.1312 0.2144 2.25% Amr::bldFineLevels() 1 0.04744 0.1531 0.212 2.22% MLMG::addInterpCorrection() 390 0.03231 0.1536 0.1886 1.98% MLCellLinOp::applyBC() 3630 0.1408 0.1585 0.1882 1.97% Amr::grid_places() 2 0.002229 0.1089 0.1686 1.77% AmrMesh::MakeNewGrids() 2 0.00212 0.1088 0.1685 1.77% TagBoxArray::mapPRD 2 0.0007038 0.1041 0.1653 1.73% FabArray::sum() 2633 0.1454 0.1469 0.1516 1.59% Gravity::fill_multipole_BCs() 10 0.1055 0.1279 0.1464 1.54% FabArray::FillBoundary() 3240 0.1115 0.1263 0.1457 1.53% Castro::expand_state() 11 0.1006 0.1213 0.1434 1.50% Castro::post_init() 1 0.1379 0.1393 0.141 1.48% MLMG::mgVcycle_up::2 65 0.005554 0.05887 0.1401 1.47% MLMG::mgVcycle_down::0 65 0.1355 0.1377 0.139 1.46% Amr::writePlotFile() 2 0.105 0.1051 0.1051 1.10% Gravity::multilevel_solve_for_new_phi() 1 0.1012 0.1012 0.1013 1.06% Gravity::actual_multilevel_solve() 1 0.1011 0.1011 0.1012 1.06% MLCellLinOp::apply() 501 0.07723 0.08256 0.0995 1.04% MLMG::mgVcycle_up::0 65 0.09851 0.09891 0.09938 1.04% FabArray::ParallelCopy_nowait() 734 0.08319 0.08794 0.09693 1.02% MLMG::mgVcycle_up::3 65 0.005502 0.01209 0.09131 0.96% Castro::post_regrid() 2 0.08631 0.08808 0.08943 0.94% Castro::derive() 18 0.07097 0.07191 0.07607 0.80% AmrLevel::derive() 18 0.07093 0.07187 0.07602 0.80% MLMG:computeResOfCorrection() 390 0.05904 0.06362 0.07487 0.78% MLCellLinOp::correctionResidual() 390 0.05832 0.06288 0.07409 0.78% FillBoundary_nowait() 3240 0.04506 0.05821 0.07257 0.76% FillBoundary_finish() 3240 0.04759 0.06093 0.06875 0.72% Castro::post_timestep() 1 0.0585 0.06115 0.06335 0.66% MLMG::mgVcycle_down::1 65 0.05754 0.05991 0.06299 0.66% Castro::construct_new_source() 48 0.05842 0.05898 0.06071 0.64% MLPoisson::Fapply() 501 0.04892 0.05114 0.0596 0.62% Castro::initialize_advance() 1 0.04937 0.05083 0.05563 0.58% MLCellLinOp::solutionResidual() 75 0.04459 0.04711 0.04943 0.52% Castro::construct_old_source() 66 0.04385 0.04429 0.04626 0.49% Amr::regrid() 1 0.04438 0.04524 0.04625 0.48% Castro::create_source_corrector() 8 0.01578 0.02886 0.04323 0.45% FabArray::Saxpy() 133 0.03538 0.03596 0.04096 0.43% MLMG::computeResidual() 65 0.03848 0.03867 0.03902 0.41% Castro::construct_old_gravity() 11 0.02306 0.03024 0.03666 0.38% MLMG::mgVcycle_bottom 65 0.0001848 0.002466 0.03663 0.38% MLMG::actualBottomSolve() 65 6.825e-05 0.002348 0.03651 0.38% Gravity::get_old_grav_vector() 11 0.02226 0.02869 0.03637 0.38% MLMG::mgVcycle_down::2 65 0.006414 0.0207 0.03603 0.38% MLCGSolver::bicgstab 4 0 0.002232 0.03571 0.37% Castro::construct_new_gravity_source() 8 0.03265 0.03293 0.03334 0.35% FabArray::Xpay() 480 0.0259 0.02706 0.03092 0.32% MLMG::mgVcycle_down::3 65 0.006161 0.009262 0.03051 0.32% VisMF::Write(FabArray) 4 0.01399 0.01562 0.02787 0.29% amrex::average_down 390 0.01672 0.01982 0.0278 0.29% Castro::construct_old_gravity_source() 11 0.02546 0.02566 0.02582 0.27% Castro::sum_integrated_quantities() 2 0.02545 0.02563 0.02564 0.27% Castro::retry_advance_ctu() 11 0.02315 0.02362 0.02508 0.26% Castro::check_for_nan() 19 0.02324 0.02339 0.02376 0.25% MultiFab::contains_nan() 19 0.02318 0.02333 0.0237 0.25% MLPoisson::define() 10 0.01555 0.01774 0.02084 0.22% TagBoxArray::collate() 2 9.587e-05 0.003441 0.02081 0.22% Gravity::get_new_grav_vector() 9 0.01776 0.01814 0.01856 0.19% MLMG::mgVcycle_down::4 65 0.006091 0.006951 0.0179 0.19% Castro::react_state() 19 0.01533 0.01567 0.01604 0.17% Castro::apply_source_to_state() 19 0.0147 0.01484 0.01513 0.16% StateDataPhysBCFunct::() 82 0.0003459 0.007779 0.01458 0.15% StateData::FillBoundary(geom) 62 0 0.007357 0.01413 0.15% Gravity::test_composite_phi() 1 0.01406 0.01408 0.01411 0.15% MLMG::mgVcycle_up::4 65 0.005548 0.006137 0.01397 0.15% Amr::writeSmallPlotFile() 2 0.01391 0.01392 0.01392 0.15% Castro::problem_post_timestep() 2 0.009111 0.01183 0.01376 0.14% Castro::wd_update() 2 0.009013 0.01181 0.01375 0.14% Castro::save_data_for_retry() 3 0.01168 0.01197 0.01329 0.14% MLCellLinOp::defineAuxData() 10 0.01001 0.01076 0.01294 0.14% MLMG::prepareForSolve() 9 0.01032 0.01074 0.01201 0.13% MLMG::mgVcycle_down::5 65 0.006103 0.006538 0.01195 0.13% MLMG::mgVcycle_up::5 65 0.005581 0.005937 0.01013 0.11% MLMG::computeMLResidual() 9 0.005472 0.00785 0.01005 0.11% FabArray::mult() 40 0.008399 0.008494 0.008727 0.09% MLMG::ResNormInf() 74 0.006823 0.00755 0.008255 0.09% FabArray::norminf() 143 0.005633 0.005884 0.007978 0.08% Castro::volWgtSum() 46 0.00738 0.007459 0.007922 0.08% FabArrayBase::getCPC() 1022 0.005021 0.005688 0.007331 0.08% FabArrayBase::getFB() 3240 0.005203 0.005592 0.007231 0.08% Other 1448 0.04547 0.05974 0.08912 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 18 B 20 B 23 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.04-22-g64d2360b209c) finalized