source: inundation/ga/storm_surge/analytical solutions/Malpasset_new.py @ 668

"""Validation study of Merimbula lake using Pyvolution.
Example of shallow water wave equation applied to
Malpasset dam break simulation.

   Copyright 2004
   Christopher Zoppou, Stephen Roberts
   Australian National University
   
Specific methods pertaining to the 2D shallow water equation
are imported from shallow_water
for use with the generic finite volume framework

Conserved quantities are h, uh and vh stored as elements 0, 1 and 2 in the
numerical vector named conserved_quantities.

Existence of file 'Malpasset_26000.tsh' is assumed.
"""

###############################
# Setup Path and import modules
import sys
from os import sep, path
sys.path.append('..'+sep+'pyvolution')

from shallow_water import Domain, Reflective_boundary, File_boundary,\
     Dirichlet_boundary, Transmissive_boundary
from pmesh2domain import pmesh_to_domain_instance

######################
# Domain
filename = 'Malpasset_26000_merged.tsh'
yieldstep = 1
finaltime = 10
    
print 'Creating domain from', filename
domain = pmesh_to_domain_instance(filename, Domain)
print "Number of triangles = ", len(domain)

domain.default_order = 1
domain.filename = filename
domain.smooth = True
domain.store = True

domain.set_quantity('friction', 0.07)

#------------------------------
# Boundary Conditions
tags = {}
tags['external'] = Reflective_boundary(domain) 
tags['open'] = Transmissive_boundary(domain)
domain.set_boundary(tags)

#-----------------
#Initial condition
domain.set_quantity('level', 20.)
         
######################
#Evolution
import time
t0 = time.time()
for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
    domain.write_time()
    
print 'That took %.2f seconds' %(time.time()-t0)