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

Last change on this file since 4003 was 4003, checked in by sexton, 17 years ago

including data and supporting files for Cairns demo + updating scripts

File size: 1.5 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
17polyAll = read_polygon('extent.csv')
18
19print 'Area of bounding polygon', polygon_area(polyAll)/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
33plot_polygons([polyAll,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 = 356300
43eastingmax = 398500
44northingmin = 8107850
45northingmax = 8148700
46
47
48slump_origin = [451871, 8128376] # move onto the continental shelf, depth = 500
49slump_depth = 500.
50
51gauge_filename = 'gauges.csv'
Note: See TracBrowser for help on using the repository browser.