source: development/momentum_sink/compare_sww.py @ 2378

"""Read in sww file, interpolate at specified locations and  times 
"""

from os import sep
import project
from pyvolution.util import file_function


#swwfile = project.outputname + '.sww'
swwfile = project.outputdir + sep  + 'friction_12330.sww'
gauge_locations = [[10,10],[20,20]]

quantities = ['stage']
gauges = file_function(swwfile,
                  quantities = quantities,
                  interpolation_points = gauge_locations,
                  verbose = True,
                  use_cache = True)

times = [0,1,2,3,4]

for gauge, _ in enumerate(gauge_locations):
    print "gauge", gauge
    for time in times:
        print "results ", gauges(time, gauge)