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

Timestamp:

Jul 25, 2008, 4:17:26 PM (15 years ago)

Author:

kristy

Message:

Addition of Polyline information with sts files

File:

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

anuga_work/production/busselton/project.py

@@ -10 +10 @@
 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
 from anuga.utilities.system_tools import get_user_name, get_host_name
+from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
 
 # file and system info
@@ -16 +17 @@
 
 home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir   
+muxhome = getenv('MUXHOME')
 user = get_user_name()
 host = get_host_name()
@@ -33 +35 @@
 scenario = 'busselton_tsunami_scenario'
 
-tide = 0
+tide = 0.6
 
 alpha = 0.1 
 friction=0.01
 starttime=0
-finaltime=50000 # changed for rev_number trial
+finaltime=80000 
 export_cellsize=25
 setup='final'
-source='exmouth'
+source='polyline'
 
 
@@ -135 +137 @@
 #boundaries_source = '1'
 
-if source=='exmouth':
-    boundaries_name = 'busselton_3103_30052008' #exmouth gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
-
-if source=='test':
-    boundaries_name = 'other' #exmouth gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep
-
+##if source=='exmouth':
+##    boundaries_name = 'busselton_3103_30052008' # exmouth gun
+##    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
+##
+##if source=='test':
+##    boundaries_name = 'other' #exmouth gun
+##    boundaries_in_dir = anuga_dir+'boundaries'+sep
+##
 
 #boundaries locations
-boundaries_in_dir_name = boundaries_in_dir + boundaries_name
+#boundaries_in_dir_name = boundaries_in_dir + boundaries_name
 boundaries_dir = anuga_dir+'boundaries'+sep
 boundaries_dir_name = boundaries_dir + scenario_name
+boundaries_dir_mux = muxhome
 
 #output locations
 output_dir = anuga_dir+'outputs'+sep
-#output_build_time_dir = output_dir +'test'+sep + build_time + dir_comment + sep
 output_build_time_dir = output_dir +build_time + dir_comment + sep
-#output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir = output_dir +'test'+sep + run_time + dir_comment +sep # take test out
+output_run_time_dir = output_dir + run_time + dir_comment +sep 
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 
@@ -170 +171 @@
 
 ###############################
-# Domain definitions
+# Interior region definitions
 ###############################
 
-from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
-
-# bounding polygon for study area
+# Initial bounding polygon for data clipping
 poly_all = read_polygon(polygons_dir+'poly_all_extend.csv')
 res_poly_all = 100000*res_factor
-
-
-###############################
-# Interior region definitions
-###############################
 
 #digitized polygons
@@ -209 +203 @@
 
 
-boundary_tags={'back': [3, 4, 5, 6, 7],
-               'side': [2, 8], 'ocean': [0, 1, 9]}
-
 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+print 'min number triangles', trigs_min
 
 poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
-
-print 'min number triangles', trigs_min