Changeset 5728 for anuga_work/development/anuga_1d/dam_h_elevation.py

Timestamp:

Sep 4, 2008, 7:12:44 AM (17 years ago)

Author:

steve

Message:

 

File:

• anuga_work/development/anuga_1d/dam_h_elevation.py (modified) (5 diffs)

anuga_work/development/anuga_1d/dam_h_elevation.py

@@ -1 +1 @@
 import os
 from math import sqrt, sin, cos, pi, exp
-from shallow_water_domain import *
+from shallow_water_domain_new import *
 from Numeric import zeros, Float
 #from analytic_dam_sudi import AnalyticDam
@@ -138 +138 @@
         #pass
     
-        from pylab import plot,title,xlabel,ylabel,legend,savefig,show,hold,subplot,ion
+        from pylab import clf,plot,title,xlabel,ylabel,legend,savefig,show,hold,subplot,ion
         #print 'Test1'
-        hold(False)
+        #hold(False)
+        clf()
         #print 'test 2'
         plot1 = subplot(211)
@@ -150 +151 @@
         xlabel('Position')
         ylabel('Stage')
-        legend( ('Bed Elevation', 'Numerical Solution'),
-           'upper right', shadow=False)
+        #legend( ('Bed Elevation', 'Numerical Solution'),
+        #   'upper right', shadow=False)
+
+        
         #legend(('Analytical Solution', 'Numerical Solution'),
         #   'lower right', shadow=False)
@@ -157 +160 @@
         plot(X,MomentumQ)
         #plot2.set_ylim([-5,35])
-        legend( ('Numerical Solution', 'for momentum'),
-           'upper right', shadow=False)
+        #legend( ('Numerical Solution', 'for momentum'),
+        #   'upper right', shadow=False)
     
         xlabel('Position')
@@ -168 +171 @@
         #savefig(file)
         show()
-
+        
 print 'That took %.2f seconds'%(time.time()-t0)