metpy
import metpy
from metpy.units import units
ds = ds.metpy.parse_cf()
g = units.Quantity(9.80665, "m / s^2")
Calculations
import metpy.calc as mpcalc
from metpy.cbook import example_data
ds = example_data()
wind_speed = mpcalc.wind_speed(ds["uwind"], ds["vwind"])
wind_direction = mpcalc.wind_direction(ds["uwind"], ds["vwind"], convention="from")
wind_direction_categorical = mpcalc.angle_to_direction(225. * units.deg)
import pint_xarray
from metpy.calc import saturation_vapor_pressure, specific_humidity_from_dewpoint, wind_direction
wdir = wind_direction(ds["u10"], ds["v10"])
wdir = wdir.pint.dequantify(format="~P")
wdir["attrs"].long_name = "10 metre wind direction from"
specific_humidity_from_dewpoint
Dev setup
$ conda create -n devel
$ conda activate devel
$ conda install --file ci/requirements.txt --file ci/extra_requirements.txt --file ci/test_requirements.txt --y
$ pip install -e .