Changeset 5489


Ignore:
Timestamp:
Jul 11, 2008, 3:38:03 PM (11 years ago)
Author:
kristy
Message:

changed data input

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/project.py

    r5467 r5489  
    6767island_name = 'rott_dli_ext' # original
    6868island_name1 = 'gard_dli_ext'
    69 island_name2 = 'carnac_island_dted'
    70 island_name3 = 'penguin_dted'
     69island_name2 = 'carnac_dli_ext'
     70island_name3 = 'penguin_dli_ext'
    7171
    7272# AHO + DPI data + colin French coastline
    73 coast_name = 'coastline'
    74 offshore_name = 'perth_bathymetry'
     73coast_name = 'coastline_perthP'
     74offshore_name = 'Perth_clip'
    7575
    7676
     
    158158from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    159159
    160 poly_all = read_polygon(polygons_dir+'bounding_area.csv')
     160poly_all = read_polygon(polygons_dir+'poly_all.csv')
    161161res_poly_all = 100000*res_factor
    162162
     
    171171res_internal_20_3 = 25000*res_factor
    172172
     173#Polygon designed to cut out the rottnest island land.
     174poly_rottnest_in = read_polygon(polygons_dir+'rottnest_internal.csv')
     175res_rottnest_in = 25000*res_factor
     176
    173177#Polygon designed to incorporate Garden Island and sand bank infront of Rockingham
    174178poly_garden_rockingham = read_polygon(polygons_dir+'garden_rockingham.csv')
    175 res_garden_rockingham = 1500*res_factor
     179res_garden_rockingham = 2500*res_factor
    176180
    177181#Polygon designed to incorporate coastline of rottnest
     
    179183res_rottnest_ex = 1000*res_factor
    180184
    181 #Polygon designed to cut out the rottnest island land.
    182 poly_rottnest_in = read_polygon(polygons_dir+'rottnest_internal.csv')
    183 res_rottnest_in = 25000*res_factor
    184 
    185185#Polygon designed to incorporate perth and Fremantle CBD
    186186poly_cbd = read_polygon(polygons_dir+'CBD_coastal.csv')
     
    188188
    189189#Polygon designed to incorporate rockingham and penguin island
    190 poly_penguin = read_polygon(polygons_dir+'rockingham_penguin.csv')
    191 res_penguin = 500*res_factor
    192 
    193 poly_geordie_bay = read_polygon(polygons_dir+'Geordie_Bay.csv')
     190poly_rockingham = read_polygon(polygons_dir+'rockingham_penguin.csv')
     191res_rockingham = 500*res_factor
     192
     193#Polygon designed to incorporate bottom of Garden Island for image verification
     194poly_garden = read_polygon(polygons_dir+'garden.csv')
     195res_garden = 500*res_factor
     196
     197poly_geordie_bay = read_polygon(polygons_dir+'geordie_bay.csv')
    194198res_geordie_bay = 500*res_factor
    195199
    196 poly_sorrento_gauge = read_polygon(polygons_dir+'Sorrento_gauge.csv')
     200poly_sorrento_gauge = read_polygon(polygons_dir+'sorrento_gauge.csv')
    197201res_sorrento_gauge = 500*res_factor
     202
     203
    198204#assert zone == refzone
    199205
    200206interior_regions = [[poly_internal_20_3,res_internal_20_3],[poly_cbd,res_cbd]
    201207                     ,[poly_garden_rockingham,res_garden_rockingham]
    202                      ,[poly_penguin,res_penguin],[poly_geordie_bay,res_geordie_bay]
     208                     ,[poly_rockingham,res_rockingham],[poly_geordie_bay,res_geordie_bay]
    203209                     ,[poly_sorrento_gauge,res_sorrento_gauge],[poly_rottnest_in, res_rottnest_in]
    204                      ,[poly_rottnest_ex, res_rottnest_ex]]
     210                     ,[poly_rottnest_ex, res_rottnest_ex], [poly_garden, res_garden]]
    205211
    206212boundary_tags={'back': [0,1,2], 'side': [3,7],'ocean': [4, 5, 6]}
     
    210216print 'min number triangles', trigs_min
    211217
    212 poly_mainland = read_polygon(polygons_dir+'Initial_Condition.csv')
     218poly_mainland = read_polygon(polygons_dir+'initial_condition.csv')
    213219
    214220###################################################################
Note: See TracChangeset for help on using the changeset viewer.