Changeset 4091 for anuga_work/production/perth_2006/project.py

Timestamp:

Dec 19, 2006, 9:19:27 AM (17 years ago)

Author:

nick

Message:

updates to dampier and perth

File:

• anuga_work/production/perth_2006/project.py (modified) (9 diffs)

anuga_work/production/perth_2006/project.py

@@ -46 +46 @@
 coast_name = 'waterline'
 offshore_name = 'perth_bathymetry'
+offshore1_name = 'missing_fairsheets'
 
 #final topo name
 combined_name ='perth_combined_elevation'
+combined_smaller_name = 'perth_combined_elevation_smaller'
+
 
 topographies_in_dir = home+sep+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
@@ -63 +66 @@
 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
@@ -76 +80 @@
 combined_dir_name = topographies_dir + combined_name
 combined_time_dir_name = topographies_time_dir + combined_name
+combined_smaller_name_dir = topographies_dir + combined_smaller_name
 #combined_time_dir_final_name = topographies_time_dir + combined_final_name
 
 
-#Derive subdirectories and filenames
 
 meshes_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'meshes'+sep
@@ -88 +92 @@
 
 
-boundaries_source = 'mag_9_corrected'
+boundaries_source = '????'
 #boundaries locations
 boundaries_in_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep+boundaries_source+sep
@@ -96 +100 @@
 #boundaries_time_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+build_time+sep
 #boundaries_time_dir_name = boundaries_time_dir + boundaries_name  #Used by post processing
-#ideas
-#boundaries_time_dir = boundaries_in_dir+'urs'+sep+boundaries_source+sep
 
 #output locations
@@ -111 +113 @@
 
 
-
-
-
-
-
-
 ###############################
 # Domain definitions
@@ -123 +119 @@
 
 bounding_polygon = read_polygon(polygons_dir+'bounding_poly.csv')
+res_bounding = 1000000
 
 refzone = 50 
 
 
-#Interior regions
+###############################
+# Interior region definitions
+###############################
 
-#cipma_polygon, zone = convert_from_latlon_to_utm([k0, k1, k2, k3])
+poly_pos20_neg20 = read_polygon(polygons_dir+'pos20_neg20pts.csv')
+res_pos20_neg20 = 500000
+
+poly_cbd = read_polygon(polygons_dir+'cbd_pts.csv')
+res_cbd = 50000
+
+poly_penguin = read_polygon(polygons_dir+'penguin_pts.csv')
+res_penguin = 50000
 #assert zone == refzone
 
-poly_pos20_neg20 = read_polygon(polygons_dir+'pos20_neg20pts.csv')
+trigs_bound = polygon_area(bounding_polygon)/res_bounding
+trigs_pos = polygon_area(poly_pos20_neg20)/res_pos20_neg20
+trigs_cbd = polygon_area(poly_cbd)/res_cbd
+trigs_penguin = polygon_area(poly_penguin)/res_penguin
+trigs_min = trigs_bound + trigs_pos + trigs_cbd + trigs_penguin
 
-poly_cbd = read_polygon(polygons_dir+'cbd_pts.csv')
-
-poly_penguin = read_polygon(polygons_dir+'penguin_pts.csv')
-#assert zone == refzone
-
-print "bounding_polygon", bounding_polygon
-print 'Area of bounding poly', polygon_area(bounding_polygon)/1000000.0
-print 'Area of pos20_neg20pts', polygon_area(poly_pos20_neg20)/1000000.0
-print 'Area of poly_cbd', polygon_area(poly_cbd)/1000000.0
-print 'Area of poly_penguin', polygon_area(poly_penguin)/1000000.0
+print 'Area of bounding poly', trigs_bound
+print 'Area of pos20_neg20pts', trigs_pos
+print 'Area of poly_cbd', trigs_cbd
+print 'Area of poly_penguin', trigs_penguin
+print 'min number triangles', trigs_min
 
 
@@ -149 +154 @@
 # Clipping regions for export to asc and regions for clipping data
 ###################################################################
-'''
+
 # exporting asc grid
 eastingmin = 406215.87
@@ -156 +161 @@
 northingmax = 8032834.52
 
-###############################
-# Interior region definitions
-###############################
 
-# perth digitized polygons
-poly_perth1 = read_polygon(polygondir+'perth_Local_Polygon_update.csv')
-poly_perth2 = read_polygon(polygondir+'perth_Close2_update.csv')
-poly_perth3 = read_polygon(polygondir+'perth_Coast_update.csv')
-#poly_perth4 = read_polygon(polygondir+'Cable_Beach_revised.csv')
 
-plot_polygons([polyAll,poly_perth1,poly_perth2,poly_perth3],'boundingpoly2',verbose=False)
-print 'Area of local polygon', polygon_area(poly_perth1)/1000000.0
-print 'Area of close polygon', polygon_area(poly_perth2)/1000000.0
-print 'Area of coastal polygon', polygon_area(poly_perth3)/1000000.0
-#print 'Area of cable beach polygon', polygon_area(poly_perth4)/1000000.0
-'''
-