Changeset 4625

Timestamp:

Jul 16, 2007, 5:53:42 PM (17 years ago)

Author:

ole

Message:

Workin' on it

Location:

anuga_work/production/shark_bay_2007

Files:

• project.py (modified) (8 diffs)
• run_shark_bay.py (modified) (3 diffs)
• run_shark_bay_frequency_sweep.py (added)

anuga_work/production/shark_bay_2007/project.py

@@ -32 +32 @@
 
 #tide = 0.0 # MSL
-tide = 0.85 # HAT average between Denham and Canarvon
+tide = 0.79 # HAT average between Denham and Canarvon
+#tide = -0.80 # Estimated LAT
+
+# For frequency response study
+amplitude = 0.5
+period = 6
+
+#momentum_scale=50 # Experiment
+momentum_scale=1 # Experiment                        
 
 #Maybe will try to make project a class to allow these parameters to be passed in.
@@ -38 +46 @@
 friction=0.01
 finaltime=40000
-starttime=6000 # Action starts at 9000
+starttime=6000 # Action starts around 9000
 setup='final'
 #setup='trial'
 source='shark_bay'
 
-#boundary_event = 'july2006'
-boundary_event = '10000'
+boundary_event = 'july2006'
+#boundary_event = '10000'
+#boundary_event = 'experimental'
 
 if setup =='trial':
@@ -64 +73 @@
 
 
+#dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
+#             str(user)+'_'+boundary_event+'_A'+str(amplitude)+'_T'+str(period)
+
 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
-             str(user)+'_'+boundary_event+'_'+'tide'+str(tide)+'_coastpolys'
-
-#dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\
-#             str(user)+'_'+boundary_event
+             str(user)+'_'+boundary_event+'_X'+str(momentum_scale)+'_fixedbathy'
 
 # elevation data filenames
-ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem']
+ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem',
+                        'bathysteeppt', 'bathyleft', 'bathyright']
 point_filenames = ['field_survey_north', 'field_survey_south',
                    'clipped_bathymetry_final', 'coast_points_final']
@@ -115 +125 @@
                             'shark_bay', 'july_2006_event')
     else:
-        raise Exception, 'Unknown boundary event specified'
-
-boundary_name = join(boundary_dir, boundary_file_name)
+        pass
+        #raise Exception, 'Unknown boundary event specified'
+        
+if boundary_event != 'experimental':
+    boundary_name = join(boundary_dir, boundary_file_name)
+else:
+    boundary_name = 'nil'
 
 #output locations
@@ -169 +183 @@
 res_bounding_polygon = 500000*res_factor
 
-#Interior regions
+# This is a restricted boundary used for experiments.
+# Cannot use URS boundary on this!
+small_bounding_polygon = read_polygon(join(polygons_dir, 'small_boundary.csv'))
+res_bounding_polygon = 500000*res_factor
+
+
+# Interior regions
 interior_regions_and_resolutions = {
     'study_area': 300000,
@@ -209 +229 @@
                  'ocean': [3, 4, 5, 6]}
 
+#Very small bounding poly (experimental boundary)
+boundary_tags_small = {'tide': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13],
+                       'ocean': [11]}
+
 
 
@@ -225 +249 @@
 
 
+
 trigs_min = number_mesh_triangles(interior_regions,
                                   bounding_polygon,
@@ -230 +255 @@
 
 
+
+# For use with resetting IC (This is temporary)
+poly_tsunami_approach_area = read_polygon(
+    join(polygons_dir, 'tsunami_approach_area.csv'))
+
+
+
 onshore_polygon = read_polygon(join(topographies_in_dir,
                                     'initial_condition.txt'))

anuga_work/production/shark_bay_2007/run_shark_bay.py

@@ -26 +26 @@
 from anuga.shallow_water import Domain
 from anuga.shallow_water import Dirichlet_boundary
-from anuga.shallow_water import File_boundary
 from anuga.shallow_water import Reflective_boundary
 from anuga.shallow_water import Field_boundary
@@ -184 +183 @@
                         domain,
                         time_thinning=time_thinning,
-                        mean_stage=tide, 
+                        mean_stage=tide,
+                        momentum_scale=project.momentum_scale,                        
                         use_cache=False,
                         verbose=True)
@@ -201 +201 @@
     #                     'ocean': Bf})
     domain.set_boundary({'tide': Bd,
+    #                     'ocean': Dirichlet_boundary([3,0,0])})
                          'ocean': Bf})     
     print'finish set boundary'