/[MITgcm]/MITgcm_contrib/submesoscale/code/GMREDI.h
ViewVC logotype

Diff of /MITgcm_contrib/submesoscale/code/GMREDI.h

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

revision 1.3 by dimitri, Fri May 30 22:24:25 2008 UTC revision 1.4 by zhc, Fri Mar 12 18:31:00 2010 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.5 _d 0 )
12          PARAMETER( op25 = 0.25 _d 0 )
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 28  C     GM_slopeSqCutoff :: slope^2 cut-of Line 49  C     GM_slopeSqCutoff :: slope^2 cut-of
49        _RL GM_Visbeck_alpha        _RL GM_Visbeck_alpha
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_minDepth
53        CHARACTER*(40) GM_taper_scheme        _RL GM_Visbeck_maxSlope
54          _RL GM_Visbeck_minVal_K
55          _RL GM_Visbeck_maxVal_K
56          _RL GM_facTrL2dz
57          _RL GM_facTrL2ML
58          _RL GM_maxTransLay
59        _RL GM_Scrit        _RL GM_Scrit
60        _RL GM_Sd        _RL GM_Sd
61    C SM (4)
62        _RL GM_SM_Ce        _RL GM_SM_Ce
63        _RL GM_SM_Lf        _RL GM_SM_Lf
64        _RL GM_SM_Lmax        _RL GM_SM_Lmax
65        _RL GM_SM_tau        _RL GM_SM_tau
66        COMMON /GM_PARAMS/ GM_isopycK, GM_background_K,                          
67          COMMON /GM_PARAMS_R/
68         &                   GM_isopycK, GM_background_K,
69       &                   GM_maxSlope,       &                   GM_maxSlope,
70       &                   GM_Kmin_horiz,       &                   GM_Kmin_horiz,
71       &                   GM_Small_Number, GM_slopeSqCutoff,       &                   GM_Small_Number, GM_slopeSqCutoff,
72       &                   GM_Visbeck_alpha,       &                   GM_Visbeck_alpha, GM_Visbeck_length,
      &                   GM_Visbeck_length,  
73       &                   GM_Visbeck_depth,       &                   GM_Visbeck_depth,
74       &                   GM_Visbeck_maxval_K,       &                   GM_Visbeck_minDepth, GM_Visbeck_maxSlope,
75       &                   GM_taper_scheme,       &                   GM_Visbeck_minVal_K, GM_Visbeck_maxVal_K,
76       &                   GM_Scrit,       &                   GM_facTrL2dz, GM_facTrL2ML, GM_maxTransLay,
77       &                   GM_Sd,       &                   GM_Scrit, GM_Sd
78       &                   GM_SM_Ce,  C SM (4)
79       &                   GM_SM_Lf,       &                   ,GM_SM_Ce,
80       &                   GM_SM_Lmax,       &                   GM_SM_Lf,
81       &                   GM_SM_tau       &                   GM_SM_Lmax,
82         &                   GM_SM_tau
83        _RL op5                      
       _RL op25  
       PARAMETER( op5 = 0.5D0 )  
       PARAMETER( op25 = 0.25D0 )  
84    
85  C     More GM/Redi parameters diagnosed from previous block  C--   More GM/Redi parameters derived from previous block
86  C     (not directly user configured)  C     (not directly user configured)
87        _RL GM_rMaxSlope        _RL GM_rMaxSlope
88        _RL GM_skewflx        _RL GM_skewflx
89        _RL GM_advect        _RL GM_advect
90        COMMON /GM_PARAMS2/ GM_rMaxSlope,        COMMON /GM_DERIVED_PAR/
91       &                    GM_skewflx, GM_advect       &                   GM_rMaxSlope,
92         &                   GM_skewflx, GM_advect
93    
94  C     GM/Redi tensor elements  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
95    C---  GM/Redi tensor elements
96    
97  C     Bottom row of tensor corresponds to W points  C     Bottom row of tensor corresponds to W points
98  C     Kwx is K_31 element, X direction at W point  C     Kwx is K_31 element, X direction at W point
# Line 77  C     Kwz is K_33 element, Z direction a Line 104  C     Kwz is K_33 element, Z direction a
104        COMMON /GM_Wtensor/ Kwx,Kwy,Kwz        COMMON /GM_Wtensor/ Kwx,Kwy,Kwz
105    
106  #ifdef GM_NON_UNITY_DIAGONAL  #ifdef GM_NON_UNITY_DIAGONAL
107  C     Horizontal part of the tensor  C     Horizontal part of the tensor
108  C     Kux is K_11 element, X direction at U point  C     Kux is K_11 element, X direction at U point
109  C     Kvy is K_22 element, Y direction at V point  C     Kvy is K_22 element, Y direction at V point
110        _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 101  C     Kvz is K_23 element, Z direction a Line 128  C     Kvz is K_23 element, Z direction a
128  #endif  #endif
129    
130  #ifdef GM_BOLUS_ADVEC  #ifdef GM_BOLUS_ADVEC
131  C     GM advection formulation: bolus velocities are derived from 2  C     GM advection formulation: bolus velocities are derived from 2
132  C        streamfunctions PsiX and PsiY :  C        streamfunctions PsiX and PsiY :
133        _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)
134        _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)
135        COMMON /GM_BOLUS/ GM_PsiX,GM_PsiY        COMMON /GM_BOLUS/ GM_PsiX,GM_PsiY
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22