HRRR

This tutorial shows you how to download HRRR (High Resolution Rapid Refresh) data from Amazon (and other sources) with the Herbie Python package.

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

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

H = Herbie(
    "2021-07-19",
    model="hrrr",
    product="sfc",
    fxx=0,
)

✅ Found ┊ model=hrrr ┊ product=sfc ┊ 2021-Jul-19 00:00 UTC F00 ┊ GRIB2 @ aws ┊ IDX @ aws

# Show available products
H.PRODUCTS

{'sfc': '2D surface level fields; 3-km resolution',
 'prs': '3D pressure level fields; 3-km resolution',
 'nat': 'Native level fields; 3-km resolution',
 'subh': 'Subhourly grids; 3-km resolution'}

# Show actual URL sources Herbie searches
H.SOURCES

{'aws': 'https://noaa-hrrr-bdp-pds.s3.amazonaws.com/hrrr.20210719/conus/hrrr.t00z.wrfsfcf00.grib2',
 'nomads': 'https://nomads.ncep.noaa.gov/pub/data/nccf/com/hrrr/prod/hrrr.20210719/conus/hrrr.t00z.wrfsfcf00.grib2',
 'google': 'https://storage.googleapis.com/high-resolution-rapid-refresh/hrrr.20210719/conus/hrrr.t00z.wrfsfcf00.grib2',
 'azure': 'https://noaahrrr.blob.core.windows.net/hrrr/hrrr.20210719/conus/hrrr.t00z.wrfsfcf00.grib2',
 'pando': 'https://pando-rgw01.chpc.utah.edu/hrrr/sfc/20210719/hrrr.t00z.wrfsfcf00.grib2',
 'pando2': 'https://pando-rgw02.chpc.utah.edu/hrrr/sfc/20210719/hrrr.t00z.wrfsfcf00.grib2'}

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

Downloading inventory file from self.idx='https://noaa-hrrr-bdp-pds.s3.amazonaws.com/hrrr.20210719/conus/hrrr.t00z.wrfsfcf00.grib2.idx'

<xarray.Dataset> Size: 38MB
Dimensions:              (y: 1059, x: 1799)
Coordinates:
    time                 datetime64[ns] 8B 2021-07-19
    step                 timedelta64[ns] 8B 00:00:00
    heightAboveGround    float64 8B 2.0
    latitude             (y, x) float64 15MB 21.14 21.15 21.15 ... 47.85 47.84
    longitude            (y, x) float64 15MB 237.3 237.3 237.3 ... 299.0 299.1
    valid_time           datetime64[ns] 8B 2021-07-19
    gribfile_projection  object 8B None
Dimensions without coordinates: y, x
Data variables:
    t2m                  (y, x) float32 8MB 295.3 295.3 295.3 ... 287.7 287.6
Attributes:
    GRIB_edition:            2
    GRIB_centre:             kwbc
    GRIB_centreDescription:  US National Weather Service - NCEP
    GRIB_subCentre:          0
    Conventions:             CF-1.7
    institution:             US National Weather Service - NCEP
    model:                   hrrr
    product:                 sfc
    description:             High-Resolution Rapid Refresh - CONUS
    remote_grib:             https://noaa-hrrr-bdp-pds.s3.amazonaws.com/hrrr....
    local_grib:              /home/blaylock/data/hrrr/20210719/subset_6eeffd5...
    search:                  TMP:2 m above

xarray.Dataset

• Dimensions:
  • y: 1059
  • x: 1799
• Coordinates: (7)
  • time
    ()
    datetime64[ns]
    2021-07-19

    long_name :

    initial time of forecast

    standard_name :

    forecast_reference_time

    array('2021-07-19T00:00:00.000000000', dtype='datetime64[ns]')

  • step
    ()
    timedelta64[ns]
    00:00:00

    long_name :

    time since forecast_reference_time

    standard_name :

    forecast_period

    array(0, dtype='timedelta64[ns]')

  • heightAboveGround
    ()
    float64
    2.0

    long_name :

    height above the surface

    units :

    m

    positive :

    up

    standard_name :

    height

    array(2.)

  • latitude
    (y, x)
    float64
    21.14 21.15 21.15 ... 47.85 47.84

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    array([[21.138123  , 21.14511004, 21.1520901 , ..., 21.1545089 ,
            21.14753125, 21.14054663],
           [21.16299459, 21.1699845 , 21.17696744, ..., 21.17938723,
            21.1724067 , 21.16541921],
           [21.18786863, 21.19486142, 21.20184723, ..., 21.20426802,
            21.19728462, 21.19029425],
           ...,
           [47.78955926, 47.799849  , 47.81012868, ..., 47.81369093,
            47.80341474, 47.79312849],
           [47.81409316, 47.82438621, 47.8346692 , ..., 47.83823259,
            47.8279531 , 47.81766354],
           [47.8386235 , 47.84891986, 47.85920615, ..., 47.86277069,
            47.85248789, 47.84219502]], shape=(1059, 1799))

  • longitude
    (y, x)
    float64
    237.3 237.3 237.3 ... 299.0 299.1

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    array([[237.280472  , 237.30713868, 237.3338097 , ..., 287.6569408 ,
            287.68361332, 287.71028151],
           [237.27297501, 237.29964881, 237.32632695, ..., 287.66442108,
            287.69110073, 287.71777603],
           [237.26547368, 237.2921546 , 237.31883986, ..., 287.6719057 ,
            287.69859247, 287.7252749 ],
           ...,
           [225.9351904 , 225.97171577, 226.00825329, ..., 298.97907406,
            299.0156158 , 299.05214538],
           [225.91986142, 225.95639874, 225.99294822, ..., 298.99437498,
            299.03092868, 299.06747022],
           [225.90452027, 225.94106954, 225.97763099, ..., 299.00968806,
            299.04625373, 299.08280723]], shape=(1059, 1799))

  • valid_time
    ()
    datetime64[ns]
    2021-07-19

    standard_name :

    time

    long_name :

    time

    array('2021-07-19T00:00:00.000000000', dtype='datetime64[ns]')

  • gribfile_projection
    ()
    object
    None

    crs_wkt :

    PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6371229,0,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Lambert Conic Conformal (2SP)",ID["EPSG",9802]],PARAMETER["Latitude of false origin",38.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",262.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",38.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",38.5,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

    semi_major_axis :

    6371229.0

    semi_minor_axis :

    6371229.0

    inverse_flattening :

    0.0

    reference_ellipsoid_name :

    unknown

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    unknown

    horizontal_datum_name :

    unknown

    projected_crs_name :

    unknown

    grid_mapping_name :

    lambert_conformal_conic

    standard_parallel :

    (38.5, 38.5)

    latitude_of_projection_origin :

    38.5

    longitude_of_central_meridian :

    262.5

    false_easting :

    0.0

    false_northing :

    0.0

    long_name :

    HRRR model grid projection

    array(None, dtype=object)

