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).

[1]:

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

import matplotlib.pyplot as plt
from 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")

/home/blaylock/GITHUB/Carpenter_Workshop/paint/radar2.py:86: SyntaxWarning: invalid escape sequence '\m'
  self.units = "m s$\mathregular{^{-1}}$"
/home/blaylock/GITHUB/Carpenter_Workshop/paint/standard2.py:153: SyntaxWarning: invalid escape sequence '\d'
  self.units = f"$\degree${units}"
/home/blaylock/GITHUB/Carpenter_Workshop/paint/standard2.py:160: SyntaxWarning: invalid escape sequence '\d'
  self.units = f"$\degree${units}"
/home/blaylock/GITHUB/Carpenter_Workshop/paint/standard2.py:309: SyntaxWarning: invalid escape sequence '\d'
  self.units = f"$\degree${units}"
/home/blaylock/GITHUB/Carpenter_Workshop/paint/standard2.py:316: SyntaxWarning: invalid escape sequence '\d'
  self.units = f"$\degree${units}"
/home/blaylock/GITHUB/Carpenter_Workshop/toolbox/stock.py:126: SyntaxWarning: invalid escape sequence '\*'
  """
/home/blaylock/GITHUB/Carpenter_Workshop/toolbox/stock.py:879: SyntaxWarning: invalid escape sequence '\*'
  '''

[2]:

# Some examples

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

✅ Found ┊ model=hrdps ┊ product=continental/2.5km ┊ 2023-Dec-29 00:00 UTC F32 ┊ GRIB2 @ msc ┊ IDX @ None

[2]:

'https://dd.weather.gc.ca/model_hrdps/continental/2.5km/00/032/20231229T00Z_MSC_HRDPS_TMP_AGL-2m_RLatLon0.0225_PT032H.grib2'

[3]:

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

✅ Found ┊ model=hrdps ┊ product=continental/2.5km ┊ 2023-Dec-29 00:00 UTC F12 ┊ GRIB2 @ msc ┊ IDX @ None

[3]:

'https://dd.weather.gc.ca/model_hrdps/continental/2.5km/00/012/20231229T00Z_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 0x7f9e3bdaa480>

../../../_images/user_guide_tutorial_model_notebooks_hrdps_9_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 0x7f9e3b961490>

../../../_images/user_guide_tutorial_model_notebooks_hrdps_11_1.png

HRDPS North domain (experimental)#

[13]:

H = Herbie(
    recent.floor("12H"),  # only run every 00 anc 12 UTC
    model="hrdps_north",
    fxx=0,
    product="north/grib2",
    variable="TMP",
    level="TGL_2",
)
ds = H.xarray()
ds

✅ Found ┊ model=hrdps_north ┊ product=north/grib2 ┊ 2023-Dec-29 00:00 UTC F00 ┊ GRIB2 @ msc ┊ IDX @ None
👨🏻‍🏭 Created directory: [/home/blaylock/data/hrdps_north/20231229]

/home/blaylock/GITHUB/Herbie/herbie/core.py:1065: UserWarning: Will not remove GRIB file because Herbie will only remove subsetted files (not full files).
  warnings.warn(

[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/user_guide_tutorial_model_notebooks_hrdps_22_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)#