Changeset 4601 for anuga_work

Timestamp:

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

Author:

ole

Message:

Work on shark bay

Location:

anuga_work/production/shark_bay_2007

Files:

• project.py (modified) (5 diffs)
• run_shark_bay.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]}
 
 

anuga_work/production/shark_bay_2007/run_shark_bay.py

@@ -74 +74 @@
     if myid == 0:
         copy_code_files(kwargs['output_dir'],__file__, 
-                 dirname(project.__file__)+sep+ project.__name__+'.py' )
+                        dirname(project.__file__)+sep+project.__name__+'.py' )
 
         store_parameters(**kwargs)
@@ -97 +97 @@
         print 'start create mesh from regions'
         
-        # FIXME (Ole): What if tags are wrong?
         create_mesh_from_regions(project.bounding_polygon,
-                                 boundary_tags={'back': [1, 2, 3, 4, 5],
-                                                'side': [0, 6],
-                                                'ocean': [7, 8, 9, 10, 11]},
+                                 boundary_tags=project.boundary_tags,
                                  maximum_triangle_area=project.res_bounding_polygon,
                                  interior_regions=project.interior_regions,
@@ -185 +182 @@
                         time_thinning=time_thinning,
                         mean_stage=tide, 
-                        use_cache=True,
+                        use_cache=False,
                         verbose=True)
                     
@@ -197 +194 @@
     print'set_boundary'
 
-    domain.set_boundary({'back': Br,
-                         'side': Bd,
-                         'ocean': Bf}) 
+    #domain.set_boundary({'back': Br,
+    #                     'side': Bd,
+    #                     'ocean': Bf})
+    domain.set_boundary({'tide': Bd,
+                         'ocean': Bf})     
     print'finish set boundary'
 
@@ -211 +210 @@
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')     
+
+
+    # Reset stage in study area
+    print 'Reset initial condition in study area'
+    domain.set_quantity('stage', 0.0,
+                        polygon=project.poly_tsunami_approach_area)
+    
 
     for t in domain.evolve(yieldstep = kwargs['yieldstep'], finaltime = 21600