Changeset 8991 for trunk/anuga_work/development/gareth/tests/parabolic/parabolaplot.py

Timestamp:

Sep 26, 2013, 1:10:23 PM (12 years ago)

Author:

davies

Message:

Updates to audusse experimental

File:

• trunk/anuga_work/development/gareth/tests/parabolic/parabolaplot.py (modified) (3 diffs)

trunk/anuga_work/development/gareth/tests/parabolic/parabolaplot.py

@@ -9 +9 @@
 #import scipy
 from matplotlib import pyplot as pyplot
-import util
+#import util
+from bal_and import plot_utils as util
 
 p=util.get_output('parabola_v2.sww', 0.001)
@@ -30 +31 @@
         pyplot.draw()
         #pyplot.plot( (0,1),(0,0), 'r' )
-        pyplot.plot(p2.x[v],p2.elev[v]) 
+        pyplot.plot(p2.x[v],p2.elev[1,v]) 
         pyplot.title(str(i)+'/200') # : velocity does not converge to zero' )
         pyplot.xlabel('x')
@@ -62 +63 @@
 w = D0 + 2*A*D0/(L**2)*numpy.cos(omega*p2.time)*( (p2.x[ppp]-4.*L/2.) -A/2.*numpy.cos(omega*p2.time))
 w2 = D0 + 2*A*D0/(L**2)*numpy.cos(omega*p2.time)*( (p2.x[ppp2]-4.*L/2.) -A/2.*numpy.cos(omega*p2.time))
-w2 = w2*(w2>p2.elev[ppp2])+p2.elev[ppp2]*(w2<=p2.elev[ppp2])
+w2 = w2*(w2>p2.elev[1,ppp2])+p2.elev[1,ppp2]*(w2<=p2.elev[1,ppp2])
 
 pyplot.clf()