Ignore:
Timestamp:
May 30, 2008, 10:52:58 AM (16 years ago)
Author:
kristy
Message:

update scripts

File:
1 edited

Legend:

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

    r5160 r5381  
    7272# AHO + DPI data
    7373coast_name = 'coastline_excluding_beach_survey'
    74 coast_name1 = 'beach_survey_final'
    75 offshore_name = 'busselton'
     74coast_name1 = 'bussleton_beach_data'
     75offshore_name = 'Busselton_bathymetry'
    7676
    7777#final topo name
     
    8484topographies_dir = home+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
    8585
    86 #input topo file location
    87 onshore_in_dir_name = topographies_in_dir + onshore_name
    88 #island_in_dir_name = topographies_in_dir + island_name
     86# input topo file location
     87onshore_in_dir_name = topographies_in_dir + onshore_name #topo
    8988
    90 coast_in_dir_name = topographies_in_dir + coast_name
    91 coast_in_dir_name1 = topographies_in_dir + coast_name1
    92 offshore_in_dir_name = topographies_in_dir + offshore_name
     89coast_in_dir_name = topographies_in_dir + coast_name #coastline
     90coast_in_dir_name1 = topographies_in_dir + coast_name1 #beach survey
     91offshore_in_dir_name = topographies_in_dir + offshore_name #bathymetry
    9392
     93#output to anuga from build file
    9494onshore_dir_name = topographies_dir + onshore_name
    95 #island_dir_name = topographies_dir + island_name
     95
    9696coast_dir_name = topographies_dir + coast_name
    9797coast_dir_name1 = topographies_dir + coast_name1
     
    177177res_busselton1 = 10000*res_factor
    178178
    179 poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly_.csv')
     179poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly.csv')
    180180res_busselton2 = 500*res_factor
    181181
     
    183183
    184184interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]]
     185print 'Hello',interior_regions
    185186
    186 boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 15],
    187                'side': [10, 14], 'ocean': [11, 12, 13]}
     187boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14],
     188               'side': [10], 'ocean': [11, 12, 13]}
    188189
    189190trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
    190191
    191 poly_mainland=read_polygon(polygons_dir+'initial_condition_south.csv')
     192poly_mainland=read_polygon(topographies_in_dir+'initial_condition_south.csv')
    192193
    193194print 'min number triangles', trigs_min
Note: See TracChangeset for help on using the changeset viewer.