import netCDF4
import wrf


# 读取数据###############################################################
data_nc = netCDF4.Dataset(r"./data/WRF/wrfout_d01_2023-01-15_12%3A00%3A00")
slp = wrf.getvar(data_nc, "slp")  # 读取计算海表面气压
print(slp.shape)
# ####slp是一个DataArray，使用标签的多维数组，主要的关键属性包括：
# values、dims、coords、attrs
# 元数据属性的字典
print(slp.attrs)
# 维度
print(slp.dims)
# 包含数组坐标的类似字典的容器
print(slp.coords)
# 数组值
print(slp.values)
# print(slp[0:1])
# smooth_slp = wrf.smooth2d(slp, 3, cenweight=4)
# lats, lons = wrf.latlon_coords(slp)