source: production/pt_hedland_2006/project.py @ 3272

"""Common filenames and locations for topographic data, meshes and outputs.
Also includes origin for slump scenario.
"""

from os import sep, environ, getenv, getcwd
from os.path import expanduser, basename
from utilities.polygon import read_polygon
import sys
from pmesh.create_mesh import convert_points_from_latlon_to_utm
from coordinate_transforms.redfearn import degminsec2decimal_degrees
from time import localtime, strftime
from geospatial_data import *

#Making assumptions about the location of scenario data
state = 'western_australia'
scenario_dir_name = 'pt_hedland_tsunami_scenario_2006'

# onshore data from 30m DTED level 2
onshore_name_dted = 'pt_hedland_onshore_30m_dted' # get from Neil/Ingo (DEM or topo data)
onshore_name_dli = 'pt_hedland_onshore_20m_dli' # get from Neil/Ingo (DEM or topo data)

# offshore data from GA digitised charts
offshore_name1 = 'pt_hedland_offshore_points'

# offshore data from AHO fairsheets
offshore_name2 = 'pt_hedland_offshore_points_fairsheet'

# coastline developed from 30m DTED
coast_name = 'pt_hedland_coastline_points.xya'

boundary_basename = 'SU-AU_clip'

#swollen/ all data output
basename = 'source'

codename = 'project.py'

if sys.platform == 'win32':
    home = getenv('INUNDATIONHOME') #Sandpit's parent dir
#    python_home = getenv('PWD')     
    #user = basename(getenv('USERPROFILE'))
    #print 'USER:', user
else:
    # original
    #home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation')
    # update to perlite 2
    home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'2'+sep+'cit'+sep+'inundation'+sep+'data')     
    user = getenv('LOGNAME')

#Derive subdirectories and filenames
time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
outputtimedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep+time+sep

#print 'outputtimedir', outputtimedir
#meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
#datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
#gaugedir = home+sep+scenario_dir_name+sep+'gauges'+sep
#polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
#boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep
#outputdir = home+sep+scenario_dir_name+sep+'outputs'+sep
#tidedir = home+sep+scenario_dir_name+sep+'tide_data'+sep

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
outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep
tidedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'tide_data'+sep

gauge_filename = gaugedir + 'gauge_location_port_hedland.csv'
buildings_filename = gaugedir + 'pt_hedland_res.csv'
community_filename = gaugedir + 'CHINS_v2.csv'
community_scenario = gaugedir + 'community_pt_hedland.csv'
tidal_filename = tidedir + 'pt_hedland_tide.txt'
coast_filename = datadir + coast_name

meshname = meshdir + basename
onshore_dem_name = datadir + onshore_name_dted
offshore_dem_name1 = datadir + offshore_name1
offshore_dem_name2 = datadir + offshore_name2
combined_dem_name = datadir + 'pt_hedland_combined_elevation'
outputname = outputtimedir + basename  #Used by post processing

# clipping region to make DEM (pts file) from onshore data
eastingmin = 594000
eastingmax = 715000
northingmin = 7720000
northingmax = 7880000

# for ferret2sww
south = degminsec2decimal_degrees(-20,30,0)
north = degminsec2decimal_degrees(-17,10,0)
west = degminsec2decimal_degrees(117,00,0)
east = degminsec2decimal_degrees(120,00,0)

# region to export (used from export_results.py)
e_min_area = 633000
e_max_area = 690000
n_min_area = 7740000
n_max_area = 7761000

refzone = 50 

# bounding polygon provided by Hamish
#d0 = [818732.55, 8062768.27]
"""
old 
d0 = [755000.0, 8025000.0]
d1 = [708940.32, 7750510.33]
d2 = [656561.15, 7732615.11]
d3 = [604415.81, 7733013.56]
d4 = [517682.34, 7899310.22]
"""

d0 = [763852.0, 7934358.0]
d1 = [710987.0, 7925797.0]
d2 = [658264.0, 7926314.0]
d3 = [552686.0, 7871580.0]
#d4 = [604415.81, 7733013.56]
d4 = [638000.0, 7733013.56]
#d5 = [656561.15, 7732615.11]
d5 = [662000.0, 7732615.11]
#d6 = [708940.32, 7750510.33]
d6 = [690000.0, 7740510.33]

polyAll = [d0, d1, d2, d3, d4, d5, d6]

#Interior region - Pt Hedland town
i0 = [668000, 7757000]
i1 = [659000, 7755000]
i2 = [660000, 7749000]
i3 = [667000, 7746000]
i4 = [678000, 7751000]

poly_pt_hedland = [i0, i1, i2, i3, i4]

#Are there other significant features?
j0 = [670000, 7760000]
j1 = [633000, 7745000]
j2 = [665000, 7743000]
j3 = [690000, 7755000]

poly_region = [j0, j1, j2, j3]