HEALPix grid cells

This example demonstrates how plotting HEALPix data with grid_cells will preserve the grid geometry.

import earthkit.data as ekd
import earthkit.plots as ekp

import numpy as np

data = ekd.from_source("sample", "healpix-h128-nested-2t.grib")

chart = ekp.Map(domain=["France", "Spain"])

chart.grid_cells(
    data,
    levels=np.arange(0, 20, 0.5),
    colors="Spectral_r",
    units="celsius",
)

chart.grid_points(data)

chart.title()
chart.legend()

chart.coastlines()

chart.show()