[MITgcm-support] diagnostics: diagstats

[David Ferreira]

On 4/16/10 11:05 AM, Martin Losch wrote:
> 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?
>    
I think the volumes are accumulated in time. You would need to divide by 
the number of time-steps
in the averaging period.
david


> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>