vtk_visualiser = True
if vtk_visualiser:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
vis.render_quantity_height("height",dynamic=True)
#vis.render_quantity_height("stage", zScale =1.0, dynamic=True)
vis.colour_height_quantity('stage', (0.0, 0.0, 1.0))
vis.start()
#===============================================================================
#------------------------------------------------------------------------------
# Evolve system through time
#------------------------------------------------------------------------------
i = 0
for t in domain.evolve(yieldstep = 0.01, finaltime = 0.25):
#print domain.timestepping_statistics(track_speeds=True)
print domain.timestepping_statistics()
if vtk_visualiser:
vis.update()
fileName = 'stage_%03d' % i + '.vtk'
i = i+1
vis.store_height_quantity('height',fileName)
#test against know data
if vtk_visualiser: vis.evolveFinished()
