Changeset 5825

Timestamp:

Oct 8, 2008, 9:42:45 AM (16 years ago)

Author:

kristy

Message:

 

File:

• anuga_work/production/perth/project.py (modified) (11 diffs)

anuga_work/production/perth/project.py

@@ -38 +38 @@
 # Note, the user needs to set up the directory system accordingly
 state = 'western_australia'
-scenario_name = 'onslow'
-scenario = 'onslow_tsunami_scenario_2008'
+scenario_name = 'perth'
+scenario = 'perth_tsunami_scenario'
 
 # Model specific parameters. One or all can be changed each time the
 # run_scenario script is executed
-tide = 0                #0.6
+tide = 0.6                #0.6
 #event_number = 27255   # linked to hazard map
 event_number = 27283
@@ -82 +82 @@
 #------------------------------------------------------------------------------
 
-# elevation data used in build_onslow.py
+# elevation data used in build_perth.py
 # onshore data: format ascii grid with accompanying projection file
-onshore_name = 'revised_onshore_20m_dli' 
+onshore_name = 'perth_dli_ext' 
 # island: format ascii grid with accompanying projection file
-island_name = 'onslow_islands_dted2' 
+island_name = 'rott_dli_ext' 
+island_name1 = 'gard_dli_ext'
+island_name2 = 'carnac_island_dli_ext'
+island_name3 = 'penguin_dli_ext'
 # coastline: format x,y,elevation (with title)
-coast_name = 'revised_coastline.txt'
+coast_name = 'coastline_perthP.txt'
 # bathymetry: format x,y,elevation (with title)
-offshore_name = 'onslow_offshore_points.txt'
-offshore_name1 = 'ONS16FINAL.txt'
-offshore_name2 = 'ONS18FINAL.txt'
-offshore_name3 = 'ONS18FINAL.txt'
-offshore_name4 = 'beach_survey.txt'
+offshore_name = 'Perth.txt'
+offshore_name1 = 'Perth_Chart.txt'
+offshore_name2 = 'Fremantle_north.txt'
+offshore_name3 = 'port_swanriver.txt'
 
 # gauges - used in get_timeseries.py
-#gauge_name = 'onslow.csv'
-#gauge_name2 = 'thinned_MGA50.csv'
-
-# BOUNDING POLYGON - used in build_boundary.py and run_onslow.py respectively
+gauge_name = 'perth.csv'
+gauge_name2 = 'thinned_MGA50.csv'
+
+# BOUNDING POLYGON - used in build_boundary.py and run_perth.py respectively
 # NOTE: when files are put together the points must be in sequence - for ease go clockwise!
-# Check the run_onslow.py for boundary_tags
+# Check the run_perth.py for boundary_tags
 # thinned ordering file from Hazard Map: format is index,latitude,longitude (with title)
 order_filename = 'thinned_boundary_ordering.csv'
@@ -113 +115 @@
 # Output filename for elevation
 # this is a combination of all the data (utilisied in build_boundary)
-combined_name ='onslow_combined_elevation'
-combined_smaller_name = 'onslow_combined_elevation_smaller'
+combined_name ='perth_combined_elevation'
+combined_smaller_name = 'perth_combined_elevation_smaller'
 
 #------------------------------------------------------------------------------
@@ -135 +137 @@
 onshore_in_dir_name = topographies_in_dir + onshore_name
 island_in_dir_name = topographies_in_dir + island_name
+island_in_dir_name1 = topographies_in_dir + island_name1
+island_in_dir_name2 = topographies_in_dir + island_name2
+island_in_dir_name3 = topographies_in_dir + island_name3
 coast_in_dir_name = topographies_in_dir + coast_name
 offshore_in_dir_name = topographies_in_dir + offshore_name
@@ -140 +145 @@
 offshore_in_dir_name2 = topographies_in_dir + offshore_name2
 offshore_in_dir_name3 = topographies_in_dir + offshore_name3
