DMI - Radar reflectivity over Denmark

This notebook describes how to access and use the zarr-version of the DMI radar composite dataset. The dataset includes ten-minute reflectivity over Denmark, which is obtained by the five radars owned by the Danish Meteorologiske Institute (DMI). The v0.1.1 zarr version of thre dataset covers 2016 to the end of 2024.

The conversion of the radar composite was done by the code hosted on https://github.com/mlcast-community/mlcast-dataset-DMI/tree/main

import matplotlib.pyplot as plt
import cartopy.crs as ccrs

import mlcast_datasets

cat = mlcast_datasets.open_catalog()
list(cat.precipitation)

['radklim_hourly',
 'radklim_5_minutes',
 'dmi_10_minutes',
 'it_dpc_sri_5min',
 'uk_metoffice_5min',
 'be_rmi_radclim_mfb_5min']

DMI 10-min reflectivity¶

Mlcast catalog contains multiple datasets, hence to call specifically for DMI radar data:

ds = cat.precipitation.dmi_10_minutes.to_dask()
ds

var_name = "dbz"

crs_name = ds[var_name].grid_mapping
data_crs = ccrs.Projection(ds[crs_name].crs_wkt)

g = (
    ds[var_name]
    .sel(time="2024-12-31")
    .isel(time=slice(None, 2))
    .plot(
        transform=data_crs,
        cmap="viridis",
        add_colorbar=True,
        col="time",
        robust=True,
        subplot_kws=dict(projection=data_crs),
    )
)

for ax in g.axes.flat:
    ax.coastlines()
    ax.gridlines(draw_labels=["top", "left"])

/tmp/ipykernel_2515/1860960605.py:15: FutureWarning: self.axes is deprecated since 2022.11 in order to align with matplotlibs plt.subplots, use self.axs instead.
  for ax in g.axes.flat: