【分享】水平风矢量的高度-时间变化

pkuzhh

上图为某一格点上空各层水平风矢量随时间的演变，可用于分析垂直风切变。代码如下：

f = addfile(r'E:\ERA\200708.nc')
lat = '36.5'
lon = '102'
tidx = 24
u = f['u'][tidx*4-2:tidx*4+5,:,lat,lon]
v = f['v'][tidx*4-2:tidx*4+5,:,lat,lon]
u = u.T
v = v.T
u = u[::-1,:]
v = v[::-1,:]
lev1 = u.dimvalue(0)
lev1 = lev1[::-1]
lev2 = meteo.p2h(lev1)
lev2 = lev2[:]/1000
u.setdimvalue(0,lev2)
v.setdimvalue(0,lev2)
lev3 = []
for j in range(0, len(lev1)):
    lev3.append('%i' % lev1[j])   
#tidx = u.dimvalue(1)
speed = sqrt(u*u+v*v)
#layer1 = quiver(u, v, speed, 10, size=20,color='gray')
layer2 = barbs(u, v, speed, 5, size=15,color='gray')
#quiverkey(layer, 0.84, 0.18, 10, bbox={'edge':True, 'fill':True})
yticks(u.dimvalue(0), lev3)
#xlim('2007-08-22,2007-08-27')
ylim(2.3,9.82)
ylabel('hPa',fontsize=18)

独孤酒见

抢占沙发，楼主强人。