使用xarray将metpy计算的变量写入netcdf

import xarray as xr
from metpy.units import units

# 假设使用metpy计算了气温和降水变量
temperature = ...
precipitation = ...

ds = xr.Dataset()
ds['temperature'] = xr.DataArray(temperature, dims=('time', 'lat', 'lon'), attrs={'units': 'C'})
ds['precipitation'] = xr.DataArray(precipitation, dims=('time', 'lat', 'lon'), attrs={'units': 'mm'})

ds.to_netcdf('output.nc')