[MITgcm-support] KPP with inverse stratification

[Andrew Grace]

Hi all,


I'm trying to use the KPP package in a freshwater setup but the water temperatures are all below the density maximum, meaning that cooler water over warm water is stable. I noticed that in kpp_calc.F, the surface buoyancy is calculated using a hardcoded linear equation of state, and therefore, when then code determines the stability of the water column, it is deemed unstable when it should not be. This is how I understand it.


I just want to confirm that the KPP package requires modification before it can be used in situations where an inverse stratification is stable. I also apologize if this question has been asked before.


Thanks in advance.


Andrew P. Grace
Graduate Student
Department of Applied Mathematics
University of Waterloo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200103/f9383734/attachment.html>