/[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.4 - (hide annotations) (download)
Wed Sep 18 16:38:02 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47e_post, checkpoint46l_post, checkpoint47c_post, checkpoint48e_post, checkpoint50c_pre, checkpoint48i_post, checkpoint46l_pre, checkpoint50, checkpoint50b_pre, checkpoint48b_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint47i_post, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint46j_pre, checkpoint48h_post, checkpoint47g_post, checkpoint46j_post, checkpoint46k_post, checkpoint48a_post, checkpoint50a_post, checkpoint47j_post, branch-exfmods-tag, checkpoint48c_post, checkpoint47b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46g_post, checkpoint47f_post, checkpoint46i_post, checkpoint47, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint48g_post, checkpoint47h_post, checkpoint50b_post
Branch point for: branch-exfmods-curt
Changes since 1.3: +9 -2 lines
File MIME type: text/plain 
o Include a new diagnostic variable phiHydLow for the ocean model
  - in z-coordinates, it is the bottom pressure anomaly
  - in p-coordinates, it is the sea surface elevation
  - in both cases, these variable have global drift, reflecting the mass
    drift in z-coordinates and the volume drift in p-coordinates
  - included time averaging for phiHydLow, be aware of the drift!
o depth-dependent computation of Bo_surf for pressure coordinates
  in the ocean (buoyancyRelation='OCEANICP')
  - requires a new routine (FIND_RHO_SCALAR) to compute density with only
    Theta, Salinity, and Pressure in the parameter list. This routine is
    presently contained in find_rho.F. This routine does not give the
    correct density for 'POLY3', which would be a z-dependent reference
    density.
o cleaned up find_rho
  - removed obsolete 'eqn' from the parameter list.
o added two new verification experiments: gop and goz
  (4x4 degree global ocean, 15 layers in pressure and height coordinates)
1 mlosch 1.4 C $Header: /u/gcmpack/MITgcm/pkg/timeave/TIMEAVE_STATV.h,v 1.3 2002/01/03 16:25:44 jmc Exp $
2 jmc 1.1 C $Name: $
3 jmc 1.3
4     #ifdef ALLOW_TIMEAVE
5    
6     CBOP
7     C !ROUTINE: TIMEAVE_STATV.h
8     C !INTERFACE:
9     C include "TIMEAVE_STATV.h"
10     C !DESCRIPTION: \bw
11     C *================================================================*
12     C | TIMEAVE_STATV.h
13     C | o Time averages of model state-variables
14     C | (common block TAVE_STATEVARS)
15     C *================================================================*
16     C | Time average of state variables is (generally) centered on the
17     C | middle of the time step (time average interval = TimeAve_half)
18     C | Time average of intermediate and tandancy variables is centered
19     C | on the time step (time average interval=TimeAve_full)
20     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 mlosch 1.4 C phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential
43     C at the fixed boundary: (ocean) bottom pressure
44     C (atmos) geo- Potential
45 jmc 1.3 C ConvectCountTave :: Average number of convective adjustment event
46 jmc 1.1
47     COMMON /TAVE_TIME/ TimeAve_half,TimeAve_full
48     _RL TimeAve_half(Nr,nSx,nSy)
49     _RL TimeAve_full(Nr,nSx,nSy)
50    
51     COMMON /TAVE_STATEVARS/
52 jmc 1.3 & etatave,Eta2tave,
53     & uVeltave,vVeltave,wVeltave,
54     & thetatave,salttave,
55 adcroft 1.2 & TTtave,UUtave,VVtave,
56 jmc 1.3 & UTtave,VTtave,WTtave,
57 mlosch 1.4 & phiHydtave,
58     & phiHydLowtave,phiHydLow2Tave,
59     & ConvectCountTave
60 jmc 1.3 c & ,KEtave
61     _RL etatave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
62     _RL eta2Tave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
63 jmc 1.1 _RL uVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
64     _RL vVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
65 jmc 1.3 _RL wVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
66 jmc 1.1 _RL thetatave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
67     _RL salttave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
68 adcroft 1.2 _RL TTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
69     _RL UUtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
70     _RL VVtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
71 jmc 1.3 c _RL KEtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
72     _RL UTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
73     _RL VTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
74     _RL WTtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
75 jmc 1.1 _RL phiHydtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
76 mlosch 1.4 _RL phiHydLowtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
77     _RL phiHydLow2Tave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
78 jmc 1.1 _RL ConvectCountTave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
79 jmc 1.3
80     #ifdef NONLIN_FRSURF
81     C hUtave :: average zonal flow (=hFacW*uVel) (still in m/s !)
82     C hVtave :: average merid.flow (=hFacS*vVel) (still in m/s !)
83     C hFacCtave :: average thickness fraction of open water, Center
84     C hFacWtave :: average thickness fraction of open water, West side
85     C hFacStave :: average thickness fraction of open water, South side
86    
87     COMMON /TAVE_THICKNESS/
88     & hUtave, hVtave
89     c & , hFacCtave, hFacWtave, hFacStave
90     _RL hUtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
91     _RL hVtave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
92     c _RL hFacCtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
93     c _RL hFacWtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
94     c _RL hFacStave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
95     #endif /* NONLIN_FRSURF */
96    
97     #endif /* ALLOW_TIMEAVE */
  ViewVC Help
Powered by ViewVC 1.1.22