###########################################
# Create advection domain with direction (1,-1)
# Initial condition is zero by default
domain = Parallel_Domain(points,
vertices,
boundary,
full_send_dict,
ghost_recv_dict,
velocity=[1.0, 0.0])
# Boundaries
T = Transmissive_boundary(domain)
domain.default_order = 2
domain.set_boundary( {'left': T, 'right': T, 'bottom': T, 'top': T, \
'ghost': None} )
# 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
if myid == 0:
import time
"no vertices = ", len(vertices), "no boundaries = ", len(boundary)
###########################################
# Start the computations on each subpartion
###########################################
# Create advection domain with direction (1,-1)
# Initial condition is zero by default
domain = Parallel_Domain(points, vertices, boundary,
full_send_dict, ghost_recv_dict, velocity=[1.0, 0.0])
# Boundaries
T = Transmissive_boundary(domain)
domain.default_order = 2
domain.set_boundary( {'left': T, 'right': T, 'bottom': T, 'top': T, \
'ghost': None} )
# 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
if myid == 0:
import time
