[MITgcm-support] Issue with Diag_stats

Chris Horvat horvat at fas.harvard.edu
Fri Nov 1 09:34:08 EDT 2013


This may be the issue,

I computed the average including the OBCS points on the boundary.

Thanks!

Chris

---
Christopher Horvat -- www.chrv.at <http://www.chrv.at/blog.html>



On Fri, Nov 1, 2013 at 9:30 AM, Jean-Michel Campin <jmc at ocean.mit.edu>wrote:

> 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
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131101/2df05b94/attachment.htm>


More information about the MITgcm-support mailing list