Calculate saturated thicknesses

rmf_saturated_thickness(hed, dis, l = NULL, na_values = NULL)

Arguments

hed

3d or 4d array object containing hydraulic head values for each cell

dis

RMODFLOW dis object

l

integer used to subset the 4th dimension of hed. If not supplied, the final time step is used.

na_values

optional; specifies which hed values should be set to NA.

Value

a rmf_3d_array with saturated thicknesses for each cell.