OceanDataStore 0.3.0__py3-none-any.whl

OceanDataStore/data/WOA23/send_WOA23_monthly_climatology_to_os.py

@@ -0,0 +1,292 @@
+# =========================================================
+# send_WOA23_monthly_climatology_to_os.py
+#
+# Script to write WOA23 monthly climatologies to Icechunk
+# repository in JASMIN cloud object storage.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# =========================================================
+import glob
+import logging
+
+import numpy as np
+import pandas as pd
+import xarray as xr
+import zarr
+
+from OceanDataStore.cli import initialise_logging, send_to_icechunk
+from OceanDataStore.data.utils import (
+    compute_cell_area,
+    compute_dx,
+    compute_dy,
+)
+
+logger = logging.getLogger(__name__)
+
+
+# -- Utility Functions -- #
+def get_filepaths(file_prefix:str, file_directory:str) -> list:
+    """
+    Get filepaths for World Ocean Atlas 2023 climate normals
+    in-situ temperature and practical salinity data. 
+
+    Parameters
+    ----------
+    file_prefix : str
+        Prefix of the WOA23 filepaths to be returned.
+    file_directory : str
+        Directory containing the WOA23 files.
+
+    Returns
+    -------
+    filepaths : list
+        List of WOA23 filepaths for practical salinity & in-situ temperature.
+    """
+    # Define salinity and temperature filenames:
+    salinity_filename = f"{file_prefix}_s*.nc"
+    temperature_filename = f"{file_prefix}_t*.nc"
+
+    # Collect filepaths for salinity and temperature data:
+    filepaths_sal = sorted(glob.glob(f"{file_directory}/{salinity_filename}"))
+    filepaths_temp = sorted(glob.glob(f"{file_directory}/{temperature_filename}"))
+
+    return filepaths_sal, filepaths_temp
+
+
+def main():
+    # ========== Initialise OceanDataStore Logging ========== #
+    initialise_logging()
+    logging.info("In Progress: Sending WOA23 Monthly Climatologies to JASMIN Object Store.")
+
+    # ========== Prepare Data ========== #
+    filedir = "/dssgfs01/scratch/otooth/npd_data/observations/WOA2023"
+    salinity_paths, temperature_paths = [], []
+    for prefix in ['woa23_decav71A0', 'woa23_decav81B0', 'woa23_decav91C0']:
+        # Define file paths for WOA23 monthly climatologies:
+        filepaths_sal, filepaths_temp = get_filepaths(file_directory=filedir, file_prefix=prefix)
+        salinity_paths.append(filepaths_sal[1:13])
+        temperature_paths.append(filepaths_temp[1:13])
+
+    # Define year bounds for each WOA23 monthly climatology:
+    year_bound_start = [1971, 1981, 1991]
+    year_bound_end = [2000, 2010, 2020]
+
+    for temp_path, sal_path, start_year, end_year in zip(temperature_paths, salinity_paths, year_bound_start, year_bound_end):
+        logging.info(f"-> In Progress: Preparing WOA23 {start_year}-{end_year} monthly climatology...")
+        # Open WOA23 monthly climatologies dataset:
+        ds_s = xr.open_mfdataset(sal_path, decode_times=False, data_vars='all')
+        ds_t = xr.open_mfdataset(temp_path, decode_times=False, data_vars='all')
+
+        ds = xr.merge([ds_s, ds_t], compat='no_conflicts').squeeze(drop=True)
+        ds = ds.rename({'time': 'month'}).assign_coords({'month': np.arange(1, 13)})
+        ds['climatology_bounds'].data = np.array([[np.datetime64(f'{start_year}-{month:02d}', 'M'), np.datetime64(f'{end_year}-{month:02d}', 'M')] for month in range(1, 13)])
+
+        ds = ds.rename({"lon": "longitude",
+                        "lat": "latitude",
+                        "climatology_bounds": "time_bnds"
+                        })
+
+        # Use OceanDataStore standard variable names:
+        for var in ds.data_vars:
+            if var.startswith("t_"):
+                ds = ds.rename({var: var.replace("t_", "thetao_")})
+            elif var.startswith("s_"):
+                ds = ds.rename({var: var.replace("s_", "so_")})
+
+        # Add ancillary variables:
+        ds['dx'] = compute_dx(ds)
+        ds['dy'] = compute_dy(ds)
+        ds['cell_area'] = compute_cell_area(ds)
+
+        # Custom ancillary variables:
+        ds['cell_thickness'] = ds['depth_bnds'].isel(nbounds=1) - ds['depth_bnds'].isel(nbounds=0)
+        ds['cell_volume'] = ds['cell_thickness'] * ds['cell_area']
+
+        # Update attributes for custom ancillary variables:
+        ds['cell_thickness'].attrs.update({
+            'long_name': "Grid-Cell Thickness",
+            'standard_name': "cell_thickness",
+            'units': "m",
+        })
+        ds['cell_volume'].attrs.update({
+            'long_name': "Grid-Cell Volume",
+            'standard_name': "cell_volume",
+            'units': "m3",
+        })
+        logging.info(f"Completed: Prepared WOA23 {start_year}-{end_year} monthly climatology dataset with ancillary variables.")
+
+        # ========== Prepare Ancillary Data ========== #
+        # Open WOA23 land sea mask:
+        df_mask = pd.read_table("/dssgfs01/working/otooth/Software/OceanDataStore/OceanDataStore/data/WOA23/data/landsea_04.msk",
+                                delimiter=',',
+                                header=0
+                                )
+        # Define level of sea floor (i.e. bottom standard level) at each grid cell:
+        ds['bottom_level'] = xr.full_like(ds['cell_area'], fill_value=np.nan).squeeze(drop=True)
+        ds['bottom_level'].data = df_mask['Bottom_Standard_Level'].values.reshape(720, 1440)
+        ds['bottom_level'].name = "bottom_level"
+        ds['bottom_level'].attrs = {"standard_name": "model_level_number_at_sea_floor",
+                                    "long_name": "Model Level Number at Sea Floor",
+                                    "units": "1"
+                                    }
+        # Define land sea mask (1 for ocean grid cells, 0 for land grid cells):
+        ds['mask'] = ds['bottom_level'] > 1
+        ds['mask'].name = "sea_binary_mask"
+        ds["mask"] = ds["mask"].assign_attrs({'long_name': "Land-Sea Binary Mask",
+                                              "standard_name": "sea_binary_mask",
+                                              "comment": "1 = sea, 0 = land"
+                                              })
+        logging.info("Completed: Prepared land sea mask & bottom level variables.")
+        
+        # Open WOA23 basin mask:
+        df_basin = pd.read_table("/dssgfs01/working/otooth/Software/OceanDataStore/OceanDataStore/data/WOA23/data/basinmask_04.msk",
+                                delimiter=',',
+                                header=0
+                                )
+        # Define basin mask (integer values for each ocean basin, NaN for land grid cells):
+        ds['basin_mask'] = xr.full_like(ds['cell_area'], fill_value=np.nan).squeeze(drop=True)
+        for basin, longitude, latitude in df_basin[["Basin_0m", "Longitude", "Latitude"]].itertuples(index=False):
+            ds['basin_mask'].loc[dict(latitude=latitude, longitude=longitude)] = basin
+
+        ds['basin_mask'].attrs = {
+            "standard_name": "ocean_basin_mask",
+            "long_name": "Ocean Basin Mask",
+            "basin_name": {
+            1: "Atlantic Ocean",
+            2: "Pacific Ocean",
+            3: "Indian Ocean",
+            4: "Mediterranean Sea",
+            5: "Baltic Sea",
+            6: "Black Sea",
+            7: "Red Sea",
+            8: "Persian Gulf",
+            9: "Hudson Bay",
+            10: "Southern Ocean",
+            11: "Arctic Ocean",
+            12: "Sea of Japan",
+            13: "Kara Sea",
+            14: "Sulu Sea",
+            15: "Baffin Bay",
+            16: "East Mediterranean",
+            17: "West Mediterranean",
+            18: "Sea of Okhotsk",
+            19: "Banda Sea",
+            20: "Caribbean Sea",
+            21: "Andaman Basin",
+            22: "North Caribbean",
+            23: "Gulf of Mexico",
+            24: "Beaufort Sea",
+            25: "South China Sea",
+            26: "Barents Sea",
+            27: "Celebes Sea",
+            28: "Aleutian Basin",
+            29: "Fiji Basin",
+            30: "North American Basin",
+            31: "West European Basin",
+            32: "Southeast Indian Basin",
+            33: "Coral Sea",
+            34: "East Indian Basin",
+            35: "Central Indian Basin",
+            36: "Southwest Atlantic Basin",
+            37: "Southeast Atlantic Basin",
+            38: "Southeast Pacific Basin",
+            39: "Guatemala Basin",
+            40: "East Caroline Basin",
+            41: "Marianas Basin",
+            42: "Philippine Sea",
+            43: "Arabian Sea",
+            44: "Chile Basin",
+            45: "Somali Basin",
+            46: "Mascarene Basin",
+            47: "Crozet Basin",
+            48: "Guinea Basin",
+            49: "Brazil Basin",
+            50: "Argentine Basin",
+            51: "Tasman Sea",
+            52: "Atlantic Indian Basin",
+            53: "Caspian Sea",
+            54: "Sulu Sea II",
+            55: "Venezuela Basin",
+            56: "Bay of Bengal",
+            57: "Java Sea",
+            58: "East Indian Atlantic Basin",
+            59: "Chiloe",
+            60: "Bransfield Strait",
+            }
+        }
+        logging.info("Completed: Prepared ocean basin mask variable.")
+
+        # ========== CF Attributes ========== #
+        ds.attrs.clear()
+        ds = ds.assign_attrs({
+            "Conventions": "CF-1.6",
+            "title": f"World Ocean Atlas 2023 temperature and salinity monthly climatology ({start_year}-{end_year}).",
+            "description": f"World Ocean Atlas 2023 (WOA23) temperature and salinity monthly climatology for the global ocean from objectively analysed, quality controlled in-situ profile data ({start_year}-{end_year}).",
+            "source": "Numerical models: Objective Analysis. In-situ observations: World Ocean Database (WOD).",
+            "dataset_type": "observation",
+            "product_type": "climatology",
+            "product_version": "1.0",
+            "institution": "NOAA National Centers for Environmental Information (NCEI)",
+            "citation": "Reagan, James R.; Boyer, Tim P.; García, Hernán E.; Locarnini, Ricardo A.; Baranova, Olga K.; Bouchard, Courtney; Cross, Scott L.; Mishonov, Alexey V.; Paver, Christopher R.; Seidov, Dan; Wang, Zhankun; Dukhovskoy, Dmitry (2023). World Ocean Atlas 2023 (NCEI Accession 0270533). https://www.ncei.noaa.gov/archive/accession/0270533. In Reagan, James R.; Boyer, Tim P.; García, Hernán E.; Locarnini, Ricardo A.; Baranova, Olga K.; Bouchard, Courtney; Cross, Scott L.; Mishonov, Alexey V.; Paver, Christopher R.; Seidov, Dan; Wang, Zhankun; Dukhovskoy, Dmitry (2023). World Ocean Atlas 2023. NOAA National Centers for Environmental Information. Dataset. https://doi.org/10.25921/va26-hv25. Accessed 06-05-2026.",
+            "references": "Locarnini, R.A., A.V. Mishonov, O.K. Baranova, J.R. Reagan, T.P. Boyer, D. Seidov, Z. Wang, H.E. Garcia, C. Bouchard, S.L. Cross, C.R. Paver, and D. Dukhovskoy (2024). World Ocean Atlas 2023, Volume 1: Temperature. A. Mishonov Technical Editor, NOAA Atlas NESDIS 89. https://doi.org/10.25923/54bh-1613. Reagan, J.R., D. Seidov, Z. Wang, D. Dukhovskoy, T.P. Boyer, R.A. Locarnini, O.K. Baranova, A.V. Mishonov, H.E. Garcia, C. Bouchard, S.L. Cross, and C.R. Paver (2023). World Ocean Atlas 2023, Volume 2: Salinity. A. Mishonov Technical Editor, NOAA Atlas NESDIS 90, https://doi.org/10.25923/70qt-9574.",
+            "acknowledgement": "None",
+            "license": "World Ocean Atlas 2023 data were obtained from https://www.ncei.noaa.gov/access/world-ocean-atlas-2023/ and are provided under a Creative Commons CC0 1.0 Universal License https://creativecommons.org/publicdomain/zero/1.0/",
+            "doi": "https://doi.org/10.25921/va26-hv25",
+            "platform": "gr",
+            "horizontal_grid_type": "regular rectilinear",
+            "horizontal_grid_resolution": "0.25 degree",
+            "vertical_grid_type": "z",
+            "vertical_grid_coordinate": "depth",
+            "vertical_grid_levels": 57,
+            "aggregation": "mean",
+            "aggregation_frequency": "monthly",
+            "status": "completed",
+            "update_frequency": "None",
+            "bbox": "[-180.0, 180.0, -90.0, 90.0]",
+        })
+        logging.info(f"Completed: Added CF-compliant global attributes to WOA23 {start_year}-{end_year} monthly climatology dataset.")
+
+        # ========== Send to Icechunk Repository ========== #
+        logging.info(f"In Progress: Sending WOA23 {start_year}-{end_year} monthly climatology to Icechunk Repository...")
+        bucket = "woa23"
+        prefix = f"woa23_{start_year}_{end_year}_monthly_climatology"
+        exists = False
+        store_credentials_json = ".../credentials/jasmin_os_credentials.json"
+        branch = "main"
+        commit_message = f"Added WOA23 {start_year}-{end_year} monthly climatology."
+        variable_commits = True
+        config_kwargs = {
+                "temporary_directory":".../OceanDataStore/OceanDataStore/data/WOA23/",
+                "local_directory":".../OceanDataStore/OceanDataStore/data/WOA23/"
+            }
+        cluster_kwargs = {
+                "n_workers" : 25,
+                "threads_per_worker" : 1,
+                "memory_limit":"3GB"
+            }
+        
+        # Optimise chunk sizes for spatial analysis:
+        ds = ds.chunk({'longitude': 1440, 'latitude': 720, 'depth': 5, 'month': 1})
+
+        # Update variable encodings:
+        blosccodec = zarr.codecs.BloscCodec(cname="zstd", clevel=3, shuffle=zarr.codecs.BloscShuffle.shuffle)
+        for var in list(ds.data_vars) + list(ds.coords):
+            ds[var].encoding['compressors'] = [blosccodec]
+
+        send_to_icechunk(
+            file=ds,
+            bucket=bucket,
+            object_prefix=prefix,
+            store_credentials_json=store_credentials_json,
+            exists=exists,
+            branch=branch,
+            commit_message=commit_message,
+            variable_commits=variable_commits,
+            dask_config_kwargs=config_kwargs,
+            dask_cluster_kwargs=cluster_kwargs,
+            )
+
+        logging.info(f"Completed: Sent WOA23 {start_year}-{end_year} monthly climatology to Icechunk Repository.")
+
+if __name__ == "__main__":
+    main()

