source: anuga_work/production/wollongong_2006/project.py @ 3624

Last change on this file since 3624 was 3624, checked in by ole, 18 years ago

Got flagstaff to work with parallel API

File size: 1.9 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2"""
3
4import sys
5from os import sep as s, environ
6from os.path import expanduser
7
8# Making assumptions about the location of scenario data
9scenario_dirname = 'new_south_wales'+s+'wollongong_tsunami_scenario_2006'+s+'anuga'
10
11# Filenames
12basename = 'flagstaff'
13demname = 'wollongong_10'
14building_footprints = 'flagstaff_building_footprints'
15
16# Parameters
17base_resolution = 5
18initial_sealevel = 0.0
19   
20# Derive subdirectories and filenames
21try:
22    home = environ['INUNDATIONHOME']
23except Exception, e:
24    print 'Environment variable INUNDATIONHOME was not found: %s' %e
25    home = expanduser('~')
26
27
28meshdir = home+s+scenario_dirname+s+'meshes'+s
29datadir = home+s+scenario_dirname+s+'topographies'+s
30outputdir = home+s+scenario_dirname+s+'output'+s
31polygondir = home+s+scenario_dirname+s+'polygons'+s
32gaugedir = home+s+scenario_dirname+s+'gauges'+s
33
34demname = datadir + demname
35meshname = meshdir + basename
36mesh_filename = meshname + '.msh'
37buildings_filename = meshdir + building_footprints + '.asc'
38
39outputname = outputdir + basename  #Used by post processing
40
41
42
43# Mesh data
44
45# Lower left corner
46xllcorner = 307000
47yllcorner = 6188940
48
49# Upper right corner
50xurcorner = 307788
51yurcorner = 6189509
52
53# Bounding polygon
54point1 = [xllcorner, (yllcorner+yurcorner)/2]
55point2 = [xllcorner + (xurcorner-xllcorner)/3, yurcorner]
56point3 = [xllcorner + 2*(xurcorner-xllcorner)/3, yurcorner]   
57point4 = [xurcorner, yllcorner + 2*(yurcorner-yllcorner)/3]       
58point5 = [xurcorner, yllcorner + (yurcorner-yllcorner)/3]   
59point6 = [xllcorner + 2*(xurcorner-xllcorner)/3, yllcorner]   
60
61bounding_polygon = [[xllcorner, yllcorner],   #sw
62                    point1,
63                    point2,
64                    point3,
65                    point4,
66                    point5,
67                    point6]   
68   
69boundary_tags = {'exterior': [0,1,5,6], 'side': [2,4], 'ocean': [3]}
70
71interior_regions = []
72       
Note: See TracBrowser for help on using the repository browser.