source: production/sydney_2006/plot_integral.py @ 3330

"""Read in csv file and plot time series of integral (should be conserved
once water stops entering the boundary.
Jane Sexton, GA 2006

"""

import project
from pyvolution.util import file_function
from pylab import *

def get_integral_from_file(filename):
    fid = open(filename)
    lines = fid.readlines()
    fid.close()

    times = []
    volumes = []
    for line in lines:
        fields = line.split(',')
        time = float(fields[0])
        volume = float(fields[1])
        times.append(time)
        volumes.append(volume)
        
    return times, volumes

t, integral = get_integral_from_file(project.integraltimeseries+'.csv')

ion()
hold(False)
plot(t, integral, '.-b')
axis([0, max(t), min(integral)-0.1, min(integral) + 1])
title('Checking for water conservation')
xlabel('Time (sec)')
ylabel('Volume')    
savefig('Volume conservation')
    
# raw_input('Next')

show()