timflow.transient.linesink1d.FluxDiffLineSink1D

class timflow.transient.linesink1d.FluxDiffLineSink1D(model, xls=0, layers=0, label=None, aq=None)

Bases: LineSink1DBase, timflow.transient.equation.FluxDiffEquation

1D flux-difference linesink element.

Used to ensure continuity of flux in a cross-section model, e.g. at the boundary of an inhomogeneity.

Parameters:

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

• xls (float) – x-coordinate of the linesink

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

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

Methods

initialize()	Initialize the element.
plot([ax])	Plot 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.
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.
equation()	Mix-in class that returns matrix rows for continuity of flow.