/[MITgcm]/MITgcm/pkg/gmredi/GMREDI.h
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/GMREDI.h

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

revision 1.14 by gforget, Fri May 30 02:50:16 2008 UTC revision 1.15 by jmc, Tue Oct 7 20:15:27 2008 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #ifdef ALLOW_GMREDI  #ifdef ALLOW_GMREDI
5    
6  C     Package flag  C---  GM/Redi package parameters
       logical GMREDIisON  
       logical GM_AdvForm  
       logical GM_AdvSeparate  
       logical GM_ExtraDiag  
       logical GM_InMomAsStress  
       logical GM_MNC  
       logical GM_MDSIO  
       COMMON /GM_PACKAGE/ GMREDIisON,  
      &                    GM_AdvForm, GM_AdvSeparate,  
      &                    GM_ExtraDiag, GM_MNC, GM_MDSIO,  
      &                    GM_InMomAsStress  
7    
8  C     GM/Redi parameters  C--   Numerical Constant
9  C     GM_Small_Number :: epsilon used in computing the slope        _RL op5
10          _RL op25
11          PARAMETER( op5 = 0.5D0 )
12          PARAMETER( op25 = 0.25D0 )
13    
14    C--   GM/Redi Logical-type parameters
15          LOGICAL GM_AdvForm
16          LOGICAL GM_AdvSeparate
17          LOGICAL GM_ExtraDiag
18          LOGICAL GM_InMomAsStress
19          LOGICAL GM_MNC
20          LOGICAL GM_MDSIO
21          COMMON /GM_PARAMS_L/
22         &                   GM_AdvForm, GM_AdvSeparate,
23         &                   GM_ExtraDiag, GM_MNC, GM_MDSIO,
24         &                   GM_InMomAsStress
25    
26    C--   GM/Redi Character-type parameters
27    C     GM_taper_scheme  :: select which tapering/clipping scheme to use
28          CHARACTER*(40) GM_taper_scheme
29          COMMON /GM_PARAMS_C/
30         &                   GM_taper_scheme
31    
32    C--   GM/Redi real-type parameters
33    C     GM_isopycK       :: Isopycnal diffusivity [m^2/s] (Redi-tensor)
34    C     GM_background_K  :: Thickness diffusivity [m^2/s] (GM bolus transport)
35    C     GM_maxSlope      :: maximum slope (tapering/clipping) [-]
36    C     GM_Kmin_horiz    :: minimum horizontal diffusivity [m^2/s]
37    C     GM_Small_Number  :: epsilon used in computing the slope
38  C     GM_slopeSqCutoff :: slope^2 cut-off value  C     GM_slopeSqCutoff :: slope^2 cut-off value
39    C-    transition layer thickness definition:
40    C     GM_facTrL2dz   :: minimum Trans. Layer Thick. as a factor of local dz
41    C     GM_facTrL2ML   :: maximum Trans. Layer Thick. as a factor of Mix-Layer Depth
42    C     GM_maxTransLay :: maximum Trans. Layer Thick. [m]
43        _RL GM_isopycK        _RL GM_isopycK
44        _RL GM_background_K        _RL GM_background_K
45        _RL GM_maxSlope        _RL GM_maxSlope
# Line 29  C     GM_slopeSqCutoff :: slope^2 cut-of Line 50  C     GM_slopeSqCutoff :: slope^2 cut-of
50        _RL GM_Visbeck_length        _RL GM_Visbeck_length
51        _RL GM_Visbeck_depth        _RL GM_Visbeck_depth
52        _RL GM_Visbeck_maxval_K        _RL GM_Visbeck_maxval_K
53        CHARACTER*(40) GM_taper_scheme        _RL GM_facTrL2dz
54          _RL GM_facTrL2ML
55          _RL GM_maxTransLay
56        _RL GM_Scrit        _RL GM_Scrit
57        _RL GM_Sd        _RL GM_Sd
58        COMMON /GM_PARAMS/ GM_isopycK, GM_background_K,        COMMON /GM_PARAMS_R/
59         &                   GM_isopycK, GM_background_K,
60       &                   GM_maxSlope,       &                   GM_maxSlope,
61       &                   GM_Kmin_horiz,       &                   GM_Kmin_horiz,
62       &                   GM_Small_Number, GM_slopeSqCutoff,       &                   GM_Small_Number, GM_slopeSqCutoff,
# Line 40  C     GM_slopeSqCutoff :: slope^2 cut-of Line 64  C     GM_slopeSqCutoff :: slope^2 cut-of
64       &                   GM_Visbeck_length,       &                   GM_Visbeck_length,
65       &                   GM_Visbeck_depth,       &                   GM_Visbeck_depth,
66       &                   GM_Visbeck_maxval_K,       &                   GM_Visbeck_maxval_K,
67       &                   GM_taper_scheme,       &                   GM_facTrL2dz, GM_facTrL2ML, GM_maxTransLay,
68       &                   GM_Scrit,       &                   GM_Scrit, GM_Sd
      &                   GM_Sd  
   
       _RL op5  
       _RL op25  
       PARAMETER( op5 = 0.5D0 )  
       PARAMETER( op25 = 0.25D0 )  
69    
70  C     More GM/Redi parameters diagnosed from previous block  C--   More GM/Redi parameters derived from previous block
71  C     (not directly user configured)  C     (not directly user configured)
72        _RL GM_rMaxSlope        _RL GM_rMaxSlope
73        _RL GM_skewflx        _RL GM_skewflx
74        _RL GM_advect        _RL GM_advect
75        COMMON /GM_PARAMS2/ GM_rMaxSlope,        COMMON /GM_DERIVED_PAR/
76       &                    GM_skewflx, GM_advect       &                   GM_rMaxSlope,
77         &                   GM_skewflx, GM_advect
78    
79  C     GM/Redi tensor elements  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
80    C---  GM/Redi tensor elements
81    
82  C     Bottom row of tensor corresponds to W points  C     Bottom row of tensor corresponds to W points
83  C     Kwx is K_31 element, X direction at W point  C     Kwx is K_31 element, X direction at W point
# Line 69  C     Kwz is K_33 element, Z direction a Line 89  C     Kwz is K_33 element, Z direction a
89        COMMON /GM_Wtensor/ Kwx,Kwy,Kwz        COMMON /GM_Wtensor/ Kwx,Kwy,Kwz
90    
91  #ifdef GM_NON_UNITY_DIAGONAL  #ifdef GM_NON_UNITY_DIAGONAL
92  C     Horizontal part of the tensor  C     Horizontal part of the tensor
93  C     Kux is K_11 element, X direction at U point  C     Kux is K_11 element, X direction at U point
94  C     Kvy is K_22 element, Y direction at V point  C     Kvy is K_22 element, Y direction at V point
95        _RL Kux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL Kux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
# Line 93  C     Kvz is K_23 element, Z direction a Line 113  C     Kvz is K_23 element, Z direction a
113  #endif  #endif
114    
115  #ifdef GM_BOLUS_ADVEC  #ifdef GM_BOLUS_ADVEC
116  C     GM advection formulation: bolus velocities are derived from 2  C     GM advection formulation: bolus velocities are derived from 2
117  C        streamfunctions PsiX and PsiY :  C        streamfunctions PsiX and PsiY :
118        _RL GM_PsiX(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL GM_PsiX(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
119        _RL GM_PsiY(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL GM_PsiY(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
120        COMMON /GM_BOLUS/ GM_PsiX,GM_PsiY        COMMON /GM_BOLUS/ GM_PsiX,GM_PsiY
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22