Tools for mirroring, manipulating (exporting, extracting) and serving AgERA5 data.
The agera5tools consist of a set of commandline scripts as well as the agera5tools python package
which can be used to
- set up a mirror for AgERA5 that can automatically build a local copy and keep it up to date with the latest AgERA5 data.
- Allow operations on the downloaded NetCDF files directly such as dumping, point extraction and clipping
- Serve AgERA5 data on web API through the HTTP protocol. By providing the latitude and longitude in the URL, agera5tools can retrieve the corresponding data and return it as JSON.
The agera5 commandline tools currently have 8 options. The first four are for setting up and managing the local AgERA5 database:
- init to generate a configuration file and initialize the set up.
- build to download the relevant AgERA5 data from Copernicus Climate Data Store (CDS) and build the local database.
- mirror to update the current database with new days from the CDS.
- serve to serve AgERA5 data through a web API and return as JSON encoded data.
The other four tools operate directly on the NetCDF files downloaded from the CDS.
- extract_point: this can be used to extract a time-series of variables for a given location
- dump which can be used to dump one day of AgERA5 data to CSV, JSON or SQLite
- clip which can be used to extract a subset of AgERA5 data which will be written to a new NetCDF file.
- dump_grid which dumps the AgERA5 grid definition to CSV, JSON or SQLite.
$ agera5tools init --help
using config from /data/agera5/agera5tools.yaml
Usage: agera5tools init [OPTIONS]
Initializes AgERA5tools
This implies the following:
- Creating a template configuration file in the current directory
- Creating a $HOME/.cdsapirc file for access to the CDS
- Creating the database tables
- Filling the grid table with the reference grid.
Options:
--help Show this message and exit.$ agera5tools build --help
using config from /data/agera5/agera5tools.yaml
Usage: agera5tools build [OPTIONS]
Builds the AgERA5 database by bulk download from CDS
Options:
-d, --to_database Load AgERA5 data into the database
-c, --to_csv Write AgERA5 data to compressed CSV files.
--help Show this message and exit.$ agera5tools mirror --help
using config from /data/agera5/agera5tools.yaml
Usage: agera5tools mirror [OPTIONS]
Incrementally updates the AgERA5 database by daily downloads from the CDS.
Options:
-c, --to_csv Write AgERA5 data to compressed CSV files.
--help Show this message and exit.$ agera5tools serve --help
using config from /data/agera5/agera5tools.yaml
Usage: agera5tools serve [OPTIONS]
Starts the http server to serve AgERA5 data through HTTP
Options:
-p, --port INTEGER Port to number to start listening, default=8080.
--help Show this message and exit.$ agera5tools extract_point --help
Usage: agera5tools extract_point [OPTIONS] AGERA5_PATH LONGITUDE LATITUDE
STARTDATE ENDDATE
Extracts AgERA5 data for given location and date range.
AGERA5_PATH: path to the AgERA5 dataset
LONGITUDE: the longitude for which to extract [dd, -180:180]
LATITUDE: the latitude for which to extract [dd, -90:90]
STARTDATE: the start date (yyyy-mm-dd, >=1979-01-01)
ENDDATE: the last date (yyyy-mm-dd, <= 1 week ago)
Options:
-o, --output PATH output file to write to: .csv, .json and .db3 (SQLite)
are supported.Giving no output will write to stdout in
CSV format
--tocelsius Convert temperature values from degrees Kelvin to Celsius
--help Show this message and exit.$ agera5tools dump --help
Usage: agera5tools dump [OPTIONS] AGERA5_PATH DAY
Dump AgERA5 data for a given day to CSV, JSON or SQLite
AGERA5_PATH: Path to the AgERA5 dataset
DAY: specifies the day to be dumped (yyyy-mm-dd)
Options:
-o, --output PATH output file to write to: .csv, .json and .db3 (SQLite)
are supported. Giving no output will write to stdout in
CSV format
--tocelsius Convert temperature values from degrees Kelvin to Celsius
--add_gridid Adds a grid ID instead of latitude/longitude columns.
--bbox FLOAT... Bounding box: <lon_min> <lon_max> <lat_min< <lat max>
--help Show this message and exit.$ agera5tools clip --help
Usage: agera5tools clip [OPTIONS] AGERA5_PATH DAY
Extracts a portion of agERA5 for the given bounding box and saves to
NetCDF.
AGERA5_PATH: Path to the AgERA5 dataset
DAY: specifies the day to be dumped (yyyy-mm-dd)
Options:
--base_fname TEXT Base file name to use, otherwise will use
'agera5_clipped'
-o, --output_dir PATH Directory to write output to. If not provided, will
use current directory.
--box FLOAT... Bounding box: <lon_min> <lon_max> <lat_min< <lat max>
--help Show this message and exit.Usage: agera5tools dump_grid [OPTIONS]
Dump the agERA5 grid to a CSV, JSON or SQLite DB.
Options:
-o, --output PATH output file to write to: .csv, .json and .db3 (SQLite)
are supported.Giving no output will write to stdout in
CSV format
--help Show this message and exit.
The shell commands described above can also be used from python directly by importing the agera5tools package.
Their working is nearly identical as the shell commands. The major difference is that the python functions
return either datasets (clip) or dataframes (extract_point, dump, dump_grid). An example for the clip function:
In [1]: import datetime as dt
...: import agera5tools
...: from agera5tools.util import BoundingBox
...: day = dt.date(2018,1,1)
...: bbox = BoundingBox(lon_min=87, lon_max=90, lat_min=24, lat_max=27)
...: ds = agera5tools.clip(day, bbox)
...:
In [2]: ds
Out[2]:
<xarray.Dataset>
Dimensions: (time: 1, lon: 30, lat: 30)
Coordinates:
* time (time) datetime64[ns] 2018-01-01
* lon (lon) float64 87.1 87.2 ... 89.9 90.0
* lat (lat) float64 26.9 26.8 ... 24.1 24.0
Data variables:
Precipitation_Flux (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Solar_Radiation_Flux (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Temperature_Air_2m_Max_Day_Time (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Temperature_Air_2m_Mean_24h (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Temperature_Air_2m_Min_Night_Time (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Vapour_Pressure_Mean (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Wind_Speed_10m_Mean (time, lat, lon) float32 dask.array<chunksize=(1, 30, 30), meta=np.ndarray>
Attributes:
CDI: Climate Data Interface version 1.9.2 (http://mpimet.mpg.de/...
history: Fri Mar 12 15:04:43 2021: cdo splitday /archive/ESG/wit015/...
Conventions: CF-1.7
CDO: Climate Data Operators version 1.9.2 (http://mpimet.mpg.de/...It works in a very similar way for the extract_point function:
In[6]: from agera5tools.util import Point
In[7]: pnt = Point(latitude=26, longitude=89)
In[8]: df = agera5tools.extract_point(pnt, startday=dt.date(2018, 1, 1), endday=dt.date(2018, 1, 31)),
In [7]: df.head(5)
Out[7]:
day precipitation_flux solar_radiation_flux ... temperature_air_2m_min_night_time vapour_pressure_mean wind_speed_10m_mean
0 2018-01-01 0.31 13282992 ... 12.156799 11.809731 1.317589
1 2018-01-02 1.91 13646220 ... 12.342041 11.711860 1.416075
2 2018-01-03 0.14 14817991 ... 11.064514 11.198871 1.524268
3 2018-01-04 0.03 14131904 ... 10.861877 11.413278 1.566405
4 2018-01-05 0.07 14315206 ... 12.292969 10.984181 1.597181
[5 rows x 8 columns]Note that extracting point data for a long timeseries can be time-consuming because all netCDF files have to be opened, decompressed and the point extracted.
The agera5tools package requires python >=3.8 and has a number of dependencies:
- pandas == 1.4.1
- PyYAML >= 6.0
- Pandas >= 1.5
- SQLAlchemy >= 1.4
- PyYAML >= 6.0
- xarray >= 2022.12.0
- dask >= 2022.7.0
- click >= 8.1
- flask >= 2.2
- cdsapi >= 0.5.1
- dotmap >= 1.3
- netCDF4 >= 1.6
- requests >= 2.28
- wsgiserver >= 1.3
Lower versions of dependencies may work, but have not been tested.
Installing agera5tools can be done through the github repository to get the latest version:
pip install https://github.com/ajwdewit/agera5tools/archive/refs/heads/main.zipor directory from PyPI:
pip install agera5tools