Ocean Modeling Discussion

Hi,

I came to this routine by tracing backwards from the python version
in pyroms. Unfortunately, I do not have a working high order version.
And, perhaps due to the bug, I have never really understood what the
function is trying to do.

The point in Shapiro 1970-paper is to make higher order by ...

Hi,

It seems to be a bug in the matlab script utility/shapiro1.m

Higher order filtering does not work properly.
The supposedly most used 2nd order is OK as far as I can see.

The loop overwrites cor(2:Im-1) each time,
making cor(2:Im-1)=2.0.*Finp(2:Im-1) - Finp(1:Im-2) - Finp(3:Im)
independent of ...

Hi,

Not a big deal, but the documentation should be correct.
In the v momentum equation, -fu should changed to +fu

Bjørn Ådlandsvik