NBM

This demonstrates using data from the National Blend of Models (NBM).

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

import matplotlib.pyplot as plt

H = Herbie("2022-05-01 12:00", model="nbm", fxx=1, product="co")

✅ Found ┊ model=nbm ┊ product=co ┊ 2022-May-01 12:00 UTC F01 ┊ GRIB2 @ aws ┊ IDX @ aws

H.PRODUCTS

{'ak': 'Alaska; 13-km resolution',
 'co': 'CONUS 13-km resolution',
 'gu': 'Guam 13-km resolution',
 'hi': 'Hawaii 13-km resolution',
 'pr': 'Puerto Rico 13-km resolution'}

H.grib

'https://noaa-nbm-grib2-pds.s3.amazonaws.com/blend.20220501/12/core/blend.t12z.core.f001.co.grib2'

H.SOURCES

{'nomads': 'https://nomads.ncep.noaa.gov/pub/data/nccf/com/blend/prod/blend.20220501/12/core/blend.t12z.core.f001.co.grib2',
 'aws': 'https://noaa-nbm-grib2-pds.s3.amazonaws.com/blend.20220501/12/core/blend.t12z.core.f001.co.grib2'}

H.inventory(":TMP:2 m.*:1 hour fcst:nan")

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	?	??	search_this
53	54	67526627	68863208.0	67526627-68863208	2022-05-01 12:00:00	2022-05-01 13:00:00	TMP	2 m above ground	1 hour fcst	NaN	NaN	:TMP:2 m above ground:1 hour fcst:nan

ds = H.xarray(":TMP:2 m.*:1 hour fcst:nan")
ds

<xarray.Dataset> Size: 75MB
Dimensions:              (y: 1597, x: 2345)
Coordinates:
    time                 datetime64[ns] 8B 2022-05-01T12:00:00
    step                 timedelta64[ns] 8B 01:00:00
    heightAboveGround    float64 8B 2.0
    latitude             (y, x) float64 30MB 19.23 19.23 19.24 ... 54.38 54.37
    longitude            (y, x) float64 30MB 233.7 233.7 233.8 ... 300.9 301.0
    valid_time           datetime64[ns] 8B 2022-05-01T13:00:00
    gribfile_projection  object 8B None
Dimensions without coordinates: y, x
Data variables:
    t2m                  (y, x) float32 15MB 292.7 292.7 292.7 ... 268.7 268.7
Attributes:
    GRIB_edition:            2
    GRIB_centre:             kwbc
    GRIB_centreDescription:  US National Weather Service - NCEP
    GRIB_subCentre:          14
    Conventions:             CF-1.7
    institution:             US National Weather Service - NCEP
    model:                   nbm
    product:                 co
    description:             National Blend of Models
    remote_grib:             https://noaa-nbm-grib2-pds.s3.amazonaws.com/blen...
    local_grib:              /home/blaylock/data/nbm/20220501/subset_6de4710f...
    search:                  :TMP:2 m.*:1 hour fcst:nan

xarray.Dataset

• Dimensions:
  • y: 1597
  • x: 2345
