[MITgcm-support] PsiVEL zero in ice shelf cavities

[Naughten, Kaitlin A.]

Just in case anyone else has this problem, I figured it out: PsiVEL is a 3D variable, not a barotropic 2D variable. So you have to put it in a diagnostic stream with other 3D variables. If you put it in a stream with 2D variables, only the surface layer will be written, and hence the cavities masked.

________________________________
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of Naughten, Kaitlin A. <kaight at bas.ac.uk>
Sent: 03 August 2018 12:02:52
To: mitgcm-support at mitgcm.org
Subject: [MITgcm-support] PsiVEL zero in ice shelf cavities

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.
________________________________
________________________________
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/20180824/d68c0a3b/attachment.html>