🌡 Vapor Pressure Deficit#

[1]:
from herbie import Herbie
from herbie.toolbox import EasyMap, pc
from herbie.toolbox.pressure import vapor_pressure_deficit
[2]:
ds = Herbie("2021-04-23").xarray(r"(?:TMP|RH):2 m")
βœ… Found β”Š model=hrrr β”Š product=sfc β”Š 2021-Apr-23 00:00 UTC F00 β”Š GRIB2 @ aws β”Š IDX @ aws
πŸ‘¨πŸ»β€πŸ­ Created directory: [/home/blaylock/data/hrrr/20210423]
[3]:
vpd = vapor_pressure_deficit(ds["t2m"] - 273.15, ds["r2"])

crs = ds.herbie.crs
[4]:
ax = EasyMap("50m", crs=crs, figsize=[10, 5], dpi=150).STATES().ax
vpd.plot.pcolormesh(
    x="longitude",
    y="latitude",
    transform=pc,
    cmap="Oranges",
    ax=ax,
    cbar_kwargs=dict(label="2 m Vapor Pressure Deficit (hPa)"),
)
[4]:
<cartopy.mpl.geocollection.GeoQuadMesh at 0x7f0756549100>
../../_images/gallery_bonus_notebooks_vapor_pressure_deficit_4_1.png
[ ]: