Changeset 8166

Timestamp:

Mar 25, 2011, 9:13:31 AM (14 years ago)

Author:

steve

Message:

committing some small changes to kinematic viscosity

Location:

trunk/anuga_core/source/anuga/operators

Files:

• kinematic_viscosity_operator.py (modified) (1 diff)
• run_dam_break.py (modified) (3 diffs)

trunk/anuga_core/source/anuga/operators/kinematic_viscosity_operator.py

@@ -58 +58 @@
 
         self.dt = 0.0 #Need to set to domain.timestep
-        self.dt_apply = 1.0
+        self.dt_apply = 0.0
 
         self.boundary_len = len(self.domain.boundary)

trunk/anuga_core/source/anuga/operators/run_dam_break.py

@@ -59 +59 @@
 domain.set_use_kinematic_viscosity(True)
 
-
+domain.set_beta(1.7)
 #------------------------------------------------------------------------------
 # Setup initial conditions
 #------------------------------------------------------------------------------
-domain.set_quantity('elevation',0.0)
-domain.set_quantity('friction', 0.0)
-
 h0 = 10.0
 h1 = 0.0
+
+def elevation(x,y):
+
+    return where(x<=95000.0, 0.0, h0)
 
 def height(x,y):
@@ -73 +74 @@
     return where(x<=50000.0, h0, h1)
 
-
-domain.set_quantity('stage', height)
+domain.set_quantity('elevation',elevation)
+domain.set_quantity('friction', 0.0)
+domain.add_quantity('stage', height)
 #-----------------------------------------------------------------------------
 # Setup boundary conditions
@@ -85 +87 @@
 
 # Associate boundary tags with boundary objects
-domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br})
+domain.set_boundary({'left': Br, 'right': Bt, 'top': Br, 'bottom': Br})