Changeset 2403 for production/sydney_2006/project.py

Timestamp:

Feb 14, 2006, 3:29:11 PM (19 years ago)

Author:

sexton

Message:

Tidying up project script for Sydney scenario

File:

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

production/sydney_2006/project.py

@@ -2 +2 @@
 Also includes origin for slump scenario.
 """
-
 
 from os import sep
@@ -10 +9 @@
 #Making assumptions about the location of scenario data
 scenario_dir_name = 'sydney_tsunami_scenario_2006'
-# original 100m data
-#basename = 'sydney_2006_100' # get from Neil/Ingo (DEM or topo data)
 # revised 100m data
 coarsename = 'bathyland100' # get from Neil/Ingo (DEM or topo data)
 # revised 25m data
-#finename = 'bathyland25' # get from Neil/Ingo (DEM or topo data)
-#finename = 'bathyland100' # get from Neil/Ingo (DEM or topo data)
 finename = 'bathy_dem25' # get from Neil/Ingo (DEM or topo data) Wed 25 Jan
 
@@ -24 +19 @@
 northingmin = 6246250
 northingmax = 6264100
-#northingmax = 6280100
 
-# creating easting and northing max and min for export viz purposes
-#eminviz = 332000
-#emaxviz = 350000
-#nminviz = 6231000
-#nmaxviz = 6283000
-
-# version2
+# creating easting and northing max and min for export viz purposes 
 eminviz = 318000
 emaxviz = 351000
 nminviz = 6231000
 nmaxviz = 6283000
-# for harbour region
-#nminviz = 6248000
-#nmaxvix = 6270000
-# for botany bay region
-#nmaxviz = 6248000
-#nminvix = 6283000
 
-basename = 'slump_ole'
+basename = 'slump'
 
 if sys.platform == 'win32':
@@ -56 +38 @@
 datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
 outputdir = home+sep+scenario_dir_name+sep+'output'+sep
-#boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep
 
 meshname = meshdir + basename
@@ -62 +43 @@
 finedemname = datadir + finename
 combineddemname = datadir + 'sydneytopo'
-#boundaryname = boundarydir + boundary_basename
 outputname = outputdir + basename  #Used by post processing
-#newoutputname = outputdir+'Combined Data v1'+sep+'d400'+sep+basename
 gauge_filename = outputdir + 'sydney_gauges.xya'
 gauge_outname = outputdir + 'gauges_max_output.xya'
 
 #Georeferencing
-#from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees
 from coordinate_transforms.redfearn import degminsec2decimal_degrees
 
-#Origin of existing dem (FIXME: Temporary measure)
-#x0_origin = 314036 #original 100m data
-#y0_origin = 6224951
-#refzone = 56
-#x0_origin = 314036.58727982 #revised 100m and 25m data
-#y0_origin = 6224951.2960092
-#mesh_origin = (refzone, x0_origin, y0_origin)  # input from Neil's data
-
 # define clipping polygon
-south = degminsec2decimal_degrees(-34,05,0)
-north = degminsec2decimal_degrees(-33,33,0)
-west = degminsec2decimal_degrees(151,1,0)
-east = degminsec2decimal_degrees(151,30,0)
-p0 = [south, west]
-p1 = [south, east]
-p2 = [north, east]
-p3 = [north, west]
-    
-polygonall = [p0, p1, p2, p3]
-
 dsouth = degminsec2decimal_degrees(-34,05,0)
 dnorth = degminsec2decimal_degrees(-33,33,0)
@@ -113 +72 @@
 dp8 = [dnorth, dwest]
     
-#diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp5, dp6, dp7]
 diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp5, dp6, dp7]
-# to put chunk back in
-#diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp8] 
 
-#Interior regions - the Harbour - take 2
+#Interior regions - the Harbour
 harbour_1x = degminsec2decimal_degrees(-33,51,0)
 harbour_1y = degminsec2decimal_degrees(151,2,30)
@@ -129 +85 @@
 harbour_3x = degminsec2decimal_degrees(-33,47,0)
 harbour_3y = degminsec2decimal_degrees(151,20,30)
-#harbour_4x = degminsec2decimal_degrees(-33,48,0)
 harbour_4x = degminsec2decimal_degrees(-33,47,50)
 harbour_4y = degminsec2decimal_degrees(151,8,10)
-#harbour_5x = degminsec2decimal_degrees(-33,49,0)
 harbour_5x = degminsec2decimal_degrees(-33,48,10)
 harbour_5y = degminsec2decimal_degrees(151,8,0)
@@ -139 +93 @@
 harbour_7x = degminsec2decimal_degrees(-33,34,30)
 harbour_7y = degminsec2decimal_degrees(151,20,20)
-#harbour_8x = degminsec2decimal_degrees(-33,33,30)
-#harbour_8y = degminsec2decimal_degrees(151,15,0)
 harbour_8x = degminsec2decimal_degrees(-33,33,30)
 harbour_8y = degminsec2decimal_degrees(151,17,0)
-#harbour_9x = degminsec2decimal_degrees(-33,45,30)
-#harbour_9y = degminsec2decimal_degrees(151,15,0)
 harbour_9x = degminsec2decimal_degrees(-33,45,30)
 harbour_9y = degminsec2decimal_degrees(151,17,0)
@@ -172 +122 @@
 k142 = [harbour_15x, harbour_15y]
 
-#harbour_polygon_2 = [k02, k112, k122, k12, k22, k62, k72, k82, k92, k102, k32, k42, k52]
-#harbour_polygon_2 = [k122, k12, k22, k62, k72, k102] #worked
 harbour_polygon_2 = [k02, k112, k122, k12, k22, k62, k72, k82, k102, k42, k52] #worked
 
-#Interior region - Botany Bay - take 2
+#Interior region - Botany Bay
 bb_1x = degminsec2decimal_degrees(-34,3,0)
 bb_1y = degminsec2decimal_degrees(151,2,30)
@@ -210 +158 @@
 
 botanybay_polygon_2 = [j92, j12, j22, j62, j82, j72, j42] # worked
-#botanybay_polygon_2 = [j02, j12, j22, j32, j42, j52] # worked
 
-# this section sets up source origin for slump scenario
-# close to harbour opening (343000,6255000)
-# x0 = 28964
-# y0 = 30049
-# around 10km harbour opening (353000,6255000)
-#x0 = 38964
-#y0 = 30049
-# around 18km across from harbour(361000,6255000) 
-#x0 = 46964
-#y0 = 30049
 # around 42km across from harbour(385000,6255000) 
-x0 = 28964 + 42000
-y0 = 30049
-
-#revised 100m and 25m data
-slump_origin = [x0+314036.58727982, y0+6224951.2960092] #Absolute UTM
-
-
-# close to botany bay opening (340000,6236000)
-# x0 = 25964
-# y0 = 11049
-# around 10km from botany bay opening (350000,6236000)
-# x0 = 35964
-# y0 = 11049
-# around 21km from botany bay opening (361000,6236000)
-#x0 = 46964
-#y0 = 11049
-
-# not used for sydney scenario, original interior regions listed though
-# 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 - take 1
-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,19,0)
-
-#harbour_south1 = degminsec2decimal_degrees(-33,53,0)
-#harbour_south2 = degminsec2decimal_degrees(-33,52,0)
-#harbour_north1 = degminsec2decimal_degrees(-33,45,0)
-#harbour_north2 = degminsec2decimal_degrees(-33,48,0)
-#harbour_west = degminsec2decimal_degrees(151,5,0)
-#harbour_east = degminsec2decimal_degrees(151,19,0)
-
-k0 = [harbour_south, harbour_west]
-k1 = [harbour_south, harbour_east]
-k2 = [harbour_north, harbour_east]
-k3 = [harbour_north, harbour_west]    
-
-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 - take 1
-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]
+#x0 = 28964 + 42000 
+#y0 = 30049 
+#slump_origin = [x0+314036.58727982, y0+6224951.2960092] #Absolute UTM
+slump_origin = [385000.0, 6255000.0] #Absolute UTM