Ctrl+K

imod.prepare.c_leakage#

imod.prepare.c_leakage(kh: Union[DataArray, UgridDataArray], kv: Union[DataArray, UgridDataArray], D: Union[DataArray, UgridDataArray], c0: Union[DataArray, UgridDataArray], c1: Union[DataArray, UgridDataArray], B: Union[DataArray, UgridDataArray], length: Union[DataArray, UgridDataArray], dx: Union[DataArray, UgridDataArray], dy: Union[DataArray, UgridDataArray]) → Union[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

Show Source