zoukankan      html  css  js  c++  java
  • MeteoInfoLab脚本示例:计算垂直螺旋度

    尝试编写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)

  • 相关阅读:
    感悟.
    近期感悟
    android 新建项目无法自动生成R文件
    重装后各种碰壁
    rhapsody
    libevent
    ui android需要解决的问题
    sql对xml的解析
    将datarow转换为实体的方法
    网页中的一些我不熟的东西
  • 原文地址:https://www.cnblogs.com/yaqiang/p/4865732.html
Copyright © 2011-2022 走看看