source: anuga_work/production/test_dampier_2006/test.py @ 4347

from os import sep, environ, getenv, getcwd
import os

import project
from anuga.pmesh.mesh_interface import create_mesh_from_regions
from anuga.shallow_water import Domain
from anuga.shallow_water import File_boundary

home = getenv('INUNDATIONHOME') #Sandpit's parent dir   
#user = get_user_name()
data = 'data'
state = 'western_australia'
scenario_name = 'dampier'
scenario = 'dampier_tsunami_scenario_2006'
an = 'anuga'
bo = 'boundaries'
urs = 'urs'
gridded = '1_10000'
best_boundary_sww_file = 'o_test_8500_12000.sww'

mesh_filename = 'small.tsh' # this will be local
best_boundary_sww = os.path.join(data,state,scenario_name,
                                 scenario,an,bo,best_boundary_sww_file)

#pro_instance = project.Project([data,state,scenario_name],
#                                       outputdir_name=outputdir_name)

poly_all = [[250000,7660000],[280000,7660000],
                             [280000,7630000]]
create_mesh_from_regions(poly_all,
                             boundary_tags={'back': [0,1,2]},
                             maximum_triangle_area=100000000000000000000,
                             filename=mesh_filename,
                             use_cache=False,
                             verbose=True)
domain = Domain( mesh_filename, verbose=True)
print domain.statistics()

domain.set_quantity('stage', 0.0)
domain.set_quantity('friction', 0.01)
domain.set_quantity('elevation', 
#                    filename = project.combined_dir_name + '.pts',
# MUST USE TXT FILES FOR CACHING TO WORK!
                    filename = project.combined_smaller_dir_name + '.txt',
                    use_cache = True,
                    verbose = True,
                    alpha = 0.1)
    
Bf = File_boundary(best_boundary_sww,
                  domain,  use_cache=False, verbose=True)