[MITgcm-support] Problem with MNC masking of ice shelf diagnostics

[Naughten, Kaitlin A.]

Hello,


I'm writing to report a small bug in the MNC masking option (useMissingValue in data.diagnostics). When I switch this option on, ice shelf diagnostics (SHIfwFlx, SHIhtFlx, SHIuStar) are incorrectly masked. The entire ice shelf is masked out, with only the open ocean (where these variables are zero) unmasked.


In the subroutine diagnostics_mnc_out, it looks like maskC is applied at klev=1, since the ice shelf variables have already been compressed into a 2D array. But in fact this 2D array represents all different values of klev, none of which are 1 beneath the ice shelf.


I'm not sure of the most efficient way to solve this problem, as I'm not very familiar with the code base yet, so I thought I'd report it here before I tried to fix it myself! If someone sees an obvious solution, I'm happy to implement it.


Many thanks,

Kaitlin Naughten

________________________________
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/20180508/1de6ccdd/attachment-0001.html>