OceanDataStore/data/update_icechunk_repo_attrs.py

@@ -0,0 +1,76 @@
+# =========================================================
+# update_icechunk_repo_attrs.py
+#
+# Script to update global and variable attributes in an
+# Icechunk repository.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# =========================================================
+from OceanDataStore.data.utils import (
+        update_icechunk_global_attrs,
+        update_icechunk_variable_attrs,
+)
+
+# ========= Shared Inputs ========= #
+# Define credential to write to JASMIN OS:
+credentials_filepath = '.../credentials/jasmin_os_credentials.json'
+
+# Define time period for climatology:
+start_year = 1991
+end_year = 2020
+
+# Define bucket and prefix to Icechunk repository:
+bucket = "armor3d"
+prefix = f"armor3d_global_my_{start_year}_{end_year}_monthly_climatogy"
+
+# ========= Update variable attributes ========= #
+vars = ["month"]
+attrs = [{"units": "1", "long_name": "Month of Year"}]
+message = f"Updated ARMOR3D Monthly Climatology ({start_year}-{end_year}) variable attributes. -> ['month']"
+
+update_icechunk_variable_attrs(
+        credentials_filepath=credentials_filepath,
+        bucket=bucket,
+        prefix=prefix,
+        vars=vars,
+        attrs=attrs,
+        commit_message=message,
+        )
+
+# ========= Update global attributes ========= #
+attrs = {
+        "Conventions": "CF-1.0",
+        "title": f"Multi Observation Global Ocean 3D Temperature Salinity Height Geostrophic Current and MLD monthly climatology ({start_year}-{end_year}).",
+        "description": f"Multi Observation Global Ocean ARMOR3D multi-year reprocessed temperature salinity, sea surface height, geostrophic current and mixed layer depth climatology on 1/8 degree regular grid and 50 depth levels ({start_year}-{end_year}).",
+        "source": "Numerical models: Multiple Linear Regression, Optimal Interpolation. In-situ observations: Copernicus In Situ TAC (including Argo, XBT, CTD and moorings) Copernicus Sea Level TAC, CNES-CLS22 Mean Dynamic Topography, OSTIA Sea Surface Temperature Analysis, Copernicus MOB TAC (Sea Surface Salinity), and World Ocean Atlas 2018 (WOA18).",
+        "dataset_type": "observation",
+        "product_type": "climatology",
+        "product_version": "2.0",
+        "institution": "Copernicus Marine Service, Mercator Ocean International, France",
+        "citation": "Multi Observation Global Ocean 3D Temperature Salinity Height Geostrophic Current and MLD. E.U. Copernicus Marine Service Information (CMEMS). Marine Data Store (MDS). DOI: 10.48670/moi-00052 (Accessed on 21 04 2026).",
+        "references": "Guinehut S., A.-L. Dhomps, G. Larnicol and P.-Y. Le Traon, 2012: High resolution 3D temperature and salinity fields derived from in situ and satellite observations. Ocean Sci., 8(5):845-857. Mulet, S., M.-H. Rio, A. Mignot, S. Guinehut and R. Morrow, 2012: A new estimate of the global 3D geostrophic ocean circulation based on satellite data and in-situ measurements. Deep Sea Research Part II : Topical Studies in Oceanography, 77-80(0):70-81.",
+        "acknowledgement": "Generated using E.U. Copernicus Marine Service Information; https://doi.org/10.48670/moi-00052.",
+        "license": "ARMOR3D data were obtained from https://doi.org/10.48670/moi-00052, and are provided under the Copernicus Marine Environment Monitoring Service Service Level Agreement (SLA) https://marine.copernicus.eu/user-corner/service-commitments-and-licence?pk_vid=42ac3e352be888641780994034c3bb6e",
+        "doi": "10.48670/moi-00052",
+        "platform": "gr",
+        "horizontal_grid_type": "regular rectilinear",
+        "horizontal_grid_resolution": "0.125 degree",
+        "vertical_grid_type": "z",
+        "vertical_grid_coordinate": "depth",
+        "vertical_grid_levels": 50,
+        "aggregation": "mean",
+        "aggregation_frequency": "monthly",
+        "status": "completed",
+        "update_frequency": "None",
+        "bbox": "[-180.0, 180.0, -90.0, 90.0]",
+    }
+
+message = f"Updated ARMOR3D Monthly Climatology ({start_year}-{end_year}) -> root group attributes."
+
+update_icechunk_global_attrs(
+        credentials_filepath=credentials_filepath,
+        bucket=bucket,
+        prefix=prefix,
+        attrs=attrs,
+        commit_message=message,
+        )

