[MITgcm-devel] conservation Problem with KPP viscosity

[Jean-Michel Campin]

Hi,

Seems that the problem a conservation of T & S in the donuts
is related to KPP viscosity:
When I comment out CALL KPP_CALC_VISC in calc_viscosity.F,
I get a good conservation of T (~1.e-5 J/m2) and S (~1.e-8 g/m2),
over 1 year, which is what I expect given the machine precision.
On the contrary, with KPP-viscosity on, the drift in T & S are
still small ( ~1 J/m2 & ~1.e-5 g/m2) over the same period,
but much larger than what they should be.

Did those tests with Olx=Oly=2 first, and repeat them with
 Olx=Oly=4, but did not change the conclusion.

Jean-Michel