imod.prepare.c_leakage#

imod.prepare.c_leakage(kh: DataArray | UgridDataArray, kv: DataArray | UgridDataArray, D: DataArray | UgridDataArray, c0: DataArray | UgridDataArray, c1: DataArray | UgridDataArray, B: DataArray | UgridDataArray, length: DataArray | UgridDataArray, dx: DataArray | UgridDataArray, dy: DataArray | UgridDataArray) → DataArray | UgridDataArray[source]#

Computes the phreatic leakage resistance.

Parameters:

kh (xr.DataArray of floats) – horizontal conductivity of phreatic aquifer
kv (xr.DataArray of floats) – vertical conductivity of phreatic aquifer
c0 (xr.DataArray of floats) – hydraulic bed resistance of water feature
c1 (xr.DataArray of floats) – hydraulic resistance of the first aquitard
D (xr.DataArray of floats) – saturated thickness of the top system
B (xr.DataArray of floats) – water feature wetted perimeter
length (xr.DataArray of floats) – water feature length per cell
dx (xr.DataArray of floats) – cellsize in x
dy (xr.DataArray of floats) – cellsize in y

Returns:

c_leakage – Hydraulic resistance of water features corrected for intra-cell hydraulic resistance and surface water interaction.

Return type:

xr.DataArray of floats