polaris.tasks.ocean.seamount.init_utils.compute_target_density

polaris.tasks.ocean.seamount.init_utils.compute_target_density(config: PolarisConfigParser, z_mid: DataArray) → DataArray

Compute the target Beckmann and Haidvogel stratification, either the linear profile of their eqn 15 or the exponential profile of their eqn 16.

The profile is a potential density referenced to the surface. Under the linear equation of state, which has no pressure dependence, that is the same thing as the density; under a nonlinear equation of state it is not, and reading it as in-situ density would be unphysical – TEOS-10 in-situ density at 5000 m is near 1050 kg m-3 from compression alone, well outside the range this profile spans.

Parameters:

• config (polaris.config.PolarisConfigParser) – Configuration with the [seamount] stratification options.

• z_mid (xarray.DataArray) – The geometric height (m, negative down) of the layer midpoints.

Returns:

xarray.DataArray – The target potential density (kg m-3) at the layer midpoints.