...
Area cropping:
Local machine
Time series extraction:
We are going to extract EFAS reforecast's timeseries at locations defined by latitude and longitude coordinates from a tiny subset of the GRDC dataset.
| Info | ||
|---|---|---|
| ||
EFAS's x and y coordinates, when converted from GRIB to NetCDF, are not projected coordinates but matrix indexes (i, j), It is necessary to download the upstream area static file that contains the projected coordinates and replace it in EFAS. |
| Code Block | ||||
|---|---|---|---|---|
| ||||
import cdsapi
c = cdsapi.Client()
c.retrieve(
'efas-reforecast', | ||||
| Code Block | ||||
| ||||
import cdsapi c = cdsapi.Client() c.retrieve( 'efas-reforecast', { 'format': 'grib', 'product_type': 'ensemble_perturbed_reforecasts', 'variable': 'river_discharge_in_the_last_6_hours', 'model_levels': 'surface_level', 'hyear': '2007', 'hmonth': 'march', 'hday': [ '04', '07', ], 'leadtime_hour': [ [ '0', '12', '18', '06', '12', '18', ], }, '6', ], }, 'efas_reforecast.grib') 'efas_reforecast.grib') |
Time series extraction:
We are going to extract EFAS reforecast's timeseries at locations defined by latitude and longitude coordinates from a tiny subset of the GRDC dataset.
| Info | ||
|---|---|---|
| ||
EFAS's x and y coordinates, when converted from GRIB to NetCDF, are not projected coordinates but matrix indexes (i, j), It is necessary to download the upstream area static file that contains the projected coordinates and replace it in EFAS. |
Copy the content into an empty file named "GRDC.csv", the file should reside in the same folder of the efas_reforecast.grib file and the upstream area.
...