"""Script for fitting data points to mesh
"""

from pyvolution.pmesh2domain import pmesh_to_domain_instance 
from caching import cache
from pyvolution.shallow_water import Domain, Reflective_boundary


meshname = 'swan_mesh.msh'

domain = cache(pmesh_to_domain_instance, (meshname, Domain),
               dependencies = [meshname],                     
               verbose = True)               


domain.set_quantity('elevation',
                    filename = 'cross_sections.pts',
                    use_cache = True,
                    verbose = True)

print domain.get_boundary_tags()
Br = Reflective_boundary(domain)
domain.set_boundary( {'exterior': Br} )


domain.set_name('swan_river')
domain.set_datadir('.')
domain.store = True

for t in domain.evolve(yieldstep = 1, finaltime = 1):
    domain.write_time()