OceanDataStore/data/update_noc_npd_era5v1_attrs.py

@@ -0,0 +1,172 @@
+# =========================================================
+# update_noc_npd_era5v1_attrs.py
+#
+# Script to update global and variable attributes in NOC
+# Near-Present Day ERA5v1 Icechunk repositories.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# =========================================================
+import logging
+
+from OceanDataStore.cli import initialise_logging
+from OceanDataStore.data.utils import update_icechunk_global_attrs
+
+
+def main(credentials_filepath: str,
+         bucket: str,
+         config_name: str,
+         nemo_config_name: str,
+         platform: str,
+         agg: str,
+         prefix_list: list
+         ) -> None:
+        # ========== Initialise OceanDataStore Logging ========== #
+        initialise_logging()
+
+        # ========= Update global attributes ========= #
+        for prefix in prefix_list:
+                logging.info(f"In Progress: Updating global attributes for {config_name} {prefix}...")
+
+                # Define aggregation frequency from prefix:
+                if "1y" in prefix:
+                        agg_freq = "annual"
+                elif "1m" in prefix:
+                        agg_freq = "monthly"
+                elif "5d" in prefix:
+                        agg_freq = "5-daily"
+                else:
+                        raise ValueError(f"Unable to determine aggregation frequency from prefix: {prefix}")
+                
+                # Define dimensionality from prefix:
+                if "_3d" in prefix:
+                        dimensionality = "3-dimensional"
+                elif "_4d" in prefix:
+                        dimensionality = "4-dimensional"
+                else:
+                        dimensionality = ""
+
+                # Define grid type from prefix:
+                if "T" in prefix:
+                        grid = "T-grid"
+                        variable_type = "scalar variables"
+                elif "U" in prefix:
+                        grid = "U-grid"
+                        variable_type = "vector variables"
+                elif "V" in prefix:
+                        grid = "V-grid"
+                        variable_type = "vector variables"
+                elif "W" in prefix:
+                        grid = "W-grid"
+                        variable_type = "vector variables"
+                elif "S" in prefix:
+                        grid = ""
+                        variable_type = "scalar variables"
+                elif "I" in prefix:
+                        grid = "T-grid"
+                        variable_type = "sea-ice variables"
+                else:
+                        raise ValueError(f"Unable to determine grid type from prefix: {prefix}")
+
+                # Define resolution from nemo_config_name:
+                if "eORCA12" in nemo_config_name:
+                        horizontal_grid_resolution = "1/12 degree"
+                elif "eORCA025" in nemo_config_name:
+                        horizontal_grid_resolution = "1/4 degree"
+                elif "eORCA1" in nemo_config_name:
+                        horizontal_grid_resolution = "1 degree"
+                else:
+                        raise ValueError(f"Unable to determine horizontal grid resolution from NEMO configuration name: {nemo_config_name}")
+
+                attrs = {
+                        "Conventions": "CF-1.6",
+                        "title": f"National Oceanography Centre Near-Present Day (NPD) {horizontal_grid_resolution} global ocean physics & sea-ice hindcast.", 
+                        "description": f"NOC Near-Present Day {agg_freq} {agg} global ocean physics & sea-ice hindcast forced using bias-corrected ERA5 atmospheric reanalysis {dimensionality} {variable_type} stored on the native {nemo_config_name} curvilinear NEMO model {grid}.",
+                        "dataset_type": "model",
+                        "product_type": "timeseries",
+                        "product_version": "1.0",
+                        "institution": "National Oceanography Centre, UK",
+                        "citation": "Blaker, A. T., Tooth, O. J., Palmiéri, J., Coward, A. C., and Mecking, J. (2025). NOC-MSM/NOC_Near_Present_Day: v0.9.0 (v0.9.0). Zenodo. https://doi.org/10.5281/zenodo.15310354.",
+                        "references": "Blaker, A.T., Tooth, O.J., Palmiéri, J., Coward, A.C., & Mecking, J. (2025). NOC-MSM/NOC_Near_Present_Day: v0.9.0 (v0.9.0). Zenodo. https://doi.org/10.5281/zenodo.15310354. Guiavarc'h, C., Storkey, D., Blaker, A. T., Blockley, E., Megann, A., Hewitt, H., Bell, M. J., Calvert, D., Copsey, D., Sinha, B., Moreton, S., Mathiot, P., and An, B.: GOSI9: UK Global Ocean and Sea Ice configurations, Geosci. Model Dev., 18, 377-403, https://doi.org/10.5194/gmd-18-377-2025, 2025.",
+                        "acknowledgement": "NOC Near-Present Day Documentation available at: https://noc-msm.github.io/NOC_Near_Present_Day/",
+                        "license": "UK Open Government License v3.0",
+                        "doi": "pending",
+                        "platform": platform,
+                        "horizontal_grid_type": "curvilinear",
+                        "horizontal_grid_resolution": horizontal_grid_resolution,
+                        "vertical_grid_type": "zps",
+                        "vertical_grid_coordinate": "depth with partial step topography",
+                        "vertical_grid_levels": 75,
+                        "aggregation": agg,
+                        "aggregation_frequency": agg_freq,
+                        "status": "ongoing",
+                        "update_frequency": "quarterly",
+                        "bbox": "[-180.0, 180.0, -90.0, 90.0]",
+                        "ocean_component": "NEMO v4.2.2",
+                        "sea_ice_component": "SI3 v4.0",
+                        "biogeochemistry_component": "None",
+                        "atmospheric_component": "None",
+                        "atmospheric_forcing": "ERA5 v1",
+                        "variant": "r1i1c1f1",
+                }
+
+                message = f"Updated {config_name} {agg_freq} {agg} -> root group attributes."
+
+                update_icechunk_global_attrs(
+                        credentials_filepath=credentials_filepath,
+                        bucket=bucket,
+                        prefix=prefix,
+                        attrs=attrs,
+                        commit_message=message,
+                        )
+
+                logging.info(f"Completed: Updated global attributes for {config_name} {prefix}.")
+        
+
+if __name__ == "__main__":
+        # ========= Define Shared Inputs ========= #
+        # Define credential to write to JASMIN OS:
+        credentials_filepath = '.../credentials/jasmin_os_credentials.json'
+
+        # Define NPD configuration propeties:
+        bucket = "npd-eorca12-era5v1"
+        config_name = "NPD eORCA12 ERA5v1"
+        nemo_config_name = "eORCA12"
+        agg = "mean"
+        platform = "gn"
+
+        # -- eORCA1 --- #
+        # prefix_list = ["I1m", "I1y",
+        #                "S1m", "S1y",
+        #                "T1m", "T1y",
+        #                "U1m", "U1y",
+        #                "V1m", "V1y",
+        #                "W1m", "W1y",
+        #                ]
+        
+        # -- eORCA025 --- #
+        # prefix_list = ["I1m_3d", "I1y_3d", "I5d_3d",
+        #                "S1m_1d", "S1y_1d", "S5d_1d",
+        #                "T1m_3d", "T1m_4d", "T1y_3d", "T1y_4d", "T5d_3d", "T5d_4d",
+        #                "U1m_3d", "U1m_4d", "U1y_3d", "U1y_4d", "U5d_3d", "U5d_4d",
+        #                "V1m_3d", "V1m_4d", "V1y_3d", "V1y_4d", "V5d_3d", "V5d_4d",
+        #                "W1m_4d", "W1y_4d", "W5d_4d"
+        #                ]
+
+        # -- eORCA12 --- #
+        prefix_list = ["I1m_3d", "I1y_3d",
+                       "S1m_1d", "S1y_1d",
+                       "T1m_3d", "T1m_4d", "T1y_3d", "T1y_4d",
+                       "U1m_3d", "U1m_4d", "U1y_3d", "U1y_4d",
+                       "V1m_3d", "V1m_4d", "V1y_3d", "V1y_4d",
+                       "W1m_4d", "W1y_4d",
+                       ]
+
+        # ========= Run Main Function ========= #
+        main(credentials_filepath=credentials_filepath,
+             bucket=bucket,
+             config_name=config_name,
+             nemo_config_name=nemo_config_name,
+             platform=platform,
+             agg=agg,
+             prefix_list=prefix_list
+             )