ERA5 atmospheric stability for usage in WAsP 12.8

The dataset "ERA5-meso.nc" is obtained by loading the variables needed for calculating the temperature scale from hourly ERA5 files. These files are available in grib format that have been obtained from the Copernicus Data Store (CDS). They are opened using xarray and cfgrib and processed using the functions stability_histogram from the python package PyWAsP. Because a conditional mean based on the 50% highest wind speeds must be calculated, all values are binned according to wind speed at 100 m and this histogram is then used to calculate the mean and root-mean-square of the temperature scale. The boundary layer height scale is calculated in a similar fashion. For more documentation see the accompanying paper. A validation of the WAsP model using these data is available in the references.

The file "ERA5-baro.nc" contains the geostrophic wind shear. The mean magnitude and direction is obtained sector-wise in similar fashion as described above. The geostrophic wind shear can be calculated from the pressure level geopotential height. The way to do this is described here:

https://orbit.dtu.dk/en/publications/implementation-of-large-scale-average-geostrophic-wind-shear-in-w

Please, also refer to the item "ERA5 data for air density calculations in WAsP" https://doi.org/10.11583/DTU.13286294

v1: Final version that corresponds to the published version of the paper:
https://doi.org/10.1007/s10546-023-00803-3
This is slightly different then the first version due to Eq. 9