[MITgcm-support] Issue with Diag_stats

[Chris Horvat]

Hi all,

I've been working with the CheapAML package and wrote a diagnostic CH_SW
which diagnoses SW radiation. This part is unimportant. Additionally, I
treat CH_SW as a diagnostic in diag_stat, in order to obtain the time
series of downward SW.

Here's the issue... if I download the full CH_SW fields, the mean between
the two is the same, as it should be. However, for what appear to be
arbitrary different configurations of albedo, even if CH_SW averaged by
"hand" is the same, there is a nonzero (~.9 W/m^2) difference in the two
diag_stat outputted timeseries.

This has got to be (right?) something to do with how diagstats_calc is
doing averaging, but I'm afraid it is a little too complex for me to find
out. Any thoughts?

Chris


---
Christopher Horvat -- www.chrv.at <http://www.chrv.at/blog.html>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131030/a34cc0ee/attachment.htm>