CFS

NCEP’s Climate Forecast System (version 2) has been operational since 2011.

https://cfs.ncep.noaa.gov/

There are three products available

• `product="time_series" <#time-series-product>`__

• `product="6_hourly" <#6-hourly-product>`__

• `product="monthly" <#monthly-product>`__

Some examples below

H = Herbie(
    date,
    model="cfs",
    product="time_series",
    member=1,          # {1, 2, 3, 4}
    variable="tmp2m",  # the file's variable short name.
)

H = Herbie(
    date,
    model="cfs",
    product="6_hourly",
    kind="flxf",  # {"flxf", "pgbf", "ocnf", "ipvf"}
    member=1,     # {1, 2, 3, 4}
    fxx=12,       # {0 ... (what's the max forecast hour??}
)

H = Herbie(
    date,
    model="cfs",
    product="monthly_means",
    kind="pgbf",  # {"flxf", "pgbf", "ocnh", "ocnf", "ipvf"}
    member=1,     # {1, 2, 3, 4}
    month=1,      # {0 through 9 (I think)}
    hour=None,    # None for daily mean, or {0, 6, 12, 18}
)

from herbie import Herbie
from herbie.toolbox import EasyMap, pc, ccrs
import matplotlib.pyplot as plt

Time Series product

For reference of what variables are avaialble: https://www.nco.ncep.noaa.gov/pmb/products/cfs/#DAILY

Herbie expects the following arguments for CFS times series product:

• member, an integer 1, 2, 3, or 4

• variable, a file’s variable short name like tmp2m, or z500, etc.

H = Herbie(
    "2025-01-01",
    model="cfs",
    product="time_series",
    member=1,
    variable="tmp2m",
)

H

