Changeset 5587 for anuga_work/development/anuga_1d/dam_h_sudi.py

Timestamp:

Jul 30, 2008, 5:03:47 PM (15 years ago)

Author:

steve

Message:

Added in minimum height

File:

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

anuga_work/development/anuga_1d/dam_h_sudi.py

@@ -1 +1 @@
 import os
 from math import sqrt
-from shallow_water_h import *
+#from shallow_water_h import *
+from shallow_water_domain import *
 from Numeric import zeros, Float
 from analytic_dam_sudi import AnalyticDam
@@ -31 +32 @@
 
 L=2000.0
-N=100
+N=400
 
 cell_len=L/N
@@ -41 +42 @@
 domain=Domain(points)
 
-domain.default_order = 1
+domain.default_order = 2
 domain.default_time_order = 1
-#domain.cfl = 1.0
-#domain.limiter = "minmod"
+domain.cfl = 1.0
+domain.limiter = "vanleer"
 
 
 
-
-
-    
-    
 def height(x):
     y=zeros(len(x), Float)
@@ -63 +60 @@
     return y
 
-domain.set_quantity('height', height)
+domain.set_quantity('stage',height) #('height', height)
 domain.order=domain.default_order
 print "domain order", domain.order
@@ -78 +75 @@
 t0=time.time()
 yieldstep=30.0
-finaltime=30.0
-print "integral", domain.quantities['height'].get_integral()
+finaltime=20.0
+print "integral", domain.quantities['stage'].get_integral() #['height'].get_integral()
 for t in domain.evolve(yieldstep=yieldstep, finaltime=finaltime):
     domain.write_time()
-    print "integral", domain.quantities['height'].get_integral()
+    print "integral", domain.quantities['stage'].get_integral() #['height'].get_integral()
     if t>0.0:
-        HeightQ=domain.quantities['height'].vertex_values
+        HeightQ=domain.quantities['stage'].vertex_values #['height'].vertex_values
         MomentumQ=domain.quantities['xmomentum'].vertex_values
         h, uh=analytical_sol(X.flat, domain.time)