MeteoinfoLab脚本示例：3维水平、经向剖面图

醉后不知天在水

MeteoInfo 发表于 2019-7-16 16:44
已加入脚本汇总帖。

王老师，三维的风场矢量图可以用quiver函数，画类似上面的切面图形吗？我尝试着写了脚本，但是还是报错f=addfile(r'E:\UVW.nc')
udata = f['u'][:,::-1,'40:60','120:140']
vdata = f['v'][:,::-1,'40:60','120:140']
wdata=f['w'][:,::-1,'40:60','120:140']
lev1=f['level'][::-1]
lat=f['latitude'][::-1]
lon=f['longitude'][:]
u=udata[16,'850','40:60','120:140']
v=vdata[16,'850','40:60','120:140']
w=wdata[16,'850','40:60','120:140']
u=u[::2,::2]
v=v[::2,::2]
w=w[::2,::2]
speed = sqrt(u*u+v*v+w*w)
ax = axes3d(position=[0.11, 0.11, 0.775, 0.88],tickfontsize=24,bbox=True)
#china = shaperead('china')
layer=ax.quiver3(lat,lon,lev1,u,v,w,size=20)

liushui513

楼主赞一个，谢谢分享