ttim.inhom1d#

Classes#

Xsection

Base class for a cross-section inhomogeneity.

XsectionMaq

Cross-section inhomogeneity consisting of stacked aquifer layers.

Xsection3D

Cross-section inhomogeneity consisting of stacked aquifer layers.

Module Contents#

class ttim.inhom1d.Xsection(model, x1, x2, kaq, z, Haq, Hll, c, Saq, Sll, poraq, porll, ltype, topboundary, phreatictop, tsandhstar, tsandN, kzoverkh=None, model3d=False, name=None)[source]#

Bases: ttim.aquifer.AquiferData

Base class for a cross-section inhomogeneity.

Parameters:

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

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

  • x2 (float) – 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.

  • Haq (array) – Thicknesses of the aquifers.

  • Hll (array) – Thicknesses of the leaky layers.

  • c (array) – Resistance of the leaky layers.

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

  • Sll (array) – Specific storage of the leaky layers.

  • poraq (array) – Porosities of the aquifers.

  • porll (array) – Porosities of the leaky layers.

  • ltype (array) – Type of each layer. ‘a’ for aquifer, ‘l’ for leaky layer.

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

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

  • kzoverkh (float, optional,) – anisotropy factor for vertical resistance, kzoverkh = kz / kh. Default is 1.

  • model3d (bool, optional) – Boolean indicating whether model is Model3D-type.

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

tiny = 1e-12#
x1#
x2#
tsandhstar#
tsandN#
inhom_number#
addlinesinks = True#
__repr__()[source]#
is_inside(x, _)[source]#

Check if a point is inside the cross-section.

Parameters:

x (float) – x-coordinate of the point.

Returns:

True if the point is inside the cross-section, False otherwise.

Return type:

bool

initialize()[source]#

Initialize the aquifer data.

eigval[naq, npval]#

Array with eigenvalues

lab[naq, npval]#

Array with lambda values

lab2[naq, nint, npint]#

Array with lambda values reorganized per interval

eigvec[naq, naq, npval]#

Array with eigenvector matrices

coef[naq, naq, npval]#

Array with coefficients coef[ilayers, :, np] are the coefficients if the element is in ilayers belonging to Laplace parameter number np.

create_elements()[source]#

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

plot(ax=None, labels=False, params=False, names=False, fmt=None, **kwargs)[source]#

Plot the cross-section.

Parameters:

  • ax (plt.Axes, optional) – Axis to plot the cross-section on. If None, a new axis will be created.

  • labels (bool, optional) – If True, add layer-name labels.

  • params (bool, optional) – If True, add parameter labels.

  • names (bool, optional) – If True, add inhomogeneity names.

  • fmt (str, optional) – format string for parameter values, e.g. ‘.2f’ for 2 decimals.

class ttim.inhom1d.XsectionMaq(model, x1, x2, kaq=1, z=(1, 0), c=(), Saq=0.001, Sll=0, poraq=0.3, porll=0.3, topboundary='conf', phreatictop=False, tsandhstar=None, tsandN=None, name=None)[source]#

Bases: Xsection

Cross-section inhomogeneity consisting of stacked aquifer layers.

Parameters:

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

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

  • x2 (float) – 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.

  • c (array) – Resistance of the leaky layers.

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

  • Sll (array) – Specific storage of the leaky layers.

  • poraq (array) – Porosities of the aquifers.

  • porll (array) – Porosities of the leaky layers.

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

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

class ttim.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)[source]#

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.