数据
本文使用数据来自:https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.derived.surface.html
北半球冬季平均经向风
难点:各种库的下载与安装,主要难点在xarray上。
因为不需要设置地图,此图还是比较简单的。结果和代码如下。
fig4.2
import matplotlib.pyplot as plt
import xarray as xr
import numpy as np
import matplotlib.ticker as ticker
dir_loc = r"D:\meteo_data\ncar\vwnd.mon.mean.nc"
def DrawVWind():
# 读取数据
vwind = xr.open_dataset(dir_loc)["vwnd"]
# 选择使用的年份
year_label = (vwind['time'].dt.year >= 1990) & (vwind['time'].dt.year <= 2019)
# 选择年份做季节平均
vwindM = vwind[year_label].groupby("time.season").mean(dim="time")
# 选择DJF,做 纬向平均
vwindM2 = vwindM.loc[dict(season="DJF")].mean("lon")