from anuga.shallow_water.data_manager import csv2array
from Numeric import array
        
        
culvert_quantities = 'Culvert No. 2_136649740_timeseries.csv'       


X = csv2array(culvert_quantities)

print X.keys()


E1 = X['E1']
E0 = X['E0']
Q = X['Discharge']

delta_E = E0 - E1


from pylab import ion, plot, show, xlabel, ylabel
ion()

plot(Q, delta_E)
xlabel('Q')
ylabel('\delta E')
show()

raw_input('done')