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. Code Block |
---|
language | bash |
---|
title | Soil depth netCDF metadata (for soil level 3, and v4.0) |
---|
collapse | true |
---|
| 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\"]}" ;
} |
Code Block |
---|
language | bash |
---|
title | Wilting point netCDF metadata (for soil level 3, and v4.0) |
---|
| 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)" ;
} |
Code Block |
---|
language | bash |
---|
title | Field capacity netCDF metadata (for soil level 3, and v4.0) |
---|
collapse | true |
---|
| 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)" ;
} |
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.
|