OceanDataStore 0.3.0__py3-none-any.whl

OceanDataStore/data/ARMOR3D/update_ARMOR3D_P1m-m_monthly_to_os.py

@@ -0,0 +1,143 @@
+# =========================================================
+# update_EN4.2.2_analyses_g10_to_os.py
+#
+# Script to update EN.4.2.2 analyses in Icechunk repository
+# in JASMIN cloud object storage.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# =========================================================
+import logging
+
+import xarray as xr
+import zarr
+
+from OceanDataStore.cli import initialise_logging, update_icechunk
+from OceanDataStore.data.utils import (
+    compute_dx,
+    compute_dy,
+    compute_cell_area,
+    compute_cell_thickness,
+    compute_land_sea_mask,
+)
+
+logger = logging.getLogger(__name__)
+
+
+def main():
+    # ========== Initialise OceanDataStore Logging ========== #
+    initialise_logging()
+
+    # ========== Prepare Data ========== #
+    # Open complete ARMOR3D REP monthly climatology dataset:
+    filepath = "/dssgfs01/scratch/otooth/npd_data/observations/ARMOR3D/armor-3d_rep_monthly_NA_*.zarr"
+    ds = xr.open_mfdataset(filepath, compat="no_conflicts", data_vars="all", engine="zarr")
+    logging.info("-> Completed: Opened ARMOR-3D REP monthly climatology dataset from Zarr stores.")
+
+    # Rename variables to standard names:
+    ds = ds.rename({"to": "thetao",
+                    "so": "so",
+                    })
+
+    # Update global attributes:
+    ds.attrs.clear()
+
+    ds = ds.assign_attrs({
+        "Conventions": "CF-1.0",
+        "title": "Multi Observation Global Ocean 3D Temperature Salinity Height Geostrophic Current and MLD.",
+        "description": "Multi Observation Global Ocean ARMOR3D multi-year reprocessed temperature salinity, sea surface height, geostrophic current and mixed layer depth monthly timeseries on 1/8 degree regular grid and 50 depth levels.",
+        "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": "timeseries",
+        "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": "ongoing",
+        "update_frequency": "quarterly",
+        "bbox": "[-180.0, 180.0, -90.0, 90.0]",
+    })
+
+    # Add ancillary variables:
+    ds['mask'] = compute_land_sea_mask(ds['thetao'].isel(time=0, depth=0))
+    ds['dx'] = compute_dx(ds)
+    ds['dy'] = compute_dy(ds)
+    ds['cell_area'] = compute_cell_area(ds)
+    # Custom ancillary variables:
+    ds['cell_thickness'] = compute_cell_thickness(ds)
+    ds['cell_volume'] = ds['cell_thickness'] * ds['cell_area']
+
+    # Update attributes for custom ancillary variables:
+    ds['cell_volume'].attrs.update({
+        'long_name': "Grid-Cell Volume",
+        'standard_name': "cell_volume",
+        'units': "m3",
+    })
+
+    # ========== Send to Icechunk Repository ========== #
+    bucket = "armor3d"
+    store_credentials_json = ".../credentials/jasmin_os_credentials.json"
+    branch = "main"
+    config_kwargs = {
+            "temporary_directory":".../OceanDataStore/OceanDataStore/data/ARMOR3D/",
+            "local_directory":".../OceanDataStore/OceanDataStore/data/ARMOR3D/"
+        }
+    cluster_kwargs = {
+            "n_workers" : 25,
+            "threads_per_worker" : 1,
+            "memory_limit":"4GB"
+        }
+    
+    # Optimise chunk sizes for spatial analysis:
+    for var in ds.data_vars:
+        if ds[var].ndim == 4:
+            ds[var] = ds[var].chunk({'time': 1, 'depth': 3, 'latitude': 689, 'longitude': 1440})
+            ds[var].encoding['chunks'] = (1, 3, 689, 1440)
+        elif ds[var].ndim == 3:
+            if "time" in ds[var].dims:
+                ds[var] = ds[var].chunk({'time': 1, 'latitude': 1378, 'longitude': 2880})
+                ds[var].encoding['chunks'] = (1, 1378, 2880)
+            elif "depth" in ds[var].dims:
+                ds[var] = ds[var].chunk({'depth': 10, 'latitude': 1378, 'longitude': 2880})
+                ds[var].encoding['chunks'] = (10, 1378, 2880)
+        elif (ds[var].ndim == 2):
+            if "latitude" in ds[var].dims and "longitude" in ds[var].dims:
+                ds[var] = ds[var].chunk({'latitude': 1378, 'longitude': 2880})
+                ds[var].encoding['chunks'] = (1378, 2880)
+        elif ds[var].ndim == 1:
+            ds[var] = ds[var].chunk({'depth': 50})
+            ds[var].encoding['chunks'] = (50,)
+
+    # Update variable encodings:
+    blosccodec = zarr.codecs.BloscCodec(cname="zstd", clevel=5, shuffle=zarr.codecs.BloscShuffle.shuffle)
+    for var in list(ds.data_vars) + list(ds.coords):
+        ds[var].encoding['compressors'] = [blosccodec]
+
+    # Define prefix and commit message based on period:
+    prefix = "armor3d_global_my_monthly"
+    commit_message = "Added ARMOR3D Global Monthly (2000-2024)."
+
+    update_icechunk(
+        file=ds.sel(time=slice("2000-01-01", "2024-12-31")),
+        bucket=bucket,
+        object_prefix=prefix,
+        store_credentials_json=store_credentials_json,
+        append_dim='time',
+        branch=branch,
+        commit_message=commit_message,
+        dask_config_kwargs=config_kwargs,
+        dask_cluster_kwargs=cluster_kwargs,
+        )
+
+if __name__ == "__main__":
+    main()

