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