polaris.tasks.ocean.mixed_layer_eddies.diagnostics.compute_mu

polaris.tasks.ocean.mixed_layer_eddies.diagnostics.compute_mu(z, ml_depth)[source]

Compute the vertical structure function of Fox-Kemper et al. (2008), zero below the mixed layer

Parameters:

z (numpy.ndarray or xarray.DataArray) – Heights (m, negative) relative to the sea surface
ml_depth (float or numpy.ndarray or xarray.DataArray) – The mixed-layer depth H (m, positive)

Returns:

mu (numpy.ndarray or xarray.DataArray) – The structure function, with a maximum of 1 at z = -H/2