mpas_tools.ocean.depth.compute_depth

mpas_tools.ocean.depth.compute_depth(refBottomDepth)[source]

Computes depth and depth bounds given refBottomDepth

Note

This produces a 1D depth coordinate from refBottomDepth, a single reference depth profile that is independent of horizontal location. This assumes a z-level/z-star vertical coordinate and is not supported for Omega or other general vertical coordinates, which do not provide refBottomDepth. Use add_zmid() for a 3D zMid coordinate instead.

Parameters:

refBottomDepth (xarray.DataArray) – the depth of the bottom of each vertical layer in the initial state (perfect z-level coordinate)

Returns:

  • depth (numpy.ndarray) – the vertical coordinate defining the middle of each layer

  • depth_bnds (numpy.ndarray) – the vertical coordinate defining the top and bottom of each layer