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, poraq=0.3, topboundary='conf', phreatictop=False, 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.

  • poraq (array) – Porosities of the aquifers.

  • topboundary (str) – Type of top boundary. Can be ‘conf’ for confined, ‘semi’ for semi-confined or “leaky” for a leaky top boundary.

  • phreatictop (bool) – If true, interpret the first specific storage coefficient as specific yield., i.e. it is not multiplied by aquifer thickness.

  • 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.