Precomputed: regridding HEALPix GRIB fieldlist

This example shows how to interpolate GRIB data defined on a HEALPix nested grid using the precomputed backend. We will also see how to inspect and plot the resulting data and how to convert it to xarray.

To make this notebook work earthkit-data and earthkit-plots have to be installed. The data will be represented as an earthkit-data GRIB FieldList.

Regridding

We perform the regridding with the regrid() method.

[1]:

from earthkit.regrid import regrid
from earthkit.data import from_source

# Get HEALPix nested GRIB data containing two fields.
ds = from_source("sample", "H8_nested_t2.grib2")

# the target grid is a global 5x5 degree regular latitude-longitude grid
out_grid = {"grid": [5,5]}

# perform interpolation for each field and add results
# to a new fieldlist stored in memory
r = regrid(ds, out_grid=out_grid, interpolation="linear", backend="precomputed")

d = r.data()
lat = d[0]
lon = d[1]
vals = d[2:]
lat.shape, lon.shape, vals.shape

[1]:

((37, 72), (37, 72), (2, 37, 72))

Please note that regridding with the precomputed backend only works between a predefined set of global grids. See the Precomputed weights for details.

Plotting the results

[2]:

import earthkit.plots as ekp

ekp.quickplot(r).show()

../_images/examples_precomp_healpix_fieldlist_7_0.png

Converting the results to xarray

[3]:

r.to_xarray()

[3]:

<xarray.Dataset> Size: 44kB
Dimensions:         (step_timedelta: 2, latitude: 37, longitude: 72)
Coordinates:
  * step_timedelta  (step_timedelta) timedelta64[ns] 16B 00:00:00 12:00:00
  * latitude        (latitude) float64 296B 90.0 85.0 80.0 ... -80.0 -85.0 -90.0
  * longitude       (longitude) float64 576B 0.0 5.0 10.0 ... 345.0 350.0 355.0
Data variables:
    2t              (step_timedelta, latitude, longitude) float64 43kB ...
Attributes:
    param:        2t
    paramId:      167
    class:        od
    stream:       oper
    levtype:      sfc
    type:         fc
    expver:       0001
    date:         20240323
    time:         1200
    domain:       g
    Conventions:  CF-1.8
    institution:  ECMWF

Writing the results to disk

Write the resulting fieldlist to disk:

[4]:

out_file = "_res_H8_nested_to_5x5.grib"
r.to_target("file", out_file)