Changeset 4081

Timestamp:

Dec 14, 2006, 3:57:12 PM (18 years ago)

Author:

nick

Message:

change to perth

Location:

anuga_work/production/perth_2006

Files:

• build_perth.py (modified) (2 diffs)
• project.py (modified) (3 diffs)

anuga_work/production/perth_2006/build_perth.py

@@ -57 +57 @@
 print"project.bounding_polygon",project.bounding_polygon
 print"project.combined_dir_name",project.combined_dir_name
-
+'''
 # topography directory filenames
 onshore_in_dir_name = project.onshore_in_dir_name
@@ -120 +120 @@
 G_clipped.export_points_file(project.combined_dir_name + '.xya')
 
+'''
+
+G_all=Geospatial_data(file_name = project.combined_dir_name + '.xya')
+print'split'
+G_all_1, G_all_2 = G_all.split(.20)
+print'export 1'
+G_all_1.export_points_file(project.combined_dir_name+'_small' + '.xya')
+print'export 2'
+G_all_2.export_points_file(project.combined_dir_name+'_other' + '.xya')
 
 '''

anuga_work/production/perth_2006/project.py

@@ -27 +27 @@
 run_time = time+'_run'
 print 'gtime: ', gtime
+
+tide = 0.6
 
 #Making assumptions about the location of scenario data
@@ -124 +126 @@
 refzone = 50 
 
-print "bounding_polygon", bounding_polygon
-print 'poly area', polygon_area(bounding_polygon)/1000000.0
-
 
 #Interior regions
@@ -139 +138 @@
 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