Changeset 9085 for trunk/anuga_core/source/anuga_1d/sqpipe/culvert.py

Timestamp:

Apr 11, 2014, 10:06:44 AM (11 years ago)

Author:

steve

Message:

Adding anuga_1d to test all

File:

• trunk/anuga_core/source/anuga_1d/sqpipe/culvert.py (modified) (4 diffs)

trunk/anuga_core/source/anuga_1d/sqpipe/culvert.py

@@ -36 +36 @@
 def top(x):
 
-    t = numpy.ones_like(x)
+    t = numpy.ones_like(x)*4
 
     t = numpy.where(x<-40, 20, t)
@@ -45 +45 @@
     return height(x)*width(x)
 
+def friction(x):
+    return numpy.ones_like(x)*0.01
+
 
 #===============================================================================
 
 def get_domain():
-    N = 100
+    N = 20
     print "Evaluating domain with %d cells" %N
 
@@ -57 +60 @@
 
     domain.set_spatial_order(2)
-    domain.set_timestepping_method('euler')
+    domain.set_timestepping_method('rk2')
     domain.set_CFL(0.5)
     domain.set_limiter("vanleer")
@@ -68 +71 @@
     domain.set_quantity('width',width)
     domain.set_quantity('top',top)
+    domain.set_quantity('friction',friction)
 
     Br = dom.Reflective_boundary(domain)
     Bt = dom.Transmissive_boundary(domain)
 
-    domain.set_boundary({'left': Br, 'right' : Bt})
+    domain.set_boundary({'left': Br, 'right' : Br})
 
     return domain