Initializing AMReX (24.05-3-gc3616816b46a)... MPI initialized with 16 MPI processes MPI initialized with thread support level 3 AMReX (24.05-3-gc3616816b46a) initialized Starting run at 07:09:31 UTC on 2024-05-06. Successfully read inputs file ... Castro git describe: 24.05 AMReX git describe: 24.05-3-gc3616816b 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.213642212 new finest: 0 Now regridding at level lbase = 0 grid_places() time: 0.001038544 new finest: 0 STEP = 0 TIME = 0 : REGRID with lbase = 0 Gravity::fill_multipole_BCs() time = 0.010748494 Gravity::multilevel_solve_for_new_phi() time = 0.103167341 Gravity::fill_multipole_BCs() time = 0.01388122 ... 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.061123659 seconds SMALL PLOTFILE: file = smallplt00000 Write small plotfile time = 0.007461907 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.006933151 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.06926609 on level 0 ... Entering construct_ctu_hydro_source() on level 0 Invalid density = -9.947567811e-05 at index 0, 62, 0 Invalid density = -1.942940124 at index 0, 63, 0 Invalid density = -1.164962232 at index 1, 63, 0 Invalid density = -0.767071355 at index 2, 63, 0 Invalid density = -0.08878802337 at index 3, 63, 0 Invalid density = -0.01174181879 at index 4, 63, 0 Invalid density = -0.0008017363671 at index 5, 63, 0 Invalid density = -1.942940124 at index 0, 192, 0 Invalid density = -1.164962232 at index 1, 192, 0 Invalid density = -0.767071355 at index 2, 192, 0 Invalid density = -0.08878802337 at index 3, 192, 0 Invalid density = -0.01174181879 at index 4, 192, 0 Invalid density = -0.0008017363671 at index 5, 192, 0 Invalid density = -9.947567811e-05 at index 0, 193, 0 Invalid density = -618.7284155 at index 0, 64, 0 Invalid density = -540.6137805 at index 1, 64, 0 Invalid density = -443.905673 at index 2, 64, 0 Invalid density = -302.5298758 at index 3, 64, 0 Invalid density = -167.9940382 at index 4, 64, 0 Invalid density = -55.18215035 at index 5, 64, 0 Invalid density = -1.089487306 at index 6, 64, 0 Invalid density = -625.0368072 at index 5, 65, 0 Invalid density = -691.1269281 at index 6, 65, 0 Invalid density = -318.5467376 at index 7, 65, 0 Invalid density = -11.08517908 at index 8, 65, 0 Invalid density = -820.8014889 at index 8, 66, 0 Invalid density = -456.4708539 at index 9, 66, 0 Invalid density = -6.239071149 at index 10, 66, 0 Invalid density = -907.0826581 at index 10, 67, 0 Invalid density = -241.113047 at index 11, 67, 0 Invalid density = -606.7538937 at index 12, 68, 0 Invalid density = -859.0333271 at index 13, 69, 0 Invalid density = -908.1731819 at index 14, 70, 0 Invalid density = -764.1508347 at index 15, 71, 0 Invalid density = -415.2911996 at index 16, 72, 0 Invalid density = -552.7694089 at index 17, 74, 0 Invalid density = -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.017763788 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.072055221 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.008571306 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.062936989 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.004568104 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.060795349 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.316920438 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015201876 Gravity::solve_for_phi() time = 0.087000258 on level 0 Castro::construct_new_gravity() time = 0.089189956 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.060460548 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.002544984 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.060133427 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.323486781 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01552079 Gravity::solve_for_phi() time = 0.085697345 on level 0 Castro::construct_new_gravity() time = 0.087869655 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.06021923 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.002284186 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.060449236 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.322714949 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.015553783 Gravity::solve_for_phi() time = 0.085748271 on level 0 Castro::construct_new_gravity() time = 0.087927867 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.060170626 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.002228857 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.06014978 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.323855319 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.01593508 Gravity::solve_for_phi() time = 0.086310399 on level 0 Castro::construct_new_gravity() time = 0.088426908 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.060289078 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.002275964 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.096815146 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.339101045 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016274087 Gravity::solve_for_phi() time = 0.086911778 on level 0 Castro::construct_new_gravity() time = 0.089073287 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.060342256 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.002278102 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.060722017 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.345736771 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.021527486 Gravity::solve_for_phi() time = 0.094227787 on level 0 Castro::construct_new_gravity() time = 0.096338922 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.087514262 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.002149149 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.059472209 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.345654918 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.019450351 Gravity::solve_for_phi() time = 0.140055307 on level 0 Castro::construct_new_gravity() time = 0.143782676 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.069430205 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.002334643 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.071303412 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.33017643 on level 0 ... new-time level Poisson gravity solve at level 0 Gravity::fill_multipole_BCs() time = 0.016191026 Gravity::solve_for_phi() time = 0.09139034 on level 0 Castro::construct_new_gravity() time = 0.094176595 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.077839072 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.004454307931 [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.414075383 [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.048057192 seconds SMALL PLOTFILE: file = smallplt00001 Write small plotfile time = 0.009037407 seconds Ending run at 07:09:40 UTC on 2024-05-06. Run time = 9.846834391 Run time without initialization = 7.472559835 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]: 9.846 ... 9.847 ... 9.847 -------------------------------------------------------------------------------------------- Name NCalls Excl. Min Excl. Avg Excl. Max Max % -------------------------------------------------------------------------------------------- Castro::construct_ctu_hydro_source() 11 2.8368 2.8876 3.0196 30.66% Amr::InitAmr() 1 1.422 1.422 1.422 14.44% Castro::computeTemp() 63 0.6873 0.7504 0.843 8.56% FabArray::ParallelCopy_finish() 734 0.02972 0.4332 0.6598 6.70% amrex::Copy() 523 0.5453 0.5542 0.5817 5.91% Castro::normalize_species() 62 0.45 0.4561 0.4652 4.72% MultiFab::AddProduct() 2633 0.4295 0.4358 0.4527 4.60% Castro::reset_internal_energy(Fab) 126 0.4261 0.433 0.4423 4.49% Castro::initMFs() 1 0.353 0.3649 0.4063 4.13% MultiFab::Subtract() 125 0.2628 0.2697 0.3187 3.24% Castro::reset_internal_energy(MultiFab) 63 0.01798 0.0742 0.2864 2.91% check_for_negative_density() 11 0.06671 0.2027 0.268 2.72% FabArray::setVal() 3453 0.2295 0.2349 0.2422 2.46% MLPoisson::Fsmooth() 3120 0.1595 0.1772 0.2083 2.12% Castro::initData() 1 0.003358 0.02779 0.1627 1.65% Gravity::fill_multipole_BCs() 10 0.1103 0.1348 0.158 1.60% FabArray::sum() 2633 0.1416 0.144 0.1497 1.52% FillBoundary_finish() 3240 0.07908 0.1003 0.1274 1.29% Castro::enforce_min_density() 62 0.01399 0.04091 0.1134 1.15% FabArray::ParallelCopy_nowait() 734 0.06713 0.07205 0.07924 0.80% FillBoundary_nowait() 3240 0.04165 0.0534 0.06617 0.67% MLPoisson::Fapply() 501 0.04805 0.05116 0.05765 0.59% MLCellLinOp::applyBC() 3630 0.02393 0.0322 0.05257 0.53% Castro::do_new_sources() 8 0.0001723 0.009217 0.0482 0.49% Castro::do_old_sources() 11 0.0002923 0.009034 0.04367 0.44% AmrLevel::derive() 18 0.04033 0.04088 0.04159 0.42% Castro::construct_new_gravity_source() 8 0.03252 0.0337 0.03871 0.39% VisMF::Write(FabArray) 4 0.01436 0.01661 0.03859 0.39% FabArray::Saxpy() 133 0.03441 0.03532 0.03813 0.39% Amr::writePlotFile() 2 0.008976 0.02777 0.0309 0.31% FabArray::Xpay() 480 0.0252 0.02676 0.02972 0.30% Castro::wd_update() 2 0.005032 0.02016 0.02928 0.30% Castro::construct_old_gravity_source() 11 0.02535 0.02596 0.02857 0.29% MultiFab::contains_nan() 19 0.02276 0.0232 0.02378 0.24% MLMG::addInterpCorrection() 390 0.01751 0.01915 0.02219 0.23% amrex::average_down 390 0.01249 0.01366 0.0163 0.17% StateData::FillBoundary(geom) 62 0 0.007228 0.01401 0.14% Castro::construct_old_gravity() 11 0.0002204 0.002355 0.01299 0.13% MLMG::ResNormInf() 74 0.005553 0.009152 0.01251 0.13% Amr::writeSmallPlotFile() 2 0.00749 0.01139 0.01214 0.12% Castro::evaluate_source_change() 144 0.009139 0.009411 0.00997 0.10% MLCellLinOp::defineAuxData() 10 0.007704 0.008119 0.009071 0.09% FabArray::mult() 40 0.008299 0.008503 0.008663 0.09% FabArray::FillBoundary() 3240 0.006818 0.00713 0.008462 0.09% Castro::construct_new_source() 48 0.007233 0.007514 0.00839 0.09% amrex::Add() 65 0.007311 0.007468 0.00798 0.08% MLCGSolver::bicgstab 4 0 0.0004967 0.007948 0.08% FabArray::norminf() 143 0.005534 0.005963 0.007711 0.08% Castro::volWgtSum() 46 0.007284 0.007404 0.007559 0.08% MLCellLinOp::smooth() 1560 0.005639 0.005979 0.006453 0.07% FillPatchIterator::Initialize 82 0.005633 0.005892 0.006153 0.06% Castro::construct_old_source() 66 0.005414 0.005554 0.006078 0.06% MLMG::prepareForSolve() 9 0.005055 0.005227 0.00539 0.05% MLMG::solve() 9 0.0006806 0.002757 0.004169 0.04% Gravity::actual_solve_with_mlmg() 10 0.003473 0.00367 0.004036 0.04% FabArrayBase::getCPC() 1022 0.003646 0.0038 0.004016 0.04% Gravity::get_old_grav_vector() 11 0.003558 0.003693 0.003888 0.04% Castro::sum_integrated_quantities() 2 0.001742 0.00347 0.003817 0.04% Gravity::get_new_grav_vector() 9 0.003079 0.003208 0.003509 0.04% Castro::initialize_do_advance() 11 0.002434 0.002822 0.003424 0.03% Castro::construct_new_gravity() 8 0.0001465 0.0008537 0.003082 0.03% MLCellLinOp::apply() 501 0.001011 0.001143 0.002358 0.02% Gravity::solve_for_phi() 8 0.0003767 0.0007609 0.002186 0.02% FabArray::ParallelCopy() 734 0.001797 0.001893 0.002053 0.02% main() 1 0.0009536 0.001284 0.001743 0.02% AmrLevel::FillPatch() 82 0.001091 0.001158 0.001272 0.01% MLCellLinOp::correctionResidual() 390 0.0009039 0.0009427 0.001 0.01% Castro::save_data_for_retry() 3 0.0007434 0.0008116 0.0009372 0.01% MLMG::mgVcycle() 65 0.0007814 0.0008289 0.000916 0.01% MLMG:computeResOfCorrection() 390 0.0007106 0.000746 0.0007974 0.01% Castro::subcycle_advance_ctu() 1 0.0001024 0.0001594 0.000745 0.01% Castro::Castro() 1 0.0004566 0.0005217 0.0005951 0.01% MLMG::actualBottomSolve() 65 6.704e-05 9.973e-05 0.0005334 0.01% AmrMesh::MakeNewGrids() 2 0.0004895 0.0005068 0.0005324 0.01% StateDataPhysBCFunct::() 82 0.0003503 0.0004164 0.0004663 0.00% FillPatchIterator::FillFromLevel0() 82 0.0003803 0.0003874 0.0004043 0.00% Castro::problem_post_timestep() 2 1.022e-05 3.611e-05 0.0004017 0.00% Castro::clean_state() 62 0.0003365 0.0003623 0.0003919 0.00% MLMG::mgVcycle_down::2 65 0.000302 0.0003354 0.0003812 0.00% MLMG::mgVcycle_down::0 65 0.0003275 0.0003375 0.0003561 0.00% Castro::do_advance_ctu() 11 0.0002217 0.0002576 0.0003556 0.00% MLMG::mgVcycle_down::1 65 0.0003245 0.0003336 0.0003529 0.00% MLMG::mgVcycle_down::3 65 0.0002846 0.0003118 0.0003493 0.00% MLMG::mgVcycle_down::4 65 0.0002777 0.0002928 0.0003283 0.00% MLMG::mgVcycle_down::5 65 0.0002879 0.000299 0.000326 0.00% FillPatchSingleLevel 82 0.0002555 0.0002639 0.0002765 0.00% MLMG::mgVcycle_up::1 65 0.0002305 0.000247 0.0002697 0.00% MLMG::mgVcycle_up::0 65 0.0002359 0.0002451 0.0002595 0.00% MLMG::mgVcycle_up::2 65 0.0001948 0.0002181 0.0002523 0.00% MLCellLinOp::solutionResidual() 75 0.0002163 0.0002265 0.0002402 0.00% MLMG::oneIter() 65 0.0002009 0.0002096 0.000226 0.00% MLMG::mgVcycle_up::3 65 0.0001919 0.0002023 0.0002253 0.00% MLMG::mgVcycle_up::5 65 0.0002018 0.00021 0.0002236 0.00% Castro::retry_advance_ctu() 11 5.05e-05 7.316e-05 0.0002213 0.00% MLMG::mgVcycle_up::4 65 0.0001937 0.0002017 0.0002188 0.00% Amr::coarseTimeStep() 1 0.0001647 0.0001876 0.0001941 0.00% Castro::initialize_advance() 1 0.0001269 0.0001478 0.0001698 0.00% Amr::defBaseLevel() 1 0.0001255 0.0001386 0.0001592 0.00% Gravity::solve_phi_with_mlmg() 10 0.0001145 0.0001231 0.0001372 0.00% MLMG::computeResidual() 65 0.0001194 0.0001251 0.0001359 0.00% MLMG::mgVcycle_bottom 65 0.0001156 0.0001196 0.0001287 0.00% TagBoxArray::mapPRD 2 0.0001102 0.0001161 0.0001269 0.00% Gravity::test_composite_phi() 1 9.406e-05 0.0001007 0.000123 0.00% Castro::react_state() 19 5.269e-05 6.268e-05 0.000113 0.00% Amr::grid_places() 2 7.531e-05 8.991e-05 0.0001116 0.00% Castro::check_for_nan() 19 4.871e-05 5.732e-05 9.552e-05 0.00% Gravity::multilevel_solve_for_new_phi() 1 6.231e-05 7.34e-05 9.529e-05 0.00% MLPoisson::define() 10 6.531e-05 7.273e-05 9.196e-05 0.00% Amr::FinalizeInit() 1 7.882e-06 1.428e-05 7.163e-05 0.00% Castro::apply_source_to_state() 19 4.863e-05 5.435e-05 5.937e-05 0.00% Castro::derive() 18 4.044e-05 4.338e-05 4.936e-05 0.00% Castro::create_source_corrector() 8 3.35e-05 3.884e-05 4.568e-05 0.00% Castro::expand_state() 11 3.224e-05 3.795e-05 4.415e-05 0.00% Amr::regrid() 1 1.101e-05 1.294e-05 3.194e-05 0.00% Gravity::actual_multilevel_solve() 1 2.76e-05 2.947e-05 3.114e-05 0.00% MLMG::computeMLResidual() 9 2.077e-05 2.36e-05 2.868e-05 0.00% Amr::timeStep() 1 1.206e-05 1.619e-05 2.491e-05 0.00% Castro::advance() 1 7.177e-06 1.183e-05 2.121e-05 0.00% Amr::InitializeInit() 1 6.625e-06 9.077e-06 1.739e-05 0.00% Castro::post_init() 1 1.134e-05 1.258e-05 1.549e-05 0.00% Castro::post_regrid() 2 1.039e-05 1.2e-05 1.526e-05 0.00% Amr::bldFineLevels() 1 6.87e-06 8.022e-06 1.018e-05 0.00% Castro::post_timestep() 1 5.722e-06 6.466e-06 9.161e-06 0.00% Amr::init() 1 6.019e-06 7.23e-06 8.864e-06 0.00% Amr::initialInit() 1 2.968e-06 3.359e-06 3.796e-06 0.00% Other 4625 0.02609 0.03532 0.05782 0.59% -------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------- Name NCalls Incl. Min Incl. Avg Incl. Max Max % -------------------------------------------------------------------------------------------- main() 1 9.846 9.846 9.847 100.00% Amr::coarseTimeStep() 1 7.472 7.472 7.472 75.88% Amr::timeStep() 1 7.414 7.414 7.414 75.29% Castro::advance() 1 7.333 7.342 7.356 74.71% Castro::subcycle_advance_ctu() 1 7.283 7.29 7.302 74.16% Castro::do_advance_ctu() 11 7.257 7.265 7.278 73.91% Castro::construct_ctu_hydro_source() 11 3.437 3.461 3.488 35.42% Castro::clean_state() 62 2.951 3.061 3.22 32.70% Castro::computeTemp() 63 1.907 2.016 2.249 22.84% Castro::reset_internal_energy(MultiFab) 63 1.19 1.265 1.467 14.90% Amr::InitAmr() 1 1.422 1.422 1.422 14.44% Amr::init() 1 0.9515 0.9515 0.9516 9.66% Amr::initialInit() 1 0.8822 0.8822 0.8823 8.96% Gravity::solve_phi_with_mlmg() 10 0.8247 0.8511 0.8762 8.90% Castro::initialize_do_advance() 11 0.7665 0.8108 0.8406 8.54% Castro::construct_new_gravity() 8 0.7284 0.752 0.7765 7.89% Castro::evaluate_source_change() 144 0.7349 0.7458 0.7748 7.87% Gravity::solve_for_phi() 8 0.7073 0.7325 0.7575 7.69% FabArray::ParallelCopy() 734 0.1198 0.5241 0.7481 7.60% Castro::do_old_sources() 11 0.6973 0.7103 0.7232 7.34% Castro::enforce_min_density() 62 0.6202 0.6568 0.7206 7.32% Gravity::actual_solve_with_mlmg() 10 0.7119 0.7144 0.7165 7.28% MLMG::solve() 9 0.6847 0.6866 0.6873 6.98% FabArray::ParallelCopy_finish() 734 0.02972 0.4332 0.6598 6.70% Amr::InitializeInit() 1 0.438 0.5127 0.6501 6.60% Amr::defBaseLevel() 1 0.438 0.5127 0.6501 6.60% MLMG::oneIter() 65 0.6062 0.6092 0.6124 6.22% MLMG::mgVcycle() 65 0.5986 0.6015 0.6042 6.14% amrex::Copy() 523 0.5453 0.5542 0.5817 5.91% Castro::do_new_sources() 8 0.4426 0.4792 0.5347 5.43% Castro::normalize_species() 62 0.45 0.4561 0.4652 4.72% MultiFab::AddProduct() 2633 0.4295 0.4358 0.4527 4.60% AmrLevel::FillPatch() 82 0.1926 0.3309 0.4446 4.51% Amr::FinalizeInit() 1 0.2322 0.3695 0.4442 4.51% Castro::reset_internal_energy(Fab) 126 0.4261 0.433 0.4423 4.49% MLCellLinOp::smooth() 1560 0.3134 0.3455 0.4141 4.21% Castro::Castro() 1 0.3573 0.3693 0.4117 4.18% Castro::initMFs() 1 0.3531 0.365 0.4064 4.13% FillPatchIterator::Initialize 82 0.1209 0.2589 0.374 3.80% FillPatchIterator::FillFromLevel0() 82 0.1099 0.2501 0.3661 3.72% FillPatchSingleLevel 82 0.1095 0.2497 0.3657 3.71% MultiFab::Subtract() 125 0.2628 0.2697 0.3187 3.24% Castro::initData() 1 0.076 0.142 0.2912 2.96% check_for_negative_density() 11 0.06671 0.2027 0.268 2.72% Amr::bldFineLevels() 1 0.04653 0.1837 0.2586 2.63% MLCellLinOp::applyBC() 3630 0.1739 0.1987 0.2458 2.50% FabArray::setVal() 3453 0.2295 0.2349 0.2422 2.46% MLMG::mgVcycle_up::1 65 0.05621 0.1429 0.2296 2.33% Amr::grid_places() 2 0.002569 0.14 0.2147 2.18% AmrMesh::MakeNewGrids() 2 0.002457 0.1399 0.2146 2.18% TagBoxArray::mapPRD 2 0.000806 0.1383 0.2131 2.16% MLPoisson::Fsmooth() 3120 0.1595 0.1772 0.2083 2.12% MLMG::addInterpCorrection() 390 0.0325 0.1566 0.1963 1.99% FabArray::FillBoundary() 3240 0.1441 0.1665 0.1932 1.96% Castro::expand_state() 11 0.09901 0.1508 0.1906 1.94% Gravity::fill_multipole_BCs() 10 0.1114 0.1359 0.159 1.62% FabArray::sum() 2633 0.1416 0.144 0.1497 1.52% MLMG::mgVcycle_down::0 65 0.1408 0.1432 0.1468 1.49% Castro::post_init() 1 0.1411 0.1424 0.1446 1.47% MLMG::mgVcycle_up::2 65 0.005469 0.0595 0.1416 1.44% FillBoundary_finish() 3240 0.07908 0.1003 0.1274 1.29% MLMG::mgVcycle_up::0 65 0.1133 0.1162 0.1184 1.20% Amr::writePlotFile() 2 0.1097 0.1098 0.1098 1.12% MLCellLinOp::apply() 501 0.08136 0.08769 0.1055 1.07% Gravity::multilevel_solve_for_new_phi() 1 0.1031 0.1032 0.1032 1.05% Gravity::actual_multilevel_solve() 1 0.1031 0.1031 0.1031 1.05% FabArray::ParallelCopy_nowait() 734 0.08363 0.08898 0.0961 0.98% Castro::post_regrid() 2 0.08487 0.08708 0.08835 0.90% MLMG::mgVcycle_up::3 65 0.005421 0.01186 0.08827 0.90% Castro::derive() 18 0.06867 0.07031 0.08145 0.83% AmrLevel::derive() 18 0.06863 0.07026 0.0814 0.83% Castro::post_timestep() 1 0.05753 0.07238 0.08129 0.83% MLMG:computeResOfCorrection() 390 0.0592 0.06559 0.07682 0.78% MLCellLinOp::correctionResidual() 390 0.05846 0.06485 0.07604 0.77% MLMG::mgVcycle_down::1 65 0.06422 0.06757 0.07251 0.74% FillBoundary_nowait() 3240 0.04692 0.059 0.07195 0.73% Castro::construct_new_source() 48 0.05804 0.05988 0.06739 0.68% Castro::initialize_advance() 1 0.04884 0.05167 0.05996 0.61% MLPoisson::Fapply() 501 0.04805 0.05116 0.05765 0.59% MLCellLinOp::solutionResidual() 75 0.04755 0.04999 0.05215 0.53% Castro::construct_old_gravity() 11 0.02907 0.03886 0.04927 0.50% Castro::construct_old_source() 66 0.04342 0.04438 0.04843 0.49% Gravity::get_old_grav_vector() 11 0.02286 0.0365 0.04736 0.48% Amr::regrid() 1 0.0437 0.04474 0.04583 0.47% MLMG::computeResidual() 65 0.04081 0.04159 0.04249 0.43% MLMG::mgVcycle_down::2 65 0.006324 0.02243 0.04026 0.41% Castro::construct_new_gravity_source() 8 0.03252 0.0337 0.03871 0.39% VisMF::Write(FabArray) 4 0.01436 0.01661 0.03859 0.39% FabArray::Saxpy() 133 0.03441 0.03532 0.03813 0.39% Castro::create_source_corrector() 8 0.01526 0.02669 0.03799 0.39% MLMG::mgVcycle_bottom 65 0.0001836 0.002384 0.03529 0.36% MLMG::actualBottomSolve() 65 6.704e-05 0.002264 0.03517 0.36% MLCGSolver::bicgstab 4 0 0.00215 0.0344 0.35% Castro::problem_post_timestep() 2 0.009311 0.02391 0.03275 0.33% Castro::wd_update() 2 0.008909 0.02387 0.03274 0.33% MLMG::mgVcycle_down::3 65 0.006015 0.00927 0.0302 0.31% amrex::average_down 390 0.01666 0.02058 0.02981 0.30% FabArray::Xpay() 480 0.0252 0.02676 0.02972 0.30% Castro::construct_old_gravity_source() 11 0.02535 0.02596 0.02857 0.29% Castro::retry_advance_ctu() 11 0.02239 0.02336 0.02761 0.28% Castro::check_for_nan() 19 0.02281 0.02325 0.02387 0.24% MultiFab::contains_nan() 19 0.02276 0.0232 0.02378 0.24% Castro::sum_integrated_quantities() 2 0.0224 0.02269 0.02272 0.23% Gravity::get_new_grav_vector() 9 0.01814 0.01968 0.02194 0.22% MLPoisson::define() 10 0.01578 0.01847 0.02146 0.22% Gravity::test_composite_phi() 1 0.018 0.01802 0.01804 0.18% MLMG::mgVcycle_down::4 65 0.00599 0.006887 0.01751 0.18% Castro::react_state() 19 0.01554 0.0161 0.01749 0.18% MLMG::ResNormInf() 74 0.0107 0.01408 0.01742 0.18% Amr::writeSmallPlotFile() 2 0.01692 0.01692 0.01693 0.17% Castro::apply_source_to_state() 19 0.01426 0.01454 0.01581 0.16% StateDataPhysBCFunct::() 82 0.0003503 0.007644 0.01447 0.15% StateData::FillBoundary(geom) 62 0 0.007228 0.01401 0.14% MLMG::mgVcycle_up::4 65 0.005483 0.006131 0.01372 0.14% Castro::save_data_for_retry() 3 0.01094 0.01134 0.01321 0.13% MLCellLinOp::defineAuxData() 10 0.009913 0.01069 0.01258 0.13% MLMG::prepareForSolve() 9 0.01018 0.0107 0.01165 0.12% MLMG::mgVcycle_down::5 65 0.005983 0.006473 0.01158 0.12% MLMG::mgVcycle_up::5 65 0.005495 0.005916 0.009887 0.10% MLMG::computeMLResidual() 9 0.005372 0.007654 0.009444 0.10% FabArray::mult() 40 0.008299 0.008503 0.008663 0.09% amrex::Add() 65 0.007311 0.007468 0.00798 0.08% FabArray::norminf() 143 0.005534 0.005963 0.007711 0.08% Castro::volWgtSum() 46 0.007284 0.007404 0.007559 0.08% FabArrayBase::getCPC() 1022 0.005028 0.00571 0.007231 0.07% Other 4625 0.04495 0.06158 0.09568 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 16 B 19 B 21 B 8192 KiB 8192 KiB 8192 KiB --------------------------------------------------------------------------------------------------------------------- AMReX (24.05-3-gc3616816b46a) finalized