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

Timestamp:

Mar 10, 2008, 3:59:09 PM (17 years ago)

Author:

nick

Message:

update busselton

File:

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

anuga_work/production/busselton/project.py

@@ -32 +32 @@
 state = 'western_australia'
 scenario_name = 'busselton'
-scenario = 'busselton_tsunami_scenario_2006'
+scenario = 'busselton_tsunami_scenario'
 
 #Maybe will try to make project a class to allow these parameters to be passed in.
@@ -39 +39 @@
 starttime=10000
 midtime=21600
-finaltime=25000
+#finaltime=25000
+finaltime=10000
 export_cellsize=50
-setup='final'
-source='test'
+setup='trial'
+source='other'
 
 
@@ -65 +66 @@
 
 # onshore data provided by WA DLI
-onshore_name = 'DLI_orthophoto_DEM' # original
+#onshore_name = 'DLI_orthophoto_DEM' # original
+onshore_name = 'dli_dem_clipped' # original
 #islands
 
 # AHO + DPI data
-coast_name = '100k_coastline'
-offshore_name = 'Bathymetry'
+coast_name = 'coastline_excluding_beach_survey'
+coast_name1 = 'beach_survey_final'
+offshore_name = 'busselton'
 
 #final topo name
-combined_name ='busselton_combined_elevation.pts'
-combined_smaller_name = 'busselton_combined_elevation_smaller'
+combined_name ='busselton_combined_elevation'
+combined_name_small = 'busselton_combined_elevation_smaller'
 
 anuga_dir = home+state+sep+scenario+sep+'anuga'+sep
 
-topographies_in_dir = home+sep+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
-topographies_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
+topographies_in_dir = home+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
+topographies_dir = home+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
 
 #input topo file location
@@ -86 +89 @@
 
 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
 
@@ -91 +95 @@
 #island_dir_name = topographies_dir + island_name
 coast_dir_name = topographies_dir + coast_name
+coast_dir_name1 = topographies_dir + coast_name1
 offshore_dir_name = topographies_dir + offshore_name
 
 #final topo files
 combined_dir_name = topographies_dir + combined_name
-combined_smaller_dir_name = topographies_dir + combined_smaller_name
+combined_dir_name_small = topographies_dir + combined_name_small
 
 meshes_dir = anuga_dir+'meshes'+sep
@@ -115 +120 @@
     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
 
-if source=='test':
+if source=='other':
     boundaries_name = 'test' #exmouth gun
     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
@@ -134 +139 @@
 #gauges
 gauge_name = '???.csv'
-gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
+gauges_dir = home+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
 gauges_dir_name = gauges_dir + gauge_name
 
@@ -171 +176 @@
 poly_busselton1 = read_polygon(polygons_dir+'neg20_pos10_polygon.csv')
 res_busselton1 = 10000*res_factor
+
 poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly_.csv')
 res_busselton2 = 500*res_factor
@@ -178 +184 @@
 interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]]
 
-boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14],
-               'side': [10,13], 'ocean': [11, 12]}
+boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9],
+               'side': [10,15], 'ocean': [12, 13, 14]}
 
 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)