RRFS 🅱#

BETA

The Rapid Refresh Forecast System (RRFS) is undergoing rapid development. There are some prototype products and forecasts available on the cloud on AWS. Some files have an index file, others do not.

[1]:

from herbie import Herbie
from toolbox import EasyMap, pc
from paint.standard2 import cm_tmp

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

[2]:

H = Herbie("2024-01-10", model="rrfs", fxx=1, member="control", domain="conus")

✅ Found ┊ model=rrfs ┊ product=prslev ┊ 2024-Jan-10 00:00 UTC F01 ┊ GRIB2 @ aws ┊ IDX @ aws

[3]:

H.help()

╭─Herbie────────────────────────────────
│ Help for model='rrfs'
│
│ Rapid Refresh Forecast System (RRFS) Ensemble
│ {'aws product description': 'https://registry.opendata.aws/noaa-rrfs/'}
│
│ Herbie(date, model='rrfs', ...)
│
│ fxx : int
│ product : {"prs", "nat", "testbed", "ififip"}
│ member : {"control", int}
│ domain : {"conus", "alaska", "hawaii", "puerto rico", None}
│
│ If product="natlev', then domain must be None
│
╰───────────────────────────────────────

[4]:

H.SOURCES

[4]:

{'aws': 'https://noaa-rrfs-pds.s3.amazonaws.com/rrfs_a/rrfs_a.20240110/00/control/rrfs.t00z.prslev.f001.conus_3km.grib2'}

[5]:

H.PRODUCTS

[5]:

{'prslev': '', 'natlev': '', 'testbed': '', 'ififip': ''}

[6]:

ds = H.xarray("TMP:2 m above ground")
ds

[7]:

ax = EasyMap("50m", crs=ds.herbie.crs, figsize=[10, 8]).STATES().ax
p = ax.pcolormesh(
    ds.longitude, ds.latitude, ds.t2m, transform=pc, **cm_tmp(units="K").cmap_kwargs
)
plt.colorbar(
    p, ax=ax, orientation="horizontal", pad=0.05, **cm_tmp(units="K").cbar_kwargs
)

ax.set_title(
    f"{ds.t2m.GRIB_name}\n{ds.time.dt.strftime('%H:%M UTC %d %b %Y').item()} F{int(ds.step.dt.total_seconds() / 60 / 60):02d}",
    loc="right",
)
ax.set_title(
    f"{H.model.upper()}: member={H.member} {H.product_description}", loc="left"
)

[7]:

Text(0.0, 1.0, 'RRFS: member=control ')

../../../_images/user_guide_tutorial_model_notebooks_rrfs_7_1.png

Other products#

Reading natlev data into xarray causes my kernel to crash

[8]:

H = Herbie("2024-01-10", model="rrfs", product="natlev", member="control", domain=None)
H

✅ Found ┊ model=rrfs ┊ product=natlev ┊ 2024-Jan-10 00:00 UTC F00 ┊ GRIB2 @ aws ┊ IDX @ aws

[8]:

▌▌Herbie RRFS model natlev product initialized 2024-Jan-10 00:00 UTC F00 ┊ source=aws

[9]:

H.inventory("TMP:5 hybrid")

[9]:

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	?	??	search_this
87	88	500393712	510769360.0	500393712-510769360	2024-01-10	2024-01-10	TMP	5 hybrid level	anl	NaN	NaN	:TMP:5 hybrid level:anl:nan

[10]:

H.download("TMP:5 hybrid")

[10]:

PosixPath('/home/blaylock/data/rrfs/20240110/control/subset_a3ef48ac__rrfs.t00z.natlev.f000.grib2')

[11]:

ds = H.xarray("TMP:5 hybrid")
ds

/home/blaylock/GITHUB/Herbie/herbie/core.py:1082: UserWarning: Will not remove GRIB file because it previously existed.
  warnings.warn("Will not remove GRIB file because it previously existed.")

[12]:

ax = EasyMap("50m", crs=ds.herbie.crs, figsize=[10, 8]).STATES().ax
p = ax.pcolormesh(
    ds.longitude, ds.latitude, ds.t, transform=pc, **cm_tmp(units="K").cmap_kwargs
)
plt.colorbar(
    p, ax=ax, orientation="horizontal", pad=0.05, **cm_tmp(units="K").cbar_kwargs
)

ax.set_title(
    f"{ds.t.GRIB_name}\n{ds.time.dt.strftime('%H:%M UTC %d %b %Y').item()} F{int(ds.step.dt.total_seconds() / 60 / 60):02d}",
    loc="right",
)
ax.set_title(
    f"{H.model.upper()}: member={H.member} {H.product_description}", loc="left"
)

The Kernel crashed while executing code in the the current cell or a previous cell. Please review the code in the cell(s) to identify a possible cause of the failure. Click <a href='https://aka.ms/vscodeJupyterKernelCrash'>here</a> for more info. View Jupyter <a href='command:jupyter.viewOutput'>log</a> for further details.

[ ]: