Changeset 3483

Timestamp:

Aug 11, 2006, 12:55:08 PM (19 years ago)

Author:

ole

Message:

Fixed topography for harbour example.

File:

• documentation/user_manual/examples/harbour.py (modified) (1 diff)

documentation/user_manual/examples/harbour.py

@@ -66 +66 @@
 def topography(x,y):
     from Numeric import zeros, Float
-    z = zeros(len(x),len(y),Float)
+    z = zeros(len(x), Float)
+    
+    xcentre = (west+east)/2.
     ycentre = (south+north)/2.
     ycentreup = ycentre+harbour_width/2.
     ycentredown = ycentre-harbour_width/2.
     for i, xi in enumerate(x):
-        for j, yi in enumerate(y):
-           if xi >= centre: z[i,j] = shelf_depth
-           if xi > west and xi <= centre-harbour_length: 
-                if yi < ycentredown: z[i,j] = 0.0
-                if yi > ycentreup: z[i,j] = 0.0
-                if yi > ycentredown and yi < ycentreup: z[i,j] = harbour_depth
-           if xi > centre-harbour_length and xi < centre: z[i,j] = harbour_depth
+        yi = y[i]
+        if xi >= xcentre: z[i] = shelf_depth
+        if xi > west and xi <= xcentre-harbour_length: 
+            if yi < ycentredown: z[i] = 0.0
+            if yi > ycentreup: z[i] = 0.0
+            if yi > ycentredown and yi < ycentreup: z[i] = harbour_depth
+        if xi > xcentre-harbour_length and xi < xcentre: z[i] = harbour_depth
+        
     return z