def __init__(self,
coordinates,
vertices,
boundary = None,
tagged_elements = None,
geo_reference = None,
use_inscribed_circle=False,
velocity = None,
full_send_dict=None,
ghost_recv_dict=None,
processor=0,
numproc=1
):
conserved_quantities = ['stage']
other_quantities = []
Generic_Domain.__init__(self,
source=coordinates,
triangles=vertices,
boundary=boundary,
conserved_quantities=conserved_quantities,
other_quantities=other_quantities,
tagged_elements=tagged_elements,
geo_reference=geo_reference,
use_inscribed_circle=use_inscribed_circle,
full_send_dict=full_send_dict,
ghost_recv_dict=ghost_recv_dict,
processor=processor,
numproc=numproc)
if velocity is None:
self.velocity = num.array([1,0],'d')
else:
self.velocity = num.array(velocity,'d')
#Only first is implemented for advection
self.set_default_order(1)
self.set_beta(1.0)
self.smooth = True
self.max_flux_update_frequency=1
def __init__(self,
coordinates,
vertices,
boundary=None,
tagged_elements=None,
geo_reference=None,
use_inscribed_circle=False,
velocity=None,
full_send_dict=None,
ghost_recv_dict=None,
processor=0,
numproc=1):
conserved_quantities = ['stage']
other_quantities = []
Generic_Domain.__init__(self,
source=coordinates,
triangles=vertices,
boundary=boundary,
conserved_quantities=conserved_quantities,
other_quantities=other_quantities,
tagged_elements=tagged_elements,
geo_reference=geo_reference,
use_inscribed_circle=use_inscribed_circle,
full_send_dict=full_send_dict,
ghost_recv_dict=ghost_recv_dict,
processor=processor,
numproc=numproc)
if velocity is None:
self.velocity = num.array([1, 0], 'd')
else:
self.velocity = num.array(velocity, 'd')
#Only first is implemented for advection
self.set_default_order(1)
self.set_beta(1.0)
self.smooth = True
self.max_flux_update_frequency = 1
