polaris.tasks.ocean.mixed_layer_eddies.init.compute_density_anomaly

polaris.tasks.ocean.mixed_layer_eddies.init.compute_density_anomaly(config, y, z, y_period)

Compute the density anomaly of the mixed layer and its fronts, Z_S(z) [Delta rho_y Y_S(y) - Delta rho_z], from Shrestha & Manucharyan (2022), their eq. 5

Parameters:

• config (polaris.config.PolarisConfigParser) – Config options with the mixed_layer_eddies section

• y (xarray.DataArray) – The y coordinate (m)

• z (xarray.DataArray) – Geometric heights (m, negative) relative to the sea surface

• y_period (float) – The periodic length of the domain in y (m), which is the wavelength of the fronts

Returns:

density_anomaly (xarray.DataArray) – The density anomaly (kg m-3)