[MITgcm-support] diagnostics: diagstats

[Martin Losch]

Hi,

two questions about the average/min/max/std/vol

1. what I find in DIAGNAME_vol is (at least for the examples that I currently care about: EXFatemp, SIarea, ETAN, THETA, so both 2D and 3D fields) consistently larger by a factor of 10 (exactly) than what I compute from the grid parameters. Is this on purpose? What am I missing? Is this a case of RTFM?

2. I was naively assuming that the averaging for atmospheric fields, such as the forcing field EXFatemp would be global, but it appears that EXFatemp is averaged over the ocean only (as an "M" field it is masked with maskC). Is that correct? How can I change that (to get truly global averages)? I saw that for fizhi and land there are special routines, but can I tweak the code to use them?

Martin