Glover (1960) recharge basin

class glover_1960.GloverRectBasinSteady

Bases: object

Glover (1960) solution class.

The default GloverRectBasinSteady object uses the default Aq2dUnconf, aquifer class, RectBasin class and BasinGrid class. Methods include impress at a point (.impress) and grid-contoured impress at specified time (.impress_grid).

Attributes:: aq (obj) : Confined aquifer object. basin (obj) : Basin object. grid (obj) : Basin grid object.

class Aq2dUnconf(K=1, Sy=0.1, B=10, bot=0, name='Aq2dUnconf class')

Bases: Aquifer

2D unconfined aquifer subclass.

Subclass of the Aquifer class defining a horizontal, 2D unconfined aquifer with infinite lateral extent and unconfined storage.

The default Aq2dUnconf object has hydraulic conductivity K=1, specific yield Sy=0.1, static saturated thickness B=10 and aquifer bottom (datum) elevation bot=0. Exceptions occur if invalid values are provided for K, Sy or B.

The .info and .draw methods display the aquifer information and diagram.

Attributes:: Sy (float) : Aquifer specific yield (units 1, default 0.1).

property B

float : Aquifer thickness (units L, default 10.0).

Setter method checks for valid values and triggers an exception if invalid values are specified.

property K

float : Aquifer hydraulic conductivity (units L/T, default 1.0).

Setter method checks for valid values and triggers an exception if invalid values are specified.

property S: float : Aquifer storage coefficient (units 1, default 0.1).

property Sy

float : Aquifer specific yield (units 1, default 0.1).

Setter method checks for valid values and triggers an exception if invalid values are specified.

draw(dw=6)

Display a drawing of the aquifer.

Args:: dw (float) : Width of figure (default 6.0).

info(): Print the aquifer information.

is_1d = False

is_2d = True

is_confined = False

is_finite = False

is_heterogeneous = False

is_homogeneous = True

is_infinite = True

is_leaky = False

is_semifinite = False

is_unconfined = True

property swl: float : Aquifer static water table elevation (units L reduced level, default 10.0).

class RectBasin(cx=0.0, cy=0.0, lx=10, ly=10, rot=0, name='Rectangle basin')

Bases: object

Rectangular recharge basin class.

Attributes:: cx (float) : Basin center x coordinate (units L, default 0.0). cy (float) : Basin center y coordinate (units L, default 0.0). lx (float) : Basin length in x direction (default 10.0) ly (float) : Basin length in y direction (default 10.0) rot (float) : Basin rotation angle in radians (default 0.0). name (str) : Basin name (default ‘Unnamed’).

property area: float : Basin area.

draw(dw=4)

Display a drawing of the basin.

Args:: dw (float) : Width of basin drawing (default 4.0).

info(): Print the basin information.

is_circular = False

is_rectangular = True

property lx

float : Basin length in x direction.

Setter method checks for valid values and triggers an exception if invalid values are specified.

property ly

float : Basin length in y direction.

Setter method checks for valid values and triggers an exception if invalid values are specified.

property rot

float : Basin rotation angle in degrees.

Setter method checks for valid values and triggers an exception if invalid values are specified.

property rot_rad: float : Basin rotation angle in radians.

property verts

dict : x and y coordinates of basin verticies.

Vertex keys: ll - lower left, ul - upper left, lr - lower right and ur - upper right.

property verts_rot

dict : x and y rotated coordinates of basin verticies.

Vertex keys: ll - lower left, ul - upper left, lr - lower right and ur - upper right.

draw(dw=8)

Display the definition diagram.

Args:: dw (float) : Width of figure (default 8.0).

h(x, y, xL, yL, T, S, t, q): Glover impress solution.

impress(t=[1], locs=[(0, 0)], q=0.0, plot=True, csv='', xlsx='')

Calculate impress at specified locations and times.

Args:

t (float) : List of times to evaluate impress (default [1.0]). locs (float tupple) : List of (dx, dy) location tupples to

evaluate impress; dx is distance from basin center in x direction; dy is distance from basin center in y direction (default [(0.0, 0.0)]).

q (float) : Basin infiltration rate (default 0.0). plot (bool) : Display a plot of the results (default True). csv (str) : Filepath for export of results to csv file; results

are exported if the string is not empty (default ‘’).

xlsx (str)Filepath for export of result to xlsx file; results: are exported if the string is not empty (default ‘’).

Returns:

Pandas dataframe containing results, hydraulic loading.

impress_grid(t=1, q=0.0, gr=100, gd=20, plot=True, local=False, csv='', xlsx='')

Calculate impress on a regular grid of points.

Args:

t (float) : Time to evaluate impress (default [1.0]). q (float) : Basin infiltration rate (default 0.0). gr (float) : Radius defining the extent of the solution grid

(default 100.0).

gd (int)Grid density defining the number of gird rows and: columns; minimum and maximum constraints are enforced (default 21).

plot (bool) : Display a plot of the results (default True). local (bool) : Display the drawdown plot in ‘local’ coordinates

with the well at 0, 0 (Default False).

csv (str)Filepath for export of results to csv file; results: are exported if the string is not empty (default ‘’).
xlsx (str)Filepath for export of result to xlsx file; results: are exported if the string is not empty (default ‘’).

Returns:

Pandas dataframe containing results, hydraulic loading.

info()

Print the solution information.

Returns:: Screen printout of solution information.

u1(x, xL, T, S, t, tau): Glover u1 solution term; tau is an integration variable.

u2(x, xL, T, S, t, tau): Glover u2 solution term; tau is an integration variable.

u3(y, yL, T, S, t, tau): Glover u3 solution term; tau is an integration variable.

u4(y, yL, T, S, t, tau): Glover u4 solution term; tau is an integration variable.