/[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.12 by dimitri, Wed May 5 07:15:41 2004 UTC revision 1.19 by jmc, Fri Oct 19 14:43:10 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 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!!!)  C Differentiate between tracers (needed for KPP - arrgh!!!)
87  cph                              and GMRedi arrgh*arrgh!!!)  cph                              and GMRedi arrgh*arrgh!!!)
88  cph  indices are used for TAF key computations, so need to  cph  indices are used for TAF key computations, so need to
# Line 65  C GAD_TR1 :: passive tracer 1 Line 100  C GAD_TR1 :: passive tracer 1
100  CEOP  CEOP
101    
102  C--   COMMON /GAD_PARM_L/ Logical parameters for GAD pkg routines  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  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  C saltMultiDimAdvec :: set to T if using multi-dim advection for Salt
107  C tempAdamsBashforth :: set to T if using Adams-Bashforth stepping for Temp  C AdamsBashforthGt  :: apply Adams-Bashforth extrapolation on T tendency (=Gt)
108  C saltAdamsBashforth :: set to T if using Adams-Bashforth stepping for Salt  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          LOGICAL tempSOM_Advection
112          LOGICAL saltSOM_Advection
113        LOGICAL tempMultiDimAdvec        LOGICAL tempMultiDimAdvec
114        LOGICAL saltMultiDimAdvec        LOGICAL saltMultiDimAdvec
115        LOGICAL tr1_MultiDimAdvec        LOGICAL AdamsBashforthGt
116        LOGICAL tempAdamsBashforth        LOGICAL AdamsBashforthGs
117        LOGICAL saltAdamsBashforth        LOGICAL AdamsBashforth_T
118        LOGICAL tr1_AdamsBashforth        LOGICAL AdamsBashforth_S
119        COMMON /GAD_PARM_L/        COMMON /GAD_PARM_L/
120       & tempMultiDimAdvec , saltMultiDimAdvec , tr1_MultiDimAdvec,       & tempSOM_Advection, saltSOM_Advection,
121       & tempAdamsBashforth, saltAdamsBashforth, tr1_AdamsBashforth       & tempMultiDimAdvec, saltMultiDimAdvec,
122         & AdamsBashforthGt, AdamsBashforthGs,
123         & AdamsBashforth_T, AdamsBashforth_S
124    
125  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***
126  CEH3 ;;; mode:fortran ***  CEH3 ;;; mode:fortran ***
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22