Ozone - 2x2 grid layout

import earthkit.data
import earthkit.maps

import cartopy.crs as ccrs
import numpy as np

data = earthkit.data.from_source(
    "url",
    "https://get.ecmwf.int/repository/test-data/metview/gallery/era5_ozone_1999.nc",
)

style = earthkit.maps.Style(
    levels = np.arange(5e-7, 4.1e-6, 2.5e-7),
    colors="BuPu",
)

chart = earthkit.maps.Chart(
    crs=ccrs.SouthPolarStereo(central_longitude=180),
    domain=[-180, 180, -90, -25],
    domain_crs=ccrs.PlateCarree(),
    rows=2, columns=2,
)

chart.plot(data, style=style)

chart.coastlines()
chart.gridlines(draw_labels=False)

chart.legend()

chart.title("ERA5 {variable_name!l} at 70$hPa$  - {time:%Y}", fontsize=14)
chart.subplot_titles("{month}")

chart.show()