source: production/karratha_2005/project.py @ 1877

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

Karratha update

File size: 3.1 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2"""
3
4
5from os import sep
6from os.path import expanduser
7import sys
8
9#Making assumptions abot the location of scenario data
10scenario_dir_name = 'karratha_tsunami_scenario_2005'
11basename = 'karratha_100m'
12boundary_basename = 'SU-AU'
13
14
15if sys.platform == 'win32':
16    home = '..\..\..\..\..'     #Sandpit's parent dir
17else:   
18    home = expanduser('~')
19
20
21#Derive subdirectories and filenames
22meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
23datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
24outputdir = home+sep+scenario_dir_name+sep+'output'+sep
25boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep
26
27meshname = meshdir + basename
28demname = datadir + basename
29boundaryname = boundarydir + boundary_basename
30outputname = outputdir + basename
31
32#Georeferencing
33from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees
34
35#Origin of existing dem (FIXME: Temporary measure)
36#mesh_origin = (50, 421544.35127423, 7677669.5257159)  #250m
37mesh_origin = (50, 420468.31429902, 7677669.5257159)  #100m
38
39
40#south = degminsec2decimal_degrees(-20,45,0)
41#north = degminsec2decimal_degrees(-20,15,0)
42#west = degminsec2decimal_degrees(116,30,0)
43#east = degminsec2decimal_degrees(117,0,0)
44
45south = degminsec2decimal_degrees(-20,55,0)
46north = degminsec2decimal_degrees(-20,15,0)
47west = degminsec2decimal_degrees(116,17,0)
48east = degminsec2decimal_degrees(117,10,0)
49
50p0 = [south, degminsec2decimal_degrees(116,32,0)]
51p1 = [south, west]
52p2 = [degminsec2decimal_degrees(-20,23,0), west]
53p3 = [north, degminsec2decimal_degrees(116,45,0)]
54p4 = [north, degminsec2decimal_degrees(117,0,0)]
55p5 = [p2[0], degminsec2decimal_degrees(117,8,0)]
56p6 = [degminsec2decimal_degrees(-20,30,0), east]
57p7 = [degminsec2decimal_degrees(-20,38,0), east]
58p8 = [south, east]
59   
60polygon = [p0, p1, p2, p3, p4, p5, p6, p7, p8]
61refzone = 50
62
63
64#Interior regions
65karratha_south = degminsec2decimal_degrees(-20,44,0)
66karratha_north = degminsec2decimal_degrees(-20,42,0)
67karratha_west = degminsec2decimal_degrees(116,48,0)
68karratha_east = degminsec2decimal_degrees(116,53,30)
69
70k0 = [karratha_south, karratha_west]
71k1 = [karratha_south, karratha_east]
72k2 = [karratha_north, karratha_east]
73k3 = [karratha_north, karratha_west]   
74
75karratha_polygon = [k0, k1, k2, k3]
76
77
78#Interior regions
79dampier_south = degminsec2decimal_degrees(-20,40,0)
80dampier_north = degminsec2decimal_degrees(-20,38,10)
81dampier_west = degminsec2decimal_degrees(116,43,0)
82dampier_east = degminsec2decimal_degrees(116,45,0)
83
84d0 = [dampier_south, dampier_west]
85d1 = [dampier_south, dampier_east]
86d2 = [dampier_north, dampier_east]
87d3 = [dampier_north, dampier_west]   
88
89dampier_polygon = [d0, d1, d2, d3]
90
91
92#Interior regions
93refinery_south = degminsec2decimal_degrees(-20,37,50)
94refinery_north = degminsec2decimal_degrees(-20,36,0)
95refinery_west = degminsec2decimal_degrees(116,44,0)
96refinery_east = degminsec2decimal_degrees(116,46,10)
97
98d0 = [refinery_south, refinery_west]
99d1 = [refinery_south, refinery_east]
100d2 = [refinery_north, refinery_east]
101d3 = [refinery_north, refinery_west]   
102
103refinery_polygon = [d0, d1, d2, d3]
Note: See TracBrowser for help on using the repository browser.