source: production/pt_hedland_2006/project.py @ 3340

"""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_orig'

# offshore data from AHO fairsheets
offshore_name2 = 'pt_hedland_offshore_points_fairsheet_orig'

# coastline developed from aerial photography and 1.5m DLI contour
#coast_name = 'pt_hedland_coastline_points_dted2_new'
coast_name = 'coast_with_extracted_z'


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 = 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'

meshname = meshdir + basename
#onshore_dem_name = datadir + onshore_name_dted
onshore_dem_name = datadir + onshore_name_dli
offshore_dem_name1 = datadir + offshore_name1
offshore_dem_name2 = datadir + offshore_name2
coast_dem_name = datadir + coast_name
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 = 648000
e_max_area = 675000
n_min_area = 7745000
n_max_area = 7761000

export_region = [[e_min_area,n_min_area],[e_min_area,n_max_area],[e_max_area,n_max_area],[e_max_area,n_min_area]]
                 
refzone = 50 

from coordinate_transforms.redfearn import redfearn
# boundary up to 50 m contour
lat1_50 = degminsec2decimal_degrees(-19,20,0)
lat2_50 = degminsec2decimal_degrees(-19,30,0)
lat3_50 = degminsec2decimal_degrees(-19,45,0)
lon1_50 = degminsec2decimal_degrees(119,05,0)
lon2_50 = degminsec2decimal_degrees(118,20,0)
lon3_50 = degminsec2decimal_degrees(117,45,0)
z, easting, northing = redfearn(lat1_50, lon1_50)
d0_50 = [easting, northing]
z, easting, northing = redfearn(lat2_50, lon2_50)
d1_50 = [easting, northing]
z, easting, northing= redfearn(lat3_50, lon3_50)
d2_50 = [easting, northing]

d4_50 = [285000, 7585000]
d6_50 = [330000, 7605000]
#bounding_poly50 = [p0_50, p1_50, p2_50, d6_50, d5, d4_50]

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]
#polyAll = [d0_50, d1_50, d2_50, d4, d5, d6]
# from Hamish
h0=[629262.17, 7747205.47]
h1=[552686.00, 7871579.99] #d3
h2=[658264.00, 7926314.00] #d2
h3=[710986.99, 7925796.99] #d1
h4=[763851.99, 7934357.99] #d0
h5=[701485.21, 7770656.86]
h6=[698273.75, 7762227.38]
h7=[698194.23, 7762018.65]
h8=[691627.41, 7744781.98]
h9=[679220.75, 7743604.59]
h10=[653512.59, 7740528.56]
h11=[634777.71, 7738247.17]
h12=[629443.86, 7746910.37]
h13=[629396.84, 7746986.75]
h14=[629352.32, 7747059.06]
h15=[629276.24, 7747182.63]
h16=[629262.17, 7747205.47] #repeat of h0
# using Hamish's new bounding polygon
#polyAll = [d0_50, d1_50, d2_50, h16,h15,h14,h13,h12,h11,h10,h9,h8,h7,h6,h5]
polyAll = [d0_50, d1_50, d2_50, h16,h11,h8,h6, h5]

#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]