Ignore:
Timestamp:
Jul 25, 2008, 4:17:26 PM (14 years ago)
Author:
kristy
Message:

Addition of Polyline information with sts files

File:
1 edited

Legend:

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

    r5526 r5575  
    1010#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
    1111from anuga.utilities.system_tools import get_user_name, get_host_name
     12from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1213
    1314# file and system info
     
    1617
    1718home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir   
     19muxhome = getenv('MUXHOME')
    1820user = get_user_name()
    1921host = get_host_name()
     
    3335scenario = 'busselton_tsunami_scenario'
    3436
    35 tide = 0
     37tide = 0.6
    3638
    3739alpha = 0.1
    3840friction=0.01
    3941starttime=0
    40 finaltime=50000 # changed for rev_number trial
     42finaltime=80000
    4143export_cellsize=25
    4244setup='final'
    43 source='exmouth'
     45source='polyline'
    4446
    4547
     
    135137#boundaries_source = '1'
    136138
    137 if source=='exmouth':
    138     boundaries_name = 'busselton_3103_30052008' #exmouth gun
    139     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
    140 
    141 if source=='test':
    142     boundaries_name = 'other' #exmouth gun
    143     boundaries_in_dir = anuga_dir+'boundaries'+sep
    144 
     139##if source=='exmouth':
     140##    boundaries_name = 'busselton_3103_30052008' # exmouth gun
     141##    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
     142##
     143##if source=='test':
     144##    boundaries_name = 'other' #exmouth gun
     145##    boundaries_in_dir = anuga_dir+'boundaries'+sep
     146##
    145147
    146148#boundaries locations
    147 boundaries_in_dir_name = boundaries_in_dir + boundaries_name
     149#boundaries_in_dir_name = boundaries_in_dir + boundaries_name
    148150boundaries_dir = anuga_dir+'boundaries'+sep
    149151boundaries_dir_name = boundaries_dir + scenario_name
     152boundaries_dir_mux = muxhome
    150153
    151154#output locations
    152155output_dir = anuga_dir+'outputs'+sep
    153 #output_build_time_dir = output_dir +'test'+sep + build_time + dir_comment + sep
    154156output_build_time_dir = output_dir +build_time + dir_comment + sep
    155 #output_run_time_dir = output_dir +run_time+dir_comment+sep
    156 output_run_time_dir = output_dir +'test'+sep + run_time + dir_comment +sep # take test out
     157output_run_time_dir = output_dir + run_time + dir_comment +sep
    157158output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
    158159
     
    170171
    171172###############################
    172 # Domain definitions
     173# Interior region definitions
    173174###############################
    174175
    175 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    176 
    177 # bounding polygon for study area
     176# Initial bounding polygon for data clipping
    178177poly_all = read_polygon(polygons_dir+'poly_all_extend.csv')
    179178res_poly_all = 100000*res_factor
    180 
    181 
    182 ###############################
    183 # Interior region definitions
    184 ###############################
    185179
    186180#digitized polygons
     
    209203
    210204
    211 boundary_tags={'back': [3, 4, 5, 6, 7],
    212                'side': [2, 8], 'ocean': [0, 1, 9]}
    213 
    214205trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
     206print 'min number triangles', trigs_min
    215207
    216208poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
    217 
    218 print 'min number triangles', trigs_min
    219209
    220210
Note: See TracChangeset for help on using the changeset viewer.