You are viewing an old version of this page. View the current version.
Compare with Current
View Page History
« Previous
Version 39
Next »
In order to interpret CEMS-Flood (EFAS & GloFAS) model outputs in a correct way, you will need a number of auxiliary data, also referred to as the 'time invariant' data on the Climate Data Store (CDS). These data are part of the EFAS and GloFAS setups of the LISFLOOD model. This page summarises the auxiliary datasets made available through the CDS and how they should be used in combination with the model outputs.
The auxiliary data are provided as small (no larger than 8MB) netCDF files, and are available to download through the Climate Data Store (CDS) on the EFAS pages.
The following five datasets are available:
- Upstream area
- Elevation
- Soil depth (3 soil levels)
- Wilting point (3 soil levels)
- Field capacity (3 soil levels)
They are in the INSPIRE compliant ETRS89 Lambert Azimuthal Equal Area Coordinate Reference System (ETRS-LAEA), also referred to as EPSG:3035. More information on the EFAS horizontal resolution is available in the model output page.
It is important for users to make sure they download the auxiliary data consistent with the dataset version of interest. This is because the auxiliary data sometimes change with a new EFAS cycle; for example, the upstream area is specific to the hydrological model configuration used which can change between cycles.
To provide an example, users accessing forecast data issued between 2019-05-13 and 2020-03-04 should use v3.0 of the auxiliary data, as forecasts issued during that period were produced with EFAS v3.0. All the version dates are shown in the table below:
Versions | For forecast data issued from: |
---|
EFAS version 4.0 | 2020-10-14 12 UTC to present |
EFAS version 3.5 | 2020-03-05 12 UTC and 2020-10-13 12 UTC |
EFAS version 3.0 | 2019-05-13 12 UTC and 2020-03-04 12 UTC |
EFAS version 2.0 | before 2019-05-13 12 UTC |
The upstream area file denotes the total upstream area for each river pixel, expressed in . It is defined as the catchment area for each river segment, i.e. the total area that contributes with water to the river at the specific grid point (or river pixel). The upstream area always includes the area of the pixel it is associated to.
dimensions:
x = 1000 ;
y = 950 ;
variables:
double x(x) ;
x:units = "m" ;
x:long_name = "x coordinate of projection" ;
x:standard_name = "projection_x_coordinate" ;
double y(y) ;
y:units = "m" ;
y:long_name = "y coordinate of projection" ;
y:standard_name = "projection_y_coordinate" ;
char lambert_azimuthal_equal_area ;
lambert_azimuthal_equal_area:inverse_flattening = 298.257222101 ;
lambert_azimuthal_equal_area:spatial_ref = "PROJCS[\"ETRS89 / LAEA Europe\",GEOGCS[\"ETRS89\",DATUM[\"European_Terrestrial_Reference_System_1989\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6258\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4258\"]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_center\",52],PARAMETER[\"longitude_of_center\",10],PARAMETER[\"false_easting\",4321000],PARAMETER[\"false_northing\",3210000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AUTHORITY[\"EPSG\",\"3035\"]]" ;
lambert_azimuthal_equal_area:latitude_of_projection_origin = 52. ;
lambert_azimuthal_equal_area:semi_major_axis = 6378137. ;
lambert_azimuthal_equal_area:GeoTransform = "2500000 5000 0 5500000 0 -5000 " ;
lambert_azimuthal_equal_area:longitude_of_projection_origin = 10. ;
lambert_azimuthal_equal_area:long_name = "CRS definition" ;
lambert_azimuthal_equal_area:false_northing = 3210000. ;
lambert_azimuthal_equal_area:grid_mapping_name = "lambert_azimuthal_equal_area" ;
lambert_azimuthal_equal_area:false_easting = 4321000. ;
lambert_azimuthal_equal_area:longitude_of_prime_meridian = 0. ;
double ec_upArea(y, x) ;
ec_upArea:_FillValue = NaN ;
ec_upArea:long_name = "ec_upArea" ;
ec_upArea:missing_value = NaN ;
ec_upArea:grid_mapping = "lambert_azimuthal_equal_area" ;
// global attributes:
:version = "EFAS Model Version - 4.0" ;
}
The elevation (sometimes referred to as orography) file denotes the mean elevation for each pixel in the EFAS domain, expressed in m above sea level. It is defined for each grid point of the modelled domain.
dimensions:
x = 1000 ;
y = 950 ;
variables:
double x(x) ;
x:standard_name = "projection_x_coordinate" ;
x:long_name = "x coordinate of projection" ;
x:units = "Meter" ;
double y(y) ;
y:standard_name = "projection_y_coordinate" ;
y:long_name = "y coordinate of projection" ;
y:units = "Meter" ;
int laea ;
laea:grid_mapping_name = "lambert_azimuthal_equal_area" ;
laea:false_easting = 4321000. ;
laea:false_northing = 3210000. ;
laea:longitude_of_projection_origin = 10. ;
laea:latitude_of_projection_origin = 52. ;
laea:semi_major_axis = 6378137. ;
laea:inverse_flattening = 298.257223563 ;
laea:proj4_params = "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs" ;
laea:EPSG_code = "EPSG:3035" ;
float dem_mean(y, x) ;
dem_mean:standard_name = "dem_mean" ;
dem_mean:long_name = "dem_mean" ;
dem_mean:units = "-" ;
dem_mean:grid_mapping = "lambert_azimuthal_equal_area" ;
dem_mean:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
// global attributes:
:history = "Created Mon Feb 29 10:20:21 2016" ;
:conventions = "CF-1.6" ;
:source_software = "Python netCDF3_Classic" ;
:title = "Lisflood maps for European setting Januar 2015" ;
:keywords = "Lisflood, Europe" ;
:source = "Lisflood European maps - pb2015" ;
:institition = "JRC H01" ;
}
The files for soil depth, wilting capacity, and field capacity are collectively known as the 'soil moisture index files', and contain all the relevant information to calculate the soil moisture index from the volumentric soil moisture dataset model output. They are provided separately for each of the three soil layers represented in LISFLOOD.
More detail on the datasets and how they should be used to derive the soil moisture index is provided here.
dimensions:
y = 950 ;
x = 1000 ;
variables:
double y(y) ;
y:_FillValue = NaN ;
y:units = "Meter" ;
y:long_name = "Y coordinate of projection" ;
y:standard_name = "projection_Y_coordinate" ;
double x(x) ;
x:_FillValue = NaN ;
x:units = "Meter" ;
x:long_name = "x coordinate of projection" ;
x:standard_name = "projection_x_coordinate" ;
double soilLayer ;
soilLayer:_FillValue = NaN ;
soilLayer:long_name = "original GRIB coordinate for key: level(soilLayer)" ;
soilLayer:units = "1" ;
float sod(y, x) ;
sod:_FillValue = NaNf ;
sod:GRIB_paramId = 260367LL ;
sod:GRIB_dataType = "cf" ;
sod:GRIB_numberOfPoints = 950000LL ;
sod:GRIB_typeOfLevel = "soilLayer" ;
sod:GRIB_stepUnits = 1LL ;
sod:GRIB_stepType = "instant" ;
sod:GRIB_gridType = "lambert_azimuthal_equal_area" ;
sod:GRIB_NV = 0LL ;
sod:GRIB_cfName = "unknown" ;
sod:GRIB_cfVarName = "sod" ;
sod:GRIB_gridDefinitionDescription = "Lambert azimuthal equal area projection " ;
sod:GRIB_missingValue = 9999LL ;
sod:GRIB_name = "Soil depth" ;
sod:GRIB_shortName = "sod" ;
sod:GRIB_totalNumber = 51LL ;
sod:GRIB_units = "m" ;
sod:long_name = "Soil depth" ;
sod:units = "m" ;
sod:standard_name = "unknown" ;
sod:grid_mapping = "lambert_azimuthal_equal_area" ;
sod:coordinates = "soilLayer" ;
int lambert_azimuthal_equal_area ;
lambert_azimuthal_equal_area:semi_major_axis = 6378137. ;
lambert_azimuthal_equal_area:EPSG_code = "EPSG:3035" ;
lambert_azimuthal_equal_area:latitude_of_projection_origin = 52. ;
lambert_azimuthal_equal_area:inverse_flattening = 298.257223563 ;
lambert_azimuthal_equal_area:longitude_of_projection_origin = 10. ;
lambert_azimuthal_equal_area:proj4_params = "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs" ;
lambert_azimuthal_equal_area:false_northing = 3210000. ;
lambert_azimuthal_equal_area:grid_mapping_name = "lambert_azimuthal_equal_area" ;
lambert_azimuthal_equal_area:false_easting = 4321000. ;
lambert_azimuthal_equal_area:coordinates = "soilLayer" ;
byte land_binary_mask(y, x) ;
land_binary_mask:standard_name = "land_binary_mask" ;
land_binary_mask:long_name = "masking_area" ;
land_binary_mask:units = "1" ;
land_binary_mask:grid_mapping = "lambert_azimuthal_equal_area" ;
land_binary_mask:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
land_binary_mask:coordinates = "soilLayer" ;
// global attributes:
:GRIB_edition = 2LL ;
:GRIB_centre = "ecmf" ;
:GRIB_centreDescription = "European Centre for Medium-Range Weather Forecasts" ;
:GRIB_subCentre = 0LL ;
:Conventions = "CF-1.7" ;
:institution = "European Centre for Medium-Range Weather Forecasts" ;
:history = "2021-11-10T10:31 GRIB to CDM+CF via cfgrib-0.9.9.1/ecCodes-2.23.0 with {\"source\": \"/cache/tmp/afa118b9-8b07-482c-aed7-26976ba41bb7-adaptor.mars.external-1636540185.293234-5781-3-tmp.grib\", \"filter_by_keys\": {}, \"encode_cf\": [\"parameter\", \"time\", \"geography\", \"vertical\"]}" ;
}
dimensions:
x = 1000 ;
y = 950 ;
variables:
double x(x) ;
x:long_name = "x coordinate of projection" ;
x:units = "Meter" ;
x:axis = "X" ;
double y(y) ;
y:long_name = "y coordinate of projection" ;
y:units = "Meter" ;
y:axis = "Y" ;
float wiltingpoint(y, x) ;
wiltingpoint:standard_name = "thetar3" ;
wiltingpoint:long_name = "residual_moisture3" ;
wiltingpoint:units = "mm" ;
wiltingpoint:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
// global attributes:
:CDI = "Climate Data Interface version 1.8.2 (http://mpimet.mpg.de/cdi)" ;
:Conventions = "CF-1.6" ;
:history = "Tue Dec 11 15:47:51 2018: cdo setname,wiltingpoint test6.nc thmin3.nc\n",
"Tue Dec 11 15:47:49 2018: cdo div test5.nc SDmap3.nc test6.nc\n",
"Tue Dec 11 15:47:49 2018: cdo add test2.nc test4.nc test5.nc\n",
"Tue Dec 11 15:47:49 2018: cdo mul test1.nc fracnoforest.nc test2.nc\n",
"Tue Dec 11 15:47:48 2018: cdo mul thetar3.nc soildepth3_o.nc test1.nc\n",
"Created Tue Feb 03 18:06:55 2015" ;
:source = "Lisflood European maps - pb2015" ;
:conventions = "CF-1.6" ;
:source_software = "Python netCDF3_Classic" ;
:title = "Lisflood maps for European setting Januar 2015" ;
:keywords = "Lisflood, Europe" ;
:institition = "JRC H01" ;
:CDO = "Climate Data Operators version 1.8.2 (http://mpimet.mpg.de/cdo)" ;
}
dimensions:
x = 1000 ;
y = 950 ;
variables:
double x(x) ;
x:long_name = "x coordinate of projection" ;
x:units = "Meter" ;
x:axis = "X" ;
double y(y) ;
y:long_name = "y coordinate of projection" ;
y:units = "Meter" ;
y:axis = "Y" ;
float fieldcapacity(y, x) ;
fieldcapacity:standard_name = "thetas3" ;
fieldcapacity:long_name = "saturated_moisture2" ;
fieldcapacity:units = "mm" ;
fieldcapacity:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
// global attributes:
:CDI = "Climate Data Interface version 1.8.2 (http://mpimet.mpg.de/cdi)" ;
:Conventions = "CF-1.6" ;
:history = "Tue Dec 11 15:53:48 2018: cdo setname,fieldcapacity thetas3.nc thmax3.nc\n",
"Created Tue Feb 03 18:06:56 2015" ;
:source = "Lisflood European maps - pb2015" ;
:conventions = "CF-1.6" ;
:source_software = "Python netCDF3_Classic" ;
:title = "Lisflood maps for European setting Januar 2015" ;
:keywords = "Lisflood, Europe" ;
:institition = "JRC H01" ;
:CDO = "Climate Data Operators version 1.8.2 (http://mpimet.mpg.de/cdo)" ;
}
GloFAS Auxiliary Data
The auxiliary data for GloFAS are provided as netCDF files in WGS 84 (EPSG:4326), and are available to download below.
Two datasets are available to download:
The two datasets provided here (upstream area and elevation) are applicable for all GloFAS versions, and therefore can be used for forecasts of all available dates.
The upstream area file denotes the total upstream area for each river pixel, expressed in m2. The upstream area is defined as the catchment area for each river segment, in other words, the total area that contributes water to the river at the specific grid point. The upstream area always includes the area of the pixel.
Download the upstream area file here: glofas_upArea.nc (21.6MB)
dimensions:
longitude = 3600 ;
latitude = 1500 ;
variables:
double longitude(longitude) ;
longitude:standard_name = "longitude" ;
longitude:long_name = "longitude" ;
longitude:units = "degrees_east" ;
longitude:axis = "X" ;
double latitude(latitude) ;
latitude:standard_name = "latitude" ;
latitude:long_name = "latitude" ;
latitude:units = "degrees_north" ;
latitude:axis = "Y" ;
float uparea(latitude, longitude) ;
uparea:long_name = "Upstream area of each river pixel" ;
uparea:units = "m2" ;
uparea:_FillValue = -3.402823e+38f ;
uparea:missing_value = -3.402823e+38f ;
// global attributes:
:CDI = "Climate Data Interface version 1.9.10 (https://mpimet.mpg.de/cdi)" ;
:Conventions = "CF-1.5" ;
:GDAL_PCRASTER_VALUESCALE = "VS_SCALAR" ;
:GDAL = "GDAL 3.2.1, released 2020/12/29" ;
:NCO = "netCDF Operators version 4.9.7 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)" ;
:CDO = "Climate Data Operators version 1.9.10 (https://mpimet.mpg.de/cdo)" ;
}
The elevation (sometimes referred to as orography) file denotes the mean elevation for each pixel in the EFAS domain, expressed in m above sea level. It is defined for each grid point of the modelled domain.
Download the elevation file here: glofas_dem.nc (21.7MB)
dimensions:
lat = 1500 ;
lon = 3600 ;
variables:
double lat(lat) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:axis = "Y" ;
double lon(lon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
float elv(lat, lon) ;
elv:long_name = "GDAL Band Number 1" ;
elv:grid_mapping = "crs" ;
elv:_FillValue = -3.402823e+38f ;
elv:missing_value = -3.402823e+38f ;
// global attributes:
:CDI = "Climate Data Interface version 1.9.6 (http://mpimet.mpg.de/cdi)" ;
:Conventions = "CF-1.5" ;
:GDAL_AREA_OR_POINT = "Area" ;
:GDAL = "GDAL 3.0.4, released 2020/01/28" ;
:CDO = "Climate Data Operators version 1.9.6 (http://mpimet.mpg.de/cdo)" ;
:NCO = "4.4.6" ;
:history = "Wed Oct 14 10:18:06 2020: ncrename -v zelv,elv dem_new2.nc dem_new3.nc\nWed Oct 14 10:17:46 2020: ncks -C -O -x -v crs dem_new.nc dem_new2.nc\nWed Oct 14 10:17:35 2020: ncrename -v elv,zelv dem.nc dem_new.nc\nTue Sep 22 17:14:54 2020: ncrename -v Band1,elv dem.nc dem2.nc\nTue Sep 22 16:50:15 2020: cdo -invertlat elv_01deg_edit.nc elv_01deg_edit_invert.nc\nTue Sep 22 14:00:31 2020: GDAL CreateCopy( elv_01deg_edit.nc, ... )" ;
}