Changeset 4625
- Timestamp:
- Jul 16, 2007, 5:53:42 PM (17 years ago)
- Location:
- anuga_work/production/shark_bay_2007
- Files:
-
- 1 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/shark_bay_2007/project.py
r4624 r4625 32 32 33 33 #tide = 0.0 # MSL 34 tide = 0.85 # HAT average between Denham and Canarvon 34 tide = 0.79 # HAT average between Denham and Canarvon 35 #tide = -0.80 # Estimated LAT 36 37 # For frequency response study 38 amplitude = 0.5 39 period = 6 40 41 #momentum_scale=50 # Experiment 42 momentum_scale=1 # Experiment 35 43 36 44 #Maybe will try to make project a class to allow these parameters to be passed in. … … 38 46 friction=0.01 39 47 finaltime=40000 40 starttime=6000 # Action starts a t900048 starttime=6000 # Action starts around 9000 41 49 setup='final' 42 50 #setup='trial' 43 51 source='shark_bay' 44 52 45 #boundary_event = 'july2006' 46 boundary_event = '10000' 53 boundary_event = 'july2006' 54 #boundary_event = '10000' 55 #boundary_event = 'experimental' 47 56 48 57 if setup =='trial': … … 64 73 65 74 75 #dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\ 76 # str(user)+'_'+boundary_event+'_A'+str(amplitude)+'_T'+str(period) 77 66 78 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\ 67 str(user)+'_'+boundary_event+'_'+'tide'+str(tide)+'_coastpolys' 68 69 #dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+\ 70 # str(user)+'_'+boundary_event 79 str(user)+'_'+boundary_event+'_X'+str(momentum_scale)+'_fixedbathy' 71 80 72 81 # elevation data filenames 73 ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem'] 82 ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem', 83 'bathysteeppt', 'bathyleft', 'bathyright'] 74 84 point_filenames = ['field_survey_north', 'field_survey_south', 75 85 'clipped_bathymetry_final', 'coast_points_final'] … … 115 125 'shark_bay', 'july_2006_event') 116 126 else: 117 raise Exception, 'Unknown boundary event specified' 118 119 boundary_name = join(boundary_dir, boundary_file_name) 127 pass 128 #raise Exception, 'Unknown boundary event specified' 129 130 if boundary_event != 'experimental': 131 boundary_name = join(boundary_dir, boundary_file_name) 132 else: 133 boundary_name = 'nil' 120 134 121 135 #output locations … … 169 183 res_bounding_polygon = 500000*res_factor 170 184 171 #Interior regions 185 # This is a restricted boundary used for experiments. 186 # Cannot use URS boundary on this! 187 small_bounding_polygon = read_polygon(join(polygons_dir, 'small_boundary.csv')) 188 res_bounding_polygon = 500000*res_factor 189 190 191 # Interior regions 172 192 interior_regions_and_resolutions = { 173 193 'study_area': 300000, … … 209 229 'ocean': [3, 4, 5, 6]} 210 230 231 #Very small bounding poly (experimental boundary) 232 boundary_tags_small = {'tide': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13], 233 'ocean': [11]} 234 211 235 212 236 … … 225 249 226 250 251 227 252 trigs_min = number_mesh_triangles(interior_regions, 228 253 bounding_polygon, … … 230 255 231 256 257 258 # For use with resetting IC (This is temporary) 259 poly_tsunami_approach_area = read_polygon( 260 join(polygons_dir, 'tsunami_approach_area.csv')) 261 262 263 232 264 onshore_polygon = read_polygon(join(topographies_in_dir, 233 265 'initial_condition.txt')) -
anuga_work/production/shark_bay_2007/run_shark_bay.py
r4616 r4625 26 26 from anuga.shallow_water import Domain 27 27 from anuga.shallow_water import Dirichlet_boundary 28 from anuga.shallow_water import File_boundary29 28 from anuga.shallow_water import Reflective_boundary 30 29 from anuga.shallow_water import Field_boundary … … 184 183 domain, 185 184 time_thinning=time_thinning, 186 mean_stage=tide, 185 mean_stage=tide, 186 momentum_scale=project.momentum_scale, 187 187 use_cache=False, 188 188 verbose=True) … … 201 201 # 'ocean': Bf}) 202 202 domain.set_boundary({'tide': Bd, 203 # 'ocean': Dirichlet_boundary([3,0,0])}) 203 204 'ocean': Bf}) 204 205 print'finish set boundary'
Note: See TracChangeset
for help on using the changeset viewer.