timflow.steady.model.Model3D#

class timflow.steady.model.Model3D(kaq=1, z=None, kzoverkh=1, npor=0.3, topboundary='conf', topres=0, topthick=0, hstar=0)#

Bases: Model

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

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

Parameters:

  • 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’)

Notes

For semi-confined aquifers, set top equal to ‘semi’ and provide:

  • topres: resistance of top

  • tophick: thickness of top

  • hstar: head above top

Examples

Create a 3D model with 10 layers:

ml = Model3D(kaq=10, z=np.arange(20, -1, -2), kzoverkh=0.1)

Methods#

remove_element(e)

Remove element e from model.

disvec(x, y[, aq])

Discharge vector at x, y.

normflux(x, y, theta)

Flux at point x, y in direction of angle theta.

intnormflux_segment(x1, y1, x2, y2[, method, ndeg])

Integrated normal (perpendicular) flux over specified line segment.

intnormflux(xy[, method, ndeg])

Integrated normal (perpendicular) flux over polyline.

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

Head at x, y.

headgrid(xg, yg[, layers, printrow])

Grid of heads.

headgrid2(x1, x2, nx, y1, y2, ny[, layers, printrow])

Grid of heads.

headalongline(x, y[, layers])

Head along line or curve.

disvecalongline(x, y[, layers])

Compute discharge vector along line.

solve([printmat, sendback, silent])

Compute solution.

solve_mp([nproc, printmat, sendback, silent])

Compute solution, multiprocessing implementation.

aquifer_summary()

Return DataFrame with summary of aquifer(s) parameters in model.