timflow.transient.model.ModelXsection#

class timflow.transient.model.ModelXsection(naq=1, tmin=1, tmax=10, tstart=0, M=10, steady=None)#

Bases: TimModel

Model class for cross-section models.

Parameters:

  • naq (integer) – number of aquifers

  • tmin (float) – the minimum time for which heads can be computed after any change in boundary condition.

  • tmax (float) – the maximum time for which heads can be computed.

  • tstart (float, optional) – time at start of simulation (default 0)

  • M (integer, optional) – the number of terms to be used in the numerical inversion algorithm. 10 is usually sufficient.

  • steady (timflow.steady.Model) – a timflow.steady model may be included to add a steady-state flow result to the computed solution.

Notes

A ModelXsection may consist of an arbitrary number of Xsection3D or XsectionMaq sections. The combined domain of all sections must span from $x = -infty$ to $x = +infty$, with no gaps.

Methods#

check_inhoms()

Check inhoms.

compute_laplace_parameters()

Compute the parameters for the Laplace transform inversion.

potential(x, y, t[, layers, aq, derivative, returnphi])

Returns pot[naq, ntimes] if layers=None, otherwise pot[len(layers), ntimes].

potentialone(x, y, time[, layers, aq, derivative, ...])

Returns pot[naq] if layers=None, otherwise pot[len(layers)].

disvec(x, y, t[, layers, aq, derivative])

Compute discharge vectgor.

head(x, y, t[, layers, aq, derivative, neglect_steady])

Head at x, y, t where t can be multiple times.

headalongline(x, y, t[, layers])

Head along line or curve.

disvecalongline(x, y, t[, layers])

Discharge vector along line or curve.

headgrid(xg, yg, t[, layers, printrow, show_progress, ...])

Grid of heads.

headgrid2(x1, x2, nx, y1, y2, ny, t[, layers, ...])

Grid of heads.

inverseLapTran(pot, t)

Returns array of potentials of len(t) t must be ordered and tmin<=t<=tmax.

solve([printmat, sendback, silent])

Compute solution.

aquifer_summary()

Return DataFrame with summary of aquifer(s) parameters in model.