Changeset 4616

Timestamp:

Jul 9, 2007, 3:24:23 PM (17 years ago)

Author:

ole

Message:

Shark bay work

Location:

anuga_work/production/shark_bay_2007

Files:

• build_shark_bay.py (modified) (1 diff)
• project.py (modified) (4 diffs)
• run_shark_bay.py (modified) (3 diffs)

anuga_work/production/shark_bay_2007/build_shark_bay.py

@@ -114 +114 @@
 print 'project.boundary_name', project.boundary_name
 
-urs_ungridded2sww(project.boundary_name, verbose=True);
-                  #verbose=True, mint=5000, maxt=35000, zscale=1)
+urs_ungridded2sww(project.boundary_name, verbose=True,
+                  mint=0, maxt=40000, zscale=1)
 
 

anuga_work/production/shark_bay_2007/project.py

@@ -31 +31 @@
 run_time = time+'_run'
 
-#tide = -3.9
-tide = 0.0
-#tide = 3.6
+#tide = 0.0 # MSL
+tide = 0.85 # HAT average between Denham and Canarvon
 
 #Maybe will try to make project a class to allow these parameters to be passed in.
 alpha = 0.1
 friction=0.01
-finaltime=25000
-starttime=3600
+finaltime=40000
+starttime=6000 # Action starts at 9000
 setup='final'
 #setup='trial'
 source='shark_bay'
 
-boundary_event = 'july2006'
-#boundary_event = '10000'
+#boundary_event = 'july2006'
+boundary_event = '10000'
 
 if setup =='trial':
@@ -59 +58 @@
 if setup =='final': 
     print'final'
-    res_factor=1.1
+    res_factor=1.0
     time_thinning=4
     yieldstep=60
 
 
+
 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
-             str(user)+'_'+boundary_event
+             str(user)+'_'+boundary_event+'_'+'tide'+str(tide)+'_coastpolys'
+
+#dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
+#             str(user)+'_'+boundary_event
 
 # elevation data filenames
@@ -160 +163 @@
                  
 
+
+#FIXME (Ole): It is bad to read in polygons in project.py. If a file does not exist
+#project.py cannot be imported
+
 # Bounding polygon
 bounding_polygon = read_polygon(join(polygons_dir, 'boundary_extent_small.csv'))
-res_bounding_polygon = 350000*res_factor
+res_bounding_polygon = 500000*res_factor
 
 #Interior regions
-poly_inundated_area = read_polygon(join(polygons_dir, 'inundated_area.csv'))
-res_inundated_area = 10*res_factor
-
-poly_tsunami_approach_area = read_polygon(join(polygons_dir, 'tsunami_approach_area.csv'))
-res_tsunami_approach_area = 50*res_factor
-
-poly_bay_area = read_polygon(join(polygons_dir, 'bay_area.csv'))                                   
-res_bay_area = 1000*res_factor
-
-poly_map_area = read_polygon(join(polygons_dir, 'map_area.csv'))                                   
-res_map_area = 5000*res_factor
-
-poly_model_area = read_polygon(join(polygons_dir, 'model_area.csv'))                                   
-res_model_area = 200000*res_factor                                   
-
-poly_exclusion_area1 = read_polygon(join(polygons_dir, 'exclusion_area1.csv'))
-res_exclusion_area1 = 2000000*res_factor
-
-interior_regions = [[poly_inundated_area, res_inundated_area],
-                    [poly_tsunami_approach_area, res_tsunami_approach_area],
-                    [poly_bay_area, res_bay_area],                    
-                    [poly_map_area, res_map_area]]
+interior_regions_and_resolutions = {
+    'study_area': 300000,
+    'bernier_dorre_islands': 250000,
+    'map_area': 10000,
+    'bay_area': 4000,
+    'south_bank_buffer': 2000,
+    'south_coast_polygon_1': 1000,
+    'south_coast_polygon_2': 1000,
+    'north_bay_coastline': 1000,
+    'tsunami_approach_area': 50,
+    'inundated_area': 10,
+    'exclusion_1': 10000000,
+    'exclusion_2': 10000000,
+    'exclusion_3': 10000000}
+
+
+#poly_exclusion_area1 = read_polygon(join(polygons_dir, 'exclusion_area1.csv'))
+#res_exclusion_area1 = 2000000*res_factor
+#
+#interior_regions = [[poly_inundated_area, res_inundated_area],
+#                    [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]]                                        
@@ -204 +212 @@
 
 
+
+# Read filenames etc (FIXME: move this to new file: read_structures)
+
+interior_regions = [] # Consider using dictionary for mesh interface.
+for area_name in interior_regions_and_resolutions.keys():
+    file_name = join(polygons_dir, area_name + '.csv')
+    
+    print 'Reading polygon from', file_name 
+    polygon = read_polygon(file_name)
+    base_resolution = interior_regions_and_resolutions[area_name]
+    
+    interior_regions.append([polygon, base_resolution*res_factor])
+
+
+
 trigs_min = number_mesh_triangles(interior_regions,
                                   bounding_polygon,

anuga_work/production/shark_bay_2007/run_shark_bay.py

@@ -93 +93 @@
     #IMPORTANT don't cache create_mesh_from_region and Domain(mesh....) as it
     # causes problems with the ability to cache set quantity which takes alot of times
+
+    # FIXME (Ole): Move to build_shark_bay.py
     if myid == 0:
     
@@ -125 +127 @@
     if myid == 0:
 
-        if project.tide == 0.0:
-            print 'Tide is MSL'
-            IC = 0.0
-        else:    
-            print 'Create onshore polygon'
-            from polygon import Polygon_function
-            #following sets the stage/water to be offcoast only
-            IC = Polygon_function([(project.onshore_polygon, -1.0)],
-                                  default = tide,
-                                  geo_reference = domain.geo_reference)
+        IC = project.tide
+        #if project.tide == 0.0:
+        #    print 'Tide is MSL'
+        #    IC = 0.0
+        #else:    
+        #    print 'Create onshore polygon'
+        #    from polygon import Polygon_function
+        #    #following sets the stage/water to be offcoast only
+        #    IC = Polygon_function([(project.onshore_polygon, -1.0)],
+        #                          default = tide,
+        #                          geo_reference = domain.geo_reference)
             
 
@@ -214 +217 @@
     # Reset stage in study area
     print 'Reset initial condition in study area'
-    domain.set_quantity('stage', 0.0,
+    domain.set_quantity('stage', project.tide,
                         polygon=project.poly_tsunami_approach_area)