• Data variables: (1)
  • t2m
    (y, x)
    float32
    295.3 295.3 295.3 ... 287.7 287.6

    GRIB_paramId :

    167

    GRIB_dataType :

    fc

    GRIB_numberOfPoints :

    1905141

    GRIB_typeOfLevel :

    heightAboveGround

    GRIB_stepUnits :

    1

    GRIB_stepType :

    instant

    GRIB_gridType :

    lambert

    GRIB_uvRelativeToGrid :

    1

    GRIB_DxInMetres :

    3000.0

    GRIB_DyInMetres :

    3000.0

    GRIB_LaDInDegrees :

    38.5

    GRIB_Latin1InDegrees :

    38.5

    GRIB_Latin2InDegrees :

    38.5

    GRIB_LoVInDegrees :

    262.5

    GRIB_NV :

    0

    GRIB_Nx :

    1799

    GRIB_Ny :

    1059

    GRIB_cfName :

    air_temperature

    GRIB_cfVarName :

    t2m

    GRIB_gridDefinitionDescription :

    Lambert Conformal can be secant or tangent, conical or bipolar

    GRIB_iScansNegatively :

    0

    GRIB_jPointsAreConsecutive :

    0

    GRIB_jScansPositively :

    1

    GRIB_latitudeOfFirstGridPointInDegrees :

    21.138123

    GRIB_latitudeOfSouthernPoleInDegrees :

    0.0

    GRIB_longitudeOfFirstGridPointInDegrees :

    237.280472

    GRIB_longitudeOfSouthernPoleInDegrees :

    0.0

    GRIB_missingValue :

    3.4028234663852886e+38

    GRIB_name :

    2 metre temperature

    GRIB_parameterName :

    Temperature

    GRIB_parameterUnits :

    K

    GRIB_shapeOfTheEarth :

    6

    GRIB_shortName :

    2t

    GRIB_stepRange :

    0

    GRIB_units :

    K

    long_name :

    2 metre temperature

    units :

    K

    standard_name :

    air_temperature

    grid_mapping :

    gribfile_projection

    array([[295.33215, 295.33215, 295.33215, ..., 301.64465, 301.64465,
            301.64465],
           [295.33215, 295.33215, 295.33215, ..., 301.64465, 301.64465,
            301.64465],
           [295.33215, 295.33215, 295.33215, ..., 301.64465, 301.64465,
            301.64465],
           ...,
           [287.26965, 287.33215, 287.33215, ..., 287.89465, 287.83215,
            287.76965],
           [287.26965, 287.33215, 287.33215, ..., 287.83215, 287.76965,
            287.70715],
           [287.26965, 287.26965, 287.26965, ..., 287.76965, 287.70715,
            287.64465]], shape=(1059, 1799), dtype=float32)

• Indexes: (0)
• Attributes: (12)

  GRIB_edition :

  2

  GRIB_centre :

  kwbc

  GRIB_centreDescription :

  US National Weather Service - NCEP

  GRIB_subCentre :

  0

  Conventions :

  CF-1.7

  institution :

  US National Weather Service - NCEP

  model :

  hrrr

  product :

  sfc

  description :

  High-Resolution Rapid Refresh - CONUS

  remote_grib :

  https://noaa-hrrr-bdp-pds.s3.amazonaws.com/hrrr.20210719/conus/hrrr.t00z.wrfsfcf00.grib2

  local_grib :

  /home/blaylock/data/hrrr/20210719/subset_6eeffd54__hrrr.t00z.wrfsfcf00.grib2

  search :

  TMP:2 m above

ax = EasyMap("50m", crs=ds.herbie.crs, figsize=[10, 8]).BORDERS().STATES().ax

p = ax.pcolormesh(
    ds.longitude,
    ds.latitude,
    ds.t2m - 273.15,
    transform=pc,
    **paint.NWSTemperature.kwargs2,
)
plt.colorbar(
    p,
    ax=ax,
    orientation="horizontal",
    pad=0.01,
    shrink=0.8,
    **paint.NWSTemperature.cbar_kwargs2,
)

ax.set_title(
    f"{ds.model.upper()}: {H.product_description}\nValid: {ds.valid_time.dt.strftime('%H:%M UTC %d %b %Y').item()}",
    loc="left",
)
ax.set_title(ds.t2m.GRIB_name, loc="right")

Text(1.0, 1.0, '2 metre temperature')