北極振動指数の時系列 1/2

*
*
前回の続き。前回もとめたfirst EOFの時系列が、NOAA CPCの発表している北極振動の値と一致しているか確かめる。黒がNOAAのCPCの値で、赤が自分でもとめたもの。
ちょっとだけ値がずれているが、ほとんどあっているのでいいんじゃないかな。前回よくわからなかった北極のとりあつかいとかがずれる原因かもしれない。

以下がプログラム。以前作ったreadAOモジュールでNOAA CPCの北極振動を読む。時系列オブジェクトを作るのと、プロットにはscikits.timeseriesを用いた。

# plot pc1
import scikits.timeseries as ts
import scikits.timeseries.lib.plotlib as tpl

dbeginmonth=ts.Date('M', '1979-1')
pc1_series=ts.time_series(pc1,start_date=dbeginmonth,freq='M')

from readAO import readAO
AOindex_series=readAO()
   
fig = tpl.tsfigure()
fsp = fig.add_tsplot(111)
l1=fsp.tsplot(AOindex_series, 'k-')
l2=fsp.tsplot(pc1_series, 'r-')
istart=ts.Date('M','1979-1').value
iend=ts.Date('M','2000-12').value
fsp.set_xlim(istart,iend)
plt.xlabel("YEAR")
plt.ylabel("AO index")
plt.legend((l1,l2),("NOAA","What I got"))
plt.savefig("check_AOindex.png")
plt.show()