Jinlun and Martin, to wrap up the "problem with 59c" thread, Chris Hill pointed out that advcheme=7, because it is not flux-limited, requires some explicit horizontal diffusivity to be applied. Indeed doing so appears to take care of Southern Ocean instability that was causing model to blow up. Dimitris