polaris.tasks.ocean.mixed_layer_eddies.diagnostics.fit_mu_amplitude

polaris.tasks.ocean.mixed_layer_eddies.diagnostics.fit_mu_amplitude(psi, z, ml_depth)

Fit psi(z) = amplitude * mu(z) by least squares over the mixed layer

Parameters:

• psi (numpy.ndarray) – Stream function values

• z (numpy.ndarray) – Heights (m, negative) of psi

• ml_depth (float) – The mixed-layer depth H (m)

Returns:

amplitude (float) – The fitted amplitude, the stream function at z = -H/2