polaris.tasks.ocean.mixed_layer_eddies.diagnostics.compute_eddy_overturning

polaris.tasks.ocean.mixed_layer_eddies.diagnostics.compute_eddy_overturning(lines, snapshots, bins, config)

Compute the eddy overturning along each line in the thickness-weighted, salinity-binned framework of Shrestha & Manucharyan (2022): the residual overturning in salinity classes minus the Eulerian-mean overturning at the mean depth of each class (Abernathey et al. 2011), both averaged over the snapshots

Parameters:

• lines (list of xarray.Dataset) – The lines from get_frontal_lines()

• snapshots (iterable of xarray.Dataset) – Model output at each time, with normalVelocity, salinity, temperature and a thickness

• bins (numpy.ndarray) – The boundaries of the salinity classes, increasing

• config (polaris.config.PolarisConfigParser) – Config options, for the equation of state

Returns:

ds_eddy (xarray.Dataset) – A dataset with dimensions nLines and nBinsP1 (salinity class boundaries) containing psiResidual, psiEulerian and psiEddy (m2 s-1) and zClass, the mean depth of each class boundary (m)