Changeset 4058 for anuga_work/production/wollongong_2006/project_slide.py

Timestamp:

Dec 5, 2006, 5:30:41 PM (18 years ago)

Author:

sexton

Message:

updates for all slide scenarios

File:

• anuga_work/production/wollongong_2006/project_slide.py (modified) (8 diffs)

anuga_work/production/wollongong_2006/project_slide.py

@@ -27 +27 @@
 # onshore data provided by NSW LPI, Krishna merged with existing
 # PMD data to make 100m DEM
-on_offshore_name = '' 
+on_offshore10_name = 'wollongong_10'
+nsw100 = 'dem_utm'
 
-# AHO
-offshore_name1 = 'XY100011610'
+# survey data
+offshore_name1 = 'surveyAreaA'
+
+# AHO data
+offshore_name4 = '1000003611export'
+offshore_name5 = '1000003613export'
+offshore_name6 = '1000003614export'
+offshore_name7 = '1000003627export'
+offshore_name8 = '1000003628export'
+offshore_name9 = 'AHDexport'
 
 #swollen/ all data output
@@ -55 +64 @@
 on_offshore_dem_name = datadir + on_offshore_name
 offshore_dem_name1 = datadir + offshore_name1
+offshore_dem_name4 = datadir + offshore_name4
+offshore_dem_name5 = datadir + offshore_name5
+offshore_dem_name6 = datadir + offshore_name6
+offshore_dem_name7 = datadir + offshore_name7
+offshore_dem_name8 = datadir + offshore_name8
+offshore_dem_name9 = datadir + offshore_name9
 combined_dem_name   = datadir + 'gong_combined_elevation'
 
@@ -62 +77 @@
 
 # bounding polygon for study area
-polyAll = read_polygon(polygondir+'extentA.csv')
+polyAll = read_polygon(polygondir+'extentA_up.csv')
 
 print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
@@ -70 +85 @@
 ###############################
 
-# broome digitized polygons
-poly_gong1 = read_polygon(polygondir+'.csv')
-poly_gong2 = read_polygon(polygondir+'.csv')
-poly_gong3 = read_polygon(polygondir+'.csv')
+# areaA digitized polygons
+poly_local = read_polygon(polygondir+'local.csv')
+poly_gong = read_polygon(polygondir+'gong.csv')
 
-print 'Area of local polygon', polygon_area(poly_gong1)/1000000.0
-print 'Area of close polygon', polygon_area(poly_gong2)/1000000.0
-print 'Area of coastal polygon', polygon_area(poly_gong3)/1000000.0
+print 'Area of local polygon', polygon_area(poly_local)/1000000.0
 
 ###################################################################
@@ -83 +95 @@
 ###################################################################
 
+poly_10mclip = read_polygon(polygondir+'10mclip.csv')
+poly_surveyclip = read_polygon(polygondir+'surveyclip.csv')
+
+# clipping for nsw grid for conversion
+eastingmin_nsw = 272900
+eastingmax_nsw = 463550
+northingmin_nsw = 6112950
+northingmax_nsw = 6269750
+
 # exporting asc grid
-eastingmin = 
-eastingmax = 
-northingmin = 
-northingmax = 
+#eastingmin = 
+#eastingmax = 
+#northingmin = 
+#northingmax = 
 
 ###################################################################
@@ -103 +124 @@
 bulli_density = 1.46
 bulli_slope = 4.0
-bulli_alpha = 126.0 + 90.0
+bulli_alpha = 126.0 - 90.0
 
 shovel_depth = 968.0
@@ -111 +132 @@
 shovel_density = 1.49
 shove_slope = 4.0
-shovel_alpha = 118.0 + 90.0
+shovel_alpha = 118.0 - 90.0
 
 yacaaba_depth = 1119.0
@@ -119 +140 @@
 yacaaba_density = 1.48
 yacaaba_slope = 3.0
-yacaaba_alpha = 133.0 + 90.0
+yacaaba_alpha = 133.0 - 90.0