source: production/karratha_2005/project.py @ 1987

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

Updates for Karratha case study:

Gauges for inundated buildings
Use built-in cache in file_function

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