MeteoInfoLab脚本示例：计算垂直螺旋度

MeteoInfoLab计算垂直螺旋度

最新推荐文章于 2022-03-31 13:33:37 发布

weixin_30809333

最新推荐文章于 2022-03-31 13:33:37 发布

阅读量302

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/yaqiang/p/4865732.html

本文介绍了一个使用MeteoInfoLab脚本计算垂直螺旋度的方法。该脚本加载了2011年特定日期的风速和角速度数据文件，并计算了指定纬度、经度范围内的垂直螺旋度。

尝试编写MeteoInfoLab脚本计算垂直螺旋度，结果未经验证。

脚本程序：

print 'Open data files...'
f_uwnd = addfile('D:/Temp/nc/uwnd.2011.nc')
f_vwnd = addfile('D:/Temp/nc/vwnd.2011.nc')
f_omega = addfile('D:/Temp/nc/omega.2011.nc')

print 'Calculate vertical helicity...'
tidx = 173    # Jun 23, 2011
t = f_uwnd.gettime(tidx)
level = [1000, 100]
lat = [15,55]
lon = [70,135]
uwnd = f_uwnd['uwnd'][tidx,level,lat,lon][:,::-1,:]
vwnd = f_vwnd['vwnd'][tidx,level,lat,lon][:,::-1,:]
omega = f_omega['omega'][tidx,level,lat,lon][:,::-1,:]
wd = hcurl(uwnd, vwnd)
lx = -(wd*omega*10.)/12.64*1e6
lx1 = lx[:,[40],:]
lev1 = lx1.dimvalue(0)
#lev2 = 1000 - lev1
lev2 = p2h(lev1)
levels = []
for i in range(0, len(lev1)):
    levels.append('%i' % lev1[i])
lx1.setdimvalue(0, lev2)

print 'Plot...'
layer = contourf(lx1, 20)
title('Vertical helicity (' + t.strftime('%Y-%m-%d') + ')')
yticks(lev2, levels)
xlabel('Longitude')
ylabel('Pressure (hPa)')
colorbar(layer)