source: production/hobart_2006/project.py @ 3559

# -*- coding: cp1252 -*-
"""Common filenames and locations for topographic data, meshes and outputs.
"""



from os import sep, environ, getenv, getcwd
from os.path import expanduser
import sys
from time import localtime, strftime, gmtime
                
#Making assumptions about the location of scenario data
state = 'tasmania'
scenario_dir_name = 'hobart_tsunami_scenario_2006'

# data provided by Tas SES and checked by NM&I
onshore_name = 'hobart_onshore' 
offshore_name = 'hobart_offshore_points'

# developed by NM&I
coast_name = 'hobart_coast'

boundary_basename = 'SU-AU_clip'

#swollen/ all data output
basename = 'source'

codename = 'project.py'

if sys.platform == 'win32':
    home = getenv('INUNDATIONHOME')
    user = getenv('USERPROFILE')    
else:    
    home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation')     
    user = getenv('LOGNAME')
    print 'USER:', user

#Derive subdirectories and filenames
#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
local_time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir

meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep
datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep
gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
#output dir without time
outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep
outputtimedir = outputdir + local_time + sep
polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep

gauge_filename = gaugedir + 'hobart_gauges.xya'
buildings_filename = gaugedir + 'hobart_res.csv'
buildings_filename_damage_out = 'hobart_res_modified.csv'

gaugetimeseries = gaugedir + 'hobart'

# boundary source data
#MOST_dir = 'f:'+sep+'3'+sep+'ehn'+sep+'users'+sep+'davidb'+sep+'tsunami'+sep+'WA_project'+sep+'SU-AU_90'+sep+'most_2'+sep+'detailed'+sep

codedir = getcwd()+sep
                                
codedirname = codedir + 'project.py'

meshname = outputtimedir + 'mesh_' + basename
coarsedemname = datadir + coarsename

onshore_dem_name = datadir + onshore_name
offshore_dem_name = datadir + offshore_name
coast_dem_name = datadir + coast_name
combined_dem_name = datadir + 'hobart_combined_elevation'

outputname = outputtimedir + basename  #Used by post processing

from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
from anuga.coordinate_transforms.redfearn import redfearn

# bounding box
south = degminsec2decimal_degrees(-43,45,0)
north = degminsec2decimal_degrees(-42,30,0)
west = degminsec2decimal_degrees(146,45,0)
east = degminsec2decimal_degrees(148,15,0)

'''
# region for visualisation
eminviz = 260000
emaxviz = 320000
nminviz = 7590000
nmaxviz = 7630000
'''
# region to export

e_min_area = 300000
e_max_area = 310000
n_min_area = 7600000
n_max_area = 7610000


#Main Domain of Hobart: first run JS 15/9/06
zone, e0, n0 = redfearn(south, west)
zone, e1, n1 = redfearn(south, east)
zone, e2, n2 = redfearn(north, east)
zone, e3, n3 = redfearn(north, west)
refzone = zone
d0 = [e0, n0]
d1 = [e1, n1]
d2 = [e2, n2]
d3 = [e3, n3]

polyAll = [d0, d1, d2, d3]

#Interior region - Hobart city area
#from anuga.utilities.polygon import read_polygon
polygonptsfile0 = polygondir + 'poly1’
polygonptsfile1 = polygondir + 'poly2’
northern_polygon = read_polygon(polygonptsfile0 + '.csv')
southern_polygon = read_polygon(polygonptsfile1 + '.csv')

i0 = [304000, 7607000]
i1 = [302000, 7605000]
i2 = [304000, 7603000]
i3 = [307000, 7602000]
i4 = [309000, 7603000]
i5 = [307000, 7606000]

poly_hobart = [i0, i1, i2, i3, i4, i5]

#med res around Hobart
l0 = [300000, 7610000]
l1 = [285000, 7600000]
l2 = [300000, 7597500]
l3 = [310000, 7600000]
l4 = [315000, 7610000]

poly_coast = [l0, l1, l2, l3, l4]

#general coast and local area to Hobart region
m0 = [270000, 7581000]
m1 = [300000, 7591000]
m2 = [339000, 7610000]
m3 = [330000, 7630000]
m4 = [290000, 7640000]
m5 = [260000, 7600000]

poly_region = [m0, m1, m2, m3, m4, m5]

# Bruny Island
poly_bruny = [b0, b1, b2, b3, b4]