mpas_tools.ocean.depth.compute_zmid

mpas_tools.ocean.depth.compute_zmid(bottomDepth, maxLevelCell, layerThickness, depth_dim='nVertLevels')[source]

Computes zMid given data arrays for bottomDepth, maxLevelCell and layerThickness

Parameters:

bottomDepth (xarray.DataArray) – the depth of the ocean bottom (positive)
maxLevelCell (xarray.DataArray) – the 1-based vertical index of the bottom of the ocean
layerThickness (xarray.DataArray) – the thickness of MPAS-Ocean layers (possibly as a function of time)
depth_dim (str, optional) – the name of the vertical dimension

Returns:

zMid (xarray.DataArray) – the vertical coordinate defining the middle of each layer, masked below the bathymetry