import numpy as np
import matplotlib.pyplot as plt
import time
from pylab import *
from datetime import datetime
from datetime import timedelta

ephindata=np.loadtxt('/home/pacifix/kuehl/work/data/dayly_ephin/ephin_hourly.dat')

yearephin=ephindata[:,0]
dayephin=ephindata[:,1]
hourephin=ephindata[:,2]/3600000.
dateephin=[(datetime(int(yearephin[date]),1,1,int(hourephin[date]))+timedelta(int(dayephin[date])-1)) for date in range(len(yearephin))]
e150=ephindata[:,6]
e300=ephindata[:,7]
e1300=ephindata[:,8]
e3000=ephindata[:,9]
inte=ephindata[:,18]


plt.figure()
p1=plt.subplot(4,1,1)
plt.plot_date( dateephin, e150, 'k-')
plt.ylim([0,2])
plt.ylabel("e150")
p1=plt.subplot(4,1,2)
plt.plot_date( dateephin, e300, 'k-')
plt.ylim([0,0.2])
plt.ylabel("e300")
p1=plt.subplot(4,1,3)
plt.plot_date( dateephin, e1300, 'k-')
plt.ylim([0,0.02])
plt.ylabel("e1300")
p1=plt.subplot(4,1,4)
plt.plot_date( dateephin, inte, 'k-')
plt.ylim([0,0.5])
plt.ylabel("inte")
plt.show()