Initializing AMReX (24.05)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05) initialized Starting run at 21:37:51 UTC on 2024-05-01. Successfully read inputs file ... Castro git describe: 24.05-1-g7d1a1c16d AMReX git describe: 24.05 Microphysics git describe: 24.05 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.201932273 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001076512 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010559452 Gravity::multilevel_solve_for_new_phi() time = 0.101849183 Gravity::fill_multipole_BCs() time = 0.010506073 ... 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.049141319 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.011671442 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.002350585 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.070888362 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -160.2179808 at index 18, 75, 0 Invalid density = -80.16477042 at index 18, 76, 0 Invalid density = -61.58489855 at index 19, 77, 0 Invalid density = -108.2035501 at index 19, 78, 0 Invalid density = -15.15828486 at index 20, 80, 0 Invalid density = -8.076810125 at index 20, 81, 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -8.076810125 at index 20, 174, 0 Invalid density = -15.15828486 at index 20, 175, 0 Invalid density = -108.2035501 at index 19, 177, 0 Invalid density = -61.58489855 at index 19, 178, 0 Invalid density = -80.16477042 at index 18, 179, 0 Invalid density = -160.2179808 at index 18, 180, 0 Invalid density = -552.7694089 at index 17, 181, 0 Invalid density = -415.2911996 at index 16, 183, 0 Invalid density = -764.1508347 at index 15, 184, 0 Invalid density = -908.1731819 at index 14, 185, 0 Invalid density = -859.0333271 at index 13, 186, 0 Invalid density = -606.7538937 at index 12, 187, 0 Invalid density = -907.0826581 at index 10, 188, 0 Invalid density = -241.113047 at index 11, 188, 0 Invalid density = -820.8014889 at index 8, 189, 0 Invalid density = -456.4708539 at index 9, 189, 0 Invalid density = -6.239071149 at index 10, 189, 0 Invalid density = -625.0368072 at index 5, 190, 0 Invalid density = -691.1269281 at index 6, 190, 0 Invalid density = -318.5467376 at index 7, 190, 0 Invalid density = -11.08517908 at index 8, 190, 0 Invalid density = -618.7284155 at index 0, 191, 0 Invalid density = -540.6137805 at index 1, 191, 0 Invalid density = -443.905673 at index 2, 191, 0 Invalid density = -302.5298758 at index 3, 191, 0 Invalid density = -167.9940382 at index 4, 191, 0 Invalid density = -55.18215035 at index 5, 191, 0 Invalid density = -1.089487306 at index 6, 191, 0 Advance was unsuccessful with reason: invalid density; proceeding to a retry.  Subcycle completed  Timestep 0.25 rejected at level 0. Performing a retry, with subcycled timesteps of maximum length dt = 0.125   Beginning subcycle 1 starting at time 0 with dt = 0.125 Estimated number of subcycles remaining: 2 Castro::construct_old_gravity() time = 0.006642942 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.067944689 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.007505589 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.067755017 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.006560472 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.066068844 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.325216064 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016479917 Gravity::solve_for_phi() time = 0.091190453 on level 0 Castro::construct_new_gravity() time = 0.093375019 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.062124042 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.002488607 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.061571677 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.332366684 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018174527 Gravity::solve_for_phi() time = 0.093794859 on level 0 Castro::construct_new_gravity() time = 0.095939972 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.061743755 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.002280782 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.063048146 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.332866027 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.017380496 Gravity::solve_for_phi() time = 0.092802558 on level 0 Castro::construct_new_gravity() time = 0.095066848 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.062224448 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.002862371 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.063692267 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.330217527 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018271402 Gravity::solve_for_phi() time = 0.089134549 on level 0 Castro::construct_new_gravity() time = 0.091279776 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.061202136 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.002209883 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.063124649 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.333424657 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018520991 Gravity::solve_for_phi() time = 0.094209893 on level 0 Castro::construct_new_gravity() time = 0.096346712 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.061398107 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.002678883 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.065009986 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.33063212 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.018867568 Gravity::solve_for_phi() time = 0.0926666 on level 0 Castro::construct_new_gravity() time = 0.094817269 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.061195587 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.0028959 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.062174954 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.332089455 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01913965 Gravity::solve_for_phi() time = 0.091290918 on level 0 Castro::construct_new_gravity() time = 0.0934435 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.061323088 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.00280493 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.061769153 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.333703621 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016922372 Gravity::solve_for_phi() time = 0.091677095 on level 0 Castro::construct_new_gravity() time = 0.093830036 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.063404226 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.004491638147 [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.356858747 [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.054426217 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.004109464 seconds Ending run at 21:38:01 UTC on 2024-05-01. Run time = 10.21219075 Run time without initialization = 7.416210926 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.21 ... 10.21 ... 10.21 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8189 2.8879 3.0406 29.77% Amr::InitAmr() 1 1.825 1.825 1.825 17.87% Castro::computeTemp() 63 0.6925 0.7589 0.8523 8.35% FabArray::ParallelCopy_finish() 734 0.05672 0.3864 0.6073 5.95% amrex::Copy() 523 0.5679 0.5786 0.5973 5.85% Castro::normalize_species() 62 0.4506 0.4592 0.4744 4.65% Castro::reset_internal_energy(Fab) 126 0.4322 0.4411 0.469 4.59% MultiFab::AddProduct() 2633 0.415 0.4228 0.4342 4.25% Castro::initMFs() 1 0.3541 0.3621 0.3797 3.72% MultiFab::Subtract() 125 0.273 0.2785 0.2926 2.87% Castro::reset_internal_energy(MultiFab) 63 0.01796 0.07405 0.2863 2.80% check_for_negative_density() 11 0.01381 0.1858 0.2627 2.57% FabArray::setVal() 3453 0.2303 0.2344 0.2437 2.39% MLPoisson::Fsmooth() 3120 0.1579 0.1733 0.2058 2.02% Castro::initData() 1 0.003314 0.0289 0.1746 1.71% Gravity::fill_multipole_BCs() 10 0.1067 0.1355 0.1622 1.59% FabArray::sum() 2633 0.1471 0.1504 0.1541 1.51% FillBoundary_finish() 3240 0.05857 0.07165 0.09802 0.96% Castro::enforce_min_density() 62 0.0137 0.02729 0.08408 0.82% FabArray::ParallelCopy_nowait() 734 0.06238 0.0685 0.07575 0.74% MLPoisson::Fapply() 501 0.0483 0.05281 0.07121 0.70% FillBoundary_nowait() 3240 0.04087 0.05283 0.06551 0.64% MLCellLinOp::applyBC() 3630 0.02475 0.03226 0.05258 0.51% AmrLevel::derive() 18 0.0405 0.04138 0.0428 0.42% Castro::sum_integrated_quantities() 2 0.00189 0.03937 0.04228 0.41% Castro::problem_post_timestep() 2 1.034e-05 0.00254 0.04046 0.40% FabArray::Saxpy() 133 0.03649 0.03761 0.04045 0.40% Castro::do_new_sources() 8 0.0001703 0.006979 0.0393 0.38% Castro::construct_new_gravity_source() 8 0.03239 0.03337 0.03616 0.35% FabArray::Xpay() 480 0.02584 0.02733 0.03041 0.30% Castro::construct_old_gravity_source() 11 0.02531 0.02582 0.02805 0.27% Castro::wd_update() 2 0.006748 0.01872 0.02746 0.27% VisMF::Write(FabArray) 4 0.01112 0.01254 0.02575 0.25% Castro::do_old_sources() 11 0.000307 0.00697 0.02554 0.25% Amr::writePlotFile() 2 0.01248 0.02252 0.02507 0.25% MultiFab::contains_nan() 19 0.02276 0.02319 0.02368 0.23% MLMG::addInterpCorrection() 390 0.01653 0.01801 0.02108 0.21% amrex::average_down 390 0.01244 0.01359 0.01619 0.16% MLMG::ResNormInf() 74 0.006771 0.01062 0.01522 0.15% StateData::FillBoundary(geom) 62 0 0.007058 0.01399 0.14% Amr::writeSmallPlotFile() 2 0.01152 0.01314 0.01363 0.13% Castro::evaluate_source_change() 144 0.009232 0.009594 0.01019 0.10% FabArray::mult() 40 0.00839 0.008605 0.009426 0.09% MLCellLinOp::defineAuxData() 10 0.007892 0.008242 0.009188 0.09% Castro::construct_new_source() 48 0.007289 0.007467 0.008714 0.09% FabArray::FillBoundary() 3240 0.006876 0.007227 0.008425 0.08% Castro::volWgtSum() 46 0.007365 0.007499 0.008267 0.08% MLCGSolver::bicgstab 4 0 0.0004978 0.007965 0.08% FabArray::norminf() 143 0.005549 0.005888 0.007819 0.08% Castro::construct_old_gravity() 11 0.0002125 0.001905 0.00712 0.07% Castro::construct_old_source() 66 0.005366 0.005539 0.006808 0.07% MLCellLinOp::smooth() 1560 0.005807 0.006112 0.00649 0.06% FillPatchIterator::Initialize 82 0.005483 0.005578 0.005747 0.06% MLMG::solve() 9 0.0009445 0.003877 0.005432 0.05% MLMG::prepareForSolve() 9 0.004987 0.00517 0.005391 0.05% Gravity::actual_solve_with_mlmg() 10 0.003433 0.003617 0.004006 0.04% Gravity::get_old_grav_vector() 11 0.003597 0.003671 0.003908 0.04% Gravity::get_new_grav_vector() 9 0.003093 0.003138 0.00332 0.03% Castro::initialize_do_advance() 11 0.002416 0.002527 0.002714 0.03% MLCellLinOp::apply() 501 0.001036 0.00115 0.00236 0.02% FabArray::ParallelCopy() 734 0.001768 0.001865 0.001976 0.02% AmrLevel::FillPatch() 82 0.001082 0.00114 0.001242 0.01% main() 1 0.0009384 0.001101 0.001226 0.01% MLLinOp::defineGrids() 10 0.0009253 0.00098 0.001159 0.01% MLCellLinOp::correctionResidual() 390 0.0009088 0.0009419 0.0009762 0.01% MLMG::mgVcycle() 65 0.000773 0.0008325 0.0009243 0.01% Castro::save_data_for_retry() 3 0.0007465 0.0007605 0.0007779 0.01% MLMG:computeResOfCorrection() 390 0.0007051 0.0007332 0.0007738 0.01% Gravity::solve_for_phi() 8 0.0003852 0.0004447 0.0006066 0.01% Castro::Castro() 1 0.0004635 0.0004958 0.000567 0.01% Castro::subcycle_advance_ctu() 1 9.234e-05 0.0001339 0.0005638 0.01% AmrMesh::MakeNewGrids() 2 0.0004979 0.0005212 0.000554 0.01% MLMG::actualBottomSolve() 65 6.747e-05 9.879e-05 0.0005333 0.01% StateDataPhysBCFunct::() 82 0.0003392 0.0004116 0.0004719 0.00% Castro::construct_new_gravity() 8 0.0001407 0.0001916 0.0004042 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003718 0.0003852 0.0004041 0.00% Castro::clean_state() 62 0.0003247 0.0003512 0.0003925 0.00% MLMG::mgVcycle_down::2 65 0.0003052 0.0003356 0.0003628 0.00% MLMG::mgVcycle_down::0 65 0.0003354 0.0003453 0.0003626 0.00% MLMG::mgVcycle_down::1 65 0.0003352 0.0003438 0.000362 0.00% MLMG::mgVcycle_down::3 65 0.0002825 0.0003104 0.0003534 0.00% MLMG::mgVcycle_down::4 65 0.0002762 0.0002913 0.0003418 0.00% MLMG::mgVcycle_down::5 65 0.0002844 0.0002987 0.0003307 0.00% Castro::do_advance_ctu() 11 0.0002304 0.0002442 0.0002859 0.00% FillPatchSingleLevel 82 0.000255 0.0002633 0.0002743 0.00% MLMG::mgVcycle_up::0 65 0.0002362 0.0002456 0.0002657 0.00% MLMG::mgVcycle_up::1 65 0.0002292 0.000244 0.0002569 0.00% MLMG::mgVcycle_up::2 65 0.0001932 0.0002197 0.0002437 0.00% MLCellLinOp::solutionResidual() 75 0.000218 0.0002268 0.0002406 0.00% MLMG::mgVcycle_up::3 65 0.0001911 0.000202 0.0002316 0.00% MLMG::mgVcycle_up::5 65 0.0002005 0.0002071 0.0002273 0.00% MLMG::oneIter() 65 0.0001993 0.0002072 0.0002244 0.00% MLMG::mgVcycle_up::4 65 0.0001941 0.0002024 0.0002239 0.00% Amr::coarseTimeStep() 1 0.0001501 0.000165 0.0001912 0.00% Amr::defBaseLevel() 1 0.0001148 0.0001348 0.0001523 0.00% Castro::retry_advance_ctu() 11 5.208e-05 6.386e-05 0.0001458 0.00% Gravity::solve_phi_with_mlmg() 10 0.000115 0.0001248 0.0001415 0.00% MLMG::computeResidual() 65 0.0001215 0.000127 0.0001393 0.00% Amr::grid_places() 2 7.192e-05 8.883e-05 0.000129 0.00% MLMG::mgVcycle_bottom 65 0.0001135 0.0001187 0.0001288 0.00% TagBoxArray::mapPRD 2 0.000108 0.0001131 0.0001169 0.00% Gravity::test_composite_phi() 1 9.296e-05 9.794e-05 0.0001126 0.00% Castro::initialize_advance() 1 0.0001016 0.0001048 0.0001109 0.00% Gravity::multilevel_solve_for_new_phi() 1 5.995e-05 6.814e-05 9.206e-05 0.00% MLPoisson::define() 10 6.684e-05 6.983e-05 7.59e-05 0.00% Castro::check_for_nan() 19 4.874e-05 5.777e-05 7.434e-05 0.00% Castro::react_state() 19 5.114e-05 5.446e-05 6.213e-05 0.00% Castro::apply_source_to_state() 19 4.997e-05 5.237e-05 5.85e-05 0.00% Amr::FinalizeInit() 1 8.261e-06 1.143e-05 5.001e-05 0.00% Castro::derive() 18 4.048e-05 4.275e-05 4.952e-05 0.00% MLLinOp::define() 10 3.849e-05 4.206e-05 4.781e-05 0.00% Castro::create_source_corrector() 8 3.438e-05 3.881e-05 4.4e-05 0.00% Castro::expand_state() 11 3.09e-05 3.526e-05 4.39e-05 0.00% Gravity::actual_multilevel_solve() 1 2.666e-05 2.829e-05 3.301e-05 0.00% MLMG::computeMLResidual() 9 2.055e-05 2.24e-05 3.127e-05 0.00% Amr::regrid() 1 1.042e-05 1.272e-05 2.936e-05 0.00% Castro::post_regrid() 2 1.097e-05 1.623e-05 2.363e-05 0.00% Amr::init() 1 5.637e-06 7.934e-06 2.348e-05 0.00% Amr::timeStep() 1 1.136e-05 1.587e-05 2.118e-05 0.00% Castro::post_init() 1 1.211e-05 1.315e-05 1.549e-05 0.00% Amr::InitializeInit() 1 6.704e-06 7.956e-06 1.39e-05 0.00% Castro::advance() 1 7.312e-06 7.979e-06 1.202e-05 0.00% Amr::bldFineLevels() 1 6.545e-06 7.255e-06 7.93e-06 0.00% Castro::post_timestep() 1 6.059e-06 6.526e-06 6.878e-06 0.00% Amr::initialInit() 1 2.943e-06 3.17e-06 3.363e-06 0.00% Other 5692 0.03562 0.04832 0.07056 0.69% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 10.21 10.21 10.21 100.00% Amr::coarseTimeStep() 1 7.416 7.416 7.416 72.62% Amr::timeStep() 1 7.357 7.357 7.357 72.04% Castro::advance() 1 7.277 7.285 7.297 71.45% Castro::subcycle_advance_ctu() 1 7.228 7.232 7.237 70.87% Castro::do_advance_ctu() 11 7.204 7.207 7.212 70.62% Castro::construct_ctu_hydro_source() 11 3.422 3.45 3.482 34.09% Castro::clean_state() 62 2.975 3.091 3.241 31.73% Castro::computeTemp() 63 1.924 2.046 2.264 22.16% Amr::InitAmr() 1 1.825 1.825 1.826 17.88% Castro::reset_internal_energy(MultiFab) 63 1.209 1.287 1.479 14.48% Amr::init() 1 0.97 0.97 0.97 9.50% Amr::initialInit() 1 0.9084 0.9084 0.9085 8.90% Gravity::solve_phi_with_mlmg() 10 0.7888 0.8206 0.8501 8.32% Castro::initialize_do_advance() 11 0.7747 0.7948 0.8307 8.13% Castro::evaluate_source_change() 144 0.7259 0.7401 0.7601 7.44% Castro::construct_new_gravity() 8 0.693 0.7244 0.7536 7.38% Gravity::solve_for_phi() 8 0.6751 0.7069 0.7363 7.21% Castro::do_old_sources() 11 0.6817 0.696 0.7112 6.96% Castro::enforce_min_density() 62 0.6273 0.6513 0.6972 6.83% FabArray::ParallelCopy() 734 0.1439 0.4739 0.6925 6.78% Gravity::actual_solve_with_mlmg() 10 0.68 0.6832 0.6861 6.72% MLMG::solve() 9 0.6529 0.6548 0.6555 6.42% Amr::InitializeInit() 1 0.437 0.5097 0.6374 6.24% Amr::defBaseLevel() 1 0.437 0.5097 0.6374 6.24% FabArray::ParallelCopy_finish() 734 0.05672 0.3864 0.6073 5.95% amrex::Copy() 523 0.5679 0.5786 0.5973 5.85% MLMG::oneIter() 65 0.5701 0.5735 0.5774 5.65% MLMG::mgVcycle() 65 0.5625 0.5658 0.5698 5.58% Castro::do_new_sources() 8 0.4449 0.474 0.4947 4.84% Castro::normalize_species() 62 0.4506 0.4592 0.4744 4.65% Amr::FinalizeInit() 1 0.271 0.3987 0.4714 4.62% Castro::reset_internal_energy(Fab) 126 0.4322 0.4411 0.469 4.59% MultiFab::AddProduct() 2633 0.415 0.4228 0.4342 4.25% AmrLevel::FillPatch() 82 0.2006 0.2956 0.4008 3.93% MLCellLinOp::smooth() 1560 0.2887 0.3163 0.3842 3.76% Castro::Castro() 1 0.3583 0.3662 0.3837 3.76% Castro::initMFs() 1 0.3542 0.3622 0.3798 3.72% FillPatchIterator::Initialize 82 0.1238 0.2208 0.3274 3.21% FillPatchIterator::FillFromLevel0() 82 0.1147 0.2125 0.3198 3.13% FillPatchSingleLevel 82 0.1143 0.2121 0.3194 3.13% MultiFab::Subtract() 125 0.273 0.2785 0.2926 2.87% Castro::initData() 1 0.07713 0.1421 0.2773 2.71% check_for_negative_density() 11 0.01381 0.1858 0.2627 2.57% Amr::bldFineLevels() 1 0.04911 0.1749 0.2462 2.41% FabArray::setVal() 3453 0.2303 0.2344 0.2437 2.39% MLMG::mgVcycle_up::1 65 0.05231 0.1359 0.2209 2.16% MLCellLinOp::applyBC() 3630 0.1479 0.1696 0.2155 2.11% MLPoisson::Fsmooth() 3120 0.1579 0.1733 0.2058 2.02% Amr::grid_places() 2 0.002616 0.1303 0.203 1.99% AmrMesh::MakeNewGrids() 2 0.002541 0.1302 0.2029 1.99% TagBoxArray::mapPRD 2 0.0008368 0.1285 0.2014 1.97% MLMG::addInterpCorrection() 390 0.03132 0.1528 0.1903 1.86% Castro::post_init() 1 0.1759 0.1795 0.1822 1.78% Castro::expand_state() 11 0.1059 0.1291 0.1686 1.65% Gravity::fill_multipole_BCs() 10 0.1079 0.1366 0.1634 1.60% FabArray::FillBoundary() 3240 0.1172 0.1374 0.1629 1.60% FabArray::sum() 2633 0.1471 0.1504 0.1541 1.51% MLMG::mgVcycle_up::2 65 0.005429 0.05956 0.1417 1.39% MLMG::mgVcycle_down::0 65 0.1319 0.134 0.1391 1.36% MLMG::mgVcycle_up::0 65 0.1018 0.1036 0.1076 1.05% MLCellLinOp::apply() 501 0.07919 0.08578 0.1042 1.02% Amr::writePlotFile() 2 0.1041 0.1042 0.1042 1.02% Gravity::multilevel_solve_for_new_phi() 1 0.1018 0.1019 0.1019 1.00% Gravity::actual_multilevel_solve() 1 0.1018 0.1018 0.1018 1.00% FillBoundary_finish() 3240 0.05857 0.07165 0.09802 0.96% FabArray::ParallelCopy_nowait() 734 0.07936 0.08561 0.09293 0.91% Castro::post_regrid() 2 0.08617 0.08888 0.09248 0.91% MLMG::mgVcycle_up::3 65 0.005381 0.01178 0.0875 0.86% MLMG:computeResOfCorrection() 390 0.05752 0.06323 0.07977 0.78% Castro::post_timestep() 1 0.06003 0.07206 0.07957 0.78% MLCellLinOp::correctionResidual() 390 0.05682 0.0625 0.07902 0.77% Castro::derive() 18 0.06979 0.07139 0.07285 0.71% AmrLevel::derive() 18 0.06975 0.07134 0.0728 0.71% FillBoundary_nowait() 3240 0.04624 0.0585 0.07141 0.70% MLPoisson::Fapply() 501 0.0483 0.05281 0.07121 0.70% Castro::construct_new_source() 48 0.05884 0.06053 0.06852 0.67% MLMG::mgVcycle_down::1 65 0.05964 0.06237 0.06569 0.64% Castro::sum_integrated_quantities() 2 0.02109 0.05894 0.06146 0.60% Castro::initialize_advance() 1 0.04876 0.05229 0.05925 0.58% MLCellLinOp::solutionResidual() 75 0.04845 0.051 0.05497 0.54% Castro::problem_post_timestep() 2 0.01046 0.0249 0.05116 0.50% Castro::construct_old_source() 66 0.04378 0.04463 0.04834 0.47% Amr::regrid() 1 0.04422 0.04565 0.04785 0.47% MLMG::computeResidual() 65 0.04113 0.04263 0.0478 0.47% Castro::create_source_corrector() 8 0.01556 0.03048 0.04328 0.42% Castro::construct_old_gravity() 11 0.02453 0.03473 0.04202 0.41% FabArray::Saxpy() 133 0.03649 0.03761 0.04045 0.40% Gravity::get_old_grav_vector() 11 0.02349 0.03283 0.03721 0.36% MLMG::mgVcycle_down::2 65 0.006309 0.02076 0.03692 0.36% Castro::construct_new_gravity_source() 8 0.03239 0.03337 0.03616 0.35% MLMG::mgVcycle_bottom 65 0.0001822 0.002385 0.03534 0.35% MLMG::actualBottomSolve() 65 6.747e-05 0.002266 0.03521 0.34% MLCGSolver::bicgstab 4 0 0.002153 0.03445 0.34% Castro::wd_update() 2 0.01045 0.02236 0.03081 0.30% MLMG::mgVcycle_down::3 65 0.005992 0.009262 0.03063 0.30% FabArray::Xpay() 480 0.02584 0.02733 0.03041 0.30% amrex::average_down 390 0.01653 0.0201 0.02887 0.28% Castro::construct_old_gravity_source() 11 0.02531 0.02582 0.02805 0.27% VisMF::Write(FabArray) 4 0.01112 0.01254 0.02575 0.25% Castro::retry_advance_ctu() 11 0.02307 0.02337 0.02456 0.24% Castro::check_for_nan() 19 0.02282 0.02325 0.02375 0.23% MultiFab::contains_nan() 19 0.02276 0.02319 0.02368 0.23% MLPoisson::define() 10 0.01727 0.01957 0.02298 0.23% MLMG::ResNormInf() 74 0.01167 0.0156 0.02013 0.20% Gravity::get_new_grav_vector() 9 0.01801 0.01833 0.0199 0.19% MLMG::mgVcycle_down::4 65 0.005956 0.006858 0.0174 0.17% Castro::apply_source_to_state() 19 0.01513 0.01558 0.01737 0.17% Castro::react_state() 19 0.01548 0.0158 0.0163 0.16% Amr::writeSmallPlotFile() 2 0.01619 0.0162 0.01621 0.16% StateDataPhysBCFunct::() 82 0.0003392 0.007469 0.01446 0.14% Gravity::test_composite_phi() 1 0.01398 0.014 0.01402 0.14% StateData::FillBoundary(geom) 62 0 0.007058 0.01399 0.14% MLMG::mgVcycle_up::4 65 0.00543 0.006088 0.01354 0.13% Castro::save_data_for_retry() 3 0.01149 0.01179 0.01288 0.13% MLCellLinOp::defineAuxData() 10 0.01007 0.01087 0.01272 0.12% MLMG::prepareForSolve() 9 0.0102 0.01074 0.01166 0.11% MLMG::mgVcycle_down::5 65 0.005936 0.006457 0.01151 0.11% MLMG::computeMLResidual() 9 0.00594 0.007796 0.01007 0.10% MLMG::mgVcycle_up::5 65 0.005477 0.005893 0.009829 0.10% FabArray::mult() 40 0.00839 0.008605 0.009426 0.09% Castro::volWgtSum() 46 0.007365 0.007499 0.008267 0.08% MLLinOp::define() 10 0.00222 0.00553 0.008193 0.08% MLLinOp::defineGrids() 10 0.002177 0.005488 0.00815 0.08% FabArray::norminf() 143 0.005549 0.005888 0.007819 0.08% Other 5692 0.05204 0.06997 0.09873 0.97% -------------------------------------------------------------------------------------------- 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.05) finalized