• Coordinates: (7)
  • time
    ()
    datetime64[ns]
    2022-05-01T12:00:00

    long_name :

    initial time of forecast

    standard_name :

    forecast_reference_time

    array('2022-05-01T12:00:00.000000000', dtype='datetime64[ns]')

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

    long_name :

    time since forecast_reference_time

    standard_name :

    forecast_period

    array(3600000000000, 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
    19.23 19.23 19.24 ... 54.38 54.37

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    array([[19.229     , 19.23419445, 19.23938511, ..., 20.34034804,
            20.33604919, 20.33174644],
           [19.25112363, 19.25631917, 19.26151093, ..., 20.36270393,
            20.3584042 , 20.35410056],
           [19.27324788, 19.27844452, 19.28363738, ..., 20.38506038,
            20.38075976, 20.37645523],
           ...,
           [53.01277632, 53.01903163, 53.02528215, ..., 54.34521296,
            54.3400818 , 54.3349458 ],
           [53.0316719 , 53.03792737, 53.04417805, ..., 54.36413781,
            54.35900656, 54.35387046],
           [53.0505631 , 53.05681873, 53.06306956, ..., 54.38305814,
            54.3779268 , 54.37279061]])

  • longitude
    (y, x)
    float64
    233.7 233.7 233.8 ... 300.9 301.0

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    array([[233.7234    , 233.74683145, 233.77026481, ..., 290.74409026,
            290.76793178, 290.79177169],
           [233.71789568, 233.74133096, 233.76476815, ..., 290.74867371,
            290.77251927, 290.79636321],
           [233.71238946, 233.73582857, 233.75926958, ..., 290.75325878,
            290.77710837, 290.80095636],
           ...,
           [221.64754381, 221.67895681, 221.71037462, ..., 300.87524264,
            300.90771111, 300.94017537],
           [221.63714024, 221.66855972, 221.69998402, ..., 300.88404341,
            300.91651908, 300.94899054],
           [221.62673185, 221.65815781, 221.6895886 , ..., 300.89284839,
            300.92533128, 300.95780994]])

  • valid_time
    ()
    datetime64[ns]
    2022-05-01T13:00:00

    standard_name :

    time

    long_name :

    time

    array('2022-05-01T13:00:00.000000000', dtype='datetime64[ns]')

  • gribfile_projection
    ()
    object
    None

    crs_wkt :

    PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6371200,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",25,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",265,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",25,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",25,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 :

    6371200.0

    semi_minor_axis :

    6371200.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 :

    (25.0, 25.0)

    latitude_of_projection_origin :

    25.0

    longitude_of_central_meridian :

    265.0

    false_easting :

    0.0

    false_northing :

    0.0

    long_name :

    NBM model grid projection

    array(None, dtype=object)

• Data variables: (1)
  • t2m
    (y, x)
    float32
    292.7 292.7 292.7 ... 268.7 268.7

    GRIB_paramId :

    167

    GRIB_dataType :

    fc

    GRIB_numberOfPoints :

    3744965

    GRIB_typeOfLevel :

    heightAboveGround

    GRIB_stepUnits :

    1

    GRIB_stepType :

    instant

    GRIB_gridType :

    lambert

    GRIB_DxInMetres :

    2539.703

    GRIB_DyInMetres :

    2539.703

    GRIB_LaDInDegrees :

    25.0

    GRIB_Latin1InDegrees :

    25.0

    GRIB_Latin2InDegrees :

    25.0

    GRIB_LoVInDegrees :

    265.0

    GRIB_NV :

    0

    GRIB_Nx :

    2345

    GRIB_Ny :

    1597

    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 :

    19.229

    GRIB_latitudeOfSouthernPoleInDegrees :

    -90.0

    GRIB_longitudeOfFirstGridPointInDegrees :

    233.7234

    GRIB_longitudeOfSouthernPoleInDegrees :

    0.0

    GRIB_missingValue :

    3.4028234663852886e+38

    GRIB_name :

    2 metre temperature

    GRIB_parameterName :

    Temperature

    GRIB_parameterUnits :

    K

    GRIB_shortName :

    2t

    GRIB_stepRange :

    1

    GRIB_units :

    K

    GRIB_uvRelativeToGrid :

    0

    long_name :

    2 metre temperature

    units :

    K

    standard_name :

    air_temperature

    grid_mapping :

    gribfile_projection

    array([[292.69   , 292.69   , 292.69   , ..., 298.28998, 298.28998,
            298.28998],
           [292.69   , 292.69   , 292.69   , ..., 298.28998, 298.28998,
            298.28998],
           [292.69   , 292.69   , 292.69   , ..., 298.28998, 298.28998,
            298.28998],
           ...,
           [279.49   , 279.49   , 279.49   , ..., 268.28998, 268.69   ,
            269.09   ],
           [279.49   , 279.49   , 279.49   , ..., 268.69   , 268.69   ,
            268.69   ],
           [279.49   , 279.49   , 279.49   , ..., 268.69   , 268.69   ,
            268.69   ]], dtype=float32)

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

  GRIB_edition :

  2

  GRIB_centre :

  kwbc

  GRIB_centreDescription :

  US National Weather Service - NCEP

  GRIB_subCentre :

  14

  Conventions :

  CF-1.7

  institution :

  US National Weather Service - NCEP

  model :

  nbm

  product :

  co

  description :

  National Blend of Models

  remote_grib :

  https://noaa-nbm-grib2-pds.s3.amazonaws.com/blend.20220501/12/core/blend.t12z.core.f001.co.grib2

  local_grib :

  /home/blaylock/data/nbm/20220501/subset_6de4710f__blend.t12z.core.f001.co.grib2

  search :

  :TMP:2 m.*:1 hour fcst:nan

ax = EasyMap(crs=ds.herbie.crs, figsize=[8, 8]).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.05, **paint.NWSTemperature.cbar_kwargs2
)

ax.set_title(ds.t2m.GRIB_name, loc="right")
ax.set_title(f"{ds.model.upper()}: {H.product_description}", loc="left")

Text(0.0, 1.0, 'NBM: CONUS 13-km resolution')

H.inventory(":WIND:10 m above ground:1 hour fcst:nan")

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	?	??	search_this
68	69	97035862	98393619.0	97035862-98393619	2022-05-01 12:00:00	2022-05-01 13:00:00	WIND	10 m above ground	1 hour fcst	NaN	NaN	:WIND:10 m above ground:1 hour fcst:nan

ds = H.xarray(":WIND:10 m above ground:1 hour fcst:nan", verbose=True)

📇 Download subset: ▌▌Herbie NBM model co product initialized 2022-May-01 12:00 UTC F01 ┊ source=aws
 cURL from https://noaa-nbm-grib2-pds.s3.amazonaws.com/blend.20220501/12/core/blend.t12z.core.f001.co.grib2
Found 1 grib messages.
Download subset group 1
  69  :WIND:10 m above ground:1 hour fcst:nan
curl -s --range 97035862-98393619 "https://noaa-nbm-grib2-pds.s3.amazonaws.com/blend.20220501/12/core/blend.t12z.core.f001.co.grib2" > "/home/blaylock/data/nbm/20220501/subset_6de4000a__blend.t12z.core.f001.co.grib2"
💾 Saved the subset to /home/blaylock/data/nbm/20220501/subset_6de4000a__blend.t12z.core.f001.co.grib2

ax = EasyMap(crs=ds.herbie.crs, figsize=[8, 8]).STATES().ax
p = ax.pcolormesh(
    ds.longitude, ds.latitude, ds.si10, transform=pc, **paint.NWSWindSpeed.kwargs2
)
plt.colorbar(
    p, ax=ax, orientation="horizontal", pad=0.05, **paint.NWSWindSpeed.cbar_kwargs2
)

ax.set_title(ds.si10.GRIB_name, loc="right")
ax.set_title(f"{ds.model.upper()}: {H.product_description}", loc="left")

Text(0.0, 1.0, 'NBM: CONUS 13-km resolution')