Changeset 3753 for anuga_core/documentation/user_manual/examples/channel_3.py

Timestamp:

Oct 11, 2006, 3:58:53 PM (18 years ago)

Author:

ole

Message:

Updates

File:

• anuga_core/documentation/user_manual/examples/channel_3.py (modified) (5 diffs)

anuga_core/documentation/user_manual/examples/channel_3.py

@@ -1 +1 @@
 """Simple water flow example using ANUGA
 
-Water flowing down a channel
+Water flowing down a channel with more complex topography
 """
 
@@ -20 +20 @@
 width = 5.
 dx = dy = 1           # Resolution: Length of subdivisions on both axes
+#dx = dy = .1           # Resolution: Length of subdivisions on both axes
 
-points, vertices, boundary = rectangular_cross(int(length/dx), int(width/dy), len1=length, len2=width)
+points, vertices, boundary = rectangular_cross(int(length/dx), int(width/dy),
+                                               len1=length, len2=width)
 domain = Domain(points, vertices, boundary)   
 domain.set_name('channel_3')                  # Output name
@@ -30 +32 @@
 #------------------------------------------------------------------------------
 def topography(x,y):
-    z = -x/10                             # linear bed slope
+    """Complex topography defined by a function of vectors x and y
+    """
+    
+    z = -x/10                                
 
     N = len(x)
@@ -50 +55 @@
 
 
-domain.set_quantity('elevation', topography)            # Use function for elevation
-domain.set_quantity('friction', 0.01)                   # Constant friction 
-domain.set_quantity('stage', expression='elevation')    # Dry
+domain.set_quantity('elevation', topography)  # Use function for elevation
+domain.set_quantity('friction', 0.01)         # Constant friction 
+domain.set_quantity('stage',
+                    expression='elevation')   # Dry initial condition
 
 
@@ -58 +64 @@
 # Setup boundary conditions
 #------------------------------------------------------------------------------
-Bi = Dirichlet_boundary([0.4, 0, 0])                            # Inflow
-Br = Reflective_boundary(domain)                                # Solid reflective wall
-Bo = Dirichlet_boundary([-5, 0, 0])                             # Outflow
+Bi = Dirichlet_boundary([0.4, 0, 0])          # Inflow
+Br = Reflective_boundary(domain)              # Solid reflective wall
+Bo = Dirichlet_boundary([-5, 0, 0])           # Outflow
 
 domain.set_boundary({'left': Bi, 'right': Bo, 'top': Br, 'bottom': Br})