Changeset 4366

Timestamp:

Apr 11, 2007, 3:43:40 PM (18 years ago)

Author:

nick

Message:

updates to broome

Location:

anuga_work/production/broome_2006

Files:

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

anuga_work/production/broome_2006/build_broome.py

@@ -67 +67 @@
 offshore_dir_name1 = project.offshore_dir_name1
 offshore_dir_name2 = project.offshore_dir_name2
-
+'''
 # creates DEM from asc data
 print "creates DEMs from asc data"
@@ -83 +83 @@
 dem2pts(offshore_dir_name1, use_cache=True, verbose=True)
 dem2pts(offshore_dir_name2, use_cache=True, verbose=True)
-
-print'create Geospatial data1 objects from topographies'
+'''
+print'create Geospatial data1 objects from topographies',onshore_dir_name + '.pts'
 G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
-print'create Geospatial data2 objects from coast'
+print'create Geospatial data2 objects from coast', coast_in_dir_name + '.txt'
 G2 = Geospatial_data(file_name = coast_in_dir_name + '.txt')
 #print'create Geospatial data3 objects from island'
 #G3 = Geospatial_data(file_name = island_dir_name + '.pts')
-print'create Geospatial data3 objects from offshore'
+print'create Geospatial data3 objects from offshore',offshore_in_dir_name + '.txt'
 G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt')
-print'create Geospatial data4 objects from offshore1'
+print'create Geospatial data4 objects from offshore1',offshore_dir_name1 + '.pts'
 G_off1 = Geospatial_data(file_name = offshore_dir_name1 + '.pts')
-print'create Geospatial data4 objects from offshore2'
+print'create Geospatial data4 objects from offshore2',offshore_dir_name2 + '.pts'
 G_off2 = Geospatial_data(file_name = offshore_dir_name2 + '.pts')
 
@@ -105 +105 @@
 #print'shape of clipped data', G_clipped.get_data_points().shape
 
+print'export combined DEM file'
+if access(project.topographies_dir,F_OK) == 0:
+    mkdir (project.topographies_dir)
+print'export',project.combined_dir_name+'_unclipped' + '.txt'
+G.export_points_file(project.combined_dir_name+'_unclipped' + '.txt')
+'''
 print'split'
 G_small, G_other = G_clipped.split(.10,verbose=True)
 
-print'export combined DEM file'
-if access(project.topographies_dir,F_OK) == 0:
-    mkdir (project.topographies_dir)
+
+print 'export',project.combined_dir_name + '.txt'
 G_clipped.export_points_file(project.combined_dir_name + '.txt')
+print 'export', project.combined_small_dir_name + '.txt'
 G_small.export_points_file(project.combined_small_dir_name + '.txt')
 #G_clipped.export_points_file(project.combined_dir_name + '.xya')
-
+'''
 
 '''

anuga_work/production/broome_2006/project.py

@@ -119 +119 @@
 
 poly_all = read_polygon(polygons_dir+'extent.csv')
-res_poly_all = 150000
+res_poly_all = 500000
 
 ###############################
@@ -126 +126 @@
 
 poly_0 = read_polygon(polygons_dir+'neg20_coast_contour_pts.csv')
-res_0 = 10000
+res_0 = 100000
+
 poly_1 = read_polygon(polygons_dir+'broome_north_coast_inside_extent.csv')
-res_1 = 5000
+res_1 = 50000
 
 poly_2 = read_polygon(polygons_dir+'broome_south_coast_inside_extent.csv')
-res_2 = 5000
+res_2 = 50000
 
 poly_3 = read_polygon(polygons_dir+'Broome_town_pts.csv')
-res_3 = 2000
+res_3 = 20000
 
 poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
-res_4 = 500
+res_4 = 5000
 #assert zone == refzone