timflow.transient.inhom1d.Xsection3D#

class timflow.transient.inhom1d.Xsection3D(model, x1, x2, kaq=1, z=(4, 3, 2, 1), Saq=0.001, kzoverkh=0.1, leffaq=0, poraq=0.3, topboundary='conf', phreatictop=None, topres=0, topthick=0, topSll=0, toppor=0.3, tsandhstar=None, tsandN=None, name=None)#

Bases: Xsection

Cross-section inhomogeneity consisting of stacked aquifer layers.

Vertical resistance is computed from vertical hydraulic conductivity and the anisotropy factor.

Parameters:

  • model (Model) – Model to add the cross-section to, usually an instance of ModelXsection.

  • x1 (scalar) – x-coordinate of the left boundary of the cross-section.

  • x2 (scalar) – x-coordinate of the right boundary of the cross-section.

  • kaq (array) – Hydraulic conductivities of the aquifers.

  • z (array) – Elevations of the tops and bottoms of the layers.

  • Saq (array) – Specific storage of the aquifers.

  • kzoverkh (scalar) – Ratio of vertical hydraulic conductivity to horizontal hydraulic conductivity.

  • leffaq (array) – Loading efficiency

  • poraq (array) – Porosities of the aquifers.

  • topboundary (string, 'confined', 'phreatic', or 'semi' (default is 'conf')) – indicating whether the top is confined (‘con’ is enough), phreatic (‘phr’ is enough) or semi-confined (‘sem’ is enough). When ‘phreatic’, the storage coefficient (Saq) of the top model layer is treated as phreatic storage (and not multiplied with the aquifer thickness) When ‘semi’, the topres and topthick must be specified.

  • topres (scalar) – Resistance of the top boundary. Only used if topboundary is ‘leaky’.

  • topthick (scalar) – Thickness of the top boundary. Only used if topboundary is ‘leaky’.

  • topSll (scalar) – Specific storage of the top boundary. Only used if topboundary is ‘leaky’.

  • toppor (scalar) – Porosity of the top boundary. Only used if topboundary is ‘leaky’.

  • tsandhstar (list of tuples) – list containing time and water level pairs for the hstar boundary condition.

  • tsandN (list of tuples) – list containing time and infiltration pairs for the infiltration boundary condition.

  • name (str) – Name of the cross-section.

Methods#

is_inside(x, _)

Check if a point is inside the cross-section.

initialize()

Initialize the aquifer data.

create_elements()

Create linesinks to meet the continuity conditions the at the boundaries.

plot([ax, labels, params, names, fmt, sep])

Plot the cross-section.

in_which_layer(z)

Get layer given elevation z.

findlayer(z)

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