串本と浦神の潮位差

串本と浦神の潮位差は黒潮の流路と密接な関係があることが知られている（参考:気象庁）。串本と浦神の潮位差の月平均値のデータは気象庁から手に入れることができる。http://www.data.kishou.go.jp/kaiyou/shindan/b_2/kuroshio_stream/ksur.txt

以下はそれを読むプログラム。基本的には前回と同じだが、欠損値があるので一手間が必要。

 

read_kushimoto_uragami.py

import numpy as np
import scikits.timeseries as ts

def read_kushimoto_uragami():
    day=[]
    value=[]
    ds = np.DataSource(None)
    f=ds.open('http://www.data.kishou.go.jp/kaiyou/shindan/b_2/kuroshio_stream/ksur.txt')
    for line in f:
         fields = line.strip().split()
         if len(fields)>=1: day.append(fields[0])
         if len(fields)==1:
             value.append(-999.)
         elif len(fields)==2:
             value.append(fields[1])
    f.close()
    valuea=np.asarray(value,dtype='float32')
    kushimoto_uragami=ts.time_series(valuea,start_date=day[0],freq='M',mask=(valuea<-990.))
    return kushimoto_uragami

if __name__ == '__main__':
    import matplotlib.pyplot as plt
    import scikits.timeseries.lib.plotlib as tpl

    kushimoto_uragami=read_kushimoto_uragami()

    fig = tpl.tsfigure()
    fsp = fig.add_tsplot(111)
    fsp.tsplot(kushimoto_uragami, '-')
    plt.xlabel("YEAR")
    plt.ylabel("Kushimoto-Uragami (cm)")
    plt.savefig("kushimoto_uragami.png")
    plt.show()