/[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.1 by cnh, Wed Apr 22 19:15:31 1998 UTC revision 1.18 by mlosch, Wed Sep 18 16:38:01 2002 UTC
# Line 1  Line 1 
1  C $Id$  C $Header$
2  C  C $Name$
3  C     /==========================================================\  CBOP
4  C     | DYNVARS.h                                                |  C     !ROUTINE: DYNVARS.h
5  C     | o Dynamical model variables (common block DYNVARS_R)     |  C     !INTERFACE:
6  C     |==========================================================|  C     include "DYNVARS.h"
7  C     | The value and two levels of time tendency are held for   |  C     !DESCRIPTION:
8  C     | each prognostic variable.                                |  C     \bv
9  C     \==========================================================/  C     *==========================================================*
10    C     | DYNVARS.h                                                
11    C     | o Dynamical model variables (common block DYNVARS_R)      
12    C     *==========================================================*
13    C     | The value and two levels of time tendency are held for    
14    C     | each prognostic variable.                                
15    C     *==========================================================*
16    C     \ev
17    CEOP
18  C  C
19    C     etaN  - free-surface r-anomaly (r unit) at current time level
20  C     uVel  - zonal velocity (m/s, i=1 held at western face)  C     uVel  - zonal velocity (m/s, i=1 held at western face)
21  C     vVel  - meridional velocity (m/s, j=1 held at southern face)  C     vVel  - meridional velocity (m/s, j=1 held at southern face)
22  C     theta - potential temperature (oC, held at pressure/tracer point)  C     theta - potential temperature (oC, held at pressure/tracer point)
23  C     salt  - salinity (ppt, held at pressure/tracer point)  C     salt  - salinity (ppt, held at pressure/tracer point)
24  C     rho   - density ( kg/m^3 )  C     gX, gXNM1 - Time tendencies at current and prvious time levels.
25        COMMON /DYNVARS_R/ uVel,vVel,theta,salt,rho,  C     uVelD  - D grid zonal velocity
26    C     vVelD  - D grid meridional velocity
27    C     etaH   - surface r-anomaly, advanced in time consistently
28    C              with 2.D flow divergence (Exact-Conservation):
29    C                etaH^n+1 = etaH^n - delta_t*Div.(H^n U^n)  
30    C  note: a) used with "exactConserv" but strictly necessary for NonLinFreeSurf
31    C        b) same as etaN but not necessarely at the same time, e.g.:
32    C           implicDiv2DFlow=0 => etaH=etaN ; =1 => etaH=etaNm1 ;
33    
34          COMMON /DYNVARS_R/
35         &                   etaN, etaH,
36         &                   uVel,vVel,wVel,theta,salt,
37       &                   gu,gv,gt,gs,guNm1,gvNm1,gtNm1,gsNm1       &                   gu,gv,gt,gs,guNm1,gvNm1,gtNm1,gsNm1
38  CcnhDebugStarts        _RL  etaN  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39  Cdbg &                   ,phSave        _RL  etaH  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40  CcnhDebugEnds        _RL  uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41        _RL  uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
42        _RL  vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  wVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
43        _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
44        _RL  salt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  salt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
45        _RL  rho  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gu(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
46        _RL  gu(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gv(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
47        _RL  gv(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
48        _RL  gt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
49        _RL  gs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  guNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
50        _RL  guNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gvNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
51        _RL  gvNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gtNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
52        _RL  gtNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        _RL  gsNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
53        _RL  gsNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)  
54  CcnhDebugStarts  #ifdef INCLUDE_CD_CODE
55  Cdbg  _RL  phSave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy)        COMMON /DYNVARS_CD/
56  CcnhDebugEnds       &                   uVelD, vVelD,
57         &                   etaNm1,
58         &                   uNM1,  vNM1,
59         &                   guCD, gvCD
60          _RL  uVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
61          _RL  vVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
62          _RL  etaNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
63          _RL  uNm1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
64          _RL  vNm1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
65          _RL  guCD  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
66          _RL  gvCD  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
67    #endif
68    
69    #ifdef ALLOW_NONHYDROSTATIC
70          COMMON /DYNVARS_NH/ phi_nh
71          _RL  phi_nh(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
72    #endif /* ALLOW_NONHYDROSTATIC */
73    
74    #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
75          COMMON /DYNVARS_DIFFKR/
76         &                       diffKr
77          _RL  diffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
78    #endif
79    #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
80          COMMON /DYNVARS_KAPGM/
81         &                       kapgm
82          _RL  kapgm  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
83    #endif
84    
85    C     diagnostic variable: holds phiHyd at the rigid surface
86    C     (bottom in z-coordinates, top in p-coordinates)
87          COMMON /DYNVARS_R_DIAG/ phiHydLow    
88          _RL  phiHydLow(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22