[MITgcm-support] PsiVEL zero in ice shelf cavities

[Naughten, Kaitlin A.]

Hello,

I am running with ice shelf cavities and the horizontal velocity streamfunction (PsiVEL diagnostic) is zero in the cavities (even though velocity is nonzero). I wonder if somewhere the code is masking with hFac at k=1 instead of the real land mask. But, I've been digging through the code for a while now and can't figure out if/where this is happening. Does anybody know?

Many thanks,
Kaitlin
________________________________
This message (and any attachments) is for the recipient only. NERC is subject to the Freedom of Information Act 2000 and the contents of this email and any reply you make may be disclosed by NERC unless it is exempt from release under the Act. Any material supplied to NERC may be stored in an electronic records management system.
________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20180803/1a80d4ca/attachment.html>