Changeset 2185

Timestamp:

Jan 4, 2006, 4:32:20 PM (19 years ago)

Author:

sexton

Message:

Update project file for Sydney tsunami scenario to include additional interior region (and potential to change clipping regions)

File:

• production/sydney_2006/project.py (modified) (2 diffs)

production/sydney_2006/project.py

@@ -51 +51 @@
 p3 = [north, west]
     
-polygon = [p0, p1, p2, p3]
+polygonall = [p0, p1, p2, p3]
 refzone = 56
 
-#Interior regions
+# setting up problem area for doing just around the harbour
+hsouth = degminsec2decimal_degrees(-33,54,0)
+hnorth = degminsec2decimal_degrees(-33,48,0)
+hwest = degminsec2decimal_degrees(151,0,0)
+heast = degminsec2decimal_degrees(151,30,0)
+
+hp0 = [hsouth, hwest]
+hp1 = [hsouth, heast]
+hp2 = [hnorth, heast]
+hp3 = [hnorth, hwest]
+polygon_h = [hp0, hp1, hp2, hp3]
+
+#Interior regions - the Harbour
 harbour_south = degminsec2decimal_degrees(-33,53,0)
 harbour_north = degminsec2decimal_degrees(-33,47,0)
 harbour_west = degminsec2decimal_degrees(151,5,0)
-harbour_east = degminsec2decimal_degrees(151,15,0)
+harbour_east = degminsec2decimal_degrees(151,19,0)
 
 k0 = [harbour_south, harbour_west]
@@ -67 +79 @@
 harbour_polygon = [k0, k1, k2, k3]
 
+# setting up problem area for doing just around Botany Bay
+bsouth = degminsec2decimal_degrees(-33,56,0)
+bnorth = degminsec2decimal_degrees(-34,3,0)
+bwest = degminsec2decimal_degrees(151,0,0)
+beast = degminsec2decimal_degrees(151,30,0)
 
+bp0 = [bsouth, bwest]
+bp1 = [bsouth, beast]
+bp2 = [bnorth, beast]
+bp3 = [bnorth, bwest]
+polygon_bb = [bp0, bp1, bp2, bp3]
+
+#Interior region - Botany Bay
+botanybay_south = degminsec2decimal_degrees(-33,58,0)
+botanybay_north = degminsec2decimal_degrees(-34,1,0)
+botanybay_west = degminsec2decimal_degrees(151,5,0)
+botanybay_east = degminsec2decimal_degrees(151,18,0)
+
+j0 = [botanybay_south, botanybay_west]
+j1 = [botanybay_south, botanybay_east]
+j2 = [botanybay_north, botanybay_east]
+j3 = [botanybay_north, botanybay_west]    
+
+botanybay_polygon = [j0, j1, j2, j3]
+