source: anuga_core/documentation/user_manual/examples/project.py @ 3869

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

data update for sydney demo

File size: 1.9 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2Also includes origin for slump scenario.
3"""
4
5from os import sep, environ
6from os.path import expanduser
7import sys
8
9# We'll do this for the final thing
10home = '.'
11
12if sys.platform == 'win32':
13    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
14else:   
15    home = expanduser('~')
16   
17basename = 'Sydney'
18
19#Setting the location of scenario data
20scenario_dir_name = 'sydney_tsunami_example_2006'
21# revised 100m data
22#coarsename = 'bathyland100'
23# revised 25m data
24#finename = 'bathy_dem25'
25# revised 25m data
26dem_name = 'sydney_MGA56'
27
28home += sep + 'data' + sep + 'new_south_wales'
29
30#Derive subdirectories and filenames
31meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
32datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
33outputdir = home+sep+scenario_dir_name+sep+'output'+sep
34polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
35
36meshname = meshdir + basename
37#coarsedemname = datadir + coarsename
38#finedemname = datadir + finename
39demname = datadir + dem_name
40combineddemname = datadir + 'sydneytopo'
41
42refzone = 56
43print 'refzone', refzone
44
45# demo poly
46j0 = [385000, 6280000]
47j1 = [360000, 6272500]
48j2 = [335000, 6272500]
49j3 = [330000, 6265000]
50j31 = [325000, 6260000]
51j4 = [316000, 6260000]
52j5 = [316000, 6247000]
53j6 = [350000, 6247000]
54j7 = [385000, 6238000]
55
56demopoly = [j0, j1, j2, j3, j31, j4, j5, j6, j7]
57
58from utilities.polygon import read_polygon
59polygonptsfile4 = polygondir + 'poly1'
60polygonptsfile0 = polygondir + 'poly2'
61polygonptsfile1 = polygondir + 'poly3'
62polygonptsfile2 = polygondir + 'poly4'
63polygonptsfile3 = polygondir + 'poly5'
64northern_polygon = read_polygon(polygonptsfile0 + '.csv')
65manly_polygon = read_polygon(polygonptsfile1 + '.csv')
66harbour_polygon = read_polygon(polygonptsfile2 + '.csv')
67southern_polygon = read_polygon(polygonptsfile3 + '.csv')
68top_polygon = read_polygon(polygonptsfile4 + '.csv')
69
70slump_origin = [372500.0, 6255000.0] #Absolute UTM
Note: See TracBrowser for help on using the repository browser.