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

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

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

-revision 1.7 by adcroft,
Wed Sep 19 20:45:09 2001 UTC
+revision 1.23 by jmc,
Tue Nov 16 17:39:13 2010 UTC
 Line 15 
 C Unfortunately, there is no easy way to
  C tokens in namelist input so for now we have to enter the
  C tokens value into "data" (ie. 2 for 2nd order etc.)
+ C !USES:
  C !DEFINED PARAMETERS:
+ C ENUM_UPWIND_1RST :: 1rst Order Upwind
+       INTEGER ENUM_UPWIND_1RST
+       PARAMETER(ENUM_UPWIND_1RST=1)
  C ENUM_CENTERED_2ND :: Centered 2nd order
        INTEGER ENUM_CENTERED_2ND
        PARAMETER(ENUM_CENTERED_2ND=2)
  C ENUM_UPWIND_3RD :: 3rd order upwind
        INTEGER ENUM_UPWIND_3RD
        PARAMETER(ENUM_UPWIND_3RD=3)
-Line 29 
 C ENUM_CENTERED_4TH :: Centered 4th orde
+Line 35 
 C ENUM_CENTERED_4TH :: Centered 4th orde
        INTEGER ENUM_CENTERED_4TH
        PARAMETER(ENUM_CENTERED_4TH=4)
+ C ENUM_DST2 :: 2nd Order Direct Space and Time (= Lax-Wendroff)
+       INTEGER ENUM_DST2
+       PARAMETER(ENUM_DST2=20)
  C ENUM_FLUX_LIMIT :: Non-linear flux limiter
        INTEGER ENUM_FLUX_LIMIT
        PARAMETER(ENUM_FLUX_LIMIT=77)
- C ENUM_DST3 :: 3-DST
+ C ENUM_DST3 :: 3rd Order Direst Space and Time
        INTEGER ENUM_DST3
        PARAMETER(ENUM_DST3=30)
-Line 41 
 C ENUM_DST3_FLUX_LIMIT :: 3-DST flux lim
+Line 51 
 C ENUM_DST3_FLUX_LIMIT :: 3-DST flux lim
        INTEGER ENUM_DST3_FLUX_LIMIT
        PARAMETER(ENUM_DST3_FLUX_LIMIT=33)
+ C ENUM_OS7MP :: 7th Order One Step method with Monotonicity Preserving Limiter
+       INTEGER ENUM_OS7MP
+       PARAMETER(ENUM_OS7MP=7)
+ C ENUM_SOM_PRATHER :: 2nd Order-Moment Advection Scheme, Prather, 1986
+       INTEGER ENUM_SOM_PRATHER
+       PARAMETER(ENUM_SOM_PRATHER=80)
+ C ENUM_SOM_LIMITER :: 2nd Order-Moment Advection Scheme, Prather Limiter
+       INTEGER ENUM_SOM_LIMITER
+       PARAMETER(ENUM_SOM_LIMITER=81)
+ C GAD_Scheme_MaxNum :: maximum possible number for an advection scheme
+       INTEGER GAD_Scheme_MaxNum
+       PARAMETER( GAD_Scheme_MaxNum = 100 )
+ C nSOM :: number of 1rst & 2nd Order-Moments: 1+1 (1D), 2+3 (2D), 3+6 (3D)
+       INTEGER nSOM
+       PARAMETER( nSOM = 3+6 )
  C oneSixth :: Third/fourth order interpolation factor
        _RL oneSixth
        PARAMETER(oneSixth=1.D0/6.D0)
- C Differentiate between tracers (only needed for KPP -  arrgh!!!)
+ C loop range for computing vertical advection tendency
+ C  iMinAdvR,iMaxAdvR  :: 1rst index (X-dir) loop range for vertical advection
+ C  jMinAdvR,jMaxAdvR  :: 2nd  index (Y-dir) loop range for vertical advection
+       INTEGER iMinAdvR, iMaxAdvR, jMinAdvR, jMaxAdvR
+ c     PARAMETER ( iMinAdvR = 1-OLx , iMaxAdvR = sNx+OLx )
+ c     PARAMETER ( jMinAdvR = 1-OLy , jMaxAdvR = sNy+OLy )
+ C- note: we use to compute vertical advection tracer tendency everywhere
+ C        (overlap included) as above, but really needs valid tracer tendency
+ C        in interior only (as below):
+       PARAMETER ( iMinAdvR = 1 , iMaxAdvR = sNx )
+       PARAMETER ( jMinAdvR = 1 , jMaxAdvR = sNy )
+ C Differentiate between tracers (needed for KPP - arrgh!!!)
+ cph                              and GMRedi arrgh*arrgh!!!)
+ cph  indices are used for TAF key computations, so need to
+ cph  running from 1, 2, ...
+ c
  C GAD_TEMPERATURE :: temperature
        INTEGER GAD_TEMPERATURE
