...
| Code Block | ||||||
|---|---|---|---|---|---|---|
| ||||||
import cdsapi
c = cdsapi.Client()
# Climatology river discharge
c.retrieve('efas-historical', {
"format": "netcdf",
"hday": ["15","16","17","18"],
"hmonth": "november",
"hyear": "2020",
"model_levels": "surface_level",
"system_version": "version_4_0",
"time": ["00:00","06:00","18:00"],
"variable": "river_discharge_in_the_last_6_hours"
},
'clim_2020111500.nc')
# Forecast river discharge
c.retrieve(
'efas-forecast',
{
'format': 'netcdf',
'originating_centre': 'ecmwf',
'product_type': 'ensemble_perturbed_forecasts',
'variable': 'river_discharge_in_the_last_6_hours',
'model_levels': 'surface_level',
'year': '2020',
'month': '11',
'day': '15',
'time': '00:00',
'leadtime_hour': [
'6', '12', '18',
'24', '30', '36',
'42', '48', '54',
'60', '66', '72',
],
},
'eue_2020111500.nc')
# Forecast soil moisture
c.retrieve(
'efas-forecast',
{
'format': 'netcdf',
'originating_centre': 'ecmwf',
'product_type': 'high_resolution_forecast',
'variable': [
'soil_depth', 'volumetric_soil_moisture',
],
'model_levels': 'soil_levels',
'year': '2019',
'month': '01',
'day': '30',
'time': '00:00',
'leadtime_hour': [
'6', '12', '18',
'24', '30', '36',
'42', '48', '54',
'60', '66', '72',
],
'soil_level': [
'1', '2', '3',
],
},
'eud_2019013000.nc')
# forecast snow depth water equivalent
c.retrieve(
'efas-forecast',
{
'format': 'netcdf',
'originating_centre': 'ecmwf',
'variable': 'snow_depth_water_equivalent',
'product_type': 'control_forecast',
'model_levels': 'surface_level',
'year': '2021',
'month': '01',
'day': '30',
'time': '00:00',
'leadtime_hour': [
'6', '12', '18',
'24', '30', '36',
'42', '48', '54',
'60', '66', '72',
'78', '84', '90',
'96', '102',
],
},
'esd_2021013000.nc')
|
Auxiliary data
From the CDS Static files link, download the following NetCDFs:
...
When the marsurl will be in production one can download the auxiliary data simply requesting them through a CDS API request:
Note: the soil depth is now part of the auxiliary data. Before it was downloaded as you would download a model output.
| Code Block | ||||||
|---|---|---|---|---|---|---|
| ||||||
import cdsapi
c = cdsapi.Client()
c.retrieve(
'efas-forecast,
{
'format': 'netcdf',
'variable': [
'field_capacity', 'wilting_point', 'soil_depth'
],
'soil_level': [
'1', '2', '3',
],
},
'auxiliary.zip') |
...
| Code Block | ||||||
|---|---|---|---|---|---|---|
| ||||||
import xarray as xr
ds = xr.open_dataset('esd_2021013000.nc')
ds.sd.plot(col="step",col_wrap=4, robust=True) |
GLOFAS
GloFAS Medium-range reforecast (download time series)
...