source: anuga_core/documentation/user_manual/demos/cairns/project.py @ 4856

Last change on this file since 4856 was 4856, checked in by sexton, 16 years ago

rename plot_polygons_points to original name of plot_polygons

File size: 1.6 KB
Line 
1# -*- coding: cp1252 -*-
2"""Common filenames and locations for topographic data, meshes and outputs.
3"""
4
5from os import sep, environ, getenv, getcwd
6from os.path import expanduser
7import sys
8from time import localtime, strftime, gmtime
9from anuga.utilities.polygon import read_polygon, plot_polygons, \
10                                    polygon_area, is_inside_polygon
11
12###############################
13# Domain definitions
14###############################
15
16# bounding polygon for study area
17bounding_polygon = read_polygon('extent.csv')
18
19print 'Area of bounding polygon', polygon_area(bounding_polygon)/1000000.0
20
21###############################
22# Interior region definitions
23###############################
24
25# interior polygons
26poly_cairns = read_polygon('cairns.csv')
27poly_island0 = read_polygon('islands.csv')
28poly_island1 = read_polygon('islands1.csv')
29poly_island2 = read_polygon('islands2.csv')
30poly_island3 = read_polygon('islands3.csv')
31poly_shallow = read_polygon('shallow.csv')
32
33#plot_polygons([bounding_polygon,poly_cairns,poly_island0,poly_island1,\
34#               poly_island2,poly_island3,poly_shallow],\
35#              'boundingpoly',verbose=False)
36
37###################################################################
38# Clipping regions for export to asc and regions for clipping data
39###################################################################
40
41# exporting asc grid
42eastingmin = 363000
43eastingmax = 418000
44northingmin = 8026600
45northingmax = 8145700
46
47
48slide_origin = [451871, 8128376] # move onto the continental shelf, depth = 500
49slide_depth = 500.
50
51gauge_filename = 'gauges.csv'
Note: See TracBrowser for help on using the repository browser.