Commit d95c01a5 authored by Dion Häfner's avatar Dion Häfner

[diagnostics] use potential density in overturning

parent a797ad3e
......@@ -72,9 +72,9 @@ class Overturning(VerosDiagnostic):
self._allocate(vs)
# sigma levels
self.sige = float(density.get_rho(vs, 35., -2., self.p_ref))
self.sigs = float(density.get_rho(vs, 35., 30., self.p_ref))
self.dsig = (self.sige - self.sigs) / (self.nlevel - 1)
self.sige = density.get_potential_rho(vs, 35., -2., press_ref=self.p_ref)
self.sigs = density.get_potential_rho(vs, 35., 30., press_ref=self.p_ref)
self.dsig = float(self.sige - self.sigs) / (self.nlevel - 1)
logger.debug(' Sigma ranges for overturning diagnostic:')
logger.debug(' Start sigma0 = {:.1f}'.format(self.sigs))
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment