Changeset 4164

Timestamp:

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

Author:

sexton

Message:

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

Location:

anuga_work/production/dampier_2006

Files:

• project.py (modified) (1 diff)
• run_dampier.py (modified) (1 diff)

anuga_work/production/dampier_2006/project.py

@@ -171 +171 @@
 ###################################################################
 
+poly_bathy = read_polygon(polygons_dir+'polybathy.csv')
+
 # exporting asc grid - Dampier
 e_min_area = 474000

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'