source: trunk/anuga_work/development/gong_2008/project_slide.py @ 7884

# -*- 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
from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon

if sys.platform == 'win32':
    home = getenv('INUNDATIONHOME')
    user = getenv('USERPROFILE')

else:    
    home = getenv('INUNDATIONHOME', sep+'d'+sep+'xrd'+sep+'gem'+sep+'2'+sep+'ramp'+sep+'risk_assessment_methods_project'+sep+'inundation')     
    user = getenv('LOGNAME')
    print 'USER:', user

# INUNDATIONHOME is the inundation directory, not the data directory.
home += sep +'data'

# Bridgette to vary this parameter from 1 to 25 in steps of ?
scale = 2.0 # None 

#Making assumptions about the location of scenario data
state = 'new_south_wales'
scenario_dir_name = 'wollongong_tsunami_scenario_2006'

# onshore data provided by NSW LPI, Krishna merged with existing
# PMD data to make 100m DEM
on_offshore10_name = 'wollongong_10'
nsw100 = 'dem_utm'

# survey data
offshore_name1 = 'surveyAreaA'

# AHO data
offshore_name4 = '1000003611export'
offshore_name5 = '1000003613export'
offshore_name6 = '1000003614export'
offshore_name7 = '1000003627export'
offshore_name8 = '1000003628export'
offshore_name9 = 'AHDexport'

#swollen/ all data output
basename = 'source'
codename = 'project_slide.py'

#Derive subdirectories and filenames
local_time = strftime('%Y%m%d_%H%M%S',gmtime()) 
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+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep
dir_comment = '_scale_' + str(scale)
outputtimedir = outputdir + local_time + dir_comment + sep
polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep

#gauge_filename = gaugedir + 'gong_gauges.csv'
gauge_filename = gaugedir + 'gauges.csv'
codedir = getcwd()+sep                            
codedirname = codedir + 'project_slide.py'
meshname = outputtimedir + 'mesh_' + basename

# Necessary if using point datasets, rather than grid
on_offshore10_dem_name = datadir + on_offshore10_name
nsw_dem_name = datadir + nsw100
offshore_dem_name1 = datadir + offshore_name1
offshore_dem_name4 = datadir + offshore_name4
offshore_dem_name5 = datadir + offshore_name5
offshore_dem_name6 = datadir + offshore_name6
offshore_dem_name7 = datadir + offshore_name7
offshore_dem_name8 = datadir + offshore_name8
offshore_dem_name9 = datadir + offshore_name9
combined_dem_name   = datadir + 'gong_combined_elevation'

###############################
# Domain definitions
###############################

# bounding polygon for study area
polyAll = read_polygon(polygondir+'extentA_up.csv')

print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0

###############################
# Interior region definitions
###############################

# areaA digitized polygons
poly_local = read_polygon(polygondir+'local2.csv')
poly_gong = read_polygon(polygondir+'gong3.csv')
poly_southgong = read_polygon(polygondir+'south_gong.csv')

print 'Area of local polygon', polygon_area(poly_local)/1000000.0

###################################################################
# Clipping regions for export to asc and regions for clipping data
###################################################################

poly_10mclip = read_polygon(polygondir+'10mclip.csv')
poly_surveyclip = read_polygon(polygondir+'surveyclip.csv')

# clipping for nsw grid for conversion
eastingmin_nsw = 272900
eastingmax_nsw = 463550
northingmin_nsw = 6112950
northingmax_nsw = 6269750

# exporting asc grid
##eastingmin = 300000
##eastingmax = 314000
##northingmin = 6167680
##northingmax = 6202130

eastingmin = 300000
eastingmax = 314000
northingmin = 6167680
northingmax = 6202130

# exporting asc grid for Amy
##eastingmin = 282000
##eastingmax = 306500
##northingmin = 6131750
##northingmax = 6170850

###################################################################
# Slide characteristics
###################################################################

# historical slides
slide_origin_bulli = [370110,6189489]
slide_origin_shovel = [354753,6187397]
#slide_origin_yacaaba = [463012,6334735] original work was actually Birubi
slide_origin_yacaaba = [476465, 6362424]
slide_origin_birubi = [462267,6335100]

# potential slides: these correspond approximately to modelled depths
slide_origin_bulli_b = [405617,6249852]
slide_origin_bulli_c0 = [438583,6308510]
slide_origin_bulli_c1 = [465283,6335210]
slide_origin_bulli_c2 = [483783,6358810]
slide_origin_shovel_b = [396262,6252144]
slide_origin_shovel_c0 = [431683,6311110]
slide_origin_shovel_c1 = [458694,6340352]
slide_origin_shovel_c2 = [475589,6363030]
slide_origin_yacaaba_a = [359796,6190943]
slide_origin_yacaaba_b = [397330,6251857]
slide_origin_yacaaba_c0 = [431283,6312110]
slide_origin_yacaaba_c1 = [459983,6339110]
slide_origin_birubi_a = [362811,6191183]
slide_origin_birubi_b = [403573,6250268]
slide_origin_birubi_c0 = [437183,6309410]
slide_origin_birubi_c2 = [483392,6359044]

#bulli_depth = 2087.0
bulli_depth = 1470.0
bulli_length = 16840.0
bulli_thickness = 424.0
bulli_width = 8860.0
bulli_density = 1.46
bulli_slope = 4.0
bulli_alpha = 126.0 - 90.0

#shovel_depth = 968.0
shovel_depth = 877.0
shovel_length = 13500.0
shovel_thickness = 165.0
shovel_width = 4350.0
shovel_density = 1.49
shovel_slope = 4.0
shovel_alpha = 118.0 - 90.0

#yacaaba_depth = 1119.0
yacaaba_depth = 938.
yacaaba_length = 4189.
yacaaba_thickness = 53.
yacaaba_width = 2898.
yacaaba_density = 1.48
yacaaba_slope = 2.3
yacaaba_alpha = 133.0 - 90.0
 
birubi_depth = 1320.
birubi_length = 9903.0
birubi_width = 4150.
birubi_density = 1.48
birubi_slope = 3.7
birubi_alpha = 133.0 - 90.0
birubi_thickness = 140.