[MITgcm-support] advection scheme for T/S

[samar khatiwala]

Chris, I do nothing to the diffusion. Everything is unchanged, except
the tempAdvScheme flag. If you look at the temperature fields, you will
see that switching has quite a big impact on the static stability in
the upper layers.

Samar

On Thu, 16 Sep 2004, Chris Hill wrote:

> Hi Samar,
>
>  What are you doing about diffusion when you switch?
>
> Chris
>
> > -----Original Message-----
> > From: mitgcm-support-bounces at mitgcm.org
> > [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of samar
> > khatiwala
> > Sent: Thursday, September 16, 2004 4:59 PM
> > To: mitgcm-support at mitgcm.org
> > Subject: [MITgcm-support] advection scheme for T/S
> >
> > Hello
> >
> > I am comparing a calculation in which temperature is advected
> > with the standard ABII scheme and one in which the DST3 (non
> > fluxlimited) scheme is used. I get substantially different
> > steady state solutions and was wondering if this is not
> > unreasonable or whether I should be looking for a bug in my
> > configuration.
> >
> > What I have is a sector domain (box with a ridge in the
> > middle) 10 levels, T is relaxed to a linear in y profile at
> > the surface. There is no salinity in this calculation.
> > Wind forcing is a simple cosine to drive a gyre. For the two
> > calculations all I do is change tempAdvScheme from 2 to 30.
> > There is implicit mixing.
> >
> > It would be great if someone could help me track down the
> > problem (assuming there is one). I've staged the steady state
> > solution here:
> > http://www.ldeo.columbia.edu/~spk/temperature.mat
> > (arrays are called T_AB and T_DST3)
> > Data file is here:
> > http://www.ldeo.columbia.edu/~spk/data
> >
> > Thanks
> >
> > Samar
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
> >
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>