./test1-helm.py /raid/testing/castro-gfortran/Castro/ Sod_stellar-test1-2d_plt00091
 slicing along x-direction at coarse grid (j,k)=(4,0) and output to test1-helm.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_stellar-test1-2d/./test1-helm.py", line 124, in <module>
    process(castro_exec_dir, plotfile)
    ~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/raid/testing/castro-gfortran/Castro-SBU-tests/2025-04-21/Sod_stellar-test1-2d/./test1-helm.py", line 72, 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.