Work on BGC from 2024 LOV visit

argopy/__init__.py

@@ -42,6 +42,7 @@
 from .utils import MonitoredThreadPoolExecutor  # noqa: E402, F401
 from .utils import monitor_status as status  # noqa: E402
 from .related import TopoFetcher, OceanOPSDeployments, ArgoNVSReferenceTables, ArgoDocs, ArgoDOI  # noqa: E402
+from .extensions import CanyonMED
 
 
 #
@@ -77,6 +78,10 @@
     "ArgoColors",  # Class
     "stores",
     "tutorial",
+
+    # Argo xarray accessor extensions
+    "CanyonMED",
+
     # Constants
     "__version__"
 )

argopy/data_fetchers/argovis_data.py

@@ -82,6 +82,7 @@ def __init__(
         """
         self.definition = "Argovis Argo data fetcher"
         self.dataset_id = OPTIONS["dataset"] if ds == "" else ds
+        self.user_mode =  kwargs["mode"] if "mode" in kwargs else OPTIONS["mode"]
         self.server = kwargs["server"] if "server" in kwargs else api_server
         timeout = OPTIONS["api_timeout"] if api_timeout == 0 else api_timeout
         self.store_opts = {
@@ -379,14 +380,20 @@ def to_xarray(self, errors: str = "ignore"):
         ds = ds[np.sort(ds.data_vars)]
         return ds
 
+    def transform_data_mode(self, ds: xr.Dataset, **kwargs):
+        # Argovis data are already curated !
+        if ds.argo._type == "point":
+            ds["N_POINTS"] = np.arange(0, len(ds["N_POINTS"]))
+        return ds
+
     def filter_data_mode(self, ds: xr.Dataset, **kwargs):
-        # Argovis data already curated !
+        # Argovis data are already curated !
         if ds.argo._type == "point":
             ds["N_POINTS"] = np.arange(0, len(ds["N_POINTS"]))
         return ds
 
     def filter_qc(self, ds: xr.Dataset, **kwargs):
-        # Argovis data already curated !
+        # Argovis data are already curated !
         if ds.argo._type == "point":
             ds["N_POINTS"] = np.arange(0, len(ds["N_POINTS"]))
         return ds
@@ -396,7 +403,7 @@ def filter_researchmode(self, ds: xr.Dataset, *args, **kwargs) -> xr.Dataset:
 
         This filter will select only QC=1 delayed mode data with pressure errors smaller than 20db
 
-        Use this filter instead of filter_data_mode and filter_qc
+        Use this filter instead of transform_data_mode and filter_qc
         """
         ds = ds.argo.filter_researchmode()
         if ds.argo._type == "point":