# Setup boundary conditions
#------------------------------------------------------------------------------
from math import sin, pi, exp
Br = anuga.Reflective_boundary(domain) # Solid reflective wall
Bt = anuga.Transmissive_boundary(domain) # Continue all values on boundary
Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values
# Associate boundary tags with boundary objects
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:
domain.set_beta(0.8)
# Ensure that the domain definitions make sense
domain.check_integrity()
# Set the inititial conditions
domain.set_quantity('stage', Set_Stage(0.2, 0.4, 1.0))
# Let processor 0 output some timing information
visualise = False
if visualise:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
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:
# Ensure that the domain definitions make sense
domain.check_integrity()
# Set the inititial conditions
domain.set_quantity('stage', Set_Stage(0.2,0.4,1.0))
# Let processor 0 output some timing information
visualise = False
if visualise:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
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()
#domain.use_centroid_velocities = False
print 'after evolve parameters'
import pdb
pdb.set_trace()
if myid == 0:
import time
t0 = time.time()
# Turn on the visualisation
visualise = False
if visualise:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
vis.render_quantity_height("elevation", offset=0.001, dynamic=False)
vis.render_quantity_height("stage", dynamic=True)
vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
vis.start()
import time
time.sleep(2.0)
yieldstep = 0.05
finaltime = 2.0
#Check that the boundary value gets propagated to all elements
for t in domain.evolve(yieldstep=yieldstep, finaltime=finaltime):
if myid == 0:
domain.write_time()
#print_test_stats(domain, tri_full_flag)
print 'after evolve parameters'
import pdb; pdb.set_trace()
if myid == 0:
import time
t0 = time.time()
# Turn on the visualisation
visualise = False
if visualise:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
vis.render_quantity_height("elevation", offset=0.001, dynamic=False)
vis.render_quantity_height("stage", dynamic=True)
vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
vis.start()
import time
time.sleep(2.0)
yieldstep = 0.05
finaltime = 2.0
#Check that the boundary value gets propagated to all elements
for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
if myid == 0:
