北極振動時の日本付近の冬の気圧パターン

*
*
12月の1000hpaのgeopotentail heightに前にもとめた北極振動のパターンを重ねてみる。北極振動指数が正だと西高東低が弱まり、負だと強まる。

北極振動指数0の時

 

北極振動+2の時

 北極振動-2の時

  

以下、プログラム

from mpl_toolkits.basemap import Basemap  # for Basemap toolkit

import matplotlib.pyplot as plt
import numpy as np

lat=np.load('lat.npy')  # latitude
lon=np.load('lon.npy')  # longitude
lat[-1]=lat[-1]-1.E-4  # lat=90 does not work

#average field
mhgt=np.load('mhgt.npy') # seasonal climatology
dec_average=mhgt[11,:,:] # december

#AO pattern
factor = np.sqrt(np.cos(np.pi*lat/180.))
eof1=np.load('eof1.npy')
AO_pattern=eof1/factor[:,np.newaxis]


#------------- plot geopotential height (AO index 0, 2,-2)
for index in [0.,2.,-2.]:
    data=dec_average+AO_pattern*index
    fig=plt.figure()
    ax = fig.add_axes([0.1,0.1,0.8,0.8])
    map = Basemap(width=15.e6,height=12.e6,\
            projection='gnom',lat_0=60.,lon_0=135.)
    map.drawcoastlines()

    lons,lats=np.meshgrid(lon[:],lat[:])
    x,y=map(lons,lats)
    cy=map.contourf(x,y,data,np.arange(-100,300,20),extend='both')
    cz=map.contour(x,y,data,np.arange(-100,300,20),colors='k')
    plt.clabel(cz,fmt='%3.0f',)
    plt.savefig('AO_patter_'+str(index)+'.png')
    plt.show()