/[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.24 by edhill, Thu Oct 30 12:00:41 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  C  
4  C     /==========================================================\  #include "PACKAGES_CONFIG.h"
5  C     | DYNVARS.h                                                |  
6  C     | o Dynamical model variables (common block DYNVARS_R)     |  CBOP
7  C     |==========================================================|  C     !ROUTINE: DYNVARS.h
8  C     | The value and two levels of time tendency are held for   |  C     !INTERFACE:
9  C     | each prognostic variable.                                |  C     include "DYNVARS.h"
10  C     \==========================================================/  C     !DESCRIPTION:
11    C     \bv
12    C     *==========================================================*
13    C     | DYNVARS.h                                                
14    C     | o Dynamical model variables (common block DYNVARS_R)      
15    C     *==========================================================*
16    C     | The value and two levels of time tendency are held for    
17    C     | each prognostic variable.                                
18    C     *==========================================================*
19    C     \ev
20    CEOP
21  C  C
22  C     etaN  - free-surface r-anomaly (r unit) at current time level  C     etaN  - free-surface r-anomaly (r unit) at current time level
23  C     uVel  - zonal velocity (m/s, i=1 held at western face)  C     uVel  - zonal velocity (m/s, i=1 held at western face)
# Line 15  C     vVel  - meridional velocity (m/s, Line 25  C     vVel  - meridional velocity (m/s,
25  C     theta - potential temperature (oC, held at pressure/tracer point)  C     theta - potential temperature (oC, held at pressure/tracer point)
26  C     salt  - salinity (ppt, held at pressure/tracer point)  C     salt  - salinity (ppt, held at pressure/tracer point)
27  C     gX, gXNM1 - Time tendencies at current and prvious time levels.  C     gX, gXNM1 - Time tendencies at current and prvious time levels.
28  C     uVelD  - D grid zonal velocity  C     etaH   - surface r-anomaly, advanced in time consistently
29  C     vVelD  - D grid meridional velocity  C              with 2.D flow divergence (Exact-Conservation):
30    C                etaH^n+1 = etaH^n - delta_t*Div.(H^n U^n)  
31    C  note: a) used with "exactConserv" but strictly necessary for NonLinFreeSurf
32    C        b) same as etaN but not necessarely at the same time, e.g.:
33    C           implicDiv2DFlow=0 => etaH=etaN ; =1 => etaH=etaNm1 ;
34    
35        COMMON /DYNVARS_R/        COMMON /DYNVARS_R/
36       &                   etaN,etaNm1,       &                   etaN, etaH,
37       &                   uVel,vVel,wVel,theta,salt,       &                   uVel,vVel,wVel,theta,salt,
38       &                   gu,gv,gt,gs,guNm1,gvNm1,gtNm1,gsNm1       &                   gu,gv,gt,gs,guNm1,gvNm1,gtNm1,gsNm1
39        _RL  etaN  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  etaN  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40        _RL  etaNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  etaH  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
41        _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)
42        _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)
43        _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)
# Line 37  C     vVelD  - D grid meridional velocit Line 51  C     vVelD  - D grid meridional velocit
51        _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)
52        _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)
53        _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)
54  #ifdef INCLUDE_CD_CODE  
       COMMON /DYNVARS_CD/  
      &                   uVelD, vVelD,  
      &                   uNM1,  vNM1,  
      &                   guCD, gvCD  
       _RL  uVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RL  vVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RL  uNm1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RL  vNm1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RL  guCD  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RL  gvCD  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
 #endif  
55  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
56        COMMON /DYNVARS_NH/ phi_nh        COMMON /DYNVARS_NH/ phi_nh
57        _RL  phi_nh(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  phi_nh(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
58  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
59    
60    cph(
61    cph the following block will eventually move to a separate
62    cph header file containing requires anomaly fields of control vars.
63    cph
64    #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
65          COMMON /DYNVARS_DIFFKR/
66         &                       diffKr
67          _RL  diffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
68    #endif
69    #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
70          COMMON /DYNVARS_KAPGM/
71         &                       kapgm
72          _RL  kapgm  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
73    #endif
74    #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_BOTTOMDRAG_CONTROL))
75          COMMON /DYNVARS_BOTTOMDRAG/
76         &                       bottomdragfld
77          _RL  bottomdragfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
78    #endif
79    cph
80    cph)
81    
82    C     diagnostic variables:
83    C     phiHydLow  :: Phi-Hydrostatic at r-lower boundary
84    C                  (bottom in z-coordinates, top in p-coordinates)
85    C     totPhiHyd :: total hydrostatic Potential (anomaly, for now),
86    C                  at cell center level ; includes surface contribution.
87    C                 (for diagnostic + used in Z-coord with EOS_funct_P)
88          COMMON /DYNVARS_DIAG/ phiHydLow, totPhiHyd
89          _RL  phiHydLow(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90          _RL  totPhiHyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22