source: trunk/anuga_work/development/mem_time_tests/hardware/cairns/project.py @ 8311

""" Common filenames and locations for topographic data, meshes and outputs.
    This file defines the parameters of the scenario you wish to run.
"""

import anuga
from anuga.abstract_2d_finite_volumes.util import add_directories
import os


home2 = os.getenv('INUNDATIONHOME')

scenariodir2 = add_directories(home2, ["data", "mem_time_test", "parallel", "cairns"])

h = 'CAIRNS.msh'
file_pathh = os.path.join(scenariodir2, h)
store ='store.txt'
file_path_store = os.path.join(scenariodir2, store)
#------------------------------------------------------------------------------
# Define scenario as either slide or fixed_wave. Choose one.
#------------------------------------------------------------------------------
scenario = 'fixed_wave' # Huge wave applied at the boundary
#scenario = 'slide'       # Slide wave form applied inside the domain

#------------------------------------------------------------------------------
# Filenames
#------------------------------------------------------------------------------
meshname = os.path.join(scenariodir2, 'cairns.msh')

# Filename for locations where timeseries are to be produced
gauge_filename = os.path.join(scenariodir2, 'gauges.csv')

#------------------------------------------------------------------------------
# Domain definitions
#------------------------------------------------------------------------------
# bounding polygon for study area
bounding_polygon = anuga.read_polygon(os.path.join(scenariodir2, 'extent.csv'))

A = anuga.polygon_area(bounding_polygon) / 1000000.0
print 'Area of bounding polygon = %.2f km^2' % A

#------------------------------------------------------------------------------
# Interior region definitions
#------------------------------------------------------------------------------
# Read interior polygons
poly_cairns = anuga.read_polygon(os.path.join(scenariodir2, 'cairns.csv'))
poly_island0 = anuga.read_polygon(os.path.join(scenariodir2, 'islands.csv'))
poly_island1 = anuga.read_polygon(os.path.join(scenariodir2, 'islands1.csv'))
poly_island2 = anuga.read_polygon(os.path.join(scenariodir2, 'islands2.csv'))
poly_island3 = anuga.read_polygon(os.path.join(scenariodir2, 'islands3.csv'))
poly_shallow = anuga.read_polygon(os.path.join(scenariodir2, 'shallow.csv'))

# Optionally plot points making up these polygons
#plot_polygons([bounding_polygon, poly_cairns, poly_island0, poly_island1,
#               poly_island2, poly_island3, poly_shallow],
#               style='boundingpoly', verbose=False)

# Define resolutions (max area per triangle) for each polygon
# Make these numbers larger to reduce the number of triangles in the model,
# and hence speed up the simulation
default_res = 10000000    # Background resolution
islands_res = 100000
cairns_res = 100000
shallow_res = 500000

# Define list of interior regions with associated resolutions
interior_regions = [[poly_cairns,  cairns_res],
                    [poly_island0, islands_res],
                    [poly_island1, islands_res],
                    [poly_island2, islands_res],
                    [poly_island3, islands_res],
                    [poly_shallow, shallow_res]]

#------------------------------------------------------------------------------
# Data for exporting ascii grid
#------------------------------------------------------------------------------
eastingmin = 363000
eastingmax = 418000
northingmin = 8026600
northingmax = 8145700

#------------------------------------------------------------------------------
# Data for landslide
#------------------------------------------------------------------------------
slide_origin = [451871, 8128376]   # Assume to be on continental shelf
slide_depth = 500.