OceanDataStore/data/EN.4.2.2/create_EN4.2.2_analysis_g10_climatology.py

@@ -0,0 +1,162 @@
+# =========================================================
+# create_EN4.2.2_analysis_g10_climatology.py
+#
+# Script to create EN.4.2.2 analyses climatologies for
+# climate normals (1971-2000, 1981-2010, 1991-2020) and
+# write to local netCDF files.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# =========================================================
+import logging
+
+import numpy as np
+import xarray as xr
+import zarr
+
+from OceanDataStore.cli import initialise_logging
+from OceanDataStore.data.utils import (
+    compute_cell_area,
+    compute_dx,
+    compute_dy,
+    compute_land_sea_mask,
+)
+
+logger = logging.getLogger(__name__)
+
+def main():
+    # ========== Initialize Logging and Print Banner ========== #
+    initialise_logging()
+
+    # ========== Prepare Data ========== #
+    # Open complete EN.4.2.2 analysis dataset:
+    filepath = "/dssgfs01/scratch/otooth/npd_data/observations/EN.4.2.2/EN.4.2.2.f.analysis.g10.*.nc"
+    ds = xr.open_mfdataset(filepath, combine="by_coords", data_vars="all", engine="netcdf4")
+    logging.info("-> Completed: Opened EN.4.2.2 analysis dataset from netCDF files.")
+
+    # Standardise coordinate dimension names:
+    ds = ds.rename({"lon": "longitude", "lat": "latitude"})
+
+    # Update longitude coordinates to be in the range [-180, 180]:
+    ds = ds.assign_coords(
+        longitude=((ds["longitude"] + 180) % 360) - 180
+    )
+    ds = ds.sortby("longitude")
+
+    # Rename variables to standard names:
+    ds = ds.rename({"temperature": "thetao",
+                    "salinity": "so",
+                    "temperature_uncertainty": "thetao_uncertainty",
+                    "salinity_uncertainty": "so_uncertainty",
+                    "temperature_observation_weights": "thetao_obs_weights",
+                    "salinity_observation_weights": "so_obs_weights"
+                    })
+    
+    # Update variable attributes:
+    ds["thetao"].attrs.update({
+        "long_name": "Potential Temperature",
+    })
+    ds["so"].attrs.update({
+        "long_name": "Practical Salinity",
+    })
+    ds["thetao_uncertainty"].attrs.update({
+        "long_name": "Potential Temperature Error Standard Deviation",
+    })
+    ds["so_uncertainty"].attrs.update({
+        "long_name": "Practical Salinity Error Standard Deviation",
+    })
+    ds["thetao_obs_weights"].attrs.update({
+        "long_name": "Potential Temperature Observation Weights",
+    })
+    ds["so_obs_weights"].attrs.update({
+        "long_name": "Practical Salinity Observation Weights",
+    })
+
+    # Update global attributes:
+    ds.attrs.clear()
+
+    ds = ds.assign_attrs({
+        "Conventions": "CF-1.0",
+        "title": "EN.4.2.2 ocean temperature and salinity monthly climatology.",
+        "description": "EN.4.2.2 quality controlled ocean temperature and salinity monthly climatology from objective analyses with uncertainty estimates using Gouretski and Reseghetti (2010) corrections.",
+        "source": "Numerical models: Objective Analysis. In-situ observations: Argo, Arctic Synoptic Basin-wide Oceanography (ASBO) project, Global Temperature and Salinity Profile Programme (GTSPP), and World Ocean Database 2018 (WOD18).",
+        "dataset_type": "observation",
+        "product_type": "climatology",
+        "product_version": "1.0",
+        "institution": "Met Office, UK",
+        "citation": "Good, S. A., Martin, M. J., and Rayner, N. A., 2013. EN4: quality controlled ocean temperature and salinity profiles and monthly objective analyses with uncertainty estimates, Journal of Geophysical Research: Oceans, 118, 6704-6716, doi:10.1002/2013JC009067.",
+        "references": "Gouretski, V., and Reseghetti, F., 2010: On depth and temperature biases in bathythermograph data: development of a new correction scheme based on analysis of a global ocean database. Deep-Sea Research I, 57, 6. doi:10.1016/j.dsr.2010.03.011.",
+        "acknowledgement": "None",
+        "license": "EN.4.2.2 data were obtained from https://www.metoffice.gov.uk/hadobs/en4/ and are © Crown Copyright, Met Office, [2026], provided under a Non-Commercial Government Licence http://www.nationalarchives.gov.uk/doc/non-commercial-government-licence/version/2/.",
+        "doi": "None",
+        "platform": "gr",
+        "horizontal_grid_type": "regular rectilinear",
+        "horizontal_grid_resolution": "1 degree",
+        "vertical_grid_type": "z",
+        "vertical_grid_coordinate": "depth",
+        "vertical_grid_levels": 42,
+        "aggregation": "mean",
+        "aggregation_frequency": "monthly",
+        "status": "completed",
+        "update_frequency": "None",
+        "bbox": "[-180.0, 180.0, -90.0, 90.0]",
+    })
+
+    logging.info("-> Completed: Updated EN.4.2.2 analysis CF-metadata.")
+
+    # -- Calculate climate normal monthly climatologies -- #
+    output_dir = "/dssgfs01/scratch/otooth/npd_data/observations/EN.4.2.2/climatology/"
+    start_years = [1971, 1981, 1991]
+    end_years = [2000, 2010, 2020]
+
+    for start_year, end_year in zip(start_years, end_years):
+        logging.info(f"In Progress: Calculating monthly climatology for {start_year}-{end_year} climate normal period...")
+        # Calculate monthly climatology for the specified period:
+        ds_climatology = ds.sel(time=slice(f'{start_year}-01', f'{end_year}-12')).groupby('time.month').mean()
+
+        # Add ancillary variables:
+        ds_climatology['mask'] = compute_land_sea_mask(ds['thetao'].isel(time=0, depth=0))
+        ds_climatology['dx'] = compute_dx(ds)
+        ds_climatology['dy'] = compute_dy(ds)
+        ds_climatology['cell_area'] = compute_cell_area(ds)
+        # Custom ancillary variables:
+        ds_climatology['cell_thickness'] = (ds_climatology['depth_bnds'].isel(bnds=1) - ds_climatology['depth_bnds'].isel(bnds=0)).isel(month=0)
+        ds_climatology['cell_volume'] = ds_climatology['cell_thickness'] * ds_climatology['cell_area']
+
+        # Update attributes for custom ancillary variables:
+        ds_climatology['cell_thickness'].attrs.update({
+            'long_name': "Grid-Cell Thickness",
+            'standard_name': "cell_thickness",
+            'units': "m",
+        })
+        ds_climatology['cell_volume'].attrs.update({
+            'long_name': "Grid-Cell Volume",
+            'standard_name': "cell_volume",
+            'units': "m3",
+        })
+
+        # Update time bounds to reflect climatological period:
+        ds_climatology['time_bnds'][:, 0] = np.arange(f'{start_year}-01', f'{start_year+1}-01', dtype='datetime64[M]')
+        ds_climatology['time_bnds'][:, 1] = np.arange(f'{end_year}-01', f'{end_year+1}-01', dtype='datetime64[M]')
+        logging.info(f"-> Completed: Calculated monthly climatology for {start_year}-{end_year} climate normal period.")
+
+        # Update title attribute to reflect climatological period:
+        ds_climatology.attrs['title'] = f"EN.4.2.2 ocean temperature and salinity monthly climatology ({start_year}-{end_year})."
+
+        ds_climatology.attrs['description'] = f"EN.4.2.2: quality controlled ocean temperature and salinity monthly climatology ({start_year}-{end_year}) from objective analyses with uncertainty estimates using Gouretski and Reseghetti (2010) corrections."
+
+        # Update variable encodings:
+        blosccodec = zarr.codecs.BloscCodec(cname="zstd", clevel=3, shuffle=zarr.codecs.BloscShuffle.shuffle)
+        for var in list(ds_climatology.data_vars) + list(ds_climatology.coords):
+            ds_climatology[var].encoding['compressors'] = [blosccodec]
+
+        # Write monthly climatology to netCDF:
+        output_filepath = f"{output_dir}EN.4.2.2.f.analysis.g10.{start_year}_{end_year}_monthly_climatology.nc"
+        ds_climatology.to_netcdf(output_filepath)
+        logging.info(f"-> Completed: Saved monthly climatology for {start_year}-{end_year} climate normal period to {output_filepath}.")
+
+    # -- Close EN.4.2.2 analysis datasets -- #
+    ds_climatology.close()
+    ds.close()
+
+if __name__ == "__main__":
+    main()

