AO 各年ごとにプロット および 気候学平均

続き
北極振動
北極振動　年平均 

AOを各年毎にプロット（黒線）
2009年だけ赤線
気候学的平均（各月ごと、青線）

2009年の12月は異常に値が小さかったことがわかる。

以下プログラム。

plot_AO_season.py
from readAO import readAO
#
import numpy as np
import matplotlib.pyplot as plt
import scikits.timeseries as ts
import scikits.timeseries.lib.plotlib as tpl


AOindex_series=readAO()
# to seasonal
AOindex_season=ts.extras.convert_to_annual(AOindex_series)


#plot each year
fig = tpl.tsfigure()
fsp = fig.add_tsplot(111)

for n in range(AOindex_season.shape[0]):
    series = ts.time_series(AOindex_season[n,:], start_date=ts.Date(freq="m",year=1990,month=1),length=12)
    lstyle='k-'
    if AOindex_season.dates.year[n] == 2009:
        lstyle='r-'
    fsp.tsplot(series, lstyle)

# climatological average

AOindex_average = ts.time_series(np.mean(AOindex_season,axis=0), start_date=ts.Date(freq="m",year=1990,month=1),length=12)
fsp.tsplot(AOindex_average, 'b-',linewidth=2)
fsp.xaxis.set_ticklabels(["Jan"]) # dirty trick for plot
plt.savefig("AO_season.png")
plt.show()