[MITgcm-support] Issue with Diag_stats

[Jean-Michel Campin]

Hi Chris,

With stat-diags (pkg/diagnostics), the standard way statistics are computed
in an ocean set-up is over ocean only (using maskC,W,S) and, if using OBCS, 
over the interior only.  And average account for grid-ceil area.

How do you compute the average ?

Cheers,
Jean-Michel

On Wed, Oct 30, 2013 at 02:33:58PM -0400, Chris Horvat wrote:
> 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>

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support