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

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

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

Revision 1.8 - (hide annotations) (download)
Thu Oct 23 07:14:49 2003 UTC (20 years, 8 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51n_post
Branch point for: checkpoint51n_branch
Changes since 1.7: +17 -6 lines
File MIME type: text/plain 
o modifications to make FREEZE flux visible to pkg/kpp
  - moved surfaceTendencyTice from pkg/seaice to main code
  - FREEZE & EXTERNAL_FORCING_SURF moved to FORWARD_STEP
  - subroutine FREEZE now limits only surface temperature
    (this means new output.txt for global_ocean.90x40x15,
     global_ocean.cs32x15, and global_with_exf)
o added surface flux output variables to TIMEAVE_STATVARS
1 dimitri 1.8 C $Header: /usr/local/gcmpack/MITgcm/pkg/timeave/TIMEAVE_STATV.h,v 1.7 2003/10/09 04:19:20 edhill Exp $
2 jmc 1.1 C $Name: $
3 edhill 1.7
4     #include "TIMEAVE_OPTIONS.h"
5 jmc 1.3
6     #ifdef ALLOW_TIMEAVE
7    
8     CBOP
9     C !ROUTINE: TIMEAVE_STATV.h
10     C !INTERFACE:
11     C include "TIMEAVE_STATV.h"
12     C !DESCRIPTION: \bw
13     C *================================================================*
14     C | TIMEAVE_STATV.h
15     C | o Time averages of model state-variables
16     C | (common block TAVE_STATEVARS)
17     C *================================================================*
18     C | Time average of state variables is (generally) centered on the
19     C | middle of the time step (time average interval = TimeAve_half)
20 dimitri 1.8 C | Time average of intermediate and tendency variables is centered
21 jmc 1.3 C | on the time step (time average interval=TimeAve_full)
22     C *================================================================*
23     C \ev
24     CEOP
25    
26     C TimeAve_* :: time of temporal integration (s) *** for each thread ***
27     C TimeAve_half :: half time_step multiple (used for state variables)
28     C TimeAve_full :: full time_step multiple (used for for intermediate var.)
29 dimitri 1.8 C uFluxtave :: zonal surface wind stress (N/m^2,
30     C >0 for increase in uVel, i=1 held at western face)
31     C vFluxtave :: meridional surface wind stress (N/m^2,
32     C >0 for increase in vVel, j=1 held at southern face)
33     C tFluxtave :: net surface heat flux (W/m^2, >0 for increase in theta)
34     C sFluxtave :: net surface salt flux (g/m^2/s, >0 for increase in salt)
35 jmc 1.3 C etatave :: surface displacement (r unit, i.e. ocean:z, atmos:p)
36     C uVeltave :: zonal velocity (m/s, i=1 held at western face)
37     C vVeltave :: meridional velocity (m/s, j=1 held at southern face)
38     C wVeltave :: vertical velocity ([r]/s, i.e.: ocean:m/s atmos:Pa/s)
39     C thetatave :: potential temperature (oC, held at pressure/tracer point)
40     C salttave :: salinity (ppt, held at pressure/tracer point)
41 dimitri 1.8 C Eta2tave :: eta * eta
42 jmc 1.3 C TTtave :: theta * theta
43     C UUtave :: uVel * uVel (used to compute the averaged KE)
44     C VVtave :: vVel * vVel (used to compute the averaged KE)
45 jmc 1.5 C UVtave :: uVel * vVel (at vorticity point, i.e. grid-corner)
46 jmc 1.3 C KEtave :: Kinetic Energy
47     C UTtave :: uVel * theta (* hFacW)
48     C VTtave :: vVel * theta (* hFacS)
49     C WTtave :: wVel * theta
50 jmc 1.5 C tDiffRtave :: vertical diffusion flux of Temperature (theta)
51 dimitri 1.8 C uZetatave :: uVel*Relativ_Vorticity_3 (computed at v point)
52     C vZetatave :: vVel*Relativ_Vorticity_3 (computed at u point)
53 jmc 1.3 C phiHydtave :: Hydrostatic (ocean) pressure / (atmos) geo- Potential
54 mlosch 1.4 C phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential
55     C at the fixed boundary: (ocean) bottom pressure
56     C (atmos) geo- Potential
57 jmc 1.3 C ConvectCountTave :: Average number of convective adjustment event
58 jmc 1.1
59     COMMON /TAVE_TIME/ TimeAve_half,TimeAve_full
60     _RL TimeAve_half(Nr,nSx,nSy)
61     _RL TimeAve_full(Nr,nSx,nSy)
62    
63 dimitri 1.8 COMMON /TAVE_STATEVARS/
64     & uFluxtave,vFluxtave,tFluxtave,sFluxtave,
65 jmc 1.3 & etatave,Eta2tave,
66     & uVeltave,vVeltave,wVeltave,
67     & thetatave,salttave,
68 jmc 1.5 & TTtave,UUtave,VVtave,UVtave,
69     & UTtave,VTtave,WTtave,TdiffRtave,
70 mlosch 1.4 & phiHydtave,
71     & phiHydLowtave,phiHydLow2Tave,
72     & ConvectCountTave
73 jmc 1.3 c & ,KEtave
74 dimitri 1.8 _RL uFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
75     _RL vFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
76     _RL tFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
77     _RL sFluxtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
78 jmc 1.3 _RL etatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
79     _RL eta2Tave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
80 jmc 1.1 _RL uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
81     _RL vVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
82 jmc 1.3 _RL wVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
83 jmc 1.1 _RL thetatave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
84     _RL salttave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
85 adcroft 1.2 _RL TTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
86     _RL UUtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
87     _RL VVtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
88 jmc 1.5 _RL UVtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
89 jmc 1.3 c _RL KEtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
90     _RL UTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
91     _RL VTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
92     _RL WTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
93 jmc 1.5 _RL TdiffRtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
94 jmc 1.6 #ifndef DISABLE_MOM_VECINV
95     _RL uZetatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
96     _RL vZetatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
97     #endif
98 jmc 1.1 _RL phiHydtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
99 mlosch 1.4 _RL phiHydLowtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100     _RL phiHydLow2Tave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
101 jmc 1.1 _RL ConvectCountTave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
102 jmc 1.3
103     #ifdef NONLIN_FRSURF
104     C hUtave :: average zonal flow (=hFacW*uVel) (still in m/s !)
105     C hVtave :: average merid.flow (=hFacS*vVel) (still in m/s !)
106     C hFacCtave :: average thickness fraction of open water, Center
107     C hFacWtave :: average thickness fraction of open water, West side
108     C hFacStave :: average thickness fraction of open water, South side
109    
110     COMMON /TAVE_THICKNESS/
111     & hUtave, hVtave
112     c & , hFacCtave, hFacWtave, hFacStave
113     _RL hUtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
114     _RL hVtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
115     c _RL hFacCtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
116     c _RL hFacWtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
117     c _RL hFacStave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
118     #endif /* NONLIN_FRSURF */
119    
120     #endif /* ALLOW_TIMEAVE */
  ViewVC Help
Powered by ViewVC 1.1.22