-offshore_in_dir_name4 = topographies_in_dir + offshore_name4
 
 # where the output data sits
 onshore_dir_name = topographies_dir + onshore_name
 island_dir_name = topographies_dir + island_name
+island_dir_name1 = topographies_dir + island_name1
+island_dir_name2 = topographies_dir + island_name2
+island_dir_name3 = topographies_dir + island_name3
 coast_dir_name = topographies_dir + coast_name
 offshore_dir_name = topographies_dir + offshore_name
@@ -150 +157 @@
 offshore_dir_name2 = topographies_dir + offshore_name2
 offshore_dir_name3 = topographies_dir + offshore_name3
-offshore_dir_name4 = topographies_dir + offshore_name4
 
 # where the combined elevation file sits
@@ -156 +162 @@
 combined_smaller_name_dir = topographies_dir + combined_smaller_name
 
-# where the mesh sits (this is created during the run_onslow.py)
+# where the mesh sits (this is created during the run_perth.py)
 meshes_dir_name = meshes_dir + scenario_name+'.msh'
 
@@ -162 +168 @@
 order_filename_dir = boundaries_dir + order_filename
 
-# where the landward points of boundary extent sit (this is used within run_onslow.py)
+# where the landward points of boundary extent sit (this is used within run_perth.py)
 landward_dir = boundaries_dir + landward
 
@@ -170 +176 @@
 
 # where the directory of the output filename sits
-output_build_time_dir = output_dir+build_time+dir_comment+sep   #used for build_onslow.py
-output_run_time_dir = output_dir+run_time+dir_comment+sep       #used for run_onslow.py
+output_build_time_dir = output_dir+build_time+dir_comment+sep   #used for build_perth.py
+output_run_time_dir = output_dir+run_time+dir_comment+sep       #used for run_perth.py
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 
 #w here the directory of the gauges sit
-#gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
-#gauges_dir_name2 = gauges_dir + gauge_name2     #used for get_timeseries.py
+gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
+gauges_dir_name2 = gauges_dir + gauge_name2     #used for get_timeseries.py
 
 #------------------------------------------------------------------------------
@@ -182 +188 @@
 #------------------------------------------------------------------------------
 
-<<<<<<< .mine
-#Land North, to set the initial stage/water to be offcoast only
-poly_mainland_N = read_polygon(polygons_dir+'initial_conditions_N.csv')
-=======
 #Land, to set the initial stage/water to be offcoast only
-poly_mainland = read_polygon(polygons_dir+'initial_conditions_mainland.csv')
->>>>>>> .r5816
-
-<<<<<<< .mine
-#Land South, to set the initial stage/water to be offcoast only
-poly_mainland_S = read_polygon(polygons_dir+'initial_conditions_S.csv')
-
-#Rottnest, to set the initial stage/water to be offcoast only
-poly_island1 = read_polygon(polygons_dir+'initial_conditions_rottnest.csv')
-
-#Garden, to set the initial stage/water to be offcoast only
-poly_island2 = read_polygon(polygons_dir+'initial_conditions_garden.csv')
-
-=======
-#Ocean
-poly_ocean = read_polygon(polygons_dir+'initial_conditions_ocean.csv')
-
->>>>>>> .r5816
+poly_mainland = read_polygon(polygons_dir+'initial_condition.csv')
+
 # Initial bounding polygon for data clipping 
 poly_all = read_polygon(polygons_dir+'poly_all.csv')
 res_poly_all = 100000*res_factor
 
-# Area of Interest 1 (Onslow)
-poly_aoi1 = read_polygon(polygons_dir+'aoi.csv')
+# Area of Interest 1 (Fremantle)
+poly_aoi1 = read_polygon(polygons_dir+'CBD_coastal.csv')
 res_aoi1 = 500*res_factor
 
