shoot.contours.get_closed_contours

shoot.contours.get_closed_contours#

shoot.contours.get_closed_contours(lon_center, lat_center, ssh, nlevels=50, robust=0.03)[source]#

Extract closed contours enclosing a center point

Parameters:

  • lon_center (float) – Center longitude in degrees.

  • lat_center (float) – Center latitude in degrees.

  • ssh (xarray.DataArray) – 2D field to contour (typically SSH).

  • nlevels (int, default 50) – Maximum number of contour levels.

  • robust (float, default 0.03) – Quantile threshold to exclude extreme values.

Returns:

Each dataset contains a closed contour with coordinates and metadata, including lon/lat coordinates, SSH level, and center position.

Return type:

list of xarray.Dataset

Example

>>> from shoot.contours import get_closed_contours
>>> contours = get_closed_contours(5.0, 43.0, ssh_field)
>>> len(contours)  # number of nested closed contours