北極振動 気圧場

*
*
2008年と2009年の気圧場（1000hpa Geopotential height）を比較する。
2009年は北極振動のパターンと逆の符号（負のフェーズ）の傾向が強い。

以下プログラム。

• python interface for grads
• matplotlib
• numpy
• basmap toolkit

を使用。
データはOpendapとい仕組みでNOAAからNCEP/NCARのデータを取得

from grads.ganum import GaNum  # for python interface for grads

from mpl_toolkits.basemap import Basemap  # for Basemap toolkit
import matplotlib.pyplot as plt
import numpy as np


for nyear in [2008,2009]:
# data import using python interface for grads
    ga = GaNum(Bin='grads')
    fh=ga.open("http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/ncep.reanalysis.derived/pressure/hgt.mon.mean.nc") # data open
    ga("set time dec%(nyear)s" % locals())
    ga("set lat 0 90")
    ga("q file")
    ga("set lev 1000") # 1000 hpa
    hgt=ga.exp("hgt")
    del ga

    lat=hgt.grid.lat
    lon=hgt.grid.lon

# map
    fig=plt.figure()
    ax = fig.add_axes([0.1,0.1,0.8,0.8])
  
    map = Basemap(projection='nplaea',boundinglat=10,lon_0=135,resolution='i')
    map.drawcoastlines()
    map.drawmapboundary()
    lons,lats=np.meshgrid(lon[:],lat[:])
    x,y=map(lons,lats)
    map.fillcontinents()
    cz=map.contour(x,y,hgt[:,:],np.arange(-100,300,20))
    plt.clabel(cz,fmt="%3.0f")
    plt.title("Geopotential height (m): Dec %(nyear)s " % locals())
    plt.savefig("gpt_%(nyear)s.png" % locals())
    plt.show()