WEMC TECH BLOG #3: ERA5 10m/100m Analysis
We wanted to see if we could improve on the standard power law coefficient (1.389495) for calculating 100m wind speed from 10m data. The currently available ten years of ERA5 U and V netCDF wind components from CDS were concatenated, then calculated for wind speed using CDO; a handy collection of command-line operators to manipulate and analyse climate and numerical weather prediction data:
The coefficient dataset was calculated with CDO, by dividing 100m by 10m datasets:
The mean coefficient for the whole domain is 1.448035, slightly higher than the standard coefficient, taking it further away from the newly available ERA5 100m forecast data. The average standard deviation is 0.816657.
The ratio as a function of the geographical location was then plotted in Python, using matplotlib package’s plot.imshow() method.
First load the necessary python packages:
Import the pre-calculated coefficient dataset:
Calculate the mean average using the numpy mean method:
Finally, plot the total wind speed average using matplotlib imshow() method, note interpolation is set to ‘none’, this can be changed to apply various smoothing algorithms:
Next we can check the standard deviation over this time period:
For this plot, the colour map (cmap=’terrain’) was adjusted to highlight the variation:
This indicates that a geographically-dependent coefficient would better reproduce the 100 m wind, when only 10 m is available.
by Luke Sanger (WEMC, 2018)