北極振動 12月だけを抜き出し

続き
北極振動
北極振動　年平均 
AO 各年ごとにプロット　および　気候学平均

北極振動の時系列から12月だけを抜き出す。

以下、プログラム。
pick_month.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

def pick_month(series,month):
    """ pick time series only month==month
    """
    mask=(series.dates.month==month)
    out=series[mask].asfreq("A")
    return out

if __name__ == '__main__':

    AOindex_series=readAO()

# pickup December
    AOindex_dec=pick_month(AOindex_series,month=12)

#plot
    fig = tpl.tsfigure()
    fsp = fig.add_tsplot(111)
    fsp.tsplot(AOindex_dec, 'r-')
    plt.xlabel("YEAR")
    plt.ylabel("AO index")
    plt.savefig("AOindex_dec.png")
    plt.show()