flopy.utils.postprocessing Module

flopy.utils.postprocessing.get_gradients(heads, m, nodata, per_idx=None)

Calculates the hydraulic gradients from the heads array for each stress period.

Parameters:

heads : 3 or 4-D np.ndarray

Heads array.

m : flopy.modflow.Modflow object

Must have a flopy.modflow.ModflowDis object attached.

nodata : real

HDRY value indicating dry cells.

per_idx : int or sequence of ints

stress periods to return. If None, returns all stress periods (default).

Returns

——-

grad : 3 or 4-D np.ndarray

Array of hydraulic gradients

flopy.utils.postprocessing.get_saturated_thickness(heads, m, nodata, per_idx=None)

Calculates the saturated thickness for each cell from the heads array for each stress period.

Parameters:

heads : 3 or 4-D np.ndarray

Heads array.

m : flopy.modflow.Modflow object

Must have a flopy.modflow.ModflowDis object attached.

nodata : real

HDRY value indicating dry cells.

per_idx : int or sequence of ints

stress periods to return. If None, returns all stress periods (default).

Returns

——-

sat_thickness : 3 or 4-D np.ndarray

Array of saturated thickness

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

flopy.utils.postprocessing.get_water_table(heads, nodata, per_idx=None)

Get a 2D array representing the water table elevation for each stress period in heads array.

Parameters:

heads : 3 or 4-D np.ndarray

Heads array.

nodata : real

HDRY value indicating dry cells.

per_idx : int or sequence of ints

stress periods to return. If None, returns all stress periods (default).

Returns

——-

wt : 2 or 3-D np.ndarray of water table elevations

for each stress period.