-# Area of Significance 1 (Onslow)
-poly_aos1 = read_polygon(polygons_dir+'aos.csv')
+# Area of Interest 2 (Rockingham)
+poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin.csv')
+res_aoi2 = 500*res_factor
+
+# Area of Interest 2 (garden Island)
+poly_aoi2a = read_polygon(polygons_dir+'garden.csv')
+res_aoi2a= 500*res_factor
+
+# Area of Interest 3 (geordie bay - record of tsunami impact)
+poly_aoi3 = read_polygon(polygons_dir+'geordie_bay.csv')
+res_aoi3 = 500*res_factor
+
+# Area of Interest 4 (sorrento - record of tsunami impact)
+poly_aoi4 = read_polygon(polygons_dir+'sorrento_gauge.csv')
+res_aoi4 = 500*res_factor
+
+# Area of Significance 1 (Garden Island and sand bank infront of Rockingham)
+poly_aos1 = read_polygon(polygons_dir+'garden_rockingham.csv')
 res_aos1 = 1000*res_factor
 
+# Area of Significance 2 (incorporate coastline of rottnest)
+poly_aos2 = read_polygon(polygons_dir+'rottnest_external.csv')
+res_aos2 = 1000*res_factor
+
+# Refined areas
+# Polygon designed to incorporate dredged area from Fremantle to
+# Rockingham as the steep incline was making the elevation go to 0
+poly_aos3 = read_polygon(polygons_dir+'DredgeArea.csv')
+res_aos3 = 1000*res_factor
+
 # Shallow water 1
-poly_sw1 = read_polygon(polygons_dir+'shallow_water.csv')
+poly_sw1 = read_polygon(polygons_dir+'internal_h20mORd3km.csv')
 res_sw1 = 25000*res_factor
 
+# Deep water (land of Rottnest, ANUGA does not do donuts!)
+poly_dw1 = read_polygon(polygons_dir+'rottnest_internal.csv')
+res_dw1 = 100000*res_factor
+
 # Combined all regions, must check that all are included!
-interior_regions = [[poly_aoi1,res_aoi1],[poly_aos1,res_aos1]
-                     ,[poly_sw1,res_sw1]]
+interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2]
+                     ,[poly_aoi2a,res_aoi2a],[poly_aoi3,res_aoi3]
+                     ,[poly_aoi4,res_aoi4],[poly_aos1,res_aos1]
+                     ,[poly_aos2,res_aos2],[poly_aos3,res_aos3]
+                     ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]]
 
     
@@ -234 +253 @@
 #------------------------------------------------------------------------------
 
-###Geordie Bay extract ascii grid
-##xminGeordie = 358000
-##xmaxGeordie = 362000
-##yminGeordie = 6458500
-##ymaxGeordie = 6461000
-##
-###Sorrento extract ascii grid
-##xminSorrento = 379000
-##xmaxSorrento = 382500
-##yminSorrento = 6477000
-##ymaxSorrento = 6480000
-##
-###Fremantle extract ascii grid
-##xminFremantle = 376000
-##xmaxFremantle = 388000
-##yminFremantle = 6449000
-##ymaxFremantle = 6461000
-##
-###Rockingham extract ascii grid
-##xminRockingham = 373500
-##xmaxRockingham = 385500
-##yminRockingham = 6424000
-##ymaxRockingham = 6433000
-
+#Geordie Bay extract ascii grid
+xminGeordie = 358000
+xmaxGeordie = 362000
+yminGeordie = 6458500
+ymaxGeordie = 6461000
+
+#Sorrento extract ascii grid
+xminSorrento = 379000
+xmaxSorrento = 382500
+yminSorrento = 6477000
+ymaxSorrento = 6480000
+
+#Fremantle extract ascii grid
+xminFremantle = 376000
+xmaxFremantle = 388000
+yminFremantle = 6449000
+ymaxFremantle = 6461000
+
+#Rockingham extract ascii grid
+xminRockingham = 373500
+xmaxRockingham = 385500
+yminRockingham = 6424000
+ymaxRockingham = 6433000
+