HRDPS 🅱#

BETA Requires MetPy >=1.6

This demonstrates using data from Canada’s High Resolution Deterministic Prediction System (HRDPS).

Data Sources

`prioriy=`	Data source	Archive Duration
`"mcs"`	Meteorological Service of Canada	Last 24 hours

Model Initialization

Model cyles every hour.

Forecast Hour

For the most recent version of HRRR…

`fxx=`	Forecast lead time
`0` through `48`, step=1	hourly forecasts available

Products

`product=`	Product Description
`"continental/2.5km"`	Continental domain

Variable and Level

You will need to specify the variable and level for each request.

NOTE: The organization of these files is different than other NWP products.

There are no index files provided.

Each GRIB2 file only contains one message. The variable name and level is in the file’s name.

Herbie requires you provide a keyword argument for both variable and level. Pay special attention to model description (linked above) to understand how the model data is organized. If you don’t provide input for variable or level, Herbie will give you some ideas. For example, variable=TMP and level=AGL-2m will give you the filename that contains
TMP_AGL-2m

Note: This requires MetPy version 1.6 or greater which has the capability to parse the rotated latitude longitude map projection type (see MetPy/#3123).

[ ]:

from herbie import Herbie
import xarray as xr
import numpy as np

import matplotlib.pyplot as plt
from herbie.toolbox import EasyMap, pc
import cartopy.crs as ccrs
import cartopy.feature as feature
import pandas as pd

recent = pd.Timestamp("now").floor("6h") - pd.Timedelta("6h")

[2]:

# Some examples

H = Herbie(
    recent,  # Datetime
    model="hrdps",
    fxx=32,
    product="continental",
    variable="TMP",
    level="AGL-2m",
)
H.grib

✅ Found ┊ model=hrdps ┊ product=continental ┊ 2026-Apr-15 06:00 UTC F32 ┊ GRIB2 @ msc ┊ IDX @ None

[2]:

'https://dd.weather.gc.ca/20260415/WXO-DD/model_hrdps/continental/2.5km/06/032/20260415T06Z_MSC_HRDPS_TMP_AGL-2m_RLatLon0.0225_PT032H.grib2'

[3]:

H = Herbie(
    recent,  # Datetime
    model="hrdps",
    fxx=12,
    product="continental",
    variable="HGT",
    level="ISBL_0500",
)
H.grib

✅ Found ┊ model=hrdps ┊ product=continental ┊ 2026-Apr-15 06:00 UTC F12 ┊ GRIB2 @ msc ┊ IDX @ None

[3]:

'https://dd.weather.gc.ca/20260415/WXO-DD/model_hrdps/continental/2.5km/06/012/20260415T06Z_MSC_HRDPS_HGT_ISBL_0500_RLatLon0.0225_PT012H.grib2'

Plot data on Plate Carree projection#

[6]:

ax = EasyMap("50m").BORDERS().STATES(alpha=0.5).ax
p = ax.pcolormesh(ds.longitude, ds.latitude, ds.t2m, transform=pc, cmap="Spectral_r")
plt.colorbar(p, ax=ax, orientation="horizontal", pad=0.01, shrink=0.8)
ax.set_title(
    f"2-m Temperature\nValid {ds.valid_time.dt.strftime('%Y-%m-%d %H:%M').item()} UTC",
    loc="right",
    fontsize=10,
)
ax.set_title(f"{ds.model.upper()}: {H.product_description}", loc="left")
ax.gridlines()

[6]:

<cartopy.mpl.gridliner.Gridliner at 0x783b9b0c4980>

../../_images/gallery_eccc_models_hrdps_8_1.png

Plot data on model grid projection - rotated latitude longitude#

Requires Metpy >=1.6.0

[7]:

ax = EasyMap("50m", crs=ds.herbie.crs).BORDERS().STATES(alpha=0.5).ax
p = ax.pcolormesh(ds.longitude, ds.latitude, ds.t2m, transform=pc, cmap="Spectral_r")
plt.colorbar(p, ax=ax, orientation="horizontal", pad=0.01, shrink=0.8)
ax.set_title(
    f"2-m Temperature\nValid {ds.valid_time.dt.strftime('%Y-%m-%d %H:%M').item()} UTC",
    loc="right",
    fontsize=10,
)
ax.set_title(f"{ds.model.upper()}: {H.product_description}", loc="left")
ax.gridlines()

[7]:

<cartopy.mpl.gridliner.Gridliner at 0x783b9b7363c0>

../../_images/gallery_eccc_models_hrdps_10_1.png

HRDPS North domain (experimental)#

Not available as of testing in Apr 2026

[8]:

H = Herbie(
    recent.floor("12h"),  # only run every 00 and 12 UTC
    model="hrdps",
    fxx=0,
    product="north",
    variable="TMP",
    level="AGL-2m",
)
ds = H.xarray()
ds

/home/kps5442/mapwall_dev/herbie-dev/src/herbie/core.py:1306: UserWarning: Will not remove GRIB file because Herbie will only remove subsetted files (not full files).
  warnings.warn(
🦨 GRIB2 file not found: self.model='hrdps' self.date=Timestamp('2026-04-15 00:00:00') self.fxx=0

💔 Did not find ┊ model=hrdps ┊ product=north ┊ 2026-Apr-15 00:00 UTC F00

---------------------------------------------------------------------------
FileNotFoundError                         Traceback (most recent call last)
Cell In[8], line 9
      5     product="north",
      6     variable="TMP",
      7     level="AGL-2m",
      8 )
----> 9 ds = H.xarray()
     10 ds

File ~/mapwall_dev/herbie-dev/src/herbie/core.py:1338, in Herbie.xarray(self, search, searchString, backend_kwargs, remove_grib, _use_pygrib_for_crs, **download_kwargs)
   1334 backend_kwargs.setdefault("errors", "raise")
   1336 # Use cfgrib.open_datasets, just in case there are multiple "hypercubes"
   1337 # for what we requested.
-> 1338 Hxr = cfgrib.open_datasets(
   1339     local_file,
   1340     backend_kwargs=backend_kwargs,
   1341     decode_timedelta=True,
   1342 )
   1344 for ds in Hxr:
   1345     # Need model attribute before using get_cf_crs
   1346     ds.attrs["model"] = str(self.model)

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/xarray_store.py:119, in open_datasets(path, backend_kwargs, **kwargs)
    117 backend_kwargs = backend_kwargs.copy()
    118 backend_kwargs["squeeze"] = False
--> 119 datasets = open_variable_datasets(path, backend_kwargs=backend_kwargs, **kwargs)
    121 type_of_level_datasets = {}  # type: T.Dict[str, T.List[xr.Dataset]]
    122 for ds in datasets:

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/xarray_store.py:101, in open_variable_datasets(path, backend_kwargs, **kwargs)
     99 errors = backend_kwargs.get("errors", "warn")
    100 stream = messages.FileStream(path, errors=errors)
--> 101 index = open_fileindex(stream, computed_keys=cfmessage.COMPUTED_KEYS, **fileindex_kwargs)
    102 datasets = []  # type: T.List[xr.Dataset]
    103 for param_id in sorted(index["paramId"]):

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/dataset.py:807, in open_fileindex(stream, indexpath, index_keys, ignore_keys, filter_by_keys, computed_keys)
    805 index_keys = sorted(set(index_keys) | set(filter_by_keys))
    806 index_keys = [key for key in index_keys if key not in ignore_keys]
--> 807 index = messages.FileIndex.from_indexpath_or_filestream(
    808     stream, index_keys, indexpath=indexpath, computed_keys=computed_keys
    809 )
    810 return index.subindex(filter_by_keys)

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/messages.py:533, in FileIndex.from_indexpath_or_filestream(cls, filestream, index_keys, indexpath, computed_keys, log)
    525 @classmethod
    526 def from_indexpath_or_filestream(
    527     cls, filestream, index_keys, indexpath=DEFAULT_INDEXPATH, computed_keys={}, log=LOG
   (...)    530
    531     # Reading and writing the index can be explicitly suppressed by passing indexpath==''.
    532     if not indexpath:
--> 533         return cls.from_fieldset(filestream, index_keys, computed_keys)
    535     hash = hashlib.md5(repr(index_keys).encode("utf-8")).hexdigest()
    536     indexpath = indexpath.format(path=filestream.path, hash=hash, short_hash=hash[:5])

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/messages.py:379, in FieldsetIndex.from_fieldset(cls, fieldset, index_keys, computed_keys)
    377 else:
    378     iteritems = enumerate(fieldset)
--> 379 return cls.from_fieldset_and_iteritems(fieldset, iteritems, index_keys, computed_keys)

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/messages.py:392, in FieldsetIndex.from_fieldset_and_iteritems(cls, fieldset, iteritems, index_keys, computed_keys)
    390 index_keys = list(index_keys)
    391 header_values_cache = {}  # type: T.Dict[T.Tuple[T.Any, type], T.Any]
--> 392 for field_id, raw_field in iteritems:
    393     field = ComputedKeysAdapter(raw_field, computed_keys)
    394     header_values = []

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/messages.py:292, in FileStreamItems.__iter__(self)
    290 old_offset = -1
    291 count = 0
--> 292 for message in self.itervalues():
    293     offset = message.message_get("offset", int)
    294     if offset == old_offset:

File /home/meteo/kps5442/miniconda3/envs/herbie-dev-local/lib/python3.14/site-packages/cfgrib/messages.py:268, in FileStreamItems.itervalues(self)
    266 def itervalues(self) -> T.Iterator[Message]:
    267     errors = self.filestream.errors
--> 268     with open(self.filestream.path, "rb") as file:
    269         # enable MULTI-FIELD support on sequential reads (like when building the index)
    270         with multi_enabled(file):
    271             valid_message_found = False

FileNotFoundError: [Errno 2] No such file or directory: '/home/meteo/kps5442/data/hrdps/20260415/20260415T00Z_MSC_HRDPS-North_TMP_AGL-2m_RLatLon0.03_PT000H.grib2'

[14]:

# This domain is run on a polar stereographic projection
ds.herbie.crs

[14]:

<cartopy.crs.Stereographic object at 0x7f9e3bbca060>

[15]:

ax = EasyMap("50m", crs=ds.herbie.crs).BORDERS().STATES(alpha=0.5).ax
p = ax.pcolormesh(
    ds.longitude, ds.latitude, ds.t2m - 273.15, transform=pc, cmap="Spectral_r"
)
plt.colorbar(
    p,
    ax=ax,
    orientation="horizontal",
    pad=0.01,
    shrink=0.8,
    label="2-m Temperature (C)",
)
ax.set_title(f"{ds.model.upper()}", loc="left")
ax.set_title(
    f"2-m Temperature\nValid {ds.valid_time.dt.strftime('%Y-%m-%d %H:%M').item()} UTC",
    loc="center",
    fontsize=10,
)
ax.gridlines()
ax.EasyMap.INSET_GLOBE()

[15]:

<GeoAxes: >

../../_images/gallery_eccc_models_hrdps_21_1.png

[ ]:

HRDPS 🅱#

Get the 2-metre temperature#

Plot data on Plate Carree projection#

Plot data on model grid projection - rotated latitude longitude#

Get 10-m U and 10-m V wind#

500 hPa Humidity and Geopotential Height#

HRDPS North domain (experimental)#