Changeset 8893 for trunk/anuga_work/development/gareth/tests/shallow_steep_slope

Timestamp:

Jun 10, 2013, 11:11:02 AM (12 years ago)

Author:

davies

Message:

Updates to bal_dev

File:

• trunk/anuga_work/development/gareth/tests/shallow_steep_slope/channel_SU_sparse.py (modified) (4 diffs)

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

@@ -22 +22 @@
 from anuga.shallow_water.shallow_water_domain import Domain as Domain
 from balanced_dev import *
-from balanced_dev import Domain as Domain
+#from balanced_dev import Domain as Domain
 #from anuga_tsunami import *
 #from anuga_tsunami import Domain as Domain
@@ -36 +36 @@
 #------------------------------------------------------------------------------
 def topography(x, y):
-        return -x/10.  + 1.*(numpy.sin(x/10.) +abs(y-50.)/10.) -0.*(x>80.) # linear bed slope
+        return -x/10. # + 1.*(numpy.sin(x/10.) +abs(y-50.)/10.) -0.*(x>80.) # linear bed slope
 
 def stagetopo(x,y):
@@ -44 +44 @@
 
 line1=[ [9.,0.], [9., 100.] ]
-Qin=2.
+Qin=0.5
 Inlet_operator(domain, line1,Qin)
 
@@ -73 +73 @@
 # Evolve system through time
 #------------------------------------------------------------------------------
-for t in domain.evolve(yieldstep=4.0, finaltime=400.0):
+for t in domain.evolve(yieldstep=4.0, finaltime=800.0):
     print domain.timestepping_statistics()
-    print domain.boundary_flux_integral
+    #print domain.boundary_flux_integral
     print (domain.areas*(domain.quantities['stage'].centroid_values - domain.quantities['elevation'].centroid_values)).sum()
     #print domain.quantities['stage'].centroid_values[myindex] - domain.quantities['elevation'].centroid_values[myindex]