Changeset 4601 for anuga_work/production/shark_bay_2007/project.py

Timestamp:

Jul 6, 2007, 2:40:57 PM (18 years ago)

Author:

ole

Message:

Work on shark bay

File:

• anuga_work/production/shark_bay_2007/project.py (modified) (5 diffs)

anuga_work/production/shark_bay_2007/project.py

@@ -41 +41 @@
 starttime=3600
 setup='final'
+#setup='trial'
 source='shark_bay'
+
+boundary_event = 'july2006'
+#boundary_event = '10000'
 
 if setup =='trial':
@@ -60 +64 @@
 
 
-#dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(user)+'_'+'nobdry'
-dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(user)
+dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
+             str(user)+'_'+boundary_event
 
 # elevation data filenames
@@ -100 +104 @@
 # locations for boundary conditions
 if source=='shark_bay':
-    boundary_file_name = 'shark_bay_3867_18052007'
-    boundary_dir = join(anuga_dir, 'boundaries', 'shark_bay', '1_10000')
+    if boundary_event == '10000':
+        boundary_file_name = 'shark_bay_3867_18052007'
+        boundary_dir = join(anuga_dir, 'boundaries', 'shark_bay', '1_10000')
+    elif boundary_event == 'july2006':
+        boundary_file_name = 'july2006'    
+        boundary_dir = join(anuga_dir, 'boundaries',
+                            'shark_bay', 'july_2006_event')
+    else:
+        raise Exception, 'Unknown boundary event specified'
 
 boundary_name = join(boundary_dir, boundary_file_name)
@@ -150 +161 @@
 
 # Bounding polygon
-bounding_polygon = read_polygon(join(polygons_dir, 'boundary_extent.csv'))
-res_bounding_polygon = 250000*res_factor
+bounding_polygon = read_polygon(join(polygons_dir, 'boundary_extent_small.csv'))
+res_bounding_polygon = 350000*res_factor
 
 #Interior regions
 poly_inundated_area = read_polygon(join(polygons_dir, 'inundated_area.csv'))
-res_inundated_area = 200*res_factor
+res_inundated_area = 10*res_factor
 
 poly_tsunami_approach_area = read_polygon(join(polygons_dir, 'tsunami_approach_area.csv'))
-res_tsunami_approach_area = 500*res_factor
+res_tsunami_approach_area = 50*res_factor
 
 poly_bay_area = read_polygon(join(polygons_dir, 'bay_area.csv'))                                   
@@ -175 +186 @@
                     [poly_tsunami_approach_area, res_tsunami_approach_area],
                     [poly_bay_area, res_bay_area],                    
-                    [poly_map_area, res_map_area],
-                    [poly_model_area, res_model_area],
-                    [poly_exclusion_area1, res_exclusion_area1]]                                        
+                    [poly_map_area, res_map_area]]
+                    #[poly_model_area, res_model_area],
+                    #[poly_exclusion_area1, res_exclusion_area1]]                                        
+
+# FIXME (Ole): What if tags are wrong?
+
+#Large bounding poly
+#boundary_tags = {'back': [1, 2, 3, 4, 5],   
+#                 'side': [0, 6],
+#                 'ocean': [7, 8, 9, 10, 11]}
+
+
+#Small bounding poly
+#(counter clockwise starting in upper right corner)
+boundary_tags = {'tide': [0, 1, 2],
+                 'ocean': [3, 4, 5, 6]}