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

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

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

Revision 1.4 - (show 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 C $Header: /u/gcmpack/MITgcm/pkg/timeave/TIMEAVE_STATV.h,v 1.3 2002/01/03 16:25:44 jmc Exp $
2 C $Name: $
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 C phiHydLowtave:: Hydrostatic (ocean) pressure / (atmos) geo- Potential
43 C at the fixed boundary: (ocean) bottom pressure
44 C (atmos) geo- Potential
45 C ConvectCountTave :: Average number of convective adjustment event
46
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 & etatave,Eta2tave,
53 & uVeltave,vVeltave,wVeltave,
54 & thetatave,salttave,
55 & TTtave,UUtave,VVtave,
56 & UTtave,VTtave,WTtave,
57 & phiHydtave,
58 & phiHydLowtave,phiHydLow2Tave,
59 & ConvectCountTave
60 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 _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 _RL wVeltave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
66 _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 _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 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 _RL phiHydtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
76 _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 _RL ConvectCountTave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
79
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