source: production/sydney_2006/project.py @ 2247

"""Common filenames and locations for topographic data, meshes and outputs.
"""


from os import sep
from os.path import expanduser
import sys

#Making assumptions about the location of scenario data
scenario_dir_name = 'sydney_tsunami_scenario_2006'
basename = 'sydney_2006_100' # get from Neil (DEM or topo data)

if sys.platform == 'win32':
    home = '..\..\..\..'     #Sandpit's parent dir
else:    
    home = expanduser('~')


#Derive subdirectories and filenames
meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
outputdir = home+sep+scenario_dir_name+sep+'output'+sep
#boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep

meshname = meshdir + basename
demname = datadir + basename
#boundaryname = boundarydir + boundary_basename
outputname = outputdir + basename  #Used by post processing


#Georeferencing
from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees

#Origin of existing dem (FIXME: Temporary measure)
mesh_origin = (56, 314036, 6224951)  # input from Neil's data

# define clipping polygon
#south = degminsec2decimal_degrees(-20,45,0)
#north = degminsec2decimal_degrees(-20,15,0)
#west = degminsec2decimal_degrees(116,30,0)
#east = degminsec2decimal_degrees(117,0,0)

south = degminsec2decimal_degrees(-34,05,0)
north = degminsec2decimal_degrees(-33,45,0)
west = degminsec2decimal_degrees(151,0,0)
east = degminsec2decimal_degrees(151,30,0)

p0 = [south, west]
p1 = [south, east]
p2 = [north, east]
p3 = [north, west]
    
polygonall = [p0, p1, p2, p3]
refzone = 56

# setting up problem area for doing just around the harbour
hsouth = degminsec2decimal_degrees(-33,54,0)
hnorth = degminsec2decimal_degrees(-33,48,0)
hwest = degminsec2decimal_degrees(151,0,0)
heast = degminsec2decimal_degrees(151,30,0)

hp0 = [hsouth, hwest]
hp1 = [hsouth, heast]
hp2 = [hnorth, heast]
hp3 = [hnorth, hwest]
polygon_h = [hp0, hp1, hp2, hp3]

#Interior regions - the Harbour
harbour_south = degminsec2decimal_degrees(-33,53,0)
harbour_north = degminsec2decimal_degrees(-33,47,0)
harbour_west = degminsec2decimal_degrees(151,5,0)
harbour_east = degminsec2decimal_degrees(151,19,0)

k0 = [harbour_south, harbour_west]
k1 = [harbour_south, harbour_east]
k2 = [harbour_north, harbour_east]
k3 = [harbour_north, harbour_west]    

harbour_polygon = [k0, k1, k2, k3]

# setting up problem area for doing just around Botany Bay
bsouth = degminsec2decimal_degrees(-33,56,0)
bnorth = degminsec2decimal_degrees(-34,3,0)
bwest = degminsec2decimal_degrees(151,0,0)
beast = degminsec2decimal_degrees(151,30,0)

bp0 = [bsouth, bwest]
bp1 = [bsouth, beast]
bp2 = [bnorth, beast]
bp3 = [bnorth, bwest]
polygon_bb = [bp0, bp1, bp2, bp3]

#Interior region - Botany Bay
botanybay_south = degminsec2decimal_degrees(-33,58,0)
botanybay_north = degminsec2decimal_degrees(-34,1,0)
botanybay_west = degminsec2decimal_degrees(151,5,0)
botanybay_east = degminsec2decimal_degrees(151,18,0)

j0 = [botanybay_south, botanybay_west]
j1 = [botanybay_south, botanybay_east]
j2 = [botanybay_north, botanybay_east]
j3 = [botanybay_north, botanybay_west]    

botanybay_polygon = [j0, j1, j2, j3]