CMORise CMEMS sea surface height zos

esmvaltool/cmorizers/data/cmor_config/CMEMS.yml

@@ -0,0 +1,21 @@
+---
+# Filename
+filename: 'dt_global_allsat_phy_l4.*.nc'
+
+# Common global attributes for Cmorizer output
+attributes:
+  project_id: OBS6
+  dataset_id: CMEMS
+  version: 'SEALEVEL-GLO-PHY-L4'
+  tier: 2
+  modeling_realm: sat
+  source: ''
+  reference: 'cmems'
+  comment: ''
+
+# Variables to cmorize
+variables:
+  zos:
+    mip: Omon
+    raw: adt
+    name: sea_surface_height_above_geoid

esmvaltool/cmorizers/data/datasets.yml

@@ -280,6 +280,14 @@ datasets:
       https://downloads.psl.noaa.gov/Datasets/cmap/enh/
         precip.mon.mean.nc
 
+  CMEMS:
+    tier: 2
+    source: https://doi.org/10.48670/moi-00148
+    last_access: 2025-02-20
+    info: |
+      Download daily files from:
+      https://data.marine.copernicus.eu/product/SEALEVEL_GLO_PHY_L4_MY_008_047/services
+
   CowtanWay:
     tier: 2
     source: https://www-users.york.ac.uk/~kdc3/papers/coverage2013/series.html

esmvaltool/cmorizers/data/formatters/datasets/cmems.py

@@ -0,0 +1,97 @@
+"""ESMValTool CMORizer for CMEMS - Sea Level Thematic Assembly Center.
+
+   This is the CMORizer script for Sea Surface Height Above Geoid.
+   Copernicus Marine Environment Monitoring Service (CMEMS)
+   Product: SEALEVEL_GLO_PHY_L4_MY_008_047
+
+Tier
+    Tier 2: open dataset.
+
+Source
+    https://doi.org/10.48670/moi-00148
+
+Last access
+    20250220
+
+Download and processing instructions
+    Download daily files from:
+    https://data.marine.copernicus.eu/product/SEALEVEL_GLO_PHY_L4_MY_008_047/services
+
+
+"""
+
+import logging
+import os
+import re
+
+import iris
+from esmvalcore.preprocessor import monthly_statistics
+
+from esmvaltool.cmorizers.data import utilities as utils
+
+logger = logging.getLogger(__name__)
+
+
+def _get_filepaths(in_dir, basename):
+    """Group daily files."""
+    regex = re.compile(basename)
+    return_ls = []
+
+    for root, _dir, files in os.walk(in_dir, followlinks=True):
+        if len(files) > 1:
+            return_files = []  # load and concat each list & process
+            for filename in files:
+                if regex.match(filename):
+                    return_files.append(os.path.join(root, filename))
+
+            return_ls.append(return_files)
+
+    return return_ls
+
+
+def _extract_variable(cmor_info, attrs, file_ls, out_dir):
+    """Extract variable and aggregate months."""
+    var = cmor_info.short_name
+    standard_name = cmor_info.standard_name
+
+    cube_prepls = iris.cube.CubeList()
+    for filepaths in file_ls:
+        cubels = iris.load(filepaths, standard_name)
+        # logger.info('Number of day files to concatenate: %d', len(filepaths))
+        iris.util.equalise_attributes(cubels)
+        iris.util.unify_time_units(cubels)
+        # if cmor.frequency is mon, check to save out
+        cube_mon = cubels.concatenate_cube()  # create mean for month
+        cube_mon = iris.util.squeeze(cube_mon)
+        cube_prepls.append(monthly_statistics(cube_mon, "mean"))
+
+    iris.util.equalise_attributes(cube_prepls)
+    cube = cube_prepls.concatenate_cube()
+    utils.fix_var_metadata(cube, cmor_info)
+    cube = utils.fix_coords(cube)
+
+    utils.set_global_atts(cube, attrs)
+    utils.save_variable(
+        cube, var, out_dir, attrs, unlimited_dimensions=["time"]
+    )
+
+
+def cmorization(in_dir, out_dir, cfg, cfg_user, start_date, end_date):
+    """Cmorization func call."""
+    glob_attrs = cfg["attributes"]
+    cmor_table = cfg["cmor_table"]
+
+    file_ls = _get_filepaths(in_dir, cfg["filename"])
+
+    if len(file_ls) > 0:
+        logger.info("Found %d input files in '%s'", len(file_ls), in_dir)
+    else:
+        logger.info("No files found, basename: %s", cfg["filename"])
+
+    # Run the cmorization
+    for var, var_info in cfg["variables"].items():
+        logger.info("CMORizing variable '%s'", var)
+        glob_attrs["mip"] = var_info["mip"]
+        cmor_info = cmor_table.get_variable(var_info["mip"], var)
+
+        _extract_variable(cmor_info, glob_attrs, file_ls, out_dir)

esmvaltool/recipes/examples/recipe_check_obs.yml

@@ -67,7 +67,17 @@ diagnostics:
       pr:
     additional_datasets:
       - {project: OBS6, dataset: CMAP, mip: Amon, tier: 2,
-        type: reanaly, version: v1}
+        type: sat, version: v1}
+    scripts: null
+
+
+  CMEMS:
+    description: CMEMS check
+    variables:
+      zos:
+    additional_datasets:
+      - {project: OBS6, dataset: CMEMS, mip: Omon, tier: 2,
+        type: reanaly, version: SEALEVEL-GLO-PHY-L4}
     scripts: null
 
 

esmvaltool/references/cmems.bibtex

@@ -0,0 +1,10 @@
+@misc{cmems,
+  doi = {10.48670/MOI-00148},
+  url = {https://resources.marine.copernicus.eu/product-detail/SEALEVEL_GLO_PHY_L4_MY_008_047/INFORMATION},
+  author = {{European Union-Copernicus Marine Service}},
+  keywords = {oceanography},
+  language = {en},
+  title = {GLOBAL OCEAN GRIDDED L4 SEA SURFACE HEIGHTS AND DERIVED VARIABLES REPROCESSED (1993-ONGOING)},
+  publisher = {Mercator Ocean International},
+  year = {2021}
+}