import sys
from os import sep
sys.path.append('..'+sep+'pyvolution')




alpha = 1
mesh_file = 'merimbula_10785.tsh'
point_file = 'meri0.xya'
mesh_output_file = 'merimbula_10785_%g.tsh'%alpha

print mesh_output_file

from shallow_water import Domain
from pmesh2domain import pmesh_to_domain_instance
from anuga.pyvolution.util import file_function, Polygon_function, read_polygon
from Numeric import zeros, Float, maximum, minimum

#-------
# Domain
filename = mesh_output_file
print 'Creating domain from', filename
domain = pmesh_to_domain_instance(filename, Domain)
print "Number of triangles = ", len(domain)


from realtime_visualisation_new import *

surface =  Surface(domain, scale_z = 1.0/50.0)