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()
vis.render_quantity_height("stage", zScale=5.0, dynamic=True)
vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
vis.start()
import time
time.sleep(2.0)
if myid == 0:
import time
t0 = time.time()
for t in domain.evolve(yieldstep=5.0, finaltime=30.0):
if myid == 0:
domain.write_time()
if visualise:
vis.update()
if visualise: vis.evolveFinished()
# Output some computation statistics
if myid == 0:
print 'That took %.2f seconds' % (time.time() - t0)
print 'Communication time %.2f seconds'\
%domain.communication_time
print 'Reduction Communication time %.2f seconds'\
%domain.communication_reduce_time
if visualise:
vis.join()
domain.set_boundary({'left': Bt, 'right': Bt, 'top': Br, 'bottom': Br})
#===============================================================================
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()
if myid == 0:
import time
t0 = time.time()
for t in domain.evolve(yieldstep = 5.0, finaltime = 30.0):
if myid == 0:
domain.write_time()
if visualise:
vis.update()
if visualise: vis.evolveFinished()
# Output some computation statistics
if myid == 0:
print 'That took %.2f seconds' %(time.time()-t0)
print 'Communication time %.2f seconds'\
%domain.communication_time
print 'Reduction Communication time %.2f seconds'\
%domain.communication_reduce_time
if visualise: vis.join()
