def initialize(self):
LineSinkHoBase.initialize(self)
self.xcin, self.ycin = controlpoints(self.Ncp-1, self.z1, self.z2, eps=1e-6, include_ends=True)
self.xcout, self.ycout = controlpoints(self.Ncp-1, self.z1, self.z2, eps=-1e-6, include_ends=True)
if self.aqin is None:
self.aqin = self.model.aq.find_aquifer_data(self.xcin[0], self.ycin[0])
if self.aqout is None:
self.aqout = self.model.aq.find_aquifer_data(self.xcout[0], self.ycout[0])
def __init__(self, model, x1=-1, y1=0, x2=1, y2=0, \
order=0, ndeg=3, label=None, addtomodel=True, aq=None, aqin=None, aqout=None):
LineSinkHoBase.__init__(self, model, x1, y1, x2, y2, Qls=0, \
layers=range(model.aq.Naq), order= order, name='IntFluxLineSink', label=label, \
addtomodel=addtomodel, aq=aq)
self.inhomelement = True
self.ndeg = ndeg
self.Xleg, self.wleg = np.polynomial.legendre.leggauss(self.ndeg)
self.Nunknowns = self.Nparam
self.aqin = aqin
self.aqout = aqout
