Changeset 4177

Timestamp:

Jan 15, 2007, 2:58:45 PM (18 years ago)

Author:

sexton

Message:

update dampier scripts based on delivery of data

Location:

anuga_work/production/dampier_2006

Files:

• build_dampier.py (modified) (2 diffs)
• project.py (modified) (6 diffs)

anuga_work/production/dampier_2006/build_dampier.py

@@ -55 +55 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
-print"project.bounding_polygon",project.bounding_polygon
+print"project.poly_all",project.poly_all
 print"project.combined_dir_name",project.combined_dir_name
 
@@ -96 +96 @@
 
 print'create Geospatial onshore objects from topographies'
-G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
-print'create Geospatial coastal objects from topographies'
-G2 = Geospatial_data(file_name = coast_in_dir_name + '.xya')
-print'create Geospatial nickel bay objects from topographies'
-G3 = Geospatial_data(file_name = offshore1_dir_name + '.pts')
-print'create Geospatial clipped TIN offshore objects from topographies'
-G4 = Geospatial_data(file_name = offshore2_dir_name + '.pts')
-print'create Geospatial offshore objects from topographies'
-G_off = Geospatial_data(file_name = offshore_in_dir_name + '.xya')
+#print'create Geospatial coastal objects from topographies'
+#print'create Geospatial nickel bay objects from topographies'
+#print'create Geospatial clipped TIN offshore objects from topographies'
+#print'create Geospatial offshore objects from topographies'
+G = Geospatial_data(file_name = onshore_dir_name + '.pts') +\
+     Geospatial_data(file_name = coast_in_dir_name + '.xya') +\
+     Geospatial_data(file_name = offshore1_dir_name + '.pts') +\
+     Geospatial_data(file_name = offshore2_dir_name + '.pts') +\
+     Geospatial_data(file_name = offshore_in_dir_name + '.xya')
 
 print'add all geospatial objects'
-G = G1 + G2 + G3 + G4 + G_off
+#G = G1 + G2 + G3 + G4 + G_off
 
 print'clip combined geospatial object by bounding polygon'

anuga_work/production/dampier_2006/project.py

@@ -7 +7 @@
 from time import localtime, strftime, gmtime
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
-from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
+from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
 from anuga.utilities.system_tools import get_user_name
 
@@ -35 +35 @@
 
 # onshore data provided by WA DLI
-onshore_name = 'dampier_dli_ext' # original
+onshore_name = 'dli' # original
 
 # AHO + DPI data + colin French coastline
-coast_name = 'waterline'
-offshore_name = 'dampier_bathymetry'
-offshore1_name = 'nickelbay'
-offshore2_name = 'tinclip'
+coast_name = 'coastline'
+offshore_name = 'bathymetry'
+offshore1_name = 'elev_501'
+offshore2_name = 'inferrec_e'
 
 #final topo name
@@ -47 +47 @@
 combined_smaller_name = 'dampier_combined_elevation_smaller'
 
-topographies_in_dir = home+sep+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
+topographies_in_dir = home+sep+state+sep+scenario+sep+'elevation_final'+sep+'070112'+sep+'points'+sep
 topographies_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
 topographies_time_dir = topographies_dir+build_time+sep
@@ -126 +126 @@
 # Interior region definitions
 ###############################
-
+"""
 poly_pos20_neg20 = read_polygon(polygons_dir+'pos20_neg20_pts.csv')
 res_pos20_neg20 = 20000
@@ -155 +155 @@
 
 print 'min number triangles', trigs_min
-
+"""
 ###################################################################
 # Clipping regions for export to asc and regions for clipping data
 ###################################################################
 
-poly_bathy = read_polygon(polygons_dir+'polybathy.csv')
+#poly_bathy = read_polygon(polygons_dir+'polybathy.csv')
 
 # exporting asc grid - Dampier
@@ -169 +169 @@
 
 # exporting asc grid - Karratha
-e_min_area = 
-e_max_area = 
-n_min_area = 
-n_max_area = 
+##e_min_area = 
+##e_max_area = 
+##n_min_area = 
+##n_max_area = 
 
 """