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

Timestamp:

Sep 4, 2008, 8:28:49 PM (16 years ago)

Author:

steve

Message:

Unit test for second order timestepping

File:

• anuga_work/development/anuga_1d/dam_h_elevation.py (modified) (4 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_new import *
+from shallow_water_domain import *
 from Numeric import zeros, Float
 #from analytic_dam_sudi import AnalyticDam
@@ -41 +41 @@
 domain=Domain(points)
 
-domain.default_order = 2
+domain.default_order = 1
 domain.default_time_order = 2
 domain.cfl = 1.0
@@ -140 +140 @@
         from pylab import clf,plot,title,xlabel,ylabel,legend,savefig,show,hold,subplot,ion
         #print 'Test1'
-        #hold(False)
-        clf()
+        hold(False)
+        #clf()
         #print 'test 2'
         plot1 = subplot(211)
         #print 'test 3'
     
+        
+        #print 'Test4'
+
         plot(X,ElevationQ,X,HeightQ)
-        #print 'Test4'
         plot1.set_ylim([-1,12])
         xlabel('Position')
@@ -159 +161 @@
         plot2 = subplot(212)
         plot(X,MomentumQ)
-        #plot2.set_ylim([-5,35])
+        plot2.set_ylim([-1,1])
+        
         #legend( ('Numerical Solution', 'for momentum'),
         #   'upper right', shadow=False)