/[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.4 by adcroft, Fri Feb 2 21:36:29 2001 UTC revision 1.12 by edhill, Fri Jul 29 18:24:35 2005 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #ifdef ALLOW_GMREDI  #ifdef ALLOW_GMREDI
5    
6  C     Package flag  C     Package flag
7        logical GMREDIisON        logical GMREDIisON
8        COMMON /GM_PACKAGE/ GMREDIisON        logical GM_AdvForm
9          logical GM_AdvSeparate
10          logical GM_ExtraDiag
11          logical GM_MNC
12          logical GM_MDSIO
13          COMMON /GM_PACKAGE/ GMREDIisON,
14         &                    GM_AdvForm, GM_AdvSeparate,
15         &                    GM_ExtraDiag, GM_MNC, GM_MDSIO
16    
17  C     GM/Redi parameters  C     GM/Redi parameters
18    C     GM_Small_Number :: epsilon used in computing the slope
19    C     GM_slopeSqCutoff :: slope^2 cut-off value
20          _RL GM_isopycK
21        _RL GM_background_K        _RL GM_background_K
22        _RL GM_maxSlope        _RL GM_maxSlope
23          _RL GM_Kmin_horiz
24          _RL GM_Small_Number
25          _RL GM_slopeSqCutoff
26        _RL GM_Visbeck_alpha        _RL GM_Visbeck_alpha
27        _RL GM_Visbeck_length        _RL GM_Visbeck_length
28        _RL GM_Visbeck_depth        _RL GM_Visbeck_depth
# Line 16  C     GM/Redi parameters Line 30  C     GM/Redi parameters
30        CHARACTER*(MAX_LEN_FNAM) GM_taper_scheme        CHARACTER*(MAX_LEN_FNAM) GM_taper_scheme
31        _RL GM_Scrit        _RL GM_Scrit
32        _RL GM_Sd        _RL GM_Sd
33        COMMON /GM_PARAMS/ GM_background_K,        COMMON /GM_PARAMS/ GM_isopycK, GM_background_K,
34       &                   GM_maxSlope,       &                   GM_maxSlope,
35         &                   GM_Kmin_horiz,
36         &                   GM_Small_Number, GM_slopeSqCutoff,
37       &                   GM_Visbeck_alpha,       &                   GM_Visbeck_alpha,
38       &                   GM_Visbeck_length,       &                   GM_Visbeck_length,
39       &                   GM_Visbeck_depth,       &                   GM_Visbeck_depth,
# Line 26  C     GM/Redi parameters Line 42  C     GM/Redi parameters
42       &                   GM_Scrit,       &                   GM_Scrit,
43       &                   GM_Sd       &                   GM_Sd
44    
45          _RL op5
46          _RL op25
47          PARAMETER( op5 = 0.5D0 )
48          PARAMETER( op25 = 0.25D0 )
49    
50  C     More GM/Redi parameters diagnosed from previous block  C     More GM/Redi parameters diagnosed from previous block
51  C     (not directly user configured)  C     (not directly user configured)
52        _RL GM_rMaxSlope        _RL GM_rMaxSlope
53        COMMON /GM_PARAMS2/ GM_rMaxSlope        _RL GM_skewflx
54          _RL GM_advect
55          COMMON /GM_PARAMS2/ GM_rMaxSlope,
56         &                    GM_skewflx, GM_advect
57    
58  C     GM/Redi tensor elements  C     GM/Redi tensor elements
59    
# Line 45  C     Kwz is K_33 element, Z direction a Line 67  C     Kwz is K_33 element, Z direction a
67        COMMON /GM_Wtensor/ Kwx,Kwy,Kwz        COMMON /GM_Wtensor/ Kwx,Kwy,Kwz
68    
69  #ifdef GM_NON_UNITY_DIAGONAL  #ifdef GM_NON_UNITY_DIAGONAL
70  C     First/second rows of tensor corresponds to U/V points  C     Horizontal part of the tensor
71  C     Kux is K_11 element, X direction at U point  C     Kux is K_11 element, X direction at U point
72  C     Kvy is K_22 element, Y direction at V point  C     Kvy is K_22 element, Y direction at V point
73        _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)
74        _RL Kvy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL Kvy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
75        COMMON /GM_UVtensor/ Kux,Kvy        COMMON /GM_HorTensor/ Kux,Kvy
76  #else  #else
77        _RL Kux,Kvy        _RL Kux,Kvy
78        PARAMETER(Kux=1.,Kvy=1.)        PARAMETER(Kux=1.,Kvy=1.)
79  #endif  #endif
80    
81    #ifdef GM_EXTRA_DIAGONAL
82    C     First/second rows of tensor corresponds to U/V points
83    C     Kuz is K_13 element, Z direction at U point
84    C     Kvz is K_23 element, Z direction at V point
85          _RL Kuz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
86          _RL Kvz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
87          COMMON /GM_UVtensor/ Kuz,Kvz
88    #else
89          _RL Kuz,Kvz
90          PARAMETER(Kuz=1.,Kvz=1.)
91    #endif
92    
93    #ifdef GM_BOLUS_ADVEC
94    C     GM advection formulation: bolus velocities are derived from 2
95    C        streamfunctions PsiX and PsiY :
96          _RL GM_PsiX(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
97          _RL GM_PsiY(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
98          COMMON /GM_BOLUS/ GM_PsiX,GM_PsiY
99    #endif
100    
101  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
102  C     GM mixing/stirring coefficient (spatially variable in horizontal  C     GM mixing/stirring coefficient (spatially variable in horizontal
103  C     for Visbeck et al. parameterization)  C     for Visbeck et al. parameterization)
# Line 64  C     for Visbeck et al. parameterizatio Line 106  C     for Visbeck et al. parameterizatio
106  #endif  #endif
107    
108  #endif /* ALLOW_GMREDI */  #endif /* ALLOW_GMREDI */
109    
110    CEH3 ;;; Local Variables: ***
111    CEH3 ;;; mode:fortran ***
112    CEH3 ;;; End: ***
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22