HEALPix interpolated

import earthkit.data
import earthkit.plots
import numpy as np

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

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

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

chart.grid_points(data)

chart.title()
chart.legend()

chart.coastlines()
chart.borders()

chart.show()