pkg/generic_advdiff/GAD_FLUX_LIMITER.h

C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/GAD_FLUX_LIMITER.h,v 1.5 2001/09/19 20:45:09 adcroft Exp $
C $Name:  $

CBOP
C !ROUTINE: GAD_FLUX_LIMITER.h

C !INTERFACE:
C #include "GAD_FLUX_LIMITER.h"
C       _RL Limiter
C       Limiter(Cr)

C !DESCRIPTION:
C Contains statement function defining limiter function.
C
C A trivial limit for the limiter recovers the upwind scheme:
C \begin{equation*}
C Limiter(Cr)=0
C \end{equation*}
C Lax-Wendroff is recovered with:
C \begin{equation*}
C Limiter(Cr)=1
C \end{equation*}
C The Min-Mod limiter is obtained with:
C \begin{equation*}
C Limiter(Cr)=max(0,min(1,Cr))
C \end{equation*}
C The current limiter of choice is the "Superbee" limiter:
C \begin{equation*}
C Limiter(Cr)=max(0,max(min(1,2*Cr),min(2,Cr)))
C \end{equation*}
C which is the default.
CEOP

C Statement function to describe flux limiter
      _RL Limiter
C Upwind        Limiter(Cr)=0.
C Lax-Wendroff  Limiter(Cr)=1.
C Min-Mod       Limiter(Cr)=max(0.,min(1.,Cr))
C Suberbee      Limiter(Cr)=max(0.,max(min(1.,2*Cr),min(2.,Cr)))

c     Limiter(Cr)=0.
c     Limiter(Cr)=1.
c     Limiter(Cr)=max(0.D0,min(1.D0,Cr))
      Limiter(Cr)=max(0.D0,max(min(1.D0,2.D0*Cr),
     &                         min(2.D0,Cr)))
1	C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/GAD_FLUX_LIMITER.h,v 1.5 2001/09/19 20:45:09 adcroft Exp $
2	C $Name: $
3
4	CBOP
5	C !ROUTINE: GAD_FLUX_LIMITER.h
6
7	C !INTERFACE:
8	C #include "GAD_FLUX_LIMITER.h"
9	C _RL Limiter
10	C Limiter(Cr)
11
12	C !DESCRIPTION:
13	C Contains statement function defining limiter function.
14	C
15	C A trivial limit for the limiter recovers the upwind scheme:
16	C \begin{equation*}
17	C Limiter(Cr)=0
18	C \end{equation*}
19	C Lax-Wendroff is recovered with:
20	C \begin{equation*}
21	C Limiter(Cr)=1
22	C \end{equation*}
23	C The Min-Mod limiter is obtained with:
24	C \begin{equation*}
25	C Limiter(Cr)=max(0,min(1,Cr))
26	C \end{equation*}
27	C The current limiter of choice is the "Superbee" limiter:
28	C \begin{equation*}
29	C Limiter(Cr)=max(0,max(min(1,2*Cr),min(2,Cr)))
30	C \end{equation*}
31	C which is the default.
32	CEOP
33
34	C Statement function to describe flux limiter
35	_RL Limiter
36	C Upwind Limiter(Cr)=0.
37	C Lax-Wendroff Limiter(Cr)=1.
38	C Min-Mod Limiter(Cr)=max(0.,min(1.,Cr))
39	C Suberbee Limiter(Cr)=max(0.,max(min(1.,2*Cr),min(2.,Cr)))
40
41	c Limiter(Cr)=0.
42	c Limiter(Cr)=1.
43	c Limiter(Cr)=max(0.D0,min(1.D0,Cr))
44	Limiter(Cr)=max(0.D0,max(min(1.D0,2.D0*Cr),
45	& min(2.D0,Cr)))