error about eddy_identification function

[zhanglan009]

Dear Author,
When I running the eddy_identification function, I always get only 1 anticyclone and 1 cyclone, or 0 anticyclones and cyclones entirely, regardless of running the example pet_eddy_sla_and_adt you provided or the example on https://py-eddy-tracker.readthedocs.io/en/latest/grid_identification.html. I'm wondering if the eddy_identification function has been modified and is now unable to calculate all a and c.


code:
from matplotlib import pyplot as plt
from py_eddy_tracker import start_logger
from py_eddy_tracker.dataset.grid import RegularGridDataset

start_logger().setLevel("DEBUG") # Available options: ERROR, WARNING, INFO, DEBUG

grid_name= ('dataset-duacs-nrt-global-merged-allsat-phy-l4_1699415795667.nc')
from datetime import datetime

h = RegularGridDataset(grid_name, "longitude", "latitude", nan_masking=True)
h.bessel_high_filter("adt", 500)
date = datetime(2019, 2, 23)
h.add_uv("adt")
u, v = h.grid("u").T, h.grid("v").T
a, c = h.eddy_identification("adt","u","v", # Variables used for identification
date, # Date of identification
0.002, # step between two isolines of detection (m)
pixel_limit=(5, 2000), # Min and max pixel count for valid contour
shape_error=55, # Error max (%) between ratio of circle fit and contour
force_speed_unit="m/s"
)

from datetime import datetime

h = RegularGridDataset(grid_name, "longitude", "latitude", nan_masking=True)
h.bessel_high_filter("adt", 500)
date = datetime(2019, 2, 23)
h.add_uv("adt")
u, v = h.grid("u").T, h.grid("v").T
a, c = h.eddy_identification("adt","u","v", # Variables used for identification
date, # Date of identification
0.002, # step between two isolines of detection (m)
pixel_limit=(5, 2000), # Min and max pixel count for valid contour
shape_error=55, # Error max (%) between ratio of circle fit and contour
force_speed_unit="m/s"
)

[AntSimi]

take a look at #219

[zhanglan009]

Thank you very much! I sincerely appreciate your invaluable assistance.