Changeset 4046 for anuga_work/production/sydney_2006/project_slide.py

Timestamp:

Dec 1, 2006, 2:46:11 PM (17 years ago)

Author:

sexton

Message:

slide updates

File:

• anuga_work/production/sydney_2006/project_slide.py (modified) (6 diffs)

anuga_work/production/sydney_2006/project_slide.py

@@ -66 +66 @@
 
 # bounding polygon for study area
-polyAll = read_polygon(polygondir+'extent.csv')
+polyAll = read_polygon(polygondir+'extent2.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
 
@@ -77 +75 @@
 
 # sydney digitized polygons
-interior_regions = []
-area = 0.
-coast_res = 500.
-for i in range(1,10):
-    poly = read_polygon(polygondir+'poly'+str(i)+'.csv')
-    area += polygon_area(poly)
-    #interior_regions.append([poly,coast_res])
-
+poly_coast = read_polygon(polygondir+'coast.csv')
 poly_syd1 = read_polygon(polygondir+'local.csv')
-local_res = 50000.
-interior_regions.append([poly_syd1,local_res])
 
 print 'Area of local polygon', polygon_area(poly_syd1)/1000000.0
-print 'Area of local polygon', area/1000000.0
+print 'Area of coastal polygon', polygon_area(poly_coast)/1000000.0
 
 ###################################################################
@@ -129 +118 @@
 ###################################################################
 
-#slide_origin = [,]
+slide_origin_a = [352928.37,6188303.94]
+slide_origin_b = [393301.80,6253255.08]
+slide_origin_c = [446025.35,6330307.46]
+
 bulli_depth = 2087.0
 bulli_length = 16840.0
@@ -135 +127 @@
 bulli_width = 8860.0
 bulli_density = 1.46
+bulli_slope = 4.0
 bulli_alpha = 126.0 + 90.0
 
@@ -142 +135 @@
 shovel_width = 4350.0
 shovel_density = 1.49
+shove_slope = 4.0
 shovel_alpha = 118.0 + 90.0
 
@@ -149 +143 @@
 yacaaba_width = 3080.0
 yacaaba_density = 1.48
+yacaaba_slope = 3.0
 yacaaba_alpha = 133.0 + 90.0