/[MITgcm]/MITgcm/model/inc/DYNVARS.h
ViewVC logotype

Diff of /MITgcm/model/inc/DYNVARS.h

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

revision 1.14 by adcroft, Fri Jun 29 17:14:49 2001 UTC revision 1.37 by heimbach, Fri Mar 28 18:48:05 2008 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  C  
4  C     /==========================================================\  CBOP
5  C     | DYNVARS.h                                                |  C     !ROUTINE: DYNVARS.h
6  C     | o Dynamical model variables (common block DYNVARS_R)     |  C     !INTERFACE:
7  C     |==========================================================|  C     include "DYNVARS.h"
8  C     | The value and two levels of time tendency are held for   |  C     !DESCRIPTION:
9  C     | each prognostic variable.                                |  C     \bv
10  C     \==========================================================/  C     *==========================================================*
11    C     | DYNVARS.h
12    C     | o Dynamical model variables (common block DYNVARS_R)
13    C     *==========================================================*
14    C     | The value and two levels of time tendency are held for
15    C     | each prognostic variable.
16    C     *==========================================================*
17    C     \ev
18    CEOP
19  C  C
20  C     etaN  - free-surface r-anomaly (r unit) at current time level  C     etaN  - free-surface r-anomaly (r unit) at current time level
21  C     uVel  - zonal velocity (m/s, i=1 held at western face)  C     uVel  - zonal velocity (m/s, i=1 held at western face)
22  C     vVel  - meridional velocity (m/s, j=1 held at southern face)  C     vVel  - meridional velocity (m/s, j=1 held at southern face)
23  C     theta - potential temperature (oC, held at pressure/tracer point)  C     theta - potential temperature (oC, held at pressure/tracer point)
24  C     salt  - salinity (ppt, held at pressure/tracer point)  C     salt  - salinity (ppt, held at pressure/tracer point)
25  C     gX, gXNM1 - Time tendencies at current and prvious time levels.  C     gX, gxNm1 - Time tendencies at current and previous time levels.
26  C     uVelD  - D grid zonal velocity  C     etaH   - surface r-anomaly, advanced in time consistently
27  C     vVelD  - D grid meridional velocity  C              with 2.D flow divergence (Exact-Conservation):
28    C                etaH^n+1 = etaH^n - delta_t*Div.(H^n U^n)
29    C  note: a) used with "exactConserv" but strictly necessary for NonLinFreeSurf
30    C        b) same as etaN but not necessarely at the same time, e.g.:
31    C           implicDiv2DFlow=0 => etaH=etaN ; =1 => etaH=etaNm1 ;
32    
33    #ifdef ALLOW_ADAMSBASHFORTH_3
34          COMMON /DYNVARS_R/
35         &                   etaN,
36         &                   uVel,vVel,wVel,theta,salt,
37         &                   gU,   gV,   gT,   gS,
38         &                   guNm, gvNm, gtNm, gsNm
39    #else /* ALLOW_ADAMSBASHFORTH_3 */
40        COMMON /DYNVARS_R/        COMMON /DYNVARS_R/
41       &                   etaN,etaNm1,       &                   etaN,
42       &                   uVel,vVel,wVel,theta,salt,       &                   uVel,vVel,wVel,theta,salt,
43       &                   gu,gv,gt,gs,guNm1,gvNm1,gtNm1,gsNm1       &                   gU,   gV,   gT,   gS,
44         &                   guNm1,gvNm1,gtNm1,gsNm1
45    #endif /* ALLOW_ADAMSBASHFORTH_3 */
46        _RL  etaN  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  etaN  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL  etaNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
47        _RL  uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
48        _RL  vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
49        _RL  wVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  wVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
50        _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
51        _RL  salt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  salt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
52        _RL  gu(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gU(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
53        _RL  gv(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gV(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
54        _RL  gt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
55        _RL  gs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
56    #ifdef ALLOW_ADAMSBASHFORTH_3
57          _RL  guNm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,2)
58          _RL  gvNm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,2)
59          _RL  gtNm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,2)
60          _RL  gsNm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,2)
61    #else /* ALLOW_ADAMSBASHFORTH_3 */
62        _RL  guNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  guNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
63        _RL  gvNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gvNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
64        _RL  gtNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gtNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
65        _RL  gsNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gsNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
66  #ifdef INCLUDE_CD_CODE  #endif /* ALLOW_ADAMSBASHFORTH_3 */
67        COMMON /DYNVARS_CD/  
68       &                   uVelD, vVelD,        COMMON /DYNVARS_R_2/
69       &                   uNM1,  vNM1,       &                   etaH
70       &                   guCD, gvCD        _RL  etaH  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
71        _RL  uVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
72        _RL  vVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  cph(
73        _RL  uNm1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  cph the following block will eventually move to a separate
74        _RL  vNm1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  cph header file containing requires anomaly fields of control vars.
75        _RL  guCD  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  cph
76        _RL  gvCD  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  #if (defined ALLOW_3D_DIFFKR || \
77         (defined (ALLOW_AUTODIFF) && defined (ALLOW_DIFFKR_CONTROL)))
78    C     diffKr :: full 3D specification of Laplacian diffusion coeff.
79    C               for mixing of tracers vertically ( units of r^2/s )
80          COMMON /DYNVARS_DIFFKR/
81         &                       diffKr
82          _RL  diffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
83    #endif
84    #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL))
85          COMMON /DYNVARS_KAPGM/
86         &                       kapgm
87          _RL  kapgm  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
88  #endif  #endif
89  #ifdef ALLOW_NONHYDROSTATIC  #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPREDI_CONTROL))
90        COMMON /DYNVARS_NH/ phi_nh        COMMON /DYNVARS_KAPREDI/
91        _RL  phi_nh(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)       &                       kapredi
92  #endif /* ALLOW_NONHYDROSTATIC */        _RL  kapredi  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
93    #endif
94    #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_BOTTOMDRAG_CONTROL))
95          COMMON /DYNVARS_BOTTOMDRAG/
96         &                       bottomdragfld
97          _RL  bottomdragfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
98    #endif
99    cph
100    cph)
101    
102    #ifdef ALLOW_BL79_LAT_VARY
103    C     BL79LatArray :: is used for latitudinal dependence of
104    C                     BryanLewis79 vertical diffusivity
105          COMMON /DYNVARS_BL79LatArray/ BL79LatArray
106          _RL BL79LatArray (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107    #endif
108    
109    C     diagnostic variables:
110    C     phiHydLow    :: Phi-Hydrostatic at r-lower boundary
111    C                     (bottom in z-coordinates, top in p-coordinates)
112    C     totPhiHyd    :: total hydrostatic Potential (anomaly, for now),
113    C                     at cell center level ; includes surface contribution.
114    C                     (for diagnostic + used in Z-coord with EOS_funct_P)
115    C     hMixLayer    :: Mixed layer depth [m]
116    C                     (for diagnostic + used GMRedi "fm07")
117    C     IVDConvCount :: Impl.Vert.Diffusion convection counter:
118    C                     = 0 (not convecting) or 1 (convecting)
119          COMMON /DYNVARS_DIAG/
120         &                       phiHydLow, totPhiHyd,
121         &                       hMixLayer, IVDConvCount
122          _RL  phiHydLow(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
123          _RL  totPhiHyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
124          _RL  hMixLayer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
125          _RL  IVDConvCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
126    
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.37
  ViewVC Help
Powered by ViewVC 1.1.22