polaris.ocean.vertical.compute_zint_zmid_from_layer_thickness

polaris.ocean.vertical.compute_zint_zmid_from_layer_thickness(layer_thickness: DataArray, bottom_depth: DataArray, min_level_cell: DataArray, max_level_cell: DataArray) → tuple[DataArray, DataArray][source]

Compute height z at layer interfaces and midpoints given layer thicknesses and bottom depth.

Parameters:

  • layer_thickness (xarray.DataArray) – The layer thickness of each layer.

  • bottom_depth (xarray.DataArray) – The positive-down depth of the seafloor.

  • min_level_cell (xarray.DataArray) – The zero-based minimum vertical index from each column.

  • max_level_cell (xarray.DataArray) – The zero-based maximum vertical index from each column.

Returns:

  • z_interface (xarray.DataArray) – The elevation of layer interfaces.

  • z_mid (xarray.DataArray) – The elevation of layer midpoints.