Changeset 5381 for anuga_work/production/busselton/project.py

Timestamp:

May 30, 2008, 10:52:58 AM (17 years ago)

Author:

kristy

Message:

update scripts

File:

• anuga_work/production/busselton/project.py (modified) (4 diffs)

anuga_work/production/busselton/project.py

@@ -72 +72 @@
 # AHO + DPI data
 coast_name = 'coastline_excluding_beach_survey'
-coast_name1 = 'beach_survey_final'
-offshore_name = 'busselton'
+coast_name1 = 'bussleton_beach_data'
+offshore_name = 'Busselton_bathymetry'
 
 #final topo name
@@ -84 +84 @@
 topographies_dir = home+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
 
-#input topo file location
-onshore_in_dir_name = topographies_in_dir + onshore_name
-#island_in_dir_name = topographies_in_dir + island_name
+# input topo file location
+onshore_in_dir_name = topographies_in_dir + onshore_name #topo
 
-coast_in_dir_name = topographies_in_dir + coast_name
-coast_in_dir_name1 = topographies_in_dir + coast_name1
-offshore_in_dir_name = topographies_in_dir + offshore_name
+coast_in_dir_name = topographies_in_dir + coast_name #coastline
+coast_in_dir_name1 = topographies_in_dir + coast_name1 #beach survey
+offshore_in_dir_name = topographies_in_dir + offshore_name #bathymetry
 
+#output to anuga from build file
 onshore_dir_name = topographies_dir + onshore_name
-#island_dir_name = topographies_dir + island_name
+
 coast_dir_name = topographies_dir + coast_name
 coast_dir_name1 = topographies_dir + coast_name1
@@ -177 +177 @@
 res_busselton1 = 10000*res_factor
 
-poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly_.csv')
+poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly.csv')
 res_busselton2 = 500*res_factor
 
@@ -183 +183 @@
 
 interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]]
+print 'Hello',interior_regions
 
-boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 15],
-               'side': [10, 14], 'ocean': [11, 12, 13]}
+boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14],
+               'side': [10], 'ocean': [11, 12, 13]}
 
 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
 
-poly_mainland=read_polygon(polygons_dir+'initial_condition_south.csv')
+poly_mainland=read_polygon(topographies_in_dir+'initial_condition_south.csv')
 
 print 'min number triangles', trigs_min