Changeset 4399

Timestamp:

Apr 18, 2007, 1:19:31 PM (16 years ago)

Author:

nick

Message:

updates to exmouth

Location:

anuga_work/production/exmouth_2006

Files:

• project.py (modified) (5 diffs)
• run_boundary_points.py (added)

anuga_work/production/exmouth_2006/project.py

@@ -27 +27 @@
 
 # onshore data provided by WA DLI
-#onshore_name =  # original
+onshore_name =  'first' # original
 
 # AHO + DPI data
@@ -130 +130 @@
 
 # bounding polygon for study area
-polyAll = read_polygon(polygondir+'extent.csv')
+poly_all = read_polygon(polygondir+'extent_zone50a.csv')
 
 # plot bounding polygon and make sure BC info surrounds it
-#plot_polygons([polyAll, poly_bc],'boundingpoly',verbose=False)
-print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
+#plot_polygons([poly_all, poly_bc],'boundingpoly',verbose=False)
+print 'Area of bounding polygon', polygon_area(poly_all)/1000000.0
 
 ###################################################################
@@ -149 +149 @@
 # Interior region definitions
 ###############################
-
+'''
 # exmouth digitized polygons
 poly_exmouth1 = read_polygon(polygondir+'exmouth_Local_Polygon_update.csv')
@@ -156 +156 @@
 #poly_exmouth4 = read_polygon(polygondir+'Cable_Beach_revised.csv')
 
-plot_polygons([polyAll,poly_exmouth1,poly_exmouth2,poly_exmouth3],'boundingpoly2',verbose=False)
+plot_polygons([poly_all,poly_exmouth1,poly_exmouth2,poly_exmouth3],'boundingpoly2',verbose=False)
 print 'Area of local polygon', polygon_area(poly_exmouth1)/1000000.0
 print 'Area of close polygon', polygon_area(poly_exmouth2)/1000000.0
@@ -179 +179 @@
     no_triangles += area/remainder_res
     return int(no_triangles/0.7)
+'''