imod.select.get_lower_active_layer_number

imod.select.get_lower_active_layer_number(active: DataArray | UgridDataArray) → DataArray | UgridDataArray

Returns planar grid of integers with the layer number of the lower most active cell.

Parameters:

active ({xr.DataArray, xu.UgridDataArray}) – Grid of booleans (…, layer, y, x) designating active cell.

Returns:

lower_active_layer – Planar grid of integers (…, y, x) with the layer number of the lowermost active cell.

Return type:

{xr.DataArray, xu.UgridDataArray}

Examples

To get the layer numbers of the lowermost active cells for an idomain grid:

>>> active = idomain > 0
>>> lower_active_layer = get_lower_active_layer_number(active)

To get the layer numbers of the lowermost active cells for a data grid with floats, where np.nan indicates inactive cells:

>>> active = data.notnull()
>>> lower_active_layer = get_lower_active_layer_number(active)