timflow.transient.circareasink.CircAreaSink#

class timflow.transient.circareasink.CircAreaSink(model, xc=0, yc=0, R=0.1, tsandN=[(0, 1)], name='CircAreaSink', label=None)#

Bases: timflow.transient.element.Element

Create a circular area-sink with uniform infiltration rate in aquifer layer 0.

Infiltration rate in length / time, positive for water entering the aquifer.

Parameters:

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

  • xc (float) – x-coordinate of center of area-sink

  • yc (float) – y-coordinate of center of area-sink

  • R (radius of area-sink)

  • tsandN (list of tuples) – tuples of starting time and infiltration rate after starting time

  • label (string or None (default: None)) – label of the area-sink

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,y value.

disvecinf(x, y[, aq])

Can be called with only one x,y value.

plot([ax, layer])

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.