/home/blaylock/GITHUB/Herbie/herbie/models/cfs.py:139: UserWarning: Herbie's CFS templates are and subject to major changes. PRs are welcome to improve it.
  warnings.warn(
/home/blaylock/GITHUB/Herbie/herbie/models/cfs.py:178: UserWarning: Variable tmp2m is not in the list of available time series variables. Expected one of {'wnd10m': 'Wind u and v at 10 m', 'wnd1000': 'Wind u and v at 1000 hPa', 'wnd925': 'Wind u and v at 925 hPa', 'wnd850': 'Wind u and v at 850 hPa', 'wnd500': 'Wind u and v at 500 hPa', 'wnd700': 'Wind u and v at 700 hPa', 'wnd250': 'Wind u and v at 250 hPa', 'wnd200': 'Wind u and v at 200 hPa', 'wndstrs': 'Stress at surface', 'prmsl': 'Pressure Reduced to MSL', 'pressfc': 'Surface Pressure', 'z1000': 'GeoPotential Height at 1000 hPa', 'z850': 'GeoPotential Height at 850 hPa', 'z700': 'GeoPotential Height at 700 hPa', 'z500': 'GeoPotential Height at 500 hPa', 'z200': 'GeoPotential Height at 200 hPa', 'tmpsfc': 'Temperature at the Surface', 'tmp2m': 'Temperature at 2 m', 'tmin': 'Minimum Temperature at 2 m', 'tmax': 'Maximum Temperature at 2 m', 't1000': 'Temperature at 1000 hPa', 't850': 'Temperature at 850 hPa', 't700': 'Temperature at 700 hPa', 't500': 'Temperature at 500 hPa', 't250': 'Temperature at 250 hPa', 't200': 'Temperature at 200 hPa', 't50': 'Temperature at 50 hPa', 't2': 'Temperature at 2 hPa', 'tmphy1': 'Temperature at hybrid level 1', 'q2m': 'Specific Humidity at 2 m', 'q925': 'Specific Humidity at 925 hPa', 'q850': 'Specific Humidity at 850 hPa', 'q700': 'Specific Humidity at 700 hPa', 'q500': 'Specific Humidity at 500 hPa', 'chi850': 'Velocity Potential at 850 hPa', 'chi200': 'Velocity Potential at 200 hPa', 'dlwsfc': 'Downward Longwave Radiation at the Surface', 'dswsfc': 'Downward Shortwave Radiation at the Surface', 'ulwsfc': 'Upward Longwave Radiation at the Surface', 'uswsfc': 'Upward Shortwave Radiation at the Surface', 'ulwtoa': 'Upward Longwave Radiation at the Top of the Atmosphere', 'uswtoa': 'Upward Shortwave Radiation at the Top of the Atmosphere', 'gflux': 'Ground Heat Flux', 'lhtfl': 'Latent Heat Net Flux at the Surface', 'shtfl': 'Sensible Heat Net Flux at the Surface', 'cprat': 'Convective Precipitation Rate', 'csdlf': 'Clear Sky Downward Longwave Flux', 'csdsf': 'Clear Sky Downward Solar Flux', 'csusf': 'Clear Sky Upward Solar Flux', 'icecon': 'Surface Ice Cover', 'icethk': 'Surface Ice Thickness', 'ipv450': 'Isentropic Potential Vorticity at 450 K', 'ipv550': 'Isentropic Potential Vorticity at 550 K', 'ipv650': 'Isentropic Potential Vorticity at 650 K', 'nddsf': 'Near IR Diffuse Downward Solar Flux', 'tcdcclm': 'Total Cloud Cover', 'prate': 'Total Precipitation Rate', 'psi200': 'Stream Function at 200 hPa', 'psi850': 'Stream Function at 850 hPa', 'pwat': 'Precipitable Water', 'runoff': 'Surface Runoff', 'snohf': 'Snowfall Rate Water Equivalent', 'soilm1': 'Vol. Soil Moisture Content, 0-10 cm Below Ground', 'soilm2': 'Vol. Soil Moisture Content, 10-40 cm Below Ground', 'soilm3': 'Vol. Soil Moisture Content, 40-100 cm Below Ground', 'soilm4': 'Vol. Soil Moisture Content, 100-200 cm Below Ground', 'soilt1': 'Soil Temperature, 0-0.1 m Below Ground', 'vddsf': 'Visible Diffuse Downward Solar Flux', 'vvel500': 'Vertical Velocity at 500 hPa', 'weasd': 'Water Equivalent of Accumulated Snow Depth', 'ocnmld': 'Mixed Layer Depth', 'ocnslh': 'Sea Level Height', 'ocnsild': 'Surface Isothermal Layer Depth', 'ocndt2.5c': 'Isothermal Layer Depth 2.5 Celsius', 'ocndt5c': 'Isothermal Layer Depth 5 Celsius', 'ocndt10c': 'Isothermal Layer Depth 10 Celsius', 'ocndt15c': 'Isothermal Layer Depth 15 Celsius', 'ocndt20c': 'Isothermal Layer Depth 20 Celsius', 'ocndt25c': 'Isothermal Layer Depth 25 Celsius', 'ocndt28c': 'Isothermal Layer Depth 28 Celsius', 'ocnsst': 'Temperature at depth 5 meters', 'ocnt15': 'Potential temperature 15 meters below', 'ocnheat': 'Heat Content', 'ocntchp': 'Tropical Cyclone Heat Potential', 'ocnsal5': 'Salinity at depth 5 meters', 'ocnsal15': 'Salinity at depth 15 meters', 'ocnu5': 'Zonal Current at depth 5 meters', 'ocnu15': 'Zonal Current at depth 15 meters', 'ocnv5': 'Meridional Current at depth 5 meters', 'ocnv15': 'Meridional Current at depth 15 meters', 'ocnvv55': 'Vertical Velocity at depth 55 meters'}
  warnings.warn(

✅ Found ┊ model=cfs ┊ product=time_series ┊ 2025-Jan-01 00:00 UTC F00 ┊ GRIB2 @ aws ┊ IDX @ aws

▌▌Herbie CFS model time_series product initialized 2025-Jan-01 00:00 UTC F00 ┊ source=aws

H.SOURCES

{'aws': 'https://noaa-cfs-pds.s3.amazonaws.com/cfs.20250101/00/time_grib_01/tmp2m.01.2025010100.daily.grb2',
 'nomads': 'https://nomads.ncep.noaa.gov/pub/data/nccf/com/cfs/prod/cfs.20250101/00/time_grib_01/tmp2m.01.2025010100.daily.grb2'}

# What are the available time series variables values?
# You may import the dictionary of values or refer to
# https://www.nco.ncep.noaa.gov/pmb/products/cfs/

from herbie.models.cfs import time_series_variables

time_series_variables

{'wnd10m': 'Wind u and v at 10 m',
 'wnd1000': 'Wind u and v at 1000 hPa',
 'wnd925': 'Wind u and v at 925 hPa',
 'wnd850': 'Wind u and v at 850 hPa',
 'wnd500': 'Wind u and v at 500 hPa',
 'wnd700': 'Wind u and v at 700 hPa',
 'wnd250': 'Wind u and v at 250 hPa',
 'wnd200': 'Wind u and v at 200 hPa',
 'wndstrs': 'Stress at surface',
 'prmsl': 'Pressure Reduced to MSL',
 'pressfc': 'Surface Pressure',
 'z1000': 'GeoPotential Height at 1000 hPa',
 'z850': 'GeoPotential Height at 850 hPa',
 'z700': 'GeoPotential Height at 700 hPa',
 'z500': 'GeoPotential Height at 500 hPa',
 'z200': 'GeoPotential Height at 200 hPa',
 'tmpsfc': 'Temperature at the Surface',
 'tmp2m': 'Temperature at 2 m',
 'tmin': 'Minimum Temperature at 2 m',
 'tmax': 'Maximum Temperature at 2 m',
 't1000': 'Temperature at 1000 hPa',
 't850': 'Temperature at 850 hPa',
 't700': 'Temperature at 700 hPa',
 't500': 'Temperature at 500 hPa',
 't250': 'Temperature at 250 hPa',
 't200': 'Temperature at 200 hPa',
 't50': 'Temperature at 50 hPa',
 't2': 'Temperature at 2 hPa',
 'tmphy1': 'Temperature at hybrid level 1',
 'q2m': 'Specific Humidity at 2 m',
 'q925': 'Specific Humidity at 925 hPa',
 'q850': 'Specific Humidity at 850 hPa',
 'q700': 'Specific Humidity at 700 hPa',
 'q500': 'Specific Humidity at 500 hPa',
 'chi850': 'Velocity Potential at 850 hPa',
 'chi200': 'Velocity Potential at 200 hPa',
 'dlwsfc': 'Downward Longwave Radiation at the Surface',
 'dswsfc': 'Downward Shortwave Radiation at the Surface',
 'ulwsfc': 'Upward Longwave Radiation at the Surface',
 'uswsfc': 'Upward Shortwave Radiation at the Surface',
 'ulwtoa': 'Upward Longwave Radiation at the Top of the Atmosphere',
 'uswtoa': 'Upward Shortwave Radiation at the Top of the Atmosphere',
 'gflux': 'Ground Heat Flux',
 'lhtfl': 'Latent Heat Net Flux at the Surface',
 'shtfl': 'Sensible Heat Net Flux at the Surface',
 'cprat': 'Convective Precipitation Rate',
 'csdlf': 'Clear Sky Downward Longwave Flux',
 'csdsf': 'Clear Sky Downward Solar Flux',
 'csusf': 'Clear Sky Upward Solar Flux',
 'icecon': 'Surface Ice Cover',
 'icethk': 'Surface Ice Thickness',
 'ipv450': 'Isentropic Potential Vorticity at 450 K',
 'ipv550': 'Isentropic Potential Vorticity at 550 K',
 'ipv650': 'Isentropic Potential Vorticity at 650 K',
 'nddsf': 'Near IR Diffuse Downward Solar Flux',
 'tcdcclm': 'Total Cloud Cover',
 'prate': 'Total Precipitation Rate',
 'psi200': 'Stream Function at 200 hPa',
 'psi850': 'Stream Function at 850 hPa',
 'pwat': 'Precipitable Water',
 'runoff': 'Surface Runoff',
 'snohf': 'Snowfall Rate Water Equivalent',
 'soilm1': 'Vol. Soil Moisture Content, 0-10 cm Below Ground',
 'soilm2': 'Vol. Soil Moisture Content, 10-40 cm Below Ground',
 'soilm3': 'Vol. Soil Moisture Content, 40-100 cm Below Ground',
 'soilm4': 'Vol. Soil Moisture Content, 100-200 cm Below Ground',
 'soilt1': 'Soil Temperature, 0-0.1 m Below Ground',
 'vddsf': 'Visible Diffuse Downward Solar Flux',
 'vvel500': 'Vertical Velocity at 500 hPa',
 'weasd': 'Water Equivalent of Accumulated Snow Depth',
 'ocnmld': 'Mixed Layer Depth',
 'ocnslh': 'Sea Level Height',
 'ocnsild': 'Surface Isothermal Layer Depth',
 'ocndt2.5c': 'Isothermal Layer Depth 2.5 Celsius',
 'ocndt5c': 'Isothermal Layer Depth 5 Celsius',
 'ocndt10c': 'Isothermal Layer Depth 10 Celsius',
 'ocndt15c': 'Isothermal Layer Depth 15 Celsius',
 'ocndt20c': 'Isothermal Layer Depth 20 Celsius',
 'ocndt25c': 'Isothermal Layer Depth 25 Celsius',
 'ocndt28c': 'Isothermal Layer Depth 28 Celsius',
 'ocnsst': 'Temperature at depth 5 meters',
 'ocnt15': 'Potential temperature 15 meters below',
 'ocnheat': 'Heat Content',
 'ocntchp': 'Tropical Cyclone Heat Potential',
 'ocnsal5': 'Salinity at depth 5 meters',
 'ocnsal15': 'Salinity at depth 15 meters',
 'ocnu5': 'Zonal Current at depth 5 meters',
 'ocnu15': 'Zonal Current at depth 15 meters',
 'ocnv5': 'Meridional Current at depth 5 meters',
 'ocnv15': 'Meridional Current at depth 15 meters',
 'ocnvv55': 'Vertical Velocity at depth 55 meters'}

H = Herbie(
    "2025-01-01",
    model="cfs",
    product="time_series",
    member=1,
    variable="wnd10m",
)

H.inventory()

/home/blaylock/GITHUB/Herbie/herbie/models/cfs.py:178: UserWarning: Variable wnd10m is not in the list of available time series variables. Expected one of {'wnd10m': 'Wind u and v at 10 m', 'wnd1000': 'Wind u and v at 1000 hPa', 'wnd925': 'Wind u and v at 925 hPa', 'wnd850': 'Wind u and v at 850 hPa', 'wnd500': 'Wind u and v at 500 hPa', 'wnd700': 'Wind u and v at 700 hPa', 'wnd250': 'Wind u and v at 250 hPa', 'wnd200': 'Wind u and v at 200 hPa', 'wndstrs': 'Stress at surface', 'prmsl': 'Pressure Reduced to MSL', 'pressfc': 'Surface Pressure', 'z1000': 'GeoPotential Height at 1000 hPa', 'z850': 'GeoPotential Height at 850 hPa', 'z700': 'GeoPotential Height at 700 hPa', 'z500': 'GeoPotential Height at 500 hPa', 'z200': 'GeoPotential Height at 200 hPa', 'tmpsfc': 'Temperature at the Surface', 'tmp2m': 'Temperature at 2 m', 'tmin': 'Minimum Temperature at 2 m', 'tmax': 'Maximum Temperature at 2 m', 't1000': 'Temperature at 1000 hPa', 't850': 'Temperature at 850 hPa', 't700': 'Temperature at 700 hPa', 't500': 'Temperature at 500 hPa', 't250': 'Temperature at 250 hPa', 't200': 'Temperature at 200 hPa', 't50': 'Temperature at 50 hPa', 't2': 'Temperature at 2 hPa', 'tmphy1': 'Temperature at hybrid level 1', 'q2m': 'Specific Humidity at 2 m', 'q925': 'Specific Humidity at 925 hPa', 'q850': 'Specific Humidity at 850 hPa', 'q700': 'Specific Humidity at 700 hPa', 'q500': 'Specific Humidity at 500 hPa', 'chi850': 'Velocity Potential at 850 hPa', 'chi200': 'Velocity Potential at 200 hPa', 'dlwsfc': 'Downward Longwave Radiation at the Surface', 'dswsfc': 'Downward Shortwave Radiation at the Surface', 'ulwsfc': 'Upward Longwave Radiation at the Surface', 'uswsfc': 'Upward Shortwave Radiation at the Surface', 'ulwtoa': 'Upward Longwave Radiation at the Top of the Atmosphere', 'uswtoa': 'Upward Shortwave Radiation at the Top of the Atmosphere', 'gflux': 'Ground Heat Flux', 'lhtfl': 'Latent Heat Net Flux at the Surface', 'shtfl': 'Sensible Heat Net Flux at the Surface', 'cprat': 'Convective Precipitation Rate', 'csdlf': 'Clear Sky Downward Longwave Flux', 'csdsf': 'Clear Sky Downward Solar Flux', 'csusf': 'Clear Sky Upward Solar Flux', 'icecon': 'Surface Ice Cover', 'icethk': 'Surface Ice Thickness', 'ipv450': 'Isentropic Potential Vorticity at 450 K', 'ipv550': 'Isentropic Potential Vorticity at 550 K', 'ipv650': 'Isentropic Potential Vorticity at 650 K', 'nddsf': 'Near IR Diffuse Downward Solar Flux', 'tcdcclm': 'Total Cloud Cover', 'prate': 'Total Precipitation Rate', 'psi200': 'Stream Function at 200 hPa', 'psi850': 'Stream Function at 850 hPa', 'pwat': 'Precipitable Water', 'runoff': 'Surface Runoff', 'snohf': 'Snowfall Rate Water Equivalent', 'soilm1': 'Vol. Soil Moisture Content, 0-10 cm Below Ground', 'soilm2': 'Vol. Soil Moisture Content, 10-40 cm Below Ground', 'soilm3': 'Vol. Soil Moisture Content, 40-100 cm Below Ground', 'soilm4': 'Vol. Soil Moisture Content, 100-200 cm Below Ground', 'soilt1': 'Soil Temperature, 0-0.1 m Below Ground', 'vddsf': 'Visible Diffuse Downward Solar Flux', 'vvel500': 'Vertical Velocity at 500 hPa', 'weasd': 'Water Equivalent of Accumulated Snow Depth', 'ocnmld': 'Mixed Layer Depth', 'ocnslh': 'Sea Level Height', 'ocnsild': 'Surface Isothermal Layer Depth', 'ocndt2.5c': 'Isothermal Layer Depth 2.5 Celsius', 'ocndt5c': 'Isothermal Layer Depth 5 Celsius', 'ocndt10c': 'Isothermal Layer Depth 10 Celsius', 'ocndt15c': 'Isothermal Layer Depth 15 Celsius', 'ocndt20c': 'Isothermal Layer Depth 20 Celsius', 'ocndt25c': 'Isothermal Layer Depth 25 Celsius', 'ocndt28c': 'Isothermal Layer Depth 28 Celsius', 'ocnsst': 'Temperature at depth 5 meters', 'ocnt15': 'Potential temperature 15 meters below', 'ocnheat': 'Heat Content', 'ocntchp': 'Tropical Cyclone Heat Potential', 'ocnsal5': 'Salinity at depth 5 meters', 'ocnsal15': 'Salinity at depth 15 meters', 'ocnu5': 'Zonal Current at depth 5 meters', 'ocnu15': 'Zonal Current at depth 15 meters', 'ocnv5': 'Meridional Current at depth 5 meters', 'ocnv15': 'Meridional Current at depth 15 meters', 'ocnvv55': 'Vertical Velocity at depth 55 meters'}
  warnings.warn(

✅ Found ┊ model=cfs ┊ product=time_series ┊ 2025-Jan-01 00:00 UTC F00 ┊ GRIB2 @ aws ┊ IDX @ aws

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
0	1	0	59189.0	0-59189	2025-01-01	2025-01-01	UGRD	10 m above ground	6 hour fcst	:UGRD:10 m above ground:6 hour fcst
1	2	59190	118009.0	59190-118009	2025-01-01	2025-01-01	UGRD	10 m above ground	12 hour fcst	:UGRD:10 m above ground:12 hour fcst
2	3	118010	176851.0	118010-176851	2025-01-01	2025-01-01	UGRD	10 m above ground	18 hour fcst	:UGRD:10 m above ground:18 hour fcst
3	4	176852	235807.0	176852-235807	2025-01-01	2025-01-01	UGRD	10 m above ground	24 hour fcst	:UGRD:10 m above ground:24 hour fcst
4	5	235808	295010.0	235808-295010	2025-01-01	2025-01-01	VGRD	10 m above ground	6 hour fcst	:VGRD:10 m above ground:6 hour fcst
...	...	...	...	...	...	...	...	...	...	...
2427	2428	144814943	144875086.0	144814943-144875086	2025-01-01	2025-01-01	VGRD	10 m above ground	7272 hour fcst	:VGRD:10 m above ground:7272 hour fcst
2428	2429	144875087	144934966.0	144875087-144934966	2025-01-01	2025-01-01	VGRD	10 m above ground	7278 hour fcst	:VGRD:10 m above ground:7278 hour fcst
2429	2430	144934967	144994972.0	144934967-144994972	2025-01-01	2025-01-01	VGRD	10 m above ground	7284 hour fcst	:VGRD:10 m above ground:7284 hour fcst
2430	2431	144994973	145054851.0	144994973-145054851	2025-01-01	2025-01-01	VGRD	10 m above ground	7290 hour fcst	:VGRD:10 m above ground:7290 hour fcst
2431	2432	145054852	NaN	145054852-	2025-01-01	2025-01-01	VGRD	10 m above ground	7296 hour fcst	:VGRD:10 m above ground:7296 hour fcst

2432 rows × 10 columns

# Only get 6 hour forecast
H.inventory("[U,V]GRD:10 m.*:6 hour")

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
0	1	0	59189.0	0-59189	2025-01-01	2025-01-01	UGRD	10 m above ground	6 hour fcst	:UGRD:10 m above ground:6 hour fcst
4	5	235808	295010.0	235808-295010	2025-01-01	2025-01-01	VGRD	10 m above ground	6 hour fcst	:VGRD:10 m above ground:6 hour fcst

ds = H.xarray("[U,V]GRD:10 m.*:6 hour")

# Calculate wind speed
ds["wind_speed"] = (ds.u10**2 + ds.v10**2) ** 0.5
ds

<xarray.Dataset> Size: 880kB
Dimensions:              (latitude: 190, longitude: 384)
Coordinates:
    time                 datetime64[ns] 8B 2025-01-01
    step                 timedelta64[ns] 8B 06:00:00
    heightAboveGround    float64 8B 10.0
  * latitude             (latitude) float64 2kB 89.28 88.34 ... -88.34 -89.28
  * longitude            (longitude) float64 3kB 0.0 0.9375 ... 358.1 359.1
    valid_time           datetime64[ns] 8B 2025-01-01T06:00:00
    gribfile_projection  object 8B None
Data variables:
    u10                  (latitude, longitude) float32 292kB 2.62 2.5 ... -3.2
    v10                  (latitude, longitude) float32 292kB -7.05 ... 1.77
    wind_speed           (latitude, longitude) float32 292kB 7.521 ... 3.657
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:                   cfs
    product:                 time_series
    description:             Climate Forecast System
    remote_grib:             https://noaa-cfs-pds.s3.amazonaws.com/cfs.202501...
    local_grib:              /home/blaylock/data/cfs/20250101/subset_6aef9f3c...
    search:                  [U,V]GRD:10 m.*:6 hour

xarray.Dataset

• Dimensions:
  • latitude: 190
  • longitude: 384
• Coordinates: (7)
  • time
    ()
    datetime64[ns]
    2025-01-01

    long_name :

    initial time of forecast

    standard_name :

    forecast_reference_time

    array('2025-01-01T00:00:00.000000000', dtype='datetime64[ns]')

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

    long_name :

    time since forecast_reference_time

    standard_name :

    forecast_period

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

  • heightAboveGround
    ()
    float64
    10.0

    long_name :

    height above the surface

    units :

    m

    positive :

    up

    standard_name :

    height

    array(10.)

  • latitude
    (latitude)
    float64
    89.28 88.34 87.4 ... -88.34 -89.28

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    stored_direction :

    decreasing

    array([ 89.276713,  88.339754,  87.397262,  86.45352 ,  85.509304,  84.564858,
            83.620284,  82.675631,  81.730926,  80.786186,  79.841419,  78.896633,
            77.951832,  77.007019,  76.062197,  75.117368,  74.172532,  73.227692,
            72.282847,  71.337998,  70.393147,  69.448293,  68.503436,  67.558577,
            66.613717,  65.668855,  64.723992,  63.779127,  62.834261,  61.889394,
            60.944527,  59.999658,  59.054789,  58.109919,  57.165049,  56.220177,
            55.275306,  54.330434,  53.385561,  52.440688,  51.495815,  50.550941,
            49.606067,  48.661193,  47.716318,  46.771444,  45.826569,  44.881693,
            43.936818,  42.991942,  42.047066,  41.10219 ,  40.157314,  39.212437,
            38.267561,  37.322684,  36.377807,  35.43293 ,  34.488053,  33.543176,
            32.598299,  31.653422,  30.708544,  29.763667,  28.818789,  27.873911,
            26.929034,  25.984156,  25.039278,  24.0944  ,  23.149522,  22.204644,
            21.259766,  20.314888,  19.370009,  18.425131,  17.480253,  16.535374,
            15.590496,  14.645618,  13.700739,  12.755861,  11.810982,  10.866104,
             9.921225,   8.976347,   8.031468,   7.08659 ,   6.141711,   5.196832,
             4.251954,   3.307075,   2.362197,   1.417318,   0.472439,  -0.472439,
            -1.417318,  -2.362197,  -3.307075,  -4.251954,  -5.196832,  -6.141711,
            -7.08659 ,  -8.031468,  -8.976347,  -9.921225, -10.866104, -11.810982,
           -12.755861, -13.700739, -14.645618, -15.590496, -16.535374, -17.480253,
           -18.425131, -19.370009, -20.314888, -21.259766, -22.204644, -23.149522,
           -24.0944  , -25.039278, -25.984156, -26.929034, -27.873911, -28.818789,
           -29.763667, -30.708544, -31.653422, -32.598299, -33.543176, -34.488053,
           -35.43293 , -36.377807, -37.322684, -38.267561, -39.212437, -40.157314,
           -41.10219 , -42.047066, -42.991942, -43.936818, -44.881693, -45.826569,
           -46.771444, -47.716318, -48.661193, -49.606067, -50.550941, -51.495815,
           -52.440688, -53.385561, -54.330434, -55.275306, -56.220177, -57.165049,
           -58.109919, -59.054789, -59.999658, -60.944527, -61.889394, -62.834261,
           -63.779127, -64.723992, -65.668855, -66.613717, -67.558577, -68.503436,
           -69.448293, -70.393147, -71.337998, -72.282847, -73.227692, -74.172532,
           -75.117368, -76.062197, -77.007019, -77.951832, -78.896633, -79.841419,
           -80.786186, -81.730926, -82.675631, -83.620284, -84.564858, -85.509304,
           -86.45352 , -87.397262, -88.339754, -89.276713])

  • longitude
    (longitude)
    float64
    0.0 0.9375 1.875 ... 358.1 359.1

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    array([  0.      ,   0.937499,   1.874997, ..., 357.187003, 358.124501,
           359.062   ], shape=(384,))

  • valid_time
    ()
    datetime64[ns]
    2025-01-01T06:00:00

    standard_name :

    time

    long_name :

    time

    array('2025-01-01T06:00:00.000000000', dtype='datetime64[ns]')

  • gribfile_projection
    ()
    object
    None

    long_name :

    CFS model grid projection

    array(None, dtype=object)

• Data variables: (3)
  • u10
    (latitude, longitude)
    float32
    2.62 2.5 2.37 ... -3.12 -3.16 -3.2

    GRIB_paramId :

    165

    GRIB_dataType :

    fc

    GRIB_numberOfPoints :

    72960

    GRIB_typeOfLevel :

    heightAboveGround

    GRIB_stepUnits :

    1

    GRIB_stepType :

    instant

    GRIB_gridType :

    regular_gg

    GRIB_uvRelativeToGrid :

    0

    GRIB_N :

    95

    GRIB_NV :

    0

    GRIB_Ni :

    384

    GRIB_Nj :

    190

    GRIB_cfName :

    eastward_wind

    GRIB_cfVarName :

    u10

    GRIB_gridDefinitionDescription :

    Gaussian latitude/longitude

    GRIB_iDirectionIncrementInDegrees :

    0.938

    GRIB_iScansNegatively :

    0

    GRIB_jPointsAreConsecutive :

    0

    GRIB_jScansPositively :

    0

    GRIB_latitudeOfFirstGridPointInDegrees :

    89.277

    GRIB_latitudeOfLastGridPointInDegrees :

    -89.277

    GRIB_longitudeOfFirstGridPointInDegrees :

    0.0

    GRIB_longitudeOfLastGridPointInDegrees :

    359.062

    GRIB_missingValue :

    3.4028234663852886e+38

    GRIB_name :

    10 metre U wind component

    GRIB_parameterName :

    u-component of wind

    GRIB_parameterUnits :

    m s-1

    GRIB_shortName :

    10u

    GRIB_stepRange :

    6

    GRIB_units :

    m s**-1

    long_name :

    10 metre U wind component

    units :

    m s**-1

    standard_name :

    eastward_wind

    grid_mapping :

    gribfile_projection

    array([[ 2.62,  2.5 ,  2.37, ...,  2.96,  2.85,  2.74],
           [ 2.87,  2.71,  2.56, ...,  3.25,  3.12,  2.99],
           [ 2.96,  2.78,  2.6 , ...,  3.43,  3.27,  3.12],
           ...,
           [-3.  , -3.08, -3.16, ..., -2.79, -2.86, -2.93],
           [-2.83, -2.88, -2.93, ..., -2.68, -2.73, -2.78],
           [-3.24, -3.27, -3.3 , ..., -3.12, -3.16, -3.2 ]],
          shape=(190, 384), dtype=float32)

  • v10
    (latitude, longitude)
    float32
    -7.05 -7.08 -7.11 ... 1.82 1.77

    GRIB_paramId :

    166

    GRIB_dataType :

    fc

    GRIB_numberOfPoints :

    72960

    GRIB_typeOfLevel :

    heightAboveGround

    GRIB_stepUnits :

    1

    GRIB_stepType :

    instant

    GRIB_gridType :

    regular_gg

    GRIB_uvRelativeToGrid :

    0

    GRIB_N :

    95

    GRIB_NV :

    0

    GRIB_Ni :

    384

    GRIB_Nj :

    190

    GRIB_cfName :

    northward_wind

    GRIB_cfVarName :

    v10

    GRIB_gridDefinitionDescription :

    Gaussian latitude/longitude

    GRIB_iDirectionIncrementInDegrees :

    0.938

    GRIB_iScansNegatively :

    0

    GRIB_jPointsAreConsecutive :

    0

    GRIB_jScansPositively :

    0

    GRIB_latitudeOfFirstGridPointInDegrees :

    89.277

    GRIB_latitudeOfLastGridPointInDegrees :

    -89.277

    GRIB_longitudeOfFirstGridPointInDegrees :

    0.0

    GRIB_longitudeOfLastGridPointInDegrees :

    359.062

    GRIB_missingValue :

    3.4028234663852886e+38

    GRIB_name :

    10 metre V wind component

    GRIB_parameterName :

    v-component of wind

    GRIB_parameterUnits :

    m s-1

    GRIB_shortName :

    10v

    GRIB_stepRange :

    6

    GRIB_units :

    m s**-1

    long_name :

    10 metre V wind component

    units :

    m s**-1

    standard_name :

    northward_wind

    grid_mapping :

    gribfile_projection

    array([[-7.05, -7.08, -7.11, ..., -6.86, -6.92, -6.99],
           [-8.21, -8.24, -8.28, ..., -7.96, -8.04, -8.12],
           [-9.33, -9.37, -9.41, ..., -9.01, -9.11, -9.22],
           ...,
           [ 2.26,  2.18,  2.1 , ...,  2.36,  2.33,  2.29],
           [ 1.81,  1.74,  1.67, ...,  1.93,  1.89,  1.85],
           [ 1.72,  1.67,  1.61, ...,  1.86,  1.82,  1.77]],
          shape=(190, 384), dtype=float32)

  • wind_speed
    (latitude, longitude)
    float32
    7.521 7.508 7.495 ... 3.647 3.657

    array([[7.521097 , 7.508422 , 7.494598 , ..., 7.471359 , 7.4839096,
            7.507842 ],
           [8.697184 , 8.674197 , 8.666718 , ..., 8.597913 , 8.624152 ,
            8.653005 ],
           [9.788284 , 9.773705 , 9.762587 , ..., 9.6407995, 9.679101 ,
            9.733592 ],
           ...,
           [3.7560084, 3.7734334, 3.7941532, ..., 3.6542714, 3.68897  ,
            3.7187364],
           [3.3593154, 3.364818 , 3.3725066, ..., 3.3026202, 3.3203914,
            3.3392963],
           [3.668242 , 3.671757 , 3.6717978, ..., 3.6323545, 3.6466424,
            3.6568978]], shape=(190, 384), dtype=float32)

• Indexes: (2)
  • latitude
    PandasIndex

    PandasIndex(Index([ 89.27671287810583,  88.33975425118209,   87.3972623773714,
            86.45352030735796,  85.50930400872335,  84.56485823005325,
            83.62028420514433,   82.6756313087777,  81.73092647110545,
            80.78618562235694,
           ...
           -80.78618562235694, -81.73092647110545,  -82.6756313087777,
           -83.62028420514433, -84.56485823005325, -85.50930400872335,
           -86.45352030735796,  -87.3972623773714, -88.33975425118209,
           -89.27671287810583],
          dtype='float64', name='latitude', length=190))

  • longitude
    PandasIndex

    PandasIndex(Index([               0.0, 0.9374986945169713, 1.8749973890339426,
            2.812496083550914, 3.7499947780678853, 4.6874934725848565,
            5.624992167101828,    6.5624908616188, 7.4999895561357715,
            8.437488250652743,
           ...
           350.62451174934785, 351.56201044386484,  352.4995091383818,
            353.4370078328988,  354.3745065274158,  355.3120052219328,
            356.2495039164498,  357.1870026109668, 358.12450130548376,
                      359.062],
          dtype='float64', name='longitude', length=384))

• 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 :

  cfs

  product :

  time_series

  description :

  Climate Forecast System

  remote_grib :

  https://noaa-cfs-pds.s3.amazonaws.com/cfs.20250101/00/time_grib_01/wnd10m.01.2025010100.daily.grb2

  local_grib :

  /home/blaylock/data/cfs/20250101/subset_6aef9f3c__wnd10m.01.2025010100.daily.grb2

  search :

  [U,V]GRD:10 m.*:6 hour

ds["wind_speed"].plot()

<matplotlib.collections.QuadMesh at 0x7f3d33a557f0>

6 Hourly Product

Herbie expects the following additional arguments for CFS 6-hourly product:

• kind, which is flxf, pgbf, ocnf, or ipvf

• member, an integer 1, 2, 3, or 4

• fxx, integer forecast lead time, in hours

# What are the available product kind values?
# You may import the dictionary of values or refer to
# https://www.nco.ncep.noaa.gov/pmb/products/cfs/

from herbie.models.cfs import product_kind

product_kind

{'flxf': 'CFS Surface, Radiative Fluxes',
 'pgbf': 'CFS 3D Pressure Level, 1 degree resolution',
 'ocnh': 'CFS 3D Ocean Data, 0.5 degree resolution',
 'ocnf': 'CFS 3D Ocean Data, 1.0 degree resolution',
 'ipvf': 'CFS 3D Isentropic Level, 1.0 degree resolution'}

H = Herbie(
    "2024-12-25",
    model="cfs",
    product="6_hourly",
    kind="flxf",
    member=1,
    fxx=12,
)

# Show the grib file and index file location
H.grib, H.idx

/home/blaylock/GITHUB/Herbie/herbie/models/cfs.py:139: UserWarning: Herbie's CFS templates are and subject to major changes. PRs are welcome to improve it.
  warnings.warn(

✅ Found ┊ model=cfs ┊ product=6_hourly ┊ 2024-Dec-25 00:00 UTC F12 ┊ GRIB2 @ aws ┊ IDX @ aws

('https://noaa-cfs-pds.s3.amazonaws.com/cfs.20241225/00/6hrly_grib_01/flxf2024122512.01.2024122500.grb2',
 'https://noaa-cfs-pds.s3.amazonaws.com/cfs.20241225/00/6hrly_grib_01/flxf2024122512.01.2024122500.grb2.idx')

H.inventory()

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
0	1.0	0	46418.0	0-46418	2024-12-25	2024-12-25 12:00:00	UFLX	surface	12 hour fcst	:UFLX:surface:12 hour fcst
1	2.0	46419	92097.0	46419-92097	2024-12-25	2024-12-25 12:00:00	VFLX	surface	12 hour fcst	:VFLX:surface:12 hour fcst
2	3.0	92098	129279.0	92098-129279	2024-12-25	2024-12-25 12:00:00	SHTFL	surface	12 hour fcst	:SHTFL:surface:12 hour fcst
3	4.0	129280	169431.0	129280-169431	2024-12-25	2024-12-25 12:00:00	LHTFL	surface	12 hour fcst	:LHTFL:surface:12 hour fcst
4	5.0	169432	252904.0	169432-252904	2024-12-25	2024-12-25 12:00:00	TMP	surface	12 hour fcst	:TMP:surface:12 hour fcst
...	...	...	...	...	...	...	...	...	...	...
98	97.0	3939672	3965411.0	3939672-3965411	2024-12-25	2024-12-25 12:00:00	VBDSF	surface	12 hour fcst	:VBDSF:surface:12 hour fcst
99	98.0	3965412	3990233.0	3965412-3990233	2024-12-25	2024-12-25 12:00:00	VDDSF	surface	12 hour fcst	:VDDSF:surface:12 hour fcst
100	99.0	3990234	4020806.0	3990234-4020806	2024-12-25	2024-12-25 12:00:00	NBDSF	surface	12 hour fcst	:NBDSF:surface:12 hour fcst
101	100.0	4020807	4048869.0	4020807-4048869	2024-12-25	2024-12-25 12:00:00	NDDSF	surface	12 hour fcst	:NDDSF:surface:12 hour fcst
102	101.0	4048870	NaN	4048870-	2024-12-25	2024-12-25 12:00:00	SRWEQ	surface	12 hour fcst	:SRWEQ:surface:12 hour fcst

103 rows × 10 columns

H.inventory("TMP:2 m above")

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
37	37.0	1863350	1941530.0	1863350-1941530	2024-12-25	2024-12-25 12:00:00	TMP	2 m above ground	12 hour fcst	:TMP:2 m above ground:12 hour fcst

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

ax = EasyMap(crs=ds.herbie.crs).ax
art = ax.pcolormesh(ds.longitude, ds.latitude, ds.t2m, transform=pc)

plt.colorbar(
    art,
    ax=ax,
    shrink=0.5,
    pad=0.01,
    orientation="horizontal",
    label="Mean 2-m Temperature (K)",
)

ax.set_title(f"Initialized:{H.date:%Y-%m-%d}\nMember:{H.member} Forecast Hour:{H.fxx}")

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Cell In[13], line 3
      1 ds = H.xarray("TMP:2 m above")
----> 3 ax = EasyMap(crs=ds.herbie.crs).ax
      4 art = ax.pcolormesh(ds.longitude, ds.latitude, ds.t2m, transform=pc)
      6 plt.colorbar(
      7     art,
      8     ax=ax,
   (...)
     12     label="Mean 2-m Temperature (K)",
     13 )

File ~/miniforge3/envs/herbie-dev/lib/python3.13/functools.py:1039, in cached_property.__get__(self, instance, owner)
   1037 val = cache.get(self.attrname, _NOT_FOUND)
   1038 if val is _NOT_FOUND:
-> 1039     val = self.func(instance)
   1040     try:
   1041         cache[self.attrname] = val

File ~/GITHUB/Herbie/herbie/accessors.py:150, in HerbieAccessor.crs(self)
    147 variables = [i for i in list(ds) if len(ds[i].dims) > 0]
    149 ds = ds.metpy.parse_cf(varname=variables)
--> 150 crs = ds.metpy_crs.item().to_cartopy()
    151 return crs

File ~/miniforge3/envs/herbie-dev/lib/python3.13/site-packages/metpy/plots/mapping.py:87, in CFProjection.to_cartopy(self)
     85     proj_handler = self.projection_registry[self.name]
     86 except KeyError:
---> 87     raise ValueError(f'Unhandled projection: {self.name}') from None
     89 return proj_handler(self._attrs, globe)

ValueError: Unhandled projection: unknown

Monthly Product

Herbie expects the following additional arguments for CFS Monthly product:

• kind, which is flxf, pgbf, ocnh, ocnf, or ipvf

• member, an integer 1, 2, 3, or 4

• month, an integer representing the month forecast (1 through 9)

• hour, None for daily mean, or 0, 6, 12, or 18

H = Herbie(
    "2024-05-25",
    model="cfs",
    product="monthly_means",
    kind="pgbf",
    member=1,
    month=1,
    hour=None,  # daily mean
)

/home/blaylock/GITHUB/Herbie/herbie/models/cfs.py:139: UserWarning: Herbie's CFS templates are and subject to major changes. PRs are welcome to improve it.
  warnings.warn(

✅ Found ┊ model=cfs ┊ product=monthly_means ┊ 2024-May-25 00:00 UTC F00 ┊ GRIB2 @ aws ┊ IDX @ aws

H.inventory()

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
0	1.0	0	59628.0	0-59628	2024-05-25	2024-05-25	PRES	mean sea level	1-2 month ave fcst	:PRES:mean sea level:1-2 month ave fcst
1	2.0	59629	81460.0	59629-81460	2024-05-25	2024-05-25	HGT	1 mb	1-2 month ave fcst	:HGT:1 mb:1-2 month ave fcst
2	3.0	81461	100657.0	81461-100657	2024-05-25	2024-05-25	TMP	1 mb	1-2 month ave fcst	:TMP:1 mb:1-2 month ave fcst
3	4.0	100658	104304.0	100658-104304	2024-05-25	2024-05-25	RH	1 mb	1-2 month ave fcst	:RH:1 mb:1-2 month ave fcst
4	5.0	104305	131107.0	104305-131107	2024-05-25	2024-05-25	SPFH	1 mb	1-2 month ave fcst	:SPFH:1 mb:1-2 month ave fcst
...	...	...	...	...	...	...	...	...	...	...
519	473.0	23632602	23656841.0	23632602-23656841	2024-05-25	2024-05-25	VPOT	3 mb	1-2 month ave fcst	:VPOT:3 mb:1-2 month ave fcst
520	474.0	23656842	23674971.0	23656842-23674971	2024-05-25	2024-05-25	STRM	2 mb	1-2 month ave fcst	:STRM:2 mb:1-2 month ave fcst
521	475.0	23674972	23699670.0	23674972-23699670	2024-05-25	2024-05-25	VPOT	2 mb	1-2 month ave fcst	:VPOT:2 mb:1-2 month ave fcst
522	476.0	23699671	23717971.0	23699671-23717971	2024-05-25	2024-05-25	STRM	1 mb	1-2 month ave fcst	:STRM:1 mb:1-2 month ave fcst
523	477.0	23717972	NaN	23717972-	2024-05-25	2024-05-25	VPOT	1 mb	1-2 month ave fcst	:VPOT:1 mb:1-2 month ave fcst

524 rows × 10 columns

H.inventory("TMP:500")

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
207	187.0	8023672	8050166.0	8023672-8050166	2024-05-25	2024-05-25	TMP	500 mb	1-2 month ave fcst	:TMP:500 mb:1-2 month ave fcst

ds = H.xarray("TMP:500 mb")
ds

<xarray.Dataset> Size: 265kB
Dimensions:              (latitude: 181, longitude: 360)
Coordinates:
    time                 datetime64[ns] 8B 2024-05-25
    step                 timedelta64[ns] 8B 60 days
    isobaricInhPa        float64 8B 500.0
  * latitude             (latitude) float64 1kB 90.0 89.0 88.0 ... -89.0 -90.0
  * longitude            (longitude) float64 3kB 0.0 1.0 2.0 ... 358.0 359.0
    valid_time           datetime64[ns] 8B 2024-07-24
    gribfile_projection  object 8B None
Data variables:
    avg_t                (latitude, longitude) float32 261kB 248.5 ... 228.8
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:                   cfs
    product:                 monthly_means
    description:             Climate Forecast System
    remote_grib:             https://noaa-cfs-pds.s3.amazonaws.com/cfs.202405...
    local_grib:              /home/blaylock/data/cfs/20240525/subset_07ef4109...
    search:                  TMP:500 mb

xarray.Dataset

• Dimensions:
  • latitude: 181
  • longitude: 360
• Coordinates: (7)
  • time
    ()
    datetime64[ns]
    2024-05-25

    long_name :

    initial time of forecast

    standard_name :

    forecast_reference_time

    array('2024-05-25T00:00:00.000000000', dtype='datetime64[ns]')

  • step
    ()
    timedelta64[ns]
    60 days

    long_name :

    time since forecast_reference_time

    standard_name :

    forecast_period

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

  • isobaricInhPa
    ()
    float64
    500.0

    long_name :

    pressure

    units :

    hPa

    positive :

    down

    stored_direction :

    decreasing

    standard_name :

    air_pressure

    array(500.)

  • latitude
    (latitude)
    float64
    90.0 89.0 88.0 ... -89.0 -90.0

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    stored_direction :

    decreasing

    array([ 90.,  89.,  88.,  87.,  86.,  85.,  84.,  83.,  82.,  81.,  80.,  79.,
            78.,  77.,  76.,  75.,  74.,  73.,  72.,  71.,  70.,  69.,  68.,  67.,
            66.,  65.,  64.,  63.,  62.,  61.,  60.,  59.,  58.,  57.,  56.,  55.,
            54.,  53.,  52.,  51.,  50.,  49.,  48.,  47.,  46.,  45.,  44.,  43.,
            42.,  41.,  40.,  39.,  38.,  37.,  36.,  35.,  34.,  33.,  32.,  31.,
            30.,  29.,  28.,  27.,  26.,  25.,  24.,  23.,  22.,  21.,  20.,  19.,
            18.,  17.,  16.,  15.,  14.,  13.,  12.,  11.,  10.,   9.,   8.,   7.,
             6.,   5.,   4.,   3.,   2.,   1.,   0.,  -1.,  -2.,  -3.,  -4.,  -5.,
            -6.,  -7.,  -8.,  -9., -10., -11., -12., -13., -14., -15., -16., -17.,
           -18., -19., -20., -21., -22., -23., -24., -25., -26., -27., -28., -29.,
           -30., -31., -32., -33., -34., -35., -36., -37., -38., -39., -40., -41.,
           -42., -43., -44., -45., -46., -47., -48., -49., -50., -51., -52., -53.,
           -54., -55., -56., -57., -58., -59., -60., -61., -62., -63., -64., -65.,
           -66., -67., -68., -69., -70., -71., -72., -73., -74., -75., -76., -77.,
           -78., -79., -80., -81., -82., -83., -84., -85., -86., -87., -88., -89.,
           -90.])

  • longitude
    (longitude)
    float64
    0.0 1.0 2.0 ... 357.0 358.0 359.0

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    array([  0.,   1.,   2., ..., 357., 358., 359.], shape=(360,))

  • valid_time
    ()
    datetime64[ns]
    2024-07-24

    standard_name :

    time

    long_name :

    time

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

  • gribfile_projection
    ()
    object
    None

    long_name :

    CFS model grid projection

    array(None, dtype=object)

• Data variables: (1)
  • avg_t
    (latitude, longitude)
    float32
    248.5 248.5 248.5 ... 228.8 228.8

    GRIB_paramId :

    235130

    GRIB_dataType :

    fc

    GRIB_numberOfPoints :

    65160

    GRIB_typeOfLevel :

    isobaricInhPa

    GRIB_stepUnits :

    1

    GRIB_stepType :

    avg

    GRIB_gridType :

    regular_ll

    GRIB_uvRelativeToGrid :

    0

    GRIB_NV :

    0

    GRIB_Nx :

    360

    GRIB_Ny :

    181

    GRIB_cfName :

    air_temperature

    GRIB_cfVarName :

    avg_t

    GRIB_gridDefinitionDescription :

    Latitude/longitude. Also called equidistant cylindrical, or Plate Carree

    GRIB_iDirectionIncrementInDegrees :

    1.0

    GRIB_iScansNegatively :

    0

    GRIB_jDirectionIncrementInDegrees :

    1.0

    GRIB_jPointsAreConsecutive :

    0

    GRIB_jScansPositively :

    0

    GRIB_latitudeOfFirstGridPointInDegrees :

    90.0

    GRIB_latitudeOfLastGridPointInDegrees :

    -90.0

    GRIB_longitudeOfFirstGridPointInDegrees :

    0.0

    GRIB_longitudeOfLastGridPointInDegrees :

    359.0

    GRIB_missingValue :

    3.4028234663852886e+38

    GRIB_name :

    Time-mean temperature

    GRIB_parameterName :

    Temperature

    GRIB_parameterUnits :

    K

    GRIB_shortName :

    avg_t

    GRIB_stepRange :

    720-1440

    GRIB_units :

    K

    long_name :

    Time-mean temperature

    units :

    K

    standard_name :

    air_temperature

    grid_mapping :

    gribfile_projection

    array([[248.5 , 248.5 , 248.5 , ..., 248.5 , 248.5 , 248.5 ],
           [248.56, 248.57, 248.56, ..., 248.56, 248.56, 248.56],
           [248.74, 248.75, 248.75, ..., 248.74, 248.74, 248.74],
           ...,
           [229.57, 229.56, 229.56, ..., 229.58, 229.57, 229.57],
           [229.14, 229.13, 229.13, ..., 229.15, 229.14, 229.14],
           [228.81, 228.81, 228.81, ..., 228.81, 228.81, 228.81]],
          shape=(181, 360), dtype=float32)

• Indexes: (2)
  • latitude
    PandasIndex

    PandasIndex(Index([ 90.0,  89.0,  88.0,  87.0,  86.0,  85.0,  84.0,  83.0,  82.0,  81.0,
           ...
           -81.0, -82.0, -83.0, -84.0, -85.0, -86.0, -87.0, -88.0, -89.0, -90.0],
          dtype='float64', name='latitude', length=181))

  • longitude
    PandasIndex

    PandasIndex(Index([  0.0,   1.0,   2.0,   3.0,   4.0,   5.0,   6.0,   7.0,   8.0,   9.0,
           ...
           350.0, 351.0, 352.0, 353.0, 354.0, 355.0, 356.0, 357.0, 358.0, 359.0],
          dtype='float64', name='longitude', length=360))

• 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 :

  cfs

  product :

  monthly_means

  description :

  Climate Forecast System

  remote_grib :

  https://noaa-cfs-pds.s3.amazonaws.com/cfs.20240525/00/monthly_grib_01/pgbf.01.2024052500.202406.avrg.grib.grb2

  local_grib :

  /home/blaylock/data/cfs/20240525/subset_07ef4109__pgbf.01.2024052500.202406.avrg.grib.grb2

  search :

  TMP:500 mb

ds.avg_t.plot()

<matplotlib.collections.QuadMesh at 0x7f3d29467890>

ax = EasyMap(crs=ds.herbie.crs).ax
art = ax.pcolormesh(ds.longitude, ds.latitude, ds.t, transform=pc)

plt.colorbar(
    art,
    ax=ax,
    shrink=0.5,
    pad=0.01,
    orientation="horizontal",
    label="Mean 2-m Temperature (K)",
)

ax.set_title(
    f"Initialized:{H.date:%Y-%m-%d}\nMember:{H.member}; Month Forecast:{H.month}; Hour:{H.hour}"
)

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Cell In[26], line 1
----> 1 ax = EasyMap(crs=ds.herbie.crs).ax
      2 art = ax.pcolormesh(ds.longitude, ds.latitude, ds.t, transform=pc)
      4 plt.colorbar(
      5     art,
      6     ax=ax,
   (...)
     10     label="Mean 2-m Temperature (K)",
     11 )

File ~/miniforge3/envs/herbie-dev/lib/python3.13/functools.py:1039, in cached_property.__get__(self, instance, owner)
   1037 val = cache.get(self.attrname, _NOT_FOUND)
   1038 if val is _NOT_FOUND:
-> 1039     val = self.func(instance)
   1040     try:
   1041         cache[self.attrname] = val

File ~/GITHUB/Herbie/herbie/accessors.py:150, in HerbieAccessor.crs(self)
    147 variables = [i for i in list(ds) if len(ds[i].dims) > 0]
    149 ds = ds.metpy.parse_cf(varname=variables)
--> 150 crs = ds.metpy_crs.item().to_cartopy()
    151 return crs

File ~/miniforge3/envs/herbie-dev/lib/python3.13/site-packages/metpy/plots/mapping.py:87, in CFProjection.to_cartopy(self)
     85     proj_handler = self.projection_registry[self.name]
     86 except KeyError:
---> 87     raise ValueError(f'Unhandled projection: {self.name}') from None
     89 return proj_handler(self._attrs, globe)

ValueError: Unhandled projection: unknown

Monthly means by hour (0, 6, 12, or 18)

H = Herbie(
    "2024-05-25",
    model="cfs",
    product="monthly_means",
    kind="pgbf",
    member=1,
    month=2,
    hour=6,  # 2-3 month forecast average for 06z
)
print(H.grib)
H.inventory()

/home/blaylock/GITHUB/Herbie/herbie/models/cfs.py:139: UserWarning: Herbie's CFS templates are and subject to major changes. PRs are welcome to improve it.
  warnings.warn(

✅ Found ┊ model=cfs ┊ product=monthly_means ┊ 2024-May-25 00:00 UTC F00 ┊ GRIB2 @ aws ┊ IDX @ aws
https://noaa-cfs-pds.s3.amazonaws.com/cfs.20240525/00/monthly_grib_01/pgbf.01.2024052500.202407.avrg.grib.06Z.grb2

	grib_message	start_byte	end_byte	range	reference_time	valid_time	variable	level	forecast_time	search_this
0	1.0	0	61046.0	0-61046	2024-05-25	2024-05-25	PRES	mean sea level	2-3 month ave fcst	:PRES:mean sea level:2-3 month ave fcst
1	2.0	61047	85673.0	61047-85673	2024-05-25	2024-05-25	HGT	1 mb	2-3 month ave fcst	:HGT:1 mb:2-3 month ave fcst
2	3.0	85674	110975.0	85674-110975	2024-05-25	2024-05-25	TMP	1 mb	2-3 month ave fcst	:TMP:1 mb:2-3 month ave fcst
3	4.0	110976	115617.0	110976-115617	2024-05-25	2024-05-25	RH	1 mb	2-3 month ave fcst	:RH:1 mb:2-3 month ave fcst
4	5.0	115618	148871.0	115618-148871	2024-05-25	2024-05-25	SPFH	1 mb	2-3 month ave fcst	:SPFH:1 mb:2-3 month ave fcst
...	...	...	...	...	...	...	...	...	...	...
519	473.0	25462570	25496004.0	25462570-25496004	2024-05-25	2024-05-25	VPOT	3 mb	2-3 month ave fcst	:VPOT:3 mb:2-3 month ave fcst
520	474.0	25496005	25515605.0	25496005-25515605	2024-05-25	2024-05-25	STRM	2 mb	2-3 month ave fcst	:STRM:2 mb:2-3 month ave fcst
521	475.0	25515606	25549621.0	25515606-25549621	2024-05-25	2024-05-25	VPOT	2 mb	2-3 month ave fcst	:VPOT:2 mb:2-3 month ave fcst
522	476.0	25549622	25566670.0	25549622-25566670	2024-05-25	2024-05-25	STRM	1 mb	2-3 month ave fcst	:STRM:1 mb:2-3 month ave fcst
523	477.0	25566671	NaN	25566671-	2024-05-25	2024-05-25	VPOT	1 mb	2-3 month ave fcst	:VPOT:1 mb:2-3 month ave fcst

524 rows × 10 columns