/[MITgcm]/MITgcm/pkg/timeave/TIMEAVE_STATV.h
ViewVC logotype

Diff of /MITgcm/pkg/timeave/TIMEAVE_STATV.h

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

revision 1.2 by adcroft, Tue May 29 14:01:40 2001 UTC revision 1.3 by jmc, Thu Jan 3 16:25:44 2002 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  C  
4  C--   98/06/12 - elk  #ifdef ALLOW_TIMEAVE
5  C  
6  C     /========================================================================\  CBOP
7  C     | TIMEAVE_STATV.h                                                        |  C     !ROUTINE: TIMEAVE_STATV.h
8  C     | o Time averages of model state-variables (common block TAVE_STATEVARS) |  C     !INTERFACE:
9  C     |========================================================================|  C     include "TIMEAVE_STATV.h"
10  C  C     !DESCRIPTION: \bw
11  C     TimeAve_* - time of temporal integration (s) *** for each thread ***  C     *================================================================*
12  C     TimeAve_half - half time_step multiple (used for state variables)  C     | TIMEAVE_STATV.h
13  C     TimeAve_full - full time_step multiple (used for for intermediate var.)  C     | o Time averages of model state-variables
14  C     etaTave   - surface displacement (r unit, i.e. ocean:z, atmos:p)  C     |   (common block TAVE_STATEVARS)
15  C     uVeltave  - zonal velocity (m/s, i=1 held at western face)  C     *================================================================*
16  C     vVeltave  - meridional velocity (m/s, j=1 held at southern face)  C     | Time average of state variables is (generally) centered on the
17  C     thetatave - potential temperature (oC, held at pressure/tracer point)  C     |  middle of the time step (time average interval = TimeAve_half)
18  C     salttave  - salinity (ppt, held at pressure/tracer point)  C     | Time average of intermediate and tandancy variables is centered
19  C     uttave    - uVel * theta  C     |  on the time step (time average interval=TimeAve_full)
20  C     vttave    - vVel * theta  C     *================================================================*
21    C     \ev
22    CEOP
23    
24    C     TimeAve_*    :: time of temporal integration (s) *** for each thread ***
25    C     TimeAve_half :: half time_step multiple (used for state variables)
26    C     TimeAve_full :: full time_step multiple (used for for intermediate var.)
27    C     etatave      :: surface displacement (r unit, i.e. ocean:z, atmos:p)
28    C     uVeltave     :: zonal velocity (m/s, i=1 held at western face)
29    C     vVeltave     :: meridional velocity (m/s, j=1 held at southern face)
30    C     wVeltave     :: vertical velocity ([r]/s, i.e.: ocean:m/s atmos:Pa/s)
31    C     thetatave    :: potential temperature (oC, held at pressure/tracer point)
32    C     salttave     :: salinity (ppt, held at pressure/tracer point)
33    C     Eta2tave     ::  eta * eta
34    C     TTtave       :: theta * theta
35    C     UUtave       :: uVel * uVel (used to compute the averaged KE)
36    C     VVtave       :: vVel * vVel (used to compute the averaged KE)
37    C     KEtave       :: Kinetic Energy
38    C     UTtave       :: uVel * theta (* hFacW)
39    C     VTtave       :: vVel * theta (* hFacS)
40    C     WTtave       :: wVel * theta
41    C     phiHydtave   :: Hydrostatic (ocean) pressure / (atmos) geo- Potential
42    C     ConvectCountTave :: Average number of convective adjustment event
43    
44        COMMON /TAVE_TIME/ TimeAve_half,TimeAve_full        COMMON /TAVE_TIME/ TimeAve_half,TimeAve_full
45        _RL TimeAve_half(Nr,nSx,nSy)        _RL TimeAve_half(Nr,nSx,nSy)
46        _RL TimeAve_full(Nr,nSx,nSy)        _RL TimeAve_full(Nr,nSx,nSy)
47    
48        COMMON /TAVE_STATEVARS/        COMMON /TAVE_STATEVARS/
49       &                  etaTave,uVeltave,vVeltave,wVeltave,       &                  etatave,Eta2tave,
50       &                  thetatave,salttave,uttave,vttave,       &                  uVeltave,vVeltave,wVeltave,
51         &                  thetatave,salttave,
52       &                  TTtave,UUtave,VVtave,       &                  TTtave,UUtave,VVtave,
53         &                  UTtave,VTtave,WTtave,
54       &                  phiHydtave,ConvectCountTave       &                  phiHydtave,ConvectCountTave
55        _RL  etaTave  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  c    &                 ,KEtave
56          _RL  etatave  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
57          _RL  eta2Tave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
58        _RL  uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
59        _RL  vVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  vVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
60          _RL  wVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
61        _RL  thetatave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  thetatave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
62        _RL  salttave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  salttave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
63        _RL  TTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  TTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
64        _RL  UUtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  UUtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
65        _RL  VVtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  VVtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
66        _RL  uttave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  c     _RL  KEtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
67        _RL  vttave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  UTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
68        _RL  wVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  VTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
69          _RL  WTtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
70        _RL phiHydtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL phiHydtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
71        _RL ConvectCountTave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL ConvectCountTave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
72    
73    #ifdef NONLIN_FRSURF
74    C     hUtave       :: average zonal flow (=hFacW*uVel) (still in m/s !)
75    C     hVtave       :: average merid.flow (=hFacS*vVel) (still in m/s !)
76    C     hFacCtave    :: average thickness fraction of open water, Center
77    C     hFacWtave    :: average thickness fraction of open water, West side
78    C     hFacStave    :: average thickness fraction of open water, South side
79    
80          COMMON /TAVE_THICKNESS/
81         &              hUtave, hVtave
82    c    &            , hFacCtave, hFacWtave, hFacStave
83          _RL  hUtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
84          _RL  hVtave   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
85    c     _RL  hFacCtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
86    c     _RL  hFacWtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
87    c     _RL  hFacStave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
88    #endif /* NONLIN_FRSURF */
89    
90    #endif /* ALLOW_TIMEAVE */
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22