_Re:_[MITgcm-support]_Advection_schemes

Patrick Heimbach heimbach at MIT.EDU
Wed Aug 18 15:27:19 EDT 2004


Hi Petter,

sorry I oversaw your reply
(I am CC-ing to mitgcm-support and suggest
that you always send to that list).

Re. your question I think you forgot to also specify
 tracerAdvScheme = 33
in your data file.
Can you try and see whether that gives the
desired results?

-Patrick



Quoting Petter StenstrXm <petters at kth.se>:

> Hi Patrick,
> 
> thanks. 33 is apparently a flux limited 3rd order DST. It works fine with S
> and T, but not with a passive tracer. I cannot understand why, but false
> extrema appear in the TR1 field, but not in the S and T fields, when I use
> 33 for all of them. Is there an obvious explanation for this that I don't
> see?
> 
> Petter
> 
> ------- Original message -------
> From: heimbach at MIT.EDU
> Date: Tue, 17 Aug 2004 14:35:19 -0400
> 
> Hi Petter,
> 
> they are described in sections 2.16, 2.17
> http://mitgcm.org/sealion/online_documents/node71.html
> http://mitgcm.org/sealion/online_documents/node76.html
> of the online documentation,
> but for the actual number associated with each sheme
> you'll have to look in the header file
> pkg/generic_advdiff/GAD.h
> 
> -Patrick
> 
> 
> 
> Quoting Petter StenstrXm <petters at kth.se>:
> 
> > Hi!
> > 
> > Is there a list of all the different advection schemes available and their
> > corresponding numbers? I cannot find any in the documentation.
> > 
> > I have been using advection scheme 33 for some time (which I thought was a
> > flux limiter scheme), but when I look now, I cannot see that the value 33
> > is interpreted anywhere in the code.
> > 
> > Regards,
> > 
> > Petter Stenström
> > Dept of Land and Water Res Eng
> > Royal Institute of Technology
> > SE-100 44 Stockholm
> > Phone: +46 8 7908657
> > Email: petters at kth.se
> > 
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
> > 
> 
> 
> --------------------------------------------------------
> Patrick Heimbach   Massachusetts Institute of Technology
> FON: +1/617/253-5259                  EAPS, Room 54-1518
> FAX: +1/617/253-4464             77 Massachusetts Avenue
> mailto:heimbach at mit.edu               Cambridge MA 02139
> http://www.mit.edu/~heimbach/                        USA
> 
> 
> 
> 


--------------------------------------------------------
Patrick Heimbach   Massachusetts Institute of Technology
FON: +1/617/253-5259                  EAPS, Room 54-1518
FAX: +1/617/253-4464             77 Massachusetts Avenue
mailto:heimbach at mit.edu               Cambridge MA 02139
http://www.mit.edu/~heimbach/                        USA




More information about the MITgcm-support mailing list