./sod.py /raid/testing/castro-gfortran/Castro/ Sod-mhd_plt00081
 slicing along x-direction at coarse grid (j,k)=(8,8) and output to sod_slice.out
analysis_routine = /home/zingale/bin/fextract.gnu.ex
Traceback (most recent call last):
  File "/raid/testing/castro-gfortran/Castro-SBU-tests/2025-04-21/Sod-mhd/./sod.py", line 117, in <module>
    process(castro_exec_dir, plotfile)
    ~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/raid/testing/castro-gfortran/Castro-SBU-tests/2025-04-21/Sod-mhd/./sod.py", line 70, in process
    if d >= 2 and data["ymom"].ptp() > data["xmom"].ptp(): dim = "ymom"
                  ^^^^^^^^^^^^^^^^
AttributeError: `ptp` was removed from the ndarray class in NumPy 2.0. Use np.ptp(arr, ...) instead.