Numerical error

[logvinov]

Hi All,
I consider a simple case of wind driven flow in the channel (Wall and periodic BCs) with sloping bottom (cross the channel). So there should be no variation in the flow along the channel. But the output of the run demonstrate significant variations in the flow along the channel.

I believe there is some kind of numerical error. I use domain depth x Ly(cross the channel) x Lx : 70-73m x 250m x 100m with number of points respectively 25 x 80 x 10.

The flow near the shallower part seams to demonstrate almost x-independent behavior and it is the deeper part that have strong along the channel variations. (Looks like it seats near the periodic boundaries occupying 2 grid points and then decreasing towards the area in the middle)

Could you give me any suggestions on this?

Thanks a lot,
Evgeny

[kate]

It sounds a lot like the UPWELLING case which has been invariant along-channel when I've checked. Do you have more terms turned on than in UPWELLING? Is this serial or parallel? Is UPWELLING behaving for you? One thing to do to confirm your suspicions is to turn on the STATIONS option and put a line of stations across the periodic join. Then you can start playing games with turning off various terms until it behaves itself and you can tell us where the trouble lies.

[jcwarner]

Evgeny-

I copied the setup files you sent and i see the disturbance that you are talking about. I changed
#define MIX_GEO_UV
to
#undef MIX_GEO_UV
#define MIX_S_UV

and it seems to work correctly now.
Give that a try and post a reply here.
This should allow you to continue while we dig deeper into mix_geo_uv.

[logvinov]

Thanks a lot for the suggestions!

Now it does work when switched to
#define MIX_S_UV