Changeset 4132


Ignore:
Timestamp:
Jan 5, 2007, 2:32:42 PM (17 years ago)
Author:
nick
Message:

update with new version

Location:
anuga_work/production/busselton_2006
Files:
2 added
1 edited

Legend:

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

    r4094 r4132  
    77import sys
    88from time import localtime, strftime, gmtime
    9 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
     9from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
    1010#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
    1111from anuga.utilities.system_tools import get_user_name
     
    2626build_time = time+'_build'
    2727run_time = time+'_run'
    28 print 'gtime: ', gtime
     28#print 'gtime: ', gtime
    2929
    3030tide = 0.6
     
    105105# bounding polygon for study area
    106106poly_all = read_polygon(polygons_dir+'poly_all.csv')
    107 res_poly_all = 100000
     107res_poly_all = 1000000
    108108
    109 
    110 '''
    111109###################################################################
    112110# Clipping regions for export to asc and regions for clipping data
     
    123121###############################
    124122
    125 # broome digitized polygons
    126 poly_busselton1 = read_polygon(polygons_dir+'buss_Local_Polygon_update.csv')
    127 res_busselton1 = 500000
    128 poly_busselton2 = read_polygon(polygons_dir+'buss_Close2_update.csv')
    129 res_busselton2 = 500000
    130 poly_busselton3 = read_polygon(polygons_dir+'buss_Coast_update.csv')
    131 res_busselton3 = 500000
     123#digitized polygons
     124poly_busselton1 = read_polygon(polygons_dir+'neg20_pos10_polygon.csv')
     125res_busselton1 = 100000
     126poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly_.csv')
     127res_busselton2 = 10000
    132128
    133129#plot_polygons([polyAll,poly_broome1,poly_broome2,poly_broome3],'boundingpoly2',verbose=False)
    134130
    135 interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2],
    136                     [poly_busselton3,res_busselton3]]
    137 #trigs_all = polygon_area(poly_all)/res_poly_all
    138 #trigs_bus1 = polygon_area(poly_busselton1)/res_busselton1
    139 #trigs_bus2 = polygon_area(poly_busselton2)/res_busselton2
    140 #trigs_bus3 = polygon_area(poly_busselton3)/res_busselton3
    141 #trigs_min = trigs_bound + trigs_pos + trigs_cbd + trigs_penguin
    142 
    143 #print 'Area of bounding poly', trigs_all
    144 #print 'Area of busselton1', trigs_bus1
    145 #print 'Area of busselton2', trigs_bus2
    146 #print 'Area of busselton3', trigs_bus3
    147 #print 'min number triangles', trigs_min
     131interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]]
    148132
    149133trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
     
    151135print 'min number triangles', trigs_min
    152136
    153 def number_mesh_triangles(interior_regions, bounding_poly, remainder_res):
    154     from anuga.utilities.polygon import polygon_area
    155    
    156     # TO DO check if any of the regions fall inside one another
    157     no_triangles = 0.0
    158     area = polygon_area(bounding_poly)
    159     for i,j in interior_regions:
    160         this_area = polygon_area(i)
    161         no_triangles += this_area/j
    162         area -= this_area
    163         print j, this_area/1000000., area/1000000.
    164     no_triangles += area/remainder_res
    165     return int(no_triangles/0.7)
    166 '''
    167137
Note: See TracChangeset for help on using the changeset viewer.