{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# 🏔 Plot model terrain" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "from herbie import Herbie\n", "from herbie import paint\n", "from herbie.toolbox import EasyMap, pc, ccrs\n", "\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## HRRR Model Terrain" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "✅ Found ┊ model=hrrr ┊ \u001b[3mproduct=sfc\u001b[0m ┊ \u001b[38;2;41;130;13m2022-Jan-01 00:00 UTC\u001b[92m F00\u001b[0m ┊ \u001b[38;2;255;153;0m\u001b[3mGRIB2 @ aws\u001b[0m ┊ \u001b[38;2;255;153;0m\u001b[3mIDX @ aws\u001b[0m\n" ] }, { "data": { "text/html": [ "
<xarray.Dataset> Size: 46MB\n",
"Dimensions: (y: 1059, x: 1799)\n",
"Coordinates:\n",
" time datetime64[ns] 8B 2022-01-01\n",
" step timedelta64[ns] 8B 00:00:00\n",
" surface float64 8B 0.0\n",
" latitude (y, x) float64 15MB 21.14 21.15 21.15 ... 47.85 47.84\n",
" longitude (y, x) float64 15MB 237.3 237.3 237.3 ... 299.0 299.1\n",
" valid_time datetime64[ns] 8B 2022-01-01\n",
" gribfile_projection object 8B None\n",
"Dimensions without coordinates: y, x\n",
"Data variables:\n",
" lsm (y, x) float32 8MB 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0\n",
" orog (y, x) float32 8MB -0.01786 -0.01786 ... -0.01786\n",
"Attributes:\n",
" GRIB_edition: 2\n",
" GRIB_centre: kwbc\n",
" GRIB_centreDescription: US National Weather Service - NCEP\n",
" GRIB_subCentre: 0\n",
" Conventions: CF-1.7\n",
" institution: US National Weather Service - NCEP\n",
" model: hrrr\n",
" product: sfc\n",
" description: High-Resolution Rapid Refresh - CONUS\n",
" remote_grib: https://noaa-hrrr-bdp-pds.s3.amazonaws.com/hrrr....\n",
" local_grib: /home/blaylock/data/hrrr/20220101/subset_04ef438...\n",
" search: (?:HGT|LAND):surface<xarray.Dataset> Size: 8MB\n",
"Dimensions: (latitude: 721, longitude: 1440)\n",
"Coordinates:\n",
" time datetime64[ns] 8B 2022-01-01\n",
" step timedelta64[ns] 8B 00:00:00\n",
" surface float64 8B 0.0\n",
" * latitude (latitude) float64 6kB 90.0 89.75 89.5 ... -89.75 -90.0\n",
" * longitude (longitude) float64 12kB 0.0 0.25 0.5 ... 359.5 359.8\n",
" valid_time datetime64[ns] 8B 2022-01-01\n",
" gribfile_projection object 8B None\n",
"Data variables:\n",
" lsm (latitude, longitude) float32 4MB 0.0 0.0 ... 1.0 1.0\n",
" orog (latitude, longitude) float32 4MB nan nan ... 2.775e+03\n",
"Attributes:\n",
" GRIB_edition: 2\n",
" GRIB_centre: kwbc\n",
" GRIB_centreDescription: US National Weather Service - NCEP\n",
" GRIB_subCentre: 0\n",
" Conventions: CF-1.7\n",
" institution: US National Weather Service - NCEP\n",
" model: gfs\n",
" product: pgrb2.0p25\n",
" description: Global Forecast System\n",
" remote_grib: https://noaa-gfs-bdp-pds.s3.amazonaws.com/gfs.20...\n",
" local_grib: /home/blaylock/data/gfs/20220101/subset_04ef4f92...\n",
" search: :(?:HGT|LAND):surface