timflow.steady.inhomogeneity.PolygonInhom3D#

class timflow.steady.inhomogeneity.PolygonInhom3D(model, xy, kaq=1, z=None, kzoverkh=1, npor=0.3, topboundary='conf', topres=0, topthick=0, hstar=0, N=None, order=3, ndeg=3)#

Bases: PolygonInhom

Create a multi-layer model object consisting of many aquifer layers.

The resistance between the layers is computed from the vertical hydraulic conductivity of the layers.

Parameters:

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

  • xy (array or list) – list or array of (x,y) pairs of coordinates of corners of the inhomogeneity. polygonal boundary is automatically closed (so first point is not repeated)

  • kaq (float, array or list) – hydraulic conductivity of each layer from the top down if float, hydraulic conductivity is the same in all aquifers

  • z (array or list) – elevation of top of system followed by bottoms of all layers from the top down bottom of layer is automatically equal to top of layer below it length is number of aquifer layers + 1

  • kzoverkh (float) – vertical anisotropy ratio vertical k divided by horizontal k if float, value is the same for all layers length is number of layers

  • npor (float, array or list) – porosity of all aquifer layers from the top down if float, porosity is the same for all layers if topboundary=’conf’: length is number of layers if topboundary=’semi’: length is number of layers + 1

  • topboundary (string, 'conf' or 'semi' (default is 'conf')) – indicating whether the top is confined (‘conf’) or semi-confined (‘semi’)

  • topres (float) – resistance of top semi-confining layer (read if topboundary=’semi’)

  • topthick (float) – thickness of top semi-confining layer (read if topboundary=’semi’)

  • hstar (float or None (default is None)) – head value above semi-confining top (read if topboundary=’semi’)

  • N (float or None (default is None)) – infiltration rate (L/T) inside inhomogeneity. Only possible if topboundary=’conf’

  • order (int) – polynomial order of flux along each segment

  • ndeg (int) – number of points used between two segments to numerically integrate normal discharge

Methods#

findlayer(z)

Returns layer-number, layer-type and model-layer-number.

summary()

Get summary of aquifer parameters.