Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

For sub-daily data for the HRES (stream=oper/wave) the analyses (type=an) are available hourly. The short forecasts, run twice daily from 06 and 18 UTC, provide hourly output forecast steps from 0 to 18 hours. For the EDA, the sub-daily non-wave data (stream=enda) are available every 3 hours but the sub-daily wave data (stream=ewda) are available hourly.

Wave spectra

The ERA5 wave model uses wave spectra with 24 directions and 30 frequencies (see "2D wave spectra (single)", Table 7),

...

titleDecoding 2D wave spectra

Download from ERA5

ERA5 wave spectra data is not available from the CDS disks. However, it is available in MARS and can be accessed through the CDS API. For more information see Data organisation and how to download ERA5 and How to download ERA5 (Option B: Download ERA5 family data that is NOT listed in the CDS online catalogue - SLOW ACCESS.

Decoding 2D wave spectra in GRIB

To decode wave spectra in GRIB format we recommend ecCodes. Wave spectra are encoded in a specific way that other tools might not decode correctly.

In GRIB, the parameter is called 2d wave spectra (single) because in GRIB, the data are stored as a single global field per each spectral bin (a given frequency and direction), but in NetCDF, the fields are nicely recombined to produce a 2d matrix representing the discretized spectra at each grid point.

...

Also note that it is NOT the spectral density that is encoded but rather log10 of it, so to recover the spectral density, expressed in m^2 /(radian Hz), one has to take the power 10 (10^) of the NON missing decoded values. Missing data are for all land points, but also, as part of the GRIB compression, all small values below a certain threshold have been discarded and so those missing spectral values are essentially 0. m^2 /(gradient Hz).

Decoding 2D wave spectra in NetCDF

The NetCDF wave spectra file will have the dimensions longitude, latitude, direction, frequency and time.

However, the direction and frequency bins are simply given as 1 to 24 and 1 to 30, respectively.

The direction bins start at 7.5 degree and increase by 15 degrees until 352.5, with 90 degree being towards the east (Oceanographic convention).

The frequency bins are non-linearly spaced. The first bin is 0.03453 Hz and the following bins are: f(n) = f(n-1)*1.1; n=2,30. The data provided is the log10 of spectra density. To obtain the spectral density one has to take to the power 10 (10 ** data). This will give the units 2D wave spectra as m**2 s radian**-1 . Very small values are discarded and set as missing values. These are essentially 0 m**2 s radian**-1.

This recoding can be done with the Python xarray package, for example:

Code Block
languagepy
import xarray as xr
import numpy as np
da = xr.open_dataarray('2d_spectra_201601.nc')
da = da.assign_coords(direction=np.arange(7.5, 352.5 + 15, 15))
da = da.assign_coords(frequency=np.full(30, 0.03453) * (1.1 ** np.arange(0, 30)))
da = 10 ** da
da = da.fillna(0)
da.to_netcdf(path='2d_spectra_201601_recoded.nc')

Units of 2D wave spectra

Once decoded, the units of 2D wave spectra are m2 s radian-1

Instantaneous parameters

All the analysed parameters and many of the forecast parameters are described as "instantaneous". For more information on what this means, see Parameters valid at the specified time.

Mean rates and accumulations

The accumulations (over the processing period) in the short forecasts (from 06 and 18 UTC) of ERA5 are treated differently compared with those in ERA-Interim and operational data (where the accumulations are from the beginning of the forecast to the validity date/time). In the short forecasts of ERA5, the accumulations are since the previous post processing (archiving), so for:

  • reanalysis: accumulations are over the hour (the processing period) ending at the validity date/time
  • ensemble: accumulations are over the 3 hours (the processing period) ending at the validity date/time
  • Monthly means (of daily means, stream=moda/edmo): accumulations have been scaled to have an "effective" processing period of one day, see section Monthly means

Mean rate parameters in ERA5 (e.g. Table 4) provide similar information to accumulations (Table 3), except that they are expressed as temporal averages, over the same processing periods, instead of accumulations, and so have units of "per second".

  • Mean rate parameters are easier to deal with than accumulations because the units do not vary with the processing period.
  • The mean rate hydrological parameters (e.g. the "Mean total precipitation rate") have units of "kg m-2 s-1", which are equivalent to "mm s-1". If you multiply by 86400 seconds (24 hours) you will produce the commonly used units of "mm day-1".

Note that:

  • For the CDS time, or validity time, of 00 UTC, the mean rates and accumulations are over the hour (3 hours for the EDA) ending at 00 UTC i.e. the mean or accumulation is during the previous day.
  • Mean rates and accumulations are not available from the analyses.
  • Mean rates and accumulations at step=0 have values of zero because the length of the processing period is zero.

Minimum/maximum since the previous post processing

The short forecasts of ERA5 contain some surface and single level parameters that are the minimum or maximum value since the previous post processing (archiving), see Table 5 below. So, for:

  • reanalysis: the minimum or maximum values are in the hour (the processing period) ending at the validity date/time
  • ensemble: the minimum or maximum values are in the 3 hours (the processing period) ending at the validity date/time

Wave spectra

The ERA5 wave model uses wave spectra with 24 directions and 30 frequencies (see "2D wave spectra (single)", Table 7),

Expand
titleDecoding 2D wave spectra

Download from ERA5

ERA5 wave spectra data is not available from the CDS disks. However, it is available in MARS and can be accessed through the CDS API. For more information see Data organisation and how to download ERA5 and How to download ERA5 (Option B: Download ERA5 family data that is NOT listed in the CDS online catalogue - SLOW ACCESS.

Decoding 2D wave spectra in GRIB

To decode wave spectra in GRIB format we recommend ecCodes. Wave spectra are encoded in a specific way that other tools might not decode correctly.

In GRIB, the parameter is called 2d wave spectra (single) because in GRIB, the data are stored as a single global field per each spectral bin (a given frequency and direction), but in NetCDF, the fields are nicely recombined to produce a 2d matrix representing the discretized spectra at each grid point.

The wave spectra are encoded in GRIB using a local table specific to ECMWF. Because of this, the conversion of the meta data containing the information about the frequencies and the directions are not properly converted from GRIB to NetCDF format. So rather than having the actual values of the frequencies and directions, values show index numbers (1,1) : first frequency, first direction, (1,2) first frequency, second direction, etc ....

For ERA, because there are a total of 24 directions, the direction increment is 15 degrees with the first direction given by half the increment, namely 7.5 degree, where direction 0. means going towards the north and 90 towards the east (Oceanographic convention), or more precisely, this should be expressed in gradient since the spectra are in m^2 /(Hz radian)
The first frequency is 0.03453 Hz and the following ones are : f(n) = f(n-1)*1.1, n=2,30

Also note that it is NOT the spectral density that is encoded but rather log10 of it, so to recover the spectral density, expressed in m^2 /(radian Hz), one has to take the power 10 (10^) of the NON missing decoded values. Missing data are for all land points, but also, as part of the GRIB compression, all small values below a certain threshold have been discarded and so those missing spectral values are essentially 0. m^2 /(gradient Hz).

Decoding 2D wave spectra in NetCDF

The NetCDF wave spectra file will have the dimensions longitude, latitude, direction, frequency and time.

However, the direction and frequency bins are simply given as 1 to 24 and 1 to 30, respectively.

The direction bins start at 7.5 degree and increase by 15 degrees until 352.5, with 90 degree being towards the east (Oceanographic convention).

The frequency bins are non-linearly spaced. The first bin is 0.03453 Hz and the following bins are: f(n) = f(n-1)*1.1; n=2,30. The data provided is the log10 of spectra density. To obtain the spectral density one has to take to the power 10 (10 ** data). This will give the units 2D wave spectra as m**2 s radian**-1 . Very small values are discarded and set as missing values. These are essentially 0 m**2 s radian**-1.

This recoding can be done with the Python xarray package, for example:

Code Block
languagepy
import xarray as xr
import numpy as np
da = xr.open_dataarray('2d_spectra_201601.nc')
da = da.assign_coords(direction=np.arange(7.5, 352.5 + 15, 15))
da = da.assign_coords(frequency=np.full(30, 0.03453) * (1.1 ** np.arange(0, 30)))
da = 10 ** da
da = da.fillna(0)
da.to_netcdf(path='2d_spectra_201601_recoded.nc')

Units of 2D wave spectra

Once decoded, the units of 2D wave spectra are m2 s radian-1

Instantaneous parameters

All the analysed parameters and many of the forecast parameters are described as "instantaneous". For more information on what this means, see Parameters valid at the specified time.

Mean rates and accumulations

The accumulations (over the processing period) in the short forecasts (from 06 and 18 UTC) of ERA5 are treated differently compared with those in ERA-Interim and operational data (where the accumulations are from the beginning of the forecast to the validity date/time). In the short forecasts of ERA5, the accumulations are since the previous post processing (archiving), so for:

  • reanalysis: accumulations are over the hour (the processing period) ending at the validity date/time
  • ensemble: accumulations are over the 3 hours (the processing period) ending at the validity date/time
  • Monthly means (of daily means, stream=moda/edmo): accumulations have been scaled to have an "effective" processing period of one day, see section Monthly means

Mean rate parameters in ERA5 (e.g. Table 4) provide similar information to accumulations (Table 3), except that they are expressed as temporal averages, over the same processing periods, instead of accumulations, and so have units of "per second".

  • Mean rate parameters are easier to deal with than accumulations because the units do not vary with the processing period.
  • The mean rate hydrological parameters (e.g. the "Mean total precipitation rate") have units of "kg m-2 s-1", which are equivalent to "mm s-1". If you multiply by 86400 seconds (24 hours) you will produce the commonly used units of "mm day-1".

Note that:

  • For the CDS time, or validity time, of 00 UTC, the mean rates and accumulations are over the hour (3 hours for the EDA) ending at 00 UTC i.e. the mean or accumulation is during the previous day.
  • Mean rates and accumulations are not available from the analyses.
  • Mean rates and accumulations at step=0 have values of zero because the length of the processing period is zero.

Minimum/maximum since the previous post processing

The short forecasts of ERA5 contain some surface and single level parameters that are the minimum or maximum value since the previous post processing (archiving), see Table 5 below. So, for:

...

Monthly means

In addition to the sub-daily data, most parameters are also available as monthly means. For the surface and single level parameters, there are some exceptions which are listed in Table 8.

...