Changeset 6499 for anuga_work
- Timestamp:
- Mar 12, 2009, 1:40:40 PM (16 years ago)
- Location:
- anuga_work/production/hobart_2009
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2009/project.py
r6457 r6499 59 59 # Used in build_elevation.py 60 60 # Format for ascii grids, as produced in ArcGIS + a projection file 61 ascii_grid_filenames = ['hob3_project', # Topo 62 'grid_250m_project'] # Busselton Topo 61 ascii_grid_filenames = ['hob3_pro_extract', 'grid_250m_project', 'ldr_bbay', 62 'ldr_bellerive', 'ldr_brunyislandneck', 63 'ldr_calvertsbeach', 'ldr_carltonbeach', 64 'ldr_connelly', 'ldr_cremorne', 'ldr_dodge1', 65 'ldr_dodge2', 'ldr_kingston', 'ldr_lauderdale', 66 'ldr_linersfarne', 'ldr_oppossumbay', 'ldr_primrose', 67 'ldr_rosny', 'ldr_sandybay', 'ldr_selfspoint', 68 'ldr_sevenmile', 'ldr_southarm1', 'ldr_southarm2', 69 'ldr_sullivanscove', 'ldr_tranmere', 'tasdem_mask'] 63 70 64 71 # Format for point is x,y,elevation (with header) … … 75 82 # Format for points easting,northing (no header) 76 83 bounding_polygon_filename = 'bounding_polygon.csv' 77 bounding_polygon_maxarea = 100000 84 bounding_polygon_maxarea = 1000000 78 85 79 86 # INTERIOR REGIONS - for designing the mesh 80 87 # Used in run_model.py 81 88 # Format for points easting,northing (no header) 82 interior_regions_data = [] 89 interior_regions_data = [['aos1.csv', 1500], 90 ['aos2.csv', 1500], 91 ['sw.csv', 30000]] 92 PriorityArea_filename = 'PriorityAreas.csv' 83 93 84 94 # LAND - used to set the initial stage/water to be offcoast only … … 237 247 mux_input = join(event_folder, mux_input_filename) 238 248 249 #Multiple polygons in one CSV file to make internal polygons 250 PriorityAreas = join(polygons_folder, PriorityArea_filename) 251 252 -
anuga_work/production/hobart_2009/setup_model.py
r6457 r6499 11 11 from os.path import join, exists 12 12 from anuga.utilities.polygon import read_polygon, number_mesh_triangles 13 13 from anuga.shallow_water.data_manager import csv2polygons 14 14 import project 15 15 … … 157 157 # Create list of interior polygons with scaling factor 158 158 project.interior_regions = [] 159 160 polygons, maxareas = csv2polygons(project.PriorityAreas) 161 print 'Creating %d internal polygons' % len(polygons) 162 #def create_polygon_function(polygons, geo_reference=None): 163 project.interior_regions = [] 164 for i, key in enumerate(polygons): 165 if i%100==0: print i 166 poly = polygons[key] 167 maxarea = float(maxareas[key]) 168 project.interior_regions.append([poly, 169 maxarea*project.scale_factor]) 170 159 171 for filename, maxarea in project.interior_regions_data: 160 172 polygon = read_polygon(join(project.polygons_folder, filename))
Note: See TracChangeset
for help on using the changeset viewer.