Changeset 8884 for trunk/anuga_work/development/gareth/tests

Timestamp:

Jun 7, 2013, 3:49:59 PM (12 years ago)

Author:

davies

Message:

Updates to bal_dev

Location:

trunk/anuga_work/development/gareth/tests

Files:

• runup_sinusoid/runup_sinusoid.py (modified) (4 diffs)
• shallow_steep_slope/channel_SU_sparse.py (modified) (1 diff)

trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py

@@ -14 +14 @@
 #from swb2_domain import *
 #from balanced_basic import *
-#from balanced_dev import *
+from balanced_dev import *
 #from anuga_tsunami import *
 
@@ -22 +22 @@
 points, vertices, boundary = anuga.rectangular_cross(20,20, len1=1., len2=1.)
 
-domain=anuga.Domain(points,vertices,boundary)    # Create Domain
+domain=Domain(points,vertices,boundary)    # Create Domain
 domain.set_name('runup_sinusoid_v2')                         # Output to file runup.sww
 #domain.set_timestepping_method('euler')
-domain.set_flow_algorithm('tsunami')
+#domain.set_flow_algorithm('tsunami')
 #------------------
 # Define topography
@@ -32 +32 @@
 #### Pathological
 #scale_me=100.0
-#boundary_ws=-0.2
+#boundary_ws=-0.2#1999
 #init_ws=-0.2
 #bumpiness=50. # Higher = shorter wavelength oscillations in topography
-#tstep=0.002
-#lasttime=1.1
+#tstep=0.2
+#lasttime=20.1
 
 ### Sensible
@@ -60 +60 @@
 domain.set_quantity('friction',0.00)             # Constant friction
 
+print '#-#', domain.minimum_allowed_height
 #def frict_change(x,y):
 #       return 0.2*(x>0.5)+0.1*(x<=0.5)

trunk/anuga_work/development/gareth/tests/shallow_steep_slope/channel_SU_sparse.py

@@ -43 +43 @@
     return stg#*(stg>topo) + topo*(stg<=topo)
 
-line1=[ [10.,0.], [10., 100.] ]
-Qin=20.
+line1=[ [9.,0.], [9., 100.] ]
+Qin=2.
 Inlet_operator(domain, line1,Qin)