timflow.steady.model.Model

class timflow.steady.model.Model(kaq, c, z, npor, ltype)

Create a model consisting of an arbitrary sequence of aquifers and leaky layers.

Notes

Use ModelMaq for regular sequence of aquifers and leaky layers. Use Model3D for multi-layer model of a single aquifer.

Parameters:

• kaq (array) – hydraulic conductivity of each aquifer from the top down

• z (array) – elevation tops and bottoms of all layers layers may have zero thickness

• c (array) – resistance between two consecutive aquifer layers if ltype[0]=’a’: length is number of aquifers - 1 if ltype[0]=’l’: length is number of aquifers

• npor (array) – porosity of all layers from the top down

• ltype (array of characters) – array indicating for each layer whether it is ‘a’ aquifer layer ‘l’ leaky layer

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.