[MITgcm-support] Rotating tank in Cartesian
aslavin at mun.ca
aslavin at mun.ca
Tue Jul 20 11:33:09 EDT 2010
Martin,
Thank you for your help once again. After removing external_forcing.F everything
begin to work as it must be. There are also lot of useful things inside
external_forcing.F which were not described in manual.
Alexander
MUN, NF, Canada
Quoting Martin Losch <Martin.Losch at awi.de>:
> 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>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
More information about the MITgcm-support
mailing list