Skip to content

Add HRRR virtual pressure_level skew-T notebook#37

Closed
aldenks wants to merge 3 commits into
mainfrom
claude/hrrr-skew-t-diagram-t85xbb
Closed

Add HRRR virtual pressure_level skew-T notebook#37
aldenks wants to merge 3 commits into
mainfrom
claude/hrrr-skew-t-diagram-t85xbb

Conversation

@aldenks

@aldenks aldenks commented Jul 11, 2026

Copy link
Copy Markdown
Member

Add a quickstart notebook that builds a skew-T log-P diagram from the NOAA
HRRR 48 hour forecast "virtual" dataset. It opens the dataset's
pressure_level group (39 isobaric levels) from the staging STAC catalog,
extracts a vertical profile at a lon/lat point via the Lambert Conformal
projection, and renders the sounding with MetPy (temperature/dew-point
curves, surface parcel path, CAPE/CIN, LCL, and wind barbs).

claude added 3 commits July 11, 2026 03:19
Add a quickstart notebook that builds a skew-T log-P diagram from the NOAA
HRRR 48 hour forecast "virtual" dataset. It opens the dataset's
pressure_level group (39 isobaric levels) from the staging STAC catalog,
extracts a vertical profile at a lon/lat point via the Lambert Conformal
projection, and renders the sounding with MetPy (temperature/dew-point
curves, surface parcel path, CAPE/CIN, LCL, and wind barbs).

Reading the virtual store requires the gribberish zarr codec (decodes the
virtual GRIB chunk references) and zarr >= 3.2 (built-in scale_offset
codec), plus metpy for the plot; these are added to pyproject.toml.

Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_016Gxiz1eWxK8cy4gt3LAm9n
dynamical-catalog 0.7.0 now declares gribberish~=1.5.0 and zarr~=3.2 as
runtime dependencies, so the virtual dataset's GRIB decoding and the
scale_offset codec are pulled in transitively. Bump the floor to >=0.7.0
and drop the explicit gribberish and zarr pins that are no longer needed.

Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_016Gxiz1eWxK8cy4gt3LAm9n
Use MetPy's .metpy accessor (quantify / unit_array) to pull units directly
from each variable's CF `units` attribute instead of attaching hPa/degC/m-s
by hand. The dataset already carries correct CF units on the pressure_level
group, so this removes the manual unit bookkeeping.

Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_016Gxiz1eWxK8cy4gt3LAm9n
@aldenks

aldenks commented Jul 12, 2026

Copy link
Copy Markdown
Member Author

going to wrap this into the official quickstart for this dataset.

@aldenks aldenks closed this Jul 12, 2026
@aldenks
aldenks deleted the claude/hrrr-skew-t-diagram-t85xbb branch July 12, 2026 14:02
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants