北極振動 年平均

続き
北極振動

年平均をつくる。

前回のreadAO.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 annual mean
AOindex_annual=ts.convert(AOindex_series,'A',func=np.ma.mean)

#plot
fig = tpl.tsfigure()
fsp = fig.add_tsplot(111)
fsp.tsplot(AOindex_annual, 'r-+')
istart=fsp.get_xlim()[0]
iend=ts.Date('A','2010').value
fsp.set_xlim(istart,iend)
plt.xlabel("YEAR")
plt.ylabel("AO index")
plt.savefig("AOindex_anu.png")
plt.show()