[MITgcm-support] Rotating tank in Cartesian

[Martin Losch]

Alexander,

if I had been familiar with the rotating tank experiment, I could have told you right away, why the rbcs-package does not do anything in this case: in the code directory there is a "private" version of external_forcing.F. This file contains the routines that apply --- surprise ---- the external forcing to the momentum and the tracer equations. Try a diff of this external_forcing.F and the one in model/src/external_forcing.F and you'll see that the private version has some extra forcing applied to the tracer equations that is specifically tailored to the experiment (in cylindrical coordinates), but that the call to rbcs_add_tendency is missing (among other things not revelant to you).

I recommend that  you remove this private routine from you code directory, do a 
"make makefile && make CLEAN && make depend && make" in your build-directory and try again. I am guessing that the strange stripes and the edge effects will go away and that you will see the effects of RBCS. If you rely on the specific code in the private version you'll have to (again by comparison with model/src) figure out where the calls to rbcs_add_tendency have to be included.

Martin

On Jul 19, 2010, at 8:55 PM, aslavin at mun.ca wrote:

> Hello Martin,
> 
> Thank you very much for your help. You was right, I have changed U and V and
> got good results: 
> [Ut,ITS]=rdmds('U',NaN);
> [Vt,ITS]=rdmds('V',NaN);
> U1=squeeze(Ut(:,:,1,2));
> V1=squeeze(Vt(:,:,1,2));
> quiver(V1,U1);
> But I still don't understand why it happens, why not U,V? :)
> There velocities near the tank edge I think due to discrete circle boundary
> cause I try square grid and they have vanish. But I don't know how to fix it..
> 
> I attach my code (RTCartesian.tar) to this letter so I will be very thankful if
> you could try it.
> 
> My second  problem concerning rbcs package. I tried to put print statement into
> rbcs_add_tendency.F, but got nothing.. so it is just not calculated new
> temperatures? I also attached my RBCS code (RBCS.tar) if it may help.. Thank
> you in advance.
> 
> Alexander
> MUN, NF, Canada
> 
> <RTCartesian.tar><RBCS.tar>