source: inundation/ga/storm_surge/validation/LWRU1/icplot.py @ 1627

filename = 'lwru1_IC.txt'


#x = []
#w = []
#fid = open(filename)
#for line in fid.readlines():
#    fields = line.split()
#    assert len(fields) == 2, '%s' %fields
#    x.append( float(fields[0]) )
#    w.append( float(fields[1]) )            
#
#print 'X', x, len(x)
#print 'W', w, len(w)

from pylab import plot, load, ion, ioff, draw, show, sin, hold

X = load(filename)
ion()
plot(X[:,0], X[:,1])
#draw()

hold(False)
for i in range(100000):
    print i*i
    plot(X[:,0], X[:,1] + sin(i/10.0))    



show()
print 'end'