timflow.transient.linesink1d.LineSink1DBase#

class timflow.transient.linesink1d.LineSink1DBase(model, xls=0, tsandbc=[(0, 1)], res=0, wh='H', layers=0, type='', name='LineSink1DBase', label=None, aq=None, inhomelement=False)#

Bases: timflow.transient.element.Element

LineSink1D Base Class.

All LineSink1D elements are derived from this class

Parameters:

  • model (Model object) – Model to which the element is added

  • xls (float) – x-coordinate of the line sink

  • tsandbc (list of tuples) – list of tuples of the form (time, bc) for boundary conditions

  • res (float) – resistance of the line sink

  • wh (string or float) – wetted perimeter of the linesink, “H” for aquifer height, “2H” for 2x aquifer height (two-sided flow) or specify any float value

  • layers (int, array or list) – layer (int) or layers (list or array) in which line sink is located

  • type (string) – type of element, “g” for given, “v” for variable and “z” for zero.

  • name (string) – name of the element

  • label (string, optional) – label of the element

  • aq (Aquifer object) – aquifer in which the element is located

  • inhomelement (boolean) – set to True if element is part of an inhomogeneity

Methods#

initialize()

Initialize the element.

setflowcoef()

Separate function so that this can be overloaded for other types.

potinf(x[, y, aq])

Can be called with only one x value.

disvecinf(x[, y, aq])

Can be called with only one x,y value.

plot([ax])

Plot the element.

potential(x, y[, aq])

Returns complex array of size (ngvbc, naq, npval).

unitpotential(x, y[, aq])

Returns complex array of size (naq, npval).

unitpotentialone(x, y, jtime[, aq])

Returns complex array of size (naq, npval).

disvec(x, y[, aq])

Returns 2 complex arrays of size (ngvbc, naq, npval).

unitdisvec(x, y[, aq])

Returns 2 complex arrays of size (naq, npval).

potinflayers(x, y[, layers, aq])

Layers can be scalar, list, or array.

potentiallayers(x, y[, layers, aq])

Returns complex array of size (ngvbc, len(layers),npval).

unitpotentiallayers(x, y[, layers, aq])

Returns complex array of size (len(layers), npval).

disvecinflayers(x, y[, layers, aq])

Layers can be scalar, list, or array.

disveclayers(x, y[, layers, aq])

Returns 2 complex array of size (ngvbc, len(layers), npval).

unitdisveclayers(x, y[, layers, aq])

Returns complex array of size (len(layers), npval).

discharge(t[, derivative])

The discharge in each layer.

dischargeold(t[, derivative])

The discharge in each layer.

run_after_solve()

Function to run after a solution is completed.