-
Notifications
You must be signed in to change notification settings - Fork 0
Download all atmospheric boundary conditions from ERA5 #4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: main
Are you sure you want to change the base?
Conversation
- temperature - wind components - humidity (to be checked)
- no remapping of meteocloud file - WIND computation from ERA5_instant file - temperature, humidity from ERA5 instant file - variable renaming adapted to ERA5 variable names
s-poll
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is there a possibility to download lowermost atmospheric model layer of ERA5 or everything at 10m height to download wind, temperature and specific humidity at the same reference height?
mkforcing/download_ERA5_input.py
Outdated
| "mean_total_precipitation_rate", | ||
| "10m_u_component_of_wind", | ||
| "10m_v_component_of_wind", | ||
| "2m_temperature", |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Wind, temperature, and specific humidify should be on the same height, as the reference height (ZBOT) is set for all variables.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I added a 2m->10m conversion script (with very simplified "conversions") at least as a starting point for more elaborate discussions of this issue!
mkforcing/download_ERA5_input.py
Outdated
| "10m_u_component_of_wind", | ||
| "10m_v_component_of_wind", | ||
| "2m_temperature", | ||
| "total_column_water_vapour", |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This is the water of the total atmospheric column. We need to download the specific humidity.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Now, dewpoint temperature at 2m is downloaded and converted to specific humidity at 2m by Python script dewpoint_to_specific_humidity.py.
`download_ERA5_input.py`: When the file target is not supplied as input variable `target`, the downloaded file has up to now always been called `*.zip`. This default behavior is replaced with the following, more robust alternative: 1. Download the data to a temporary target without extension/suffix. 2. Detect the file type (supports zip, netCDF and GRIB). 3. Add the extension of the detected file type.
both dimension and variable
Renaming both dimension and variable lead to erroneous values for unclear reasons.
Make possible downloading all atmospheric boundary conditions from ERA5: - https://cds.climate.copernicus.eu/datasets/reanalysis-era5-single-levels?tab=download
Implementations
download_ERA5_input.py: More flexible input next to default Copernicus API request. Custom requests possible. (This is also a preparation for #12 )download_ERA5_input.py; Additional variables to be added in custom request for ERA5 data:t2m)u10,v10)d2m)Added two Python scripts
dewpoint_to_specific_humidity.py: Compute specific humidity at 2m from dewpoint temperature at 2m.2m_to_10m_conversion.py: Convert temperature and specific humidity from 2m to 10m.Adapted
prepare_ERA5_input.shto be backward-compatible, but also able to use all-Copernicus-ERA5.Related issue: #1