timflow.steady.model.Model.intnormflux#

timflow.steady.model.Model.intnormflux(xy, method='legendre', ndeg=10)#

Integrated normal (perpendicular) flux over polyline.

Computes the flux per segment and per aquifer. Flux to the left is positive when going from (x1, y1) to (x2, y2).

Parameters:

  • xy (list [(x0, y0), (x1, y1),... , (xn, yn)] or 2D array) – if 2D-array, x in first column and y in second column

  • method (str, optional) – integration method, either “quad” (numerical integration using scipy) or “legendre” (approximate integral using Gauss-Legendre quadrature), by default “legendre”.

  • ndeg (int, optional) – degree for legendre polynomial, by default 10, only used when method=”legendre”

Returns:

Qn – integrated normal flux along specified polyline

Return type:

np.array of shape (naq, nsegments)

Examples

Total flow across polyline can be obtained using:

>>> np.sum(Qn)

Total flow across segments summed over aquifers using

>>> np.sum(Qn, axis=0)