# coding: utf-8
import numpy as np
import matplotlib.pyplot as pl
import datetime as DT

f=open('dat_TR_maebashi.csv','r')

data_xt=[]
data_xr=[]
data_yt=[]
data_yr=[]
data_du=[]
i=0
hm = DT.timedelta(days=14)
for str in f.readlines():
    if i==0:
        i=i+1
        continue
    else:
        str=str.strip()
        str=str.split(',')
        data_xt=data_xt+[DT.datetime.strptime(str[0],'%Y/%m/%d')]
        data_xr=data_xr+[DT.datetime.strptime(str[0],'%Y/%m/%d')-hm]
        data_yt=data_yt+[float(str[1])]
        data_yr=data_yr+[float(str[4])]
        data_du=data_du+[-10]
f.close()

fig = pl.figure(figsize=(10,6),dpi=100)
pl.subplot(1,1,1)
pl.grid()
pl.title('Temperaturs and Rainfall in Maebashi',fontsize=16)
pl.xlim(xmin=DT.date(2000,1,1))
pl.ylim(-30,40)
pl.xlabel('Year')
pl.ylabel('Average Temperature in month ($^\circ$C)')
pl.plot(data_xt,data_yt,'r-',linewidth=2,label='Temperature')
pl.legend(shadow=True, prop={'size':12},loc='upper right')

pl.twinx()
pl.xlim(xmin=DT.date(2000,1,1))
pl.ylim(0,700)
pl.ylabel('Rainfall (mm/month)')
pl.bar(data_xr,data_yr,color='cyan',edgecolor='blue',width=30,label='Rainfall',alpha=1.0)
pl.plot(data_xt,data_du,'r-',linewidth=2,label='Temperature')
pl.legend(shadow=True, prop={'size':12},loc='upper right')

#pl.show()
pl.savefig("fig_time_his_TR.png", dpi=100)