Changeset 6391 for anuga_work/production/perth/project.py
- Timestamp:
- Feb 24, 2009, 9:18:23 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/perth/project.py
r6342 r6391 44 44 # run_scenario script is executed 45 45 tide = 0.6 #0.6 46 #event_number = 27255 # Java 9.3 worst case for Perth46 event_number = 27255 # Java 9.3 worst case for Perth 47 47 #event_number = 68693 # Sumatra 9.2 48 event_number = 27283 # Java 9.3 original48 #event_number = 27283 # Java 9.3 original 49 49 alpha = 0.1 # smoothing parameter for mesh 50 50 friction=0.01 # manning's friction coefficient … … 77 77 # Important to distinguish each run - ensure str(user) is included! 78 78 # Note, the user is free to include as many parameters as desired 79 dir_comment='_'+setup+'_'+str(tide)+'_'+str(event_number)+'_'+ 'R erun_'+str(user)79 dir_comment='_'+setup+'_'+str(tide)+'_'+str(event_number)+'_'+ 'Rottnest_'+str(user) 80 80 81 81 #------------------------------------------------------------------------------ … … 168 168 169 169 # where the mesh sits (this is created during the run_perth.py) 170 meshes_dir_name = meshes_dir + scenario_name+' .msh'170 meshes_dir_name = meshes_dir + scenario_name+'_rottnest.msh' 171 171 172 172 # where the boundary ordering files sit (this is used within build_boundary.py) … … 213 213 ##res_aoi1 = 500*res_factor 214 214 215 # Area of Interest 1 (Fremantle)216 poly_aoi1 = read_polygon(polygons_dir+'CBD_coastal.csv')217 res_aoi1 = 500*res_factor218 219 # Area of Interest 2 (Rockingham)220 poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin_increase.csv')221 res_aoi2 = 500*res_factor222 223 # Area of Interest 2 (Rockingham)224 poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin.csv')225 res_aoi2 = 500*res_factor226 227 # Area of Interest 2 (garden Island)228 poly_aoi2a = read_polygon(polygons_dir+'garden.csv')229 res_aoi2a= 500*res_factor230 231 # Area of Interest 3 (geordie bay - record of tsunami impact)232 poly_aoi3 = read_polygon(polygons_dir+'geordie_bay.csv')233 res_aoi3 = 500*res_factor234 235 # Area of Interest 4 (sorrento - record of tsunami impact)236 poly_aoi4 = read_polygon(polygons_dir+'sorrento_gauge.csv')237 res_aoi4 = 500*res_factor238 239 # Area of Significance 1 (Garden Island and sand bank infront of Rockingham)240 poly_aos1 = read_polygon(polygons_dir+'garden_rockingham_increase.csv')241 res_aos1 = 1000*res_factor215 ### Area of Interest 1 (Fremantle) 216 ##poly_aoi1 = read_polygon(polygons_dir+'CBD_coastal.csv') 217 ##res_aoi1 = 500*res_factor 218 ## 219 ### Area of Interest 2 (Rockingham) 220 ##poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin_increase.csv') 221 ##res_aoi2 = 500*res_factor 222 ## 223 ### Area of Interest 2 (Rockingham) 224 ##poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin.csv') 225 ##res_aoi2 = 500*res_factor 226 ## 227 ### Area of Interest 2 (garden Island) 228 ##poly_aoi2a = read_polygon(polygons_dir+'garden.csv') 229 ##res_aoi2a= 500*res_factor 230 ## 231 ### Area of Interest 3 (geordie bay - record of tsunami impact) 232 ##poly_aoi3 = read_polygon(polygons_dir+'geordie_bay.csv') 233 ##res_aoi3 = 500*res_factor 234 ## 235 ### Area of Interest 4 (sorrento - record of tsunami impact) 236 ##poly_aoi4 = read_polygon(polygons_dir+'sorrento_gauge.csv') 237 ##res_aoi4 = 500*res_factor 238 ## 239 ### Area of Significance 1 (Garden Island and sand bank infront of Rockingham) 240 ##poly_aos1 = read_polygon(polygons_dir+'garden_rockingham_increase.csv') 241 ##res_aos1 = 1000*res_factor 242 242 243 243 # Area of Significance 2 (incorporate coastline of rottnest) 244 244 poly_aos2 = read_polygon(polygons_dir+'rottnest_external.csv') 245 res_aos2 = 1000*res_factor 246 247 # Refined areas 248 # Polygon designed to incorporate dredged area from Fremantle to 249 # Rockingham as the steep incline was making the elevation go to 0 250 poly_aos3 = read_polygon(polygons_dir+'DredgeArea.csv') 251 res_aos3 = 1000*res_factor 252 253 # Shallow water 1 254 poly_sw1 = read_polygon(polygons_dir+'internal_h20mORd3km.csv') 255 res_sw1 = 25000*res_factor 245 res_aos2 = 500*res_factor 246 247 ### Refined areas 248 ### Polygon designed to incorporate dredged area from Fremantle to 249 ### Rockingham as the steep incline was making the elevation go to 0 250 ##poly_aos3 = read_polygon(polygons_dir+'DredgeArea.csv') 251 ##res_aos3 = 1000*res_factor 252 ## 253 ### Shallow water 1 254 ##poly_sw1 = read_polygon(polygons_dir+'internal_h20mORd3km.csv') 255 ##res_sw1 = 25000*res_factor 256 257 ### Shallow water 1 258 poly_sw1 = read_polygon(polygons_dir+'rottnest_20m.csv') 259 res_sw1 = 1000*res_factor 256 260 257 261 # Deep water (land of Rottnest, ANUGA does not do donuts!) … … 264 268 ## ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]] 265 269 266 interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2] 267 ,[poly_aoi2a,res_aoi2a],[poly_aoi3,res_aoi3] 268 ,[poly_aoi4,res_aoi4],[poly_aos1,res_aos1] 269 ,[poly_aos2,res_aos2],[poly_aos3,res_aos3] 270 ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]] 270 #For Rottnest zoom in 271 interior_regions = [[poly_aos2,res_aos2],[poly_sw1,res_sw1] 272 ,[poly_dw1,res_dw1]] 273 274 ##interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2] 275 ## ,[poly_aoi2a,res_aoi2a],[poly_aoi3,res_aoi3] 276 ## ,[poly_aoi4,res_aoi4],[poly_aos1,res_aos1] 277 ## ,[poly_aos2,res_aos2],[poly_aos3,res_aos3] 278 ## ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]] 271 279 272 280 … … 309 317 ymaxPerth = 6467000 310 318 319 #Rottnest 320 xminRottnest = 352700 321 xmaxRottnest = 364800 322 yminRottnest = 6455000 323 ymaxRottnest = 6460500 324 325
Note: See TracChangeset
for help on using the changeset viewer.