[MITgcm-devel] Bug in KPP

[David Ferreira]

Hi,
I found a bug in KPP. The bug comes in subroutine Ri_iwmix when KPP
is used with variable vertical diffusivities (for example with Bryan-Lewis).

Is anybody concerned by this ? Dimitris, you use KPP and Bryan-Lewis
in the cube-sphere, no ?
Cheers,
David


In case you care, the problem appears because the diffusivities in KPP 
are carried
in a Nr+2 long vector with k=0 at the surface while the diffusivities in 
the normal
form of the MITgcm are in a Nr-long vector with k=1 at the surface. But, 
in Ri_iwmix,
the 2 are matched with the same k...