source: anuga_core/documentation/user_manual/demos/cairns/project.py @ 3979

# -*- 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'

#Making assumptions about the location of scenario data
state = 'queensland'
scenario_dir_name = 'cairns_tsunami_scenario_demo'

# onshore data provided by WA DLI
demname = 'cairns'

#swollen/ all data output
basename = 'source'
codename = 'project.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
outputtimedir = outputdir + local_time + sep
polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep

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

# Necessary if using point datasets, rather than grid
dem_name = datadir + demname

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

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

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

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

# exporting asc grid
eastingmin = 356300
eastingmax = 398500
northingmin = 8107850
northingmax = 8148700

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

# broome digitized polygons
poly_cairns = read_polygon(polygondir+'cairns.csv')
poly_island0 = read_polygon(polygondir+'islands.csv')
poly_island1 = read_polygon(polygondir+'islands1.csv')
poly_island2 = read_polygon(polygondir+'islands2.csv')
poly_island3 = read_polygon(polygondir+'islands3.csv')
poly_shallow = read_polygon(polygondir+'shallow.csv')

plot_polygons([polyAll,poly_cairns,poly_island0,poly_island1,\
               poly_island2,poly_island3,poly_shallow],\
              'boundingpoly',verbose=False)

#slump_origin = [633361, 8129700]
#slump_depth = 986
slump_origin = [451871, 8128376] # move onto the continental shelf, depth = 500
slump_depth = 500.