Changeset 4164 for anuga_work/production/dampier_2006/run_dampier.py

Timestamp:

Jan 10, 2007, 4:47:39 PM (17 years ago)

Author:

sexton

Message:

update Dampier script to incorporate setting stage to 0 onshore and tide offshore

File:

• anuga_work/production/dampier_2006/run_dampier.py (modified) (1 diff)

anuga_work/production/dampier_2006/run_dampier.py

@@ -129 +129 @@
 
     print 'Setup initial conditions'
-
-    domain.set_quantity('stage', tide)
+    
+    def IC(x,y):
+        
+        from anuga.utilities.polygon import is_inside_polygon
+        from Numeric import zeros, Float
+
+        w = zeros(len(x), Float)
+        for i, xi in enumerate(x):
+            yi = y[i]
+            if is_inside_polygon([xi,yi],project.poly_bathy,verbose=False) == True:
+                w[i] = tide
+            else:
+                w[i] = 0.
+
+        return w
+
+    domain.set_quantity('stage', IC)
     domain.set_quantity('friction', 0.01) 
     print 'Start Set quantity'