-       PARAMETER(GAD_TEMPERATURE=101)
+       PARAMETER(GAD_TEMPERATURE=1)
- C GAD_TEMPERATURE :: salinity
+ C GAD_SALINITY :: salinity
        INTEGER GAD_SALINITY
-       PARAMETER(GAD_SALINITY=102)
+       PARAMETER(GAD_SALINITY=2)
- C GAD_TR1 :: pssive tracer 1
+ C GAD_TR1 :: passive tracer 1
        INTEGER GAD_TR1
-       PARAMETER(GAD_TR1=103)
+       PARAMETER(GAD_TR1=3)
  CEOP
+ C--   COMMON /GAD_PARM_C/ Character parameters for GAD pkg routines
+ C      somSfx       :: 1rst & 2nd Order moment suffix
+       CHARACTER*2 somSfx(nSOM)
+       COMMON /GAD_PARM_C/
+      & somSfx
+ C--   COMMON /GAD_PARM_I/ Integer parameters for GAD pkg routines
+ C GAD_Scheme_olSize :: overlap minimum size for this advection scheme
+ C GAD_OlMinSize     :: overlap minimum size for GAD routines
+ C           1: min required; 2: to add to current min; 3: factor to apply
+       INTEGER GAD_Scheme_olSize(GAD_Scheme_MaxNum)
+       INTEGER GAD_OlMinSize(3)
+       COMMON /GAD_PARM_I/
+      &        GAD_Scheme_olSize,
+      &        GAD_OlMinSize
+ C--   COMMON /GAD_PARM_L/ Logical parameters for GAD pkg routines
+ C tempSOM_Advection :: set to T if using 2nd-Order Moment advection for Temp
+ C saltSOM_Advection :: set to T if using 2nd-Order Moment advection for Salt
+ C tempMultiDimAdvec :: set to T if using multi-dim advection for Temp
+ C saltMultiDimAdvec :: set to T if using multi-dim advection for Salt
+ C AdamsBashforthGt  :: apply Adams-Bashforth extrapolation on T tendency (=Gt)
+ C AdamsBashforthGs  :: apply Adams-Bashforth extrapolation on S tendency (=Gs)
+ C AdamsBashforth_T  :: apply Adams-Bashforth extrapolation on Pot.Temp.
+ C AdamsBashforth_S  :: apply Adams-Bashforth extrapolation on Salinity
+       LOGICAL tempSOM_Advection
+       LOGICAL saltSOM_Advection
+       LOGICAL tempMultiDimAdvec
+       LOGICAL saltMultiDimAdvec
+       LOGICAL AdamsBashforthGt
+       LOGICAL AdamsBashforthGs
+       LOGICAL AdamsBashforth_T
+       LOGICAL AdamsBashforth_S
+       COMMON /GAD_PARM_L/
+      & tempSOM_Advection, saltSOM_Advection,
+      & tempMultiDimAdvec, saltMultiDimAdvec,
+      & AdamsBashforthGt, AdamsBashforthGs,
+      & AdamsBashforth_T, AdamsBashforth_S
+       _RL SmolarkiewiczMaxFrac
+       COMMON /GAD_SMOL/ SmolarkiewiczMaxFrac
+ CEH3 ;;; Local Variables: ***
+ CEH3 ;;; mode:fortran ***
+ CEH3 ;;; End: ***

 Legend:



Removed from v.1.7
 


changed lines


 
Added in v.1.23
 Legend:



Removed from v.1.7
 


changed lines


 
Added in v.1.23
-Removed from v.1.7
+Added in v.1.23

	ViewVC Help
Powered by ViewVC 1.1.22