前言
有些时候为了研究不同高度上气象要素之间的联系,纬度-高度作为xy轴的图在一些SCI论文中比较常见。
这是我研究的CMIP6数据中ua,va,wap,ta这几个气象要素在纬度-高度图上的不同。
读取数据的地方就略去了,需要注意的是需要进行平均,温度转为摄氏度
#W,V为垂直方向和纬向的速度
W=W.wap.mean(("lon", "time",),keep_attrs=True,skipna=True)*100
V=V.va.mean(("lon", "time",),keep_attrs=True,skipna=True)
T_c=T_all.ta-274.15
T=T_c.mean(("lon", "time",),keep_attrs=True,skipna=True)
# 这时数据格式为(plevel,lat)的格式
数据的x轴为lat,y轴为plev。
理解了数据的格式后,转换一下思路就可以画图。
之后就是画图,因为垂直和纬向上风速的量级差的有点大,所以我对这两个风速进行了均一化的处理。
这样箭头就一样大了。
操作
plot=ax.contourf(V.lat, V.plev, V,levels=np.arange(-0.06,0.06+0.01,0.01),
cmap=cmaps.GMT_panoply,
extend='both'
# pivot='middle',
)
cb&