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