OceanDataStore/data/EN.4.2.2/download_EN4.2.2_analysis_g10_data.sh

@@ -0,0 +1,51 @@
+#!/bin/bash
+
+# ----------------------------------------------------------------
+# download_EN4.2.2_analyses_g10_data.sh
+#
+# Description: Download the EN.4.2.2 analyses.g10 dataset from the
+# Met Office Hadley Centre EN.4.2.2 website:
+# http://www.metoffice.gov.uk/hadobs/en4
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# Created On: 2026-05-27
+# ----------------------------------------------------------------
+set -euo pipefail
+
+# --- Inputs --- #
+output_dir="/dssgfs01/scratch/otooth/npd_data/observations/EN.4.2.2/"
+
+# -- Defaults -- #
+base_url="http://www.metoffice.gov.uk/hadobs/en4/data/en4-2-1"
+
+# --- Main Script --- #
+echo "==================================================="
+echo "          Downloading EN.4.2.2 Analyses"
+echo "                     v0.1.0"
+echo "               Oliver J. Tooth, NOC"
+echo "==================================================="
+echo "In Progress: Downloading EN.4.2.2 analyses dataset..."
+# Iterate over years:
+for yr in {1990..2026}; do
+    # Construct URL for current year:
+    if [ $yr -ge 2021 ]; then
+        url="$base_url/EN.4.2.2.analyses.g10.${yr}.zip"
+    else 
+        url="$base_url/EN.4.2.2/EN.4.2.2.analyses.g10.${yr}.zip"
+    fi
+
+    # Download and unzip file if not in output directory:
+    nc_files=("${output_dir}/EN.4.2.2.f.analysis.g10.${yr}"*.nc)
+    filepath="$output_dir/$(basename $url)"
+    if [ ${#nc_files[@]} -ne 12 ]; then
+        wget -P $output_dir $url
+        echo "-> Completed: Downloaded $filepath."
+
+        unzip "$filepath" -d $output_dir
+        echo "-> Completed: Unzipped $filepath."
+    else
+        echo "-> Skipping Download: NetCDF files for ${yr} already exist in $output_dir."
+    fi
+done
+
+echo "======================================="

OceanDataStore/data/EN.4.2.2/run_send_EN4.2.2_analysis_g10_climatology_to_os.slurm

@@ -0,0 +1,32 @@
+#!/bin/bash
+#SBATCH --job-name=en4.2.2_analysis_g10_climatology
+#SBATCH --partition=test
+#SBATCH --time=00:20:00
+#SBATCH --ntasks-per-core=1
+#SBATCH --ntasks-per-node=64
+#SBATCH --ntasks-per-socket=32
+#SBATCH --nodes=1
+
+# ==============================================================
+# run_send_EN4.2.2_analysis_g10_climatology_to_os.slurm
+#
+# Description: SLURM script to send the EN.4.2.2 analysis g10 
+# climatology datasets to Icechunk repository.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# Created On: 2026-05-29
+#
+# ==============================================================
+set -euo pipefail
+
+# -- Python Environment -- #
+# Activate miniconda environment:
+source .../miniforge3/bin/activate
+conda activate env_ods
+
+# -- Send EN.4.2.2.analysis.g10 climatology datasets to JASMIN OS -- #
+echo "In Progress: Sending EN.4.2.2.analysis.g10 climatology to Icechunk..."
+
+python3 send_EN4.2.2_analysis_g10_monthly_climatology_to_os.py
+
+echo "Completed: Sent EN.4.2.2.analysis.g10 climatology to Icechunk."

OceanDataStore/data/EN.4.2.2/run_send_EN4.2.2_analysis_g10_monthly_to_os.slurm

@@ -0,0 +1,32 @@
+#!/bin/bash
+#SBATCH --job-name=en4.2.2_analysis_g10_monthly
+#SBATCH --partition=compute
+#SBATCH --time=01:00:00
+#SBATCH --ntasks-per-core=1
+#SBATCH --ntasks-per-node=64
+#SBATCH --ntasks-per-socket=32
+#SBATCH --nodes=1
+
+# ==============================================================
+# run_send_EN4.2.2_analysis_g10_monthly_to_os.slurm
+#
+# Description: SLURM script to send the EN.4.2.2 analysis g10 
+# monthly datasets to Icechunk repository.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# Created On: 2026-05-29
+#
+# ==============================================================
+set -euo pipefail
+
+# -- Python Environment -- #
+# Activate miniconda environment:
+source .../miniforge3/bin/activate
+conda activate env_ods
+
+# -- Send EN.4.2.2.analysis.g10 monthly to JASMIN OS -- #
+echo "In Progress: Sending EN.4.2.2.analysis.g10 monthly to Icechunk..."
+
+python3 send_EN4.2.2_analysis_g10_monthly_to_os.py
+
+echo "Completed: Sent EN.4.2.2.analysis.g10 monthly to Icechunk."

OceanDataStore/data/EN.4.2.2/run_update_EN4.2.2_analysis_g10_monthly_to_os.slurm

@@ -0,0 +1,32 @@
+#!/bin/bash
+#SBATCH --job-name=en4.2.2_analysis_g10_monthly
+#SBATCH --partition=compute
+#SBATCH --time=01:00:00
+#SBATCH --ntasks-per-core=1
+#SBATCH --ntasks-per-node=64
+#SBATCH --ntasks-per-socket=32
+#SBATCH --nodes=1
+
+# ==============================================================
+# run_update_EN4.2.2_analysis_g10_monthly_to_os.slurm
+#
+# Description: SLURM script to update the EN.4.2.2 analysis g10 
+# monthly datasets in Icechunk repository.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# Created On: 2026-05-29
+#
+# ==============================================================
+set -euo pipefail
+
+# -- Python Environment -- #
+# Activate miniconda environment:
+source .../miniforge3/bin/activate
+conda activate env_ods
+
+# -- Update EN.4.2.2.analysis.g10 monthly in JASMIN OS -- #
+echo "In Progress: Updating EN.4.2.2.analysis.g10 monthly in Icechunk..."
+
+python3 update_EN4.2.2_analysis_g10_monthly_to_os.py
+
+echo "Completed: Updated EN.4.2.2.analysis.g10 monthly in Icechunk."

OceanDataStore/data/EN.4.2.2/send_EN4.2.2_analysis_g10_monthly_climatology_to_os.py

@@ -0,0 +1,76 @@
+# =========================================================
+# send_EN4.2.2_analysis_g10_climatology_to_os.py
+#
+# Script to write EN.4.2.2 analysis climatologies to
+# Icechunk repositories in JASMIN cloud object storage.
+#
+# Created By: Ollie Tooth (oliver.tooth@noc.ac.uk)
+# =========================================================
+import logging
+
+import xarray as xr
+
+from OceanDataStore.cli import send_to_icechunk, initialise_logging
+
+logger = logging.getLogger(__name__)
+
+
+def main():
+    # ========== Initialise OceanDataStore Logging ========== #
+    initialise_logging()
+
+    # ========== Prepare Data ========== #
+    # Open EN.4.2.2 analysis climatology datasets:
+    filepaths = [
+        "/dssgfs01/scratch/otooth/npd_data/observations/EN.4.2.2/climatology/EN.4.2.2.f.analysis.g10.1971_2000_monthly_climatology.nc",
+        "/dssgfs01/scratch/otooth/npd_data/observations/EN.4.2.2/climatology/EN.4.2.2.f.analysis.g10.1981_2010_monthly_climatology.nc",
+        "/dssgfs01/scratch/otooth/npd_data/observations/EN.4.2.2/climatology/EN.4.2.2.f.analysis.g10.1991_2020_monthly_climatology.nc"
+        ]
+    
+    # Define start & end years of climatology periods:
+    start_years = [1971, 1981, 1991]
+    end_years = [2000, 2010, 2020]
+
+    # ========== Send to Icechunk Repository ========== #
+    bucket = "en4.2.2"
+    exists = False
+    store_credentials_json = ".../credentials/jasmin_os_credentials.json"
+    branch = "main"
+    variable_commits = True
+    config_kwargs = {
+            "temporary_directory":".../OceanDataStore/OceanDataStore/data/EN.4.2.2/",
+            "local_directory":".../OceanDataStore/OceanDataStore/data/EN.4.2.2/"
+        }
+    cluster_kwargs = {
+            "n_workers" : 10,
+            "threads_per_worker" : 1,
+            "memory_limit":"3GB"
+        }
+    
+    for filepath, start_yr, end_yr in zip(filepaths, start_years, end_years):
+        # Open EN.4.2.2 analysis climatology dataset:
+        ds = xr.open_dataset(filepath, engine='netcdf4')
+
+        # Optimise chunk sizes for spatial analysis:
+        ds = ds.chunk({'month': 1, 'depth': 20, 'latitude': 173, 'longitude': 360})
+
+        # Define prefix and commit message based on climatology period:
+        prefix = f"en4.2.2_analysis_g10_{start_yr}_{end_yr}_monthly_climatology"
+        commit_message = f"Added EN.4.2.2.analysis.g10 climatology ({start_yr}-{end_yr})."
+
+        send_to_icechunk(
+            file=ds,
+            bucket=bucket,
+            object_prefix=prefix,
+            store_credentials_json=store_credentials_json,
+            exists=exists,
+            append_dim='month',
+            branch=branch,
+            commit_message=commit_message,
+            variable_commits=variable_commits,
+            dask_config_kwargs=config_kwargs,
+            dask_cluster_kwargs=cluster_kwargs,
+            )
+
+if __name__ == "__main__":
+    main()