Project

General

Profile

Statistics over a Column

Added by Matt Thompson over 7 years ago

All,

I might just be fatigued, but I'm wondering if anyone out there has a good solution for this.

Namely, every so often in the model I work on non-zero-diff changes will be introduced. When this happens, CDO will rightly return diffn saying, nope, not the same. What I'm looking to do is make a little script that might return a number on how not-same they are.

So, I strolled through the docs and saw 'fldcor' and thought, ooh! So when I use fldcor, I get a correlation per level, nice. Now I'm wondering, what might be the best way to do statistics over this column. I know there are 'vertmean' and 'vertstd' and when I use them, they do what I want but I see:

cdo(2) vertstd1 (Warning): Layer bounds not available, using constant vertical weights for variable U!

What variable might it be looking for? Or, for my "Matt is a dunce at statistics" skill level, is this not a big worry. Obviously, when I get:

(1605) $ ./test.bash
     0.999989 ±   1.62427e-05

for an ugly quick script, things are pretty close, but it'd be nice to do the right statistics, if not just okay ones.

Thanks,
Matt


Replies (1)

RE: Statistics over a Column - Added by Ralf Mueller over 7 years ago

Hi Matthew!
what I found useful in this context is the rightmost column of the diff output: Relative Error. the absolute error can be quite different from variable to variable, but if the max relative error is 10^-12, you might be in a save area.

hth
ralf

    (1-1/1)