ERA-Interim Land

esmvaltool/cmorizers/obs/cmor_config/ERA-Interim-Land.yml

@@ -0,0 +1,26 @@
+---
+# Global attributes of NetCDF file
+attributes:
+  dataset_id: ERA-Interim-Land
+  project_id: OBS6
+  tier: 3
+  version: '1'
+  modeling_realm: reanaly
+  source: 'https://apps.ecmwf.int/datasets/data/interim-land/type=an/'
+  reference: 'era-interim-land'
+  comment: |
+    'Contains modified Copernicus Climate Change Service Information {year}'
+
+# Variables to CMORize
+variables:
+  sm_daily:
+    short_name: sm
+    mip: CFday
+    raw: swvl1
+    file: 'eraintland_swvl1_????.nc'
+  sm_monthly:
+    short_name: sm
+    mip: Lmon
+    raw: swvl1
+    file: 'eraintland_swvl1_????.nc'
+

esmvaltool/cmorizers/obs/cmorize_obs_era_interim.py

@@ -95,7 +95,7 @@
 import numpy as np
 
 from esmvalcore.cmor.table import CMOR_TABLES
-from esmvalcore.preprocessor import daily_statistics
+from esmvalcore.preprocessor import daily_statistics, monthly_statistics
 
 from . import utilities as utils
 
@@ -145,8 +145,8 @@ def _fix_coordinates(cube, definition):
     if 'height10m' in definition.dimensions:
         utils.add_scalar_height_coord(cube, 10.)
 
-    for axis in 'T', 'X', 'Y', 'Z':
-        coord_def = definition.coordinates.get(axis)
+    for definition_axis,axis in zip(['time','longitude','latitude'],['T','X','Y']):
+        coord_def = definition.coordinates.get(definition_axis)
         if coord_def:
             coord = cube.coord(axis=axis)
             if axis == 'T':
@@ -178,6 +178,17 @@ def _fix_monthly_time_coord(cube):
     coord.points = 0.5 * (start + end)
     coord.bounds = np.column_stack([start, end])
 
+def _compute_monthly(cube):
+    """Convert various frequencies to daily frequency.
+
+    ERA-Interim-Land is in 6hr freq need to convert to monthly
+
+    """
+    cube = monthly_statistics(cube, operator='mean')
+    # Remove monthly statistics aux coordinates
+    cube.remove_coord(cube.coord('month_number'))
+    cube.remove_coord(cube.coord('year'))
+    return cube
 
 def _compute_daily(cube):
     """Convert various frequencies to daily frequency.
@@ -308,7 +319,8 @@ def _extract_variable(in_files, var, cfg, out_dir):
 
     # Set correct names
     cube.var_name = definition.short_name
-    cube.standard_name = definition.standard_name
+    if definition.standard_name:
+        cube.standard_name = definition.standard_name
     cube.long_name = definition.long_name
 
     _fix_units(cube, definition)
@@ -318,15 +330,23 @@ def _extract_variable(in_files, var, cfg, out_dir):
 
     cube = _fix_coordinates(cube, definition)
 
-    if 'mon' in var['mip']:
-        _fix_monthly_time_coord(cube)
-
-    if 'day' in var['mip']:
-        cube = _compute_daily(cube)
-
-    if 'fx' in var['mip']:
-        cube = iris.util.squeeze(cube)
-        cube.remove_coord('time')
+    if attributes['dataset_id']=='ERA-Interim':
+        if 'mon' in var['mip']:
+            _fix_monthly_time_coord(cube)
+        if 'day' in var['mip']:
+            cube = _compute_daily(cube)
+        if 'fx' in var['mip']:
+            cube = iris.util.squeeze(cube)
+            cube.remove_coord('time')
+
+    # Specific to ERA Interim Land
+    elif attributes['dataset_id']=='ERA-Interim-Land':
+        if 'mon' in var['mip']:
+            cube = _compute_monthly(cube)
+        if 'day' in var['mip']:
+            cube = _compute_daily(cube)
+    else:
+        raise ValueError("Unknown dataset_id for this script: {attributes['dataset_id']}")
 
     # Convert units if required
     cube.convert_units(definition.units)

esmvaltool/cmorizers/obs/cmorize_obs_era_interim_land.py

@@ -0,0 +1,20 @@
+"""ESMValTool CMORizer for ERA-Interim-Land data.
+
+Tier
+    Tier 3: restricted datasets (i.e., dataset which requires a registration
+ to be retrieved or provided upon request to the respective contact or PI).
+
+Source
+    https://apps.ecmwf.int/datasets/data/interim-land/type=fc/
+
+Last access
+    20191104
+
+Download and processing instructions
+     See script cmorize_obs_era_interim.py
+
+"""
+
+from .cmorize_obs_era_interim import cmorization
+
+__all__ = ['cmorization']

esmvaltool/config-references.yml

@@ -674,6 +674,7 @@ references:
   duveiller2018: "Duveiller et al., Sci. Data, doi:10.1038/sdata.2018.14, 2018."
   eppley-vgpm-modis: "Behrenfeld and Falkowski, Limnol. Oceanogr., doi:10.4319/lo.1997.42.1.0001, 1997."
   era-interim: "Dee et al., Q. J. Roy. Meteor. Soc., doi:10.1002/qj.828, 2011."
+  era-interim-land: "Balsamo et al. ERA-Interim/Land: a global land surface reanalysis data set, Hydrol. Earth Syst. Sci., 19, 389–407, https://doi.org/10.5194/hess-19-389-2015, 2015."
   era5: "Copernicus Climate Change Service (C3S) (2017): ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate. Copernicus Climate Change Service Climate Data Store (CDS)"
   esacci-aerosol: "Popp et al., Remote Sens., doi:10.3390/rs8050421, 2016."
   esacci-cloud: "Stengel et al., Earth Syst. Sci. Data, doi:10.5194/essd-9-881-2017, 2017."