π΅ 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>

[ ]: