polaris.ocean.vertical.ztilde.pressure_from_geom_thickness

polaris.ocean.vertical.ztilde.pressure_from_geom_thickness(surf_pressure: DataArray, geom_layer_thickness: DataArray, spec_vol: DataArray) → tuple[DataArray, DataArray]

Compute the pressure at layer interfaces and midpoints given surface pressure, geometric layer thicknesses, and specific volume. This calculation assumes a constant specific volume within each layer.

Parameters:

• surf_pressure (xarray.DataArray) – The surface pressure at the top of the water column.

• geom_layer_thickness (xarray.DataArray) – The geometric thickness of each layer, set to zero for invalid layers.

• spec_vol (xarray.DataArray) – The specific volume at each layer.

Returns:

• p_interface (xarray.DataArray) – The pressure at layer interfaces.

• p_mid (xarray.DataArray) – The pressure at layer midpoints.