Changeset 5357 for anuga_work/production/perth/project.py

Timestamp:

May 23, 2008, 9:00:49 AM (17 years ago)

Author:

kristy

Message:

updated perth project.py. Added and changed internal polygons

File:

• anuga_work/production/perth/project.py (modified) (4 diffs)

anuga_work/production/perth/project.py

@@ -71 +71 @@
 
 # AHO + DPI data + colin French coastline
-coast_name = 'waterline'
+coast_name = 'coastline'
 offshore_name = 'perth_bathymetry'
-offshore1_name = 'missing_fairsheets'
+
 
 #final topo name
@@ -94 +94 @@
 coast_in_dir_name = topographies_in_dir + coast_name
 offshore_in_dir_name = topographies_in_dir + offshore_name
-offshore1_in_dir_name = topographies_in_dir + offshore1_name
+
 
 onshore_dir_name = topographies_dir + onshore_name
@@ -163 +163 @@
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
 
-poly_all = read_polygon(polygons_dir+'bounding_poly.csv')
+poly_all = read_polygon(polygons_dir+'bounding_area.csv')
 res_poly_all = 100000*res_factor
 
@@ -172 +172 @@
 ###############################
 
-#poly_pos20_neg20 = read_polygon(polygons_dir+'pos20_neg20pts.csv')
-poly_pos20_neg20 = read_polygon(polygons_dir+'pos20_neg20_new_pts.csv')
-res_pos20_neg20 = 20000*res_factor
-
-#poly_cbd = read_polygon(polygons_dir+'cbd_pts.csv')
-poly_cbd = read_polygon(polygons_dir+'cbd_smaller_pts.csv')
+#Polygon designed by 20m contours, or 3km from the coastline
+poly_internal_20_3 = read_polygon(polygons_dir+'internalh20mORd3km.csv')
+res_internal_20_3 = 20000*res_factor
+
+#Polygon designed to incorporate Garden Island and sand bank infront of Rockingham
+poly_garden_rockingham = read_polygon(polygons_dir+'garden_rockingham.csv')
+res_garden_rockingham = 1000*res_factor
+
+#Polygon designed to incorporate coastline of rottnest
+poly_rottnest = read_polygon(polygons_dir+'rottnest.csv')
+res_rottnest = 1000*res_factor
+
+#Polygon designed to incorporate perth and Fremantle CBD
+poly_cbd = read_polygon(polygons_dir+'CBD_coastal.csv')
 res_cbd = 500*res_factor
 
-poly_penguin = read_polygon(polygons_dir+'penguin_pts.csv')
+#Polygon designed to incorporate rockingham and penguin island
+poly_penguin = read_polygon(polygons_dir+'rockingham_penguin.csv')
 res_penguin = 500*res_factor
 
-poly_geordie_bay = read_polygon(polygons_dir+'geordie_bay_pts.csv')
+poly_geordie_bay = read_polygon(polygons_dir+'Geordie_Bay.csv')
 res_geordie_bay = 500*res_factor
 
-poly_sorrento_gauge = read_polygon(polygons_dir+'sorrento_gauge_pts.csv')
+poly_sorrento_gauge = read_polygon(polygons_dir+'Sorrento_gauge.csv')
 res_sorrento_gauge = 500*res_factor
 #assert zone == refzone
 
-interior_regions = [[poly_pos20_neg20,res_pos20_neg20],[poly_cbd,res_cbd]
+interior_regions = [[poly_internal_20_3,res_internal_20_3],[poly_cbd,res_cbd]
+                     ,[poly_garden_rockingham,res_garden_rockingham]
                      ,[poly_penguin,res_penguin],[poly_geordie_bay,res_geordie_bay]
-                     ,[poly_sorrento_gauge,res_sorrento_gauge]]
+                     ,[poly_sorrento_gauge,res_sorrento_gauge],[poly_rottnest, res_rottnest]]
 
 boundary_tags={'back': [4], 'side': [0,3],'ocean': [1, 2]}