Ignore:
Timestamp:
Jul 25, 2008, 8:53:51 AM (14 years ago)
Author:
kristy
Message:

updated scripts for polyline boundary

File:
1 edited

Legend:

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

    r5558 r5569  
    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.shallow_water.data_manager import urs2sts,create_sts_boundary
     13from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    1214
    1315# file and system info
     
    3638tide = 0.6
    3739
    38 alpha = 0.2
     40alpha = 0.1
    3941friction=0.01
    4042starttime=0
    4143finaltime=80000
    4244export_cellsize=25
    43 setup='trial'
     45setup='final'
    4446source='polyline'
    4547
     
    160162buildings_filename_out = 'Perth_res_Project_modified.csv'
    161163
    162 ###############################
    163 # Domain definitions
    164 ###############################
    165 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
    166 
    167 poly_all = read_polygon(polygons_dir+'poly_all.csv')
    168 res_poly_all = 100000*res_factor
    169 
    170 #refzone = 50
     164
    171165
    172166###############################
    173167# Interior region definitions
    174168###############################
     169
     170#Initial bounding polygon for data clipping
     171poly_all = read_polygon(polygons_dir+'poly_all.csv')
     172res_poly_all = 100000*res_factor
    175173
    176174#Polygon designed by 20m contours, or 3km from the coastline
     
    217215                     ,[poly_rottnest_ex, res_rottnest_ex], [poly_garden, res_garden]]
    218216
    219 boundary_tags={'back': [0,1,2], 'side': [3,7],'ocean': [4, 5, 6]}
    220 
     217   
    221218trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
    222 
    223219print 'min number triangles', trigs_min
     220   
    224221
    225222poly_mainland = read_polygon(polygons_dir+'initial_condition.csv')
Note: See TracChangeset for help on using the changeset viewer.