Changeset 6015


Ignore:
Timestamp:
Nov 28, 2008, 11:26:25 AM (15 years ago)
Author:
kristy
Message:

Updated script for Pt Hedland Scenario

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/pt_hedland_2008/project_250m.py

    r5808 r6015  
    3838# Note, the user needs to set up the directory system accordingly
    3939state = 'western_australia'
    40 scenario_name = 'perth'
    41 scenario = 'perth_tsunami_scenario'
     40scenario_name = 'pt_hedland'
     41scenario = 'pt_hedland_tsunami_scenario_2008'
    4242
    4343# Model specific parameters. One or all can be changed each time the
     
    5151finaltime=80000         # final time for simulation
    5252
    53 interior_mesh = 'all' # Can have 'all' or 'none' for Phase 2 study
     53interior_mesh = 'none' # Can have 'all' or 'none' for Phase 2 study
    5454
    5555setup='final'  # Final can be replaced with trial or basic.
     
    8484#------------------------------------------------------------------------------
    8585
    86 # elevation data used in build_perth.py
     86# elevation data used in build_pt_hedland.py
    8787# onshore data: format ascii grid with accompanying projection file
    8888onshore_name = 'grid_250m_2005'
    8989
    9090# gauges - used in get_timeseries.py
    91 gauge_name = 'perth.csv'
     91gauge_name = 'pt_hedland.csv'
    9292gauge_name2 = 'thinned_MGA50.csv'
    9393
    94 # BOUNDING POLYGON - used in build_boundary.py and run_perth.py respectively
     94# BOUNDING POLYGON - used in build_boundary.py and run_pt_hedland.py respectively
    9595# NOTE: when files are put together the points must be in sequence - for ease go clockwise!
    96 # Check the run_perth.py for boundary_tags
     96# Check the run_pt_hedland.py for boundary_tags
    9797# thinned ordering file from Hazard Map: format is index,latitude,longitude (with title)
    98 order_filename = 'thinned_boundary_ordering.txt'
     98order_filename = 'thinned_boundary_ordering.csv'
    9999#landward bounding points
    100 landward = 'landward_bounding_polygon.txt'
     100landward = 'landward_bounding_polygon.csv'
    101101
    102102#------------------------------------------------------------------------------
     
    105105# Output filename for elevation
    106106# this is a combination of all the data (utilisied in build_boundary)
    107 combined_name ='perth_combined_elevation_250m'
     107combined_name ='pt_hedland_combined_elevation_250m'
    108108
    109109#------------------------------------------------------------------------------
     
    132132combined_dir_name = topographies_dir + combined_name
    133133
    134 # where the mesh sits (this is created during the run_perth.py)
     134# where the mesh sits (this is created during the run_pt_hedland.py)
    135135meshes_dir_name = meshes_dir + scenario_name+ interior_mesh +'.msh'
    136136
     
    138138order_filename_dir = boundaries_dir + order_filename
    139139
    140 # where the landward points of boundary extent sit (this is used within run_perth.py)
     140# where the landward points of boundary extent sit (this is used within run_pt_hedland.py)
    141141landward_dir = boundaries_dir + landward
    142142
     
    146146
    147147# where the directory of the output filename sits
    148 output_build_time_dir = output_dir+build_time+dir_comment+sep   #used for build_perth.py
    149 output_run_time_dir = output_dir+run_time+dir_comment+sep       #used for run_perth.py
     148output_build_time_dir = output_dir+build_time+dir_comment+sep   #used for build_pt_hedland.py
     149output_run_time_dir = output_dir+run_time+dir_comment+sep       #used for run_pt_hedland.py
    150150output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
    151151
     
    159159
    160160#Land, to set the initial stage/water to be offcoast only
    161 poly_mainland = read_polygon(polygons_dir+'initial_condition.csv')
     161poly_mainland = read_polygon(polygons_dir+'initial_conditions_mainland.csv')
     162
     163#Ocean
     164poly_ocean = read_polygon(polygons_dir+'initial_conditions_ocean.csv')
    162165
    163166# Initial bounding polygon for data clipping
     
    165168res_poly_all = 100000*res_factor
    166169
    167 # Area of Interest 1 (Fremantle)
    168 poly_aoi1 = read_polygon(polygons_dir+'CBD_coastal.csv')
     170# Area of Interest 1 (pt_hedland)
     171poly_aoi1 = read_polygon(polygons_dir+'area_of_interest.csv')
    169172res_aoi1 = 500*res_factor
    170173
    171 # Area of Interest 2 (Rockingham)
    172 poly_aoi2 = read_polygon(polygons_dir+'rockingham_penguin.csv')
    173 res_aoi2 = 500*res_factor
    174 
    175 # Area of Interest 2 (garden Island)
    176 poly_aoi2a = read_polygon(polygons_dir+'garden.csv')
    177 res_aoi2a= 500*res_factor
    178 
    179 # Area of Interest 3 (geordie bay - record of tsunami impact)
    180 poly_aoi3 = read_polygon(polygons_dir+'geordie_bay.csv')
    181 res_aoi3 = 500*res_factor
    182 
    183 # Area of Interest 4 (sorrento - record of tsunami impact)
    184 poly_aoi4 = read_polygon(polygons_dir+'sorrento_gauge.csv')
    185 res_aoi4 = 500*res_factor
    186 
    187 # Area of Significance 1 (Garden Island and sand bank infront of Rockingham)
    188 poly_aos1 = read_polygon(polygons_dir+'garden_rockingham.csv')
     174# Area of Significance 1 (pt_hedland)
     175poly_aos1 = read_polygon(polygons_dir+'area_of_significance.csv')
    189176res_aos1 = 1000*res_factor
    190177
    191 # Area of Significance 2 (incorporate coastline of rottnest)
    192 poly_aos2 = read_polygon(polygons_dir+'rottnest_external.csv')
    193 res_aos2 = 1000*res_factor
    194 
    195 # Refined areas
    196 # Polygon designed to incorporate dredged area from Fremantle to
    197 # Rockingham as the steep incline was making the elevation go to 0
    198 poly_aos3 = read_polygon(polygons_dir+'DredgeArea.csv')
    199 res_aos3 = 1000*res_factor
    200 
    201178# Shallow water 1
    202 poly_sw1 = read_polygon(polygons_dir+'internal_h20mORd3km.csv')
     179poly_sw1 = read_polygon(polygons_dir+'shallow_water.csv')
    203180res_sw1 = 25000*res_factor
    204 
    205 # Deep water (land of Rottnest, ANUGA does not do donuts!)
    206 poly_dw1 = read_polygon(polygons_dir+'rottnest_internal.csv')
    207 res_dw1 = 100000*res_factor
    208181
    209182# Combined all regions, must check that all are included!
     
    211184if interior_mesh =='all':
    212185    print'Mesh = all'
    213     interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2]
    214                      ,[poly_aoi2a,res_aoi2a],[poly_aoi3,res_aoi3]
    215                      ,[poly_aoi4,res_aoi4],[poly_aos1,res_aos1]
    216                      ,[poly_aos2,res_aos2],[poly_aos3,res_aos3]
    217                      ,[poly_sw1,res_sw1], [poly_dw1,res_dw1]]
     186    interior_regions = [[poly_aoi1,res_aoi1],[poly_aos1,res_aos1]
     187                     ,[poly_sw1,res_sw1]]
     188
    218189if interior_mesh =='none':
    219190    print'Mesh = none'
     
    223194print 'min estimated number of triangles', trigs_min
    224195   
    225 #------------------------------------------------------------------------------
    226 # Clipping regions for export to asc and regions for clipping data
    227 # Final inundation maps should only be created in regions of the finest mesh
    228 #------------------------------------------------------------------------------
    229 
    230 #Geordie Bay extract ascii grid
    231 xminGeordie = 358000
    232 xmaxGeordie = 362000
    233 yminGeordie = 6458500
    234 ymaxGeordie = 6461000
    235 
    236 #Sorrento extract ascii grid
    237 xminSorrento = 379000
    238 xmaxSorrento = 382500
    239 yminSorrento = 6477000
    240 ymaxSorrento = 6480000
    241 
    242 #Fremantle extract ascii grid
    243 xminFremantle = 376000
    244 xmaxFremantle = 388000
    245 yminFremantle = 6449000
    246 ymaxFremantle = 6461000
    247 
    248 #Rockingham extract ascii grid
    249 xminRockingham = 373500
    250 xmaxRockingham = 385500
    251 yminRockingham = 6424000
    252 ymaxRockingham = 6433000
    253 
Note: See TracChangeset for help on using the changeset viewer.