flopy.utils.postprocessing Module

flopy.utils.postprocessing.get_transmissivities(heads, m, r=None, c=None, x=None, y=None, sctop=None, scbot=None, nodata=-999)

Computes transmissivity in each model layer at specified locations and open intervals. A saturated thickness is determined for each row, column or x, y location supplied, based on the open interval (sctop, scbot), if supplied, otherwise the layer tops and bottoms and the water table are used.

Parameters:

heads : 2D array

numpy array of shape nlay by n locations

m : flopy.modflow.Modflow object

Must have dis, sr, and lpf or upw packages.

r : 1D array-like of ints, of length n locations

row indices (optional; alternately specify x, y)

c : 1D array-like of ints, of length n locations

column indicies (optional; alternately specify x, y)

x : 1D array-like of floats, of length n locations

x locations in real world coordinates (optional)

y : 1D array-like of floats, of length n locations

y locations in real world coordinates (optional)

sctop : 1D array-like of floats, of length n locations

open interval tops (optional; default is model top)

scbot : 1D array-like of floats, of length n locations

open interval bottoms (optional; default is model bottom)

nodata : numeric

optional; locations where heads=nodata will be assigned T=0

Returns:

T : 2D array of same shape as heads (nlay x n locations)

Transmissivities in each layer at each location