Changeset 4399
- Timestamp:
- Apr 18, 2007, 1:19:31 PM (17 years ago)
- Location:
- anuga_work/production/exmouth_2006
- Files:
-
- 1 added
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/exmouth_2006/project.py
r4228 r4399 27 27 28 28 # onshore data provided by WA DLI 29 #onshore_name =# original29 onshore_name = 'first' # original 30 30 31 31 # AHO + DPI data … … 130 130 131 131 # bounding polygon for study area 132 poly All = read_polygon(polygondir+'extent.csv')132 poly_all = read_polygon(polygondir+'extent_zone50a.csv') 133 133 134 134 # plot bounding polygon and make sure BC info surrounds it 135 #plot_polygons([poly All, poly_bc],'boundingpoly',verbose=False)136 print 'Area of bounding polygon', polygon_area(poly All)/1000000.0135 #plot_polygons([poly_all, poly_bc],'boundingpoly',verbose=False) 136 print 'Area of bounding polygon', polygon_area(poly_all)/1000000.0 137 137 138 138 ################################################################### … … 149 149 # Interior region definitions 150 150 ############################### 151 151 ''' 152 152 # exmouth digitized polygons 153 153 poly_exmouth1 = read_polygon(polygondir+'exmouth_Local_Polygon_update.csv') … … 156 156 #poly_exmouth4 = read_polygon(polygondir+'Cable_Beach_revised.csv') 157 157 158 plot_polygons([poly All,poly_exmouth1,poly_exmouth2,poly_exmouth3],'boundingpoly2',verbose=False)158 plot_polygons([poly_all,poly_exmouth1,poly_exmouth2,poly_exmouth3],'boundingpoly2',verbose=False) 159 159 print 'Area of local polygon', polygon_area(poly_exmouth1)/1000000.0 160 160 print 'Area of close polygon', polygon_area(poly_exmouth2)/1000000.0 … … 179 179 no_triangles += area/remainder_res 180 180 return int(no_triangles/0.7) 181 '''
Note: See TracChangeset
for help on using the changeset viewer.