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

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

    Parameters differ from example to make the code run faster, and so the data is stored and found in the correct places
"""

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.