/[MITgcm]/MITgcm/pkg/generic_advdiff/GAD.h
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/GAD.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by adcroft, Fri Sep 28 02:26:57 2001 UTC revision 1.18 by jmc, Wed Aug 15 22:35:41 2007 UTC
# Line 16  C tokens in namelist input so for now we Line 16  C tokens in namelist input so for now we
16  C tokens value into "data" (ie. 2 for 2nd order etc.)  C tokens value into "data" (ie. 2 for 2nd order etc.)
17    
18  C !USES:  C !USES:
 #include "GAD_OPTIONS.h"  
19    
20  C !DEFINED PARAMETERS:  C !DEFINED PARAMETERS:
21    
22    C ENUM_UPWIND_1RST :: 1rst Order Upwind
23          INTEGER ENUM_UPWIND_1RST
24          PARAMETER(ENUM_UPWIND_1RST=1)
25    
26  C ENUM_CENTERED_2ND :: Centered 2nd order  C ENUM_CENTERED_2ND :: Centered 2nd order
27        INTEGER ENUM_CENTERED_2ND        INTEGER ENUM_CENTERED_2ND
28        PARAMETER(ENUM_CENTERED_2ND=2)        PARAMETER(ENUM_CENTERED_2ND=2)
29    
30  C ENUM_UPWIND_3RD :: 3rd order upwind  C ENUM_UPWIND_3RD :: 3rd order upwind
31        INTEGER ENUM_UPWIND_3RD        INTEGER ENUM_UPWIND_3RD
32        PARAMETER(ENUM_UPWIND_3RD=3)        PARAMETER(ENUM_UPWIND_3RD=3)
33    
# Line 32  C ENUM_CENTERED_4TH :: Centered 4th orde Line 35  C ENUM_CENTERED_4TH :: Centered 4th orde
35        INTEGER ENUM_CENTERED_4TH        INTEGER ENUM_CENTERED_4TH
36        PARAMETER(ENUM_CENTERED_4TH=4)        PARAMETER(ENUM_CENTERED_4TH=4)
37    
38    C ENUM_DST2 :: 2nd Order Direct Space and Time (= Lax-Wendroff)
39          INTEGER ENUM_DST2
40          PARAMETER(ENUM_DST2=20)
41    
42  C ENUM_FLUX_LIMIT :: Non-linear flux limiter  C ENUM_FLUX_LIMIT :: Non-linear flux limiter
43        INTEGER ENUM_FLUX_LIMIT        INTEGER ENUM_FLUX_LIMIT
44        PARAMETER(ENUM_FLUX_LIMIT=77)        PARAMETER(ENUM_FLUX_LIMIT=77)
45    
46  C ENUM_DST3 :: 3-DST  C ENUM_DST3 :: 3rd Order Direst Space and Time
47        INTEGER ENUM_DST3        INTEGER ENUM_DST3
48        PARAMETER(ENUM_DST3=30)        PARAMETER(ENUM_DST3=30)
49    
# Line 44  C ENUM_DST3_FLUX_LIMIT :: 3-DST flux lim Line 51  C ENUM_DST3_FLUX_LIMIT :: 3-DST flux lim
51        INTEGER ENUM_DST3_FLUX_LIMIT        INTEGER ENUM_DST3_FLUX_LIMIT
52        PARAMETER(ENUM_DST3_FLUX_LIMIT=33)        PARAMETER(ENUM_DST3_FLUX_LIMIT=33)
53    
54    C ENUM_OS7MP :: 7th Order One Step method with Monotonicity Preserving Limiter
55          INTEGER ENUM_OS7MP
56          PARAMETER(ENUM_OS7MP=7)
57    
58    C ENUM_SOM_PRATHER :: 2nd Order-Moment Advection Scheme, Prather, 1986
59          INTEGER ENUM_SOM_PRATHER
60          PARAMETER(ENUM_SOM_PRATHER=80)
61    
62    C ENUM_SOM_LIMITER :: 2nd Order-Moment Advection Scheme, Prather Limiter
63          INTEGER ENUM_SOM_LIMITER
64          PARAMETER(ENUM_SOM_LIMITER=81)
65    
66    C nSOM :: number of 1rst & 2nd Order-Moments: 1+1 (1D), 2+3 (2D), 3+6 (3D)
67          INTEGER nSOM
68          PARAMETER( nSOM = 3+6 )
69    
70  C oneSixth :: Third/fourth order interpolation factor  C oneSixth :: Third/fourth order interpolation factor
71        _RL oneSixth        _RL oneSixth
72        PARAMETER(oneSixth=1.D0/6.D0)        PARAMETER(oneSixth=1.D0/6.D0)
73    
74  C Differentiate between tracers (only needed for KPP -  arrgh!!!)  C loop range for computing vertical advection tendency
75    C  iMinAdvR,iMaxAdvR  :: 1rst index (X-dir) loop range for vertical advection
76    C  jMinAdvR,jMaxAdvR  :: 2nd  index (Y-dir) loop range for vertical advection
77          INTEGER iMinAdvR, iMaxAdvR, jMinAdvR, jMaxAdvR
78    c     PARAMETER ( iMinAdvR = 1-OLx , iMaxAdvR = sNx+OLx )
79    c     PARAMETER ( jMinAdvR = 1-OLy , jMaxAdvR = sNy+OLy )
80    C- note: we use to compute vertical advection tracer tendency everywhere
81    C        (overlap included) as above, but really needs valid tracer tendency
82    C        in interior only (as below):
83          PARAMETER ( iMinAdvR = 1 , iMaxAdvR = sNx )
84          PARAMETER ( jMinAdvR = 1 , jMaxAdvR = sNy )
85    
86    C Differentiate between tracers (needed for KPP - arrgh!!!)
87    cph                              and GMRedi arrgh*arrgh!!!)
88    cph  indices are used for TAF key computations, so need to
89    cph  running from 1, 2, ...
90    c
91  C GAD_TEMPERATURE :: temperature  C GAD_TEMPERATURE :: temperature
92        INTEGER GAD_TEMPERATURE        INTEGER GAD_TEMPERATURE
93        PARAMETER(GAD_TEMPERATURE=101)        PARAMETER(GAD_TEMPERATURE=1)
94  C GAD_TEMPERATURE :: salinity  C GAD_SALINITY :: salinity
95        INTEGER GAD_SALINITY        INTEGER GAD_SALINITY
96        PARAMETER(GAD_SALINITY=102)        PARAMETER(GAD_SALINITY=2)
97  C GAD_TR1 :: pssive tracer 1  C GAD_TR1 :: passive tracer 1
98        INTEGER GAD_TR1        INTEGER GAD_TR1
99        PARAMETER(GAD_TR1=103)        PARAMETER(GAD_TR1=3)
100  CEOP  CEOP
101    
102    C--   COMMON /GAD_PARM_L/ Logical parameters for GAD pkg routines
103    C tempSOM_Advection :: set to T if using 2nd-Order Moment advection for Temp
104    C saltSOM_Advection :: set to T if using 2nd-Order Moment advection for Salt
105    C tempMultiDimAdvec :: set to T if using multi-dim advection for Temp
106    C saltMultiDimAdvec :: set to T if using multi-dim advection for Salt
107    C AdamsBashforthGt  :: apply Adams-Bashforth extrapolation on T tendency (=Gt)
108    C AdamsBashforthGs  :: apply Adams-Bashforth extrapolation on S tendency (=Gs)
109    C AdamsBashforth_T  :: apply Adams-Bashforth extrapolation on Pot.Temp.
110    C AdamsBashforth_S  :: apply Adams-Bashforth extrapolation on Salinity
111    C tempStartAB       :: start this run (or restart) without previous gT/Temp
112    C saltStartAB       :: start this run (or restart) without previous gS/Salt
113          LOGICAL tempSOM_Advection
114          LOGICAL saltSOM_Advection
115          LOGICAL tempMultiDimAdvec
116          LOGICAL saltMultiDimAdvec
117          LOGICAL AdamsBashforthGt
118          LOGICAL AdamsBashforthGs
119          LOGICAL AdamsBashforth_T
120          LOGICAL AdamsBashforth_S
121          LOGICAL tempStartAB
122          LOGICAL saltStartAB
123          COMMON /GAD_PARM_L/
124         & tempSOM_Advection, saltSOM_Advection,
125         & tempMultiDimAdvec, saltMultiDimAdvec,
126         & AdamsBashforthGt, AdamsBashforthGs,
127         & AdamsBashforth_T, AdamsBashforth_S,
128         & tempStartAB, saltStartAB
129    
130    CEH3 ;;; Local Variables: ***
131    CEH3 ;;; mode:fortran ***
132    CEH3 ;;; End: ***
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22