source: anuga_work/production/wollongong_2006/project_kembla.py @ 4853

# -*- 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_points, polygon_area, is_inside_polygon

#Making assumptions about the location of scenario data

dem_name = 'kembla'
basename = 'source'
meshname = 'mesh_' + basename
combined_dem_name = 'kembla_elevation'

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

# bounding polygon for study area
polyAll = read_polygon('kembla_domain.csv')

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

# areaA digitized polygons
poly_bay = read_polygon('bay.csv')

print 'Area of bay polygon', polygon_area(poly_bay)/1000000.0

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

# exporting asc grid
eastingmin = 305500
eastingmax = 309500
northingmin = 6181500
northingmax = 6185500