/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

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


Revision 1.38 - (hide annotations) (download)
Tue Jun 14 00:18:37 2011 UTC (12 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62z
Changes since 1.37: +20 -1 lines
first attempt to compute Velocity Potential (from UVELMASS & VVELMASS diags).

1 jmc 1.38 C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_main_init.F,v 1.37 2010/01/15 00:24:37 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "DIAG_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP 0
8     C !ROUTINE: DIAGNOSTICS_MAIN_INIT
9    
10     C !INTERFACE:
11     SUBROUTINE DIAGNOSTICS_MAIN_INIT( myThid )
12    
13     C !DESCRIPTION:
14     C Initialize available diagnostics list for variables of the main code
15     C (not part of a package): set the following attributes:
16     C name (=cdiag), parsing code (=gdiag), units (=udiag), and title (=tdiag)
17     C Notes: 1) diagnostics defined here do not require any EQUIVALENCE
18     C since they get filled-in with S/R FILL_DIAGNOSTICS
19     C 2) GDIAG is defined as character*16 and can be to character*1
20     C parse(16) with the following codes currently defined:
21    
22     C \begin{center}
23     C \begin{tabular}[h]{|c|c|}\hline
24     C \textbf{Positions} & \textbf{Characters}
25     C & \textbf{Meanings} \\\hline
26     C parse(1) & S & scalar \\
27 jmc 1.7 C & U & vector component in X direction \\
28     C & V & vector component in Y direction \\
29     C & W & vector component in vertical direction \\
30 jmc 1.1 C parse(2) & U & C-grid U-Point \\
31     C & V & C-grid V-Point \\
32     C & M & C-grid Mass Point \\
33     C & Z & C-grid Corner Point \\
34 jmc 1.36 C parse(3) & & Used for Level Integrated output: cumulate levels \\
35     C & r & same but cumulate product by model level thickness \\
36     C & R & same but cumulate product by hFac & level thickness \\
37 jmc 1.1 C parse(4) & P & positive definite \\
38     C parse(5 ) & C & with counter array \\
39 jmc 1.38 C & P & post-processed (not filled up) from other diags \\
40 jmc 1.1 C & D & disable an array for output \\
41 jmc 1.36 C parse(6--8) & '123' & retired, formerly: 3-digit mate number \\
42 jmc 1.1 C parse(9) & U & model-level plus 1/2 \\
43     C & M & model-level middle \\
44     C & L & model-level minus 1/2 \\
45     C parse(10) & 0 & levels = 0 \\
46     C & 1 & levels = 1 \\
47     C & R & levels = Nr \\
48     C & L & levels = MAX(Nr,NrPhys) \\
49     C & M & levels = MAX(Nr,NrPhys) - 1 \\
50 jmc 1.14 C & G & levels = Ground_level Number \\
51     C & I & levels = sea-Ice_level Number \\
52 jmc 1.37 C & X & free levels option (need to be set explicitly) \\
53 jmc 1.1 C \end{tabular}
54     C \end{center}
55    
56     C !USES:
57     IMPLICIT NONE
58     #include "SIZE.h"
59     #include "EEPARAMS.h"
60     #include "PARAMS.h"
61    
62     C !INPUT PARAMETERS:
63     INTEGER myThid
64     CEOP
65    
66     C !LOCAL VARIABLES:
67 jmc 1.26 C rTitle :: r-coordinate title
68     C eTitle :: free-surface title
69     C fTitle :: fixed boundary title
70     C pTitle :: "Phi" title
71     C sTitle :: "salt" title
72     INTEGER diagNum
73 jmc 1.32 INTEGER diagMate
74 jmc 1.26 CHARACTER*8 diagName
75     CHARACTER*16 diagCode
76     CHARACTER*16 diagUnits
77 jmc 1.1 CHARACTER*(80) diagTitle
78 jmc 1.26 CHARACTER*2 rUnit2c
79     CHARACTER*4 tUnit4c, sUnit4c
80     CHARACTER*(10) rTitle, eTitle, fTitle
81     CHARACTER*(20) pTitle, sTitle
82 jmc 1.8
83     CHARACTER*(16) DIAGS_MK_UNITS
84     EXTERNAL DIAGS_MK_UNITS
85 jmc 1.26 CHARACTER*(80) DIAGS_MK_TITLE
86     EXTERNAL DIAGS_MK_TITLE
87 jmc 1.1
88     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
89     C For each output variable,
90 jmc 1.26 C specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
91 jmc 1.1 C and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
92     C----------------------------------------------------------------------
93    
94 jmc 1.8 IF ( usingPCoords ) THEN
95     rUnit2c= 'Pa'
96 jmc 1.26 rTitle = ' Pressure '
97     pTitle = ' Geopotential '
98 jmc 1.8 ELSE
99     rUnit2c= 'm '
100 jmc 1.26 rTitle = ' Height '
101     pTitle = 'Pressure Pot.(p/rho)'
102 jmc 1.8 ENDIF
103     IF ( fluidIsAir ) THEN
104     tUnit4c= 'K '
105     sUnit4c= 'g/kg'
106 jmc 1.26 sTitle = ' Specific Humidity '
107 molod 1.19 #ifdef ALLOW_FIZHI
108 jmc 1.26 IF (useFIZHI) sUnit4c= 'kg/kg'
109 molod 1.19 #endif /* ALLOW_FIZHI */
110 jmc 1.8 ELSE
111     tUnit4c= 'degC'
112     sUnit4c= 'psu '
113 jmc 1.26 sTitle = ' Salinity '
114     ENDIF
115     C- free-surface (eTitle) and fixed-boundary (fTitle) position:
116     IF ( fluidIsAir ) THEN
117     eTitle = ' Surface '
118     fTitle = ' Top '
119     ELSEIF ( usingPCoords ) THEN
120     eTitle = ' Bottom '
121     fTitle = ' Surface '
122     ELSE
123     eTitle = ' Surface '
124     fTitle = ' Bottom '
125 jmc 1.8 ENDIF
126    
127 jmc 1.6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
128     C- state variables of the main code (and related quadratic var):
129    
130 jmc 1.32
131 jmc 1.6 diagName = 'ETAN '
132 jmc 1.26 diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
133     c IF ( fluidIsWater .AND. usingZCoords )
134     c &diagTitle = 'Sea Surface Elevation'
135 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
136 jmc 1.6 diagCode = 'SM M1 '
137 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
138     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
139 jmc 1.6
140     diagName = 'ETANSQ '
141 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
142     I //' Anomaly', myThid )
143 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
144 jmc 1.32 diagCode = 'SM P M1 '
145     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
146     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
147 jmc 1.6
148     diagName = 'DETADT2 '
149 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
150     I //' Anomaly Tendency', myThid )
151 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
152 jmc 1.6 diagCode = 'SM M1 '
153 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
154     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
155 jmc 1.6
156     diagName = 'THETA '
157 jmc 1.26 diagTitle = 'Potential Temperature'
158 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
159 jmc 1.36 diagCode = 'SMR MR '
160 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
161     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
162 jmc 1.6
163 jmc 1.21 c diagName = 'SST '
164     c diagTitle = 'Sea Surface Temperature (degC,K)'
165     c diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
166     c diagCode = 'SM M1 '
167 jmc 1.32 c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
168     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
169 dimitri 1.12
170 jmc 1.6 diagName = 'SALT '
171 jmc 1.26 diagTitle = DIAGS_MK_TITLE( sTitle, myThid )
172 jmc 1.8 diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
173 jmc 1.36 diagCode = 'SMR MR '
174 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
175     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
176 jmc 1.6
177 molod 1.20 diagName = 'RELHUM '
178 jmc 1.26 diagTitle = 'Relative Humidity'
179 molod 1.20 diagUnits = 'percent '
180 jmc 1.36 diagCode = 'SMR MR '
181 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
182     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
183 molod 1.20
184 jmc 1.21 c diagName = 'SSS '
185     c diagTitle = 'Sea Surface Salinity '
186     c diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
187     c diagCode = 'SM M1 '
188 jmc 1.32 c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
189     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
190 dimitri 1.12
191 jmc 1.26 IF ( fluidIsWater ) THEN
192 dimitri 1.11 diagName = 'SALTanom'
193 jmc 1.21 diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
194 dimitri 1.11 diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
195 jmc 1.36 diagCode = 'SMR MR '
196 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
197     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
198 jmc 1.26 ENDIF
199 dimitri 1.11
200 jmc 1.6 diagName = 'UVEL '
201     diagTitle = 'Zonal Component of Velocity (m/s)'
202 jmc 1.8 diagUnits = 'm/s '
203 jmc 1.36 diagCode = 'UUR MR '
204 jmc 1.32 diagMate = diagNum + 2
205     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
206     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
207 jmc 1.6
208     diagName = 'VVEL '
209     diagTitle = 'Meridional Component of Velocity (m/s)'
210 jmc 1.8 diagUnits = 'm/s '
211 jmc 1.36 diagCode = 'VVR MR '
212 jmc 1.32 diagMate = diagNum
213     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
214     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
215 jmc 1.6
216     diagName = 'WVEL '
217 jmc 1.8 diagTitle = 'Vertical Component of Velocity (r_units/s)'
218     diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
219 jmc 1.6 diagCode = 'WM LR '
220 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
221     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
222 jmc 1.6
223     diagName = 'THETASQ '
224 jmc 1.26 diagTitle = 'Square of Potential Temperature'
225 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
226 jmc 1.36 diagCode = 'SMRP MR '
227 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
229 jmc 1.6
230     diagName = 'SALTSQ '
231 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
232 jmc 1.8 diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233 jmc 1.36 diagCode = 'SMRP MR '
234 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
235     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
236 jmc 1.6
237 jmc 1.26 IF ( fluidIsWater ) THEN
238 dimitri 1.11 diagName = 'SALTSQan'
239 dimitri 1.13 diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
240 dimitri 1.11 diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
241 jmc 1.36 diagCode = 'SMRP MR '
242 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
243     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
244 jmc 1.26 ENDIF
245 dimitri 1.11
246 jmc 1.6 diagName = 'UVELSQ '
247     diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
248 jmc 1.8 diagUnits = 'm^2/s^2 '
249 jmc 1.36 diagCode = 'UURP MR '
250 jmc 1.32 diagMate = diagNum + 2
251     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
252     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
253 jmc 1.6
254     diagName = 'VVELSQ '
255     diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
256 jmc 1.8 diagUnits = 'm^2/s^2 '
257 jmc 1.36 diagCode = 'VVRP MR '
258 jmc 1.32 diagMate = diagNum
259     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
260     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
261 jmc 1.6
262     diagName = 'WVELSQ '
263 jmc 1.26 diagTitle = 'Square of Vertical Comp of Velocity'
264 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
265 jmc 1.32 diagCode = 'WM P LR '
266     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
267     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
268 jmc 1.6
269 jmc 1.34 diagName = 'UE_VEL_C'
270     diagTitle = 'Eastward Velocity (m/s) (cell center)'
271     diagUnits = 'm/s '
272 jmc 1.36 diagCode = 'UMR MR '
273 jmc 1.34 diagMate = diagNum + 2
274     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
275     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
276    
277     diagName = 'VN_VEL_C'
278     diagTitle = 'Northward Velocity (m/s) (cell center)'
279     diagUnits = 'm/s '
280 jmc 1.36 diagCode = 'VMR MR '
281 jmc 1.34 diagMate = diagNum
282     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
283     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
284    
285 jmc 1.8 diagName = 'UV_VEL_C'
286     diagTitle ='Product of horizontal Comp of velocity (cell center)'
287     diagUnits = 'm^2/s^2 '
288 jmc 1.36 diagCode = 'UMR MR '
289 jmc 1.32 diagMate = diagNum + 1
290     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
291     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
292 jmc 1.8
293     diagName = 'UV_VEL_Z'
294 jmc 1.6 diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
295 jmc 1.8 diagUnits = 'm^2/s^2 '
296 jmc 1.36 diagCode = 'UZR MR '
297 jmc 1.32 diagMate = diagNum + 1
298     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
299     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
300 jmc 1.6
301 jmc 1.9 diagName = 'WU_VEL '
302 jmc 1.26 diagTitle = 'Vertical Transport of Zonal Momentum'
303 jmc 1.9 diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
304     diagCode = 'WU LR '
305 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
306     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
307 jmc 1.9
308     diagName = 'WV_VEL '
309 jmc 1.26 diagTitle ='Vertical Transport of Meridional Momentum'
310 jmc 1.9 diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
311     diagCode = 'WV LR '
312 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
313     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
314 jmc 1.8
315 jmc 1.6 diagName = 'UVELMASS'
316     diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
317 jmc 1.8 diagUnits = 'm/s '
318 jmc 1.36 diagCode = 'UUr MR '
319 jmc 1.32 diagMate = diagNum + 2
320     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
321     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
322 jmc 1.6
323     diagName = 'VVELMASS'
324     diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
325 jmc 1.8 diagUnits = 'm/s '
326 jmc 1.36 diagCode = 'VVr MR '
327 jmc 1.32 diagMate = diagNum
328     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
329     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
330 jmc 1.6
331     diagName = 'WVELMASS'
332 jmc 1.26 diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
333 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
334 jmc 1.10 diagCode = 'WM LR '
335 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
336     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
337 jmc 1.6
338 jmc 1.38 diagName = 'PhiVEL '
339     diagTitle = 'Horizontal Velocity Potential (m^2/s)'
340     diagUnits = 'm^2/s '
341     diagCode = 'SMR P MR '
342     C- use 'UVELMASS' as mate.
343     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
344     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
345    
346     c diagName = 'PsiVEL '
347     c diagTitle = 'Horizontal Velocity Stream-Function'
348     c diagUnits = 'm^2/s '
349     c diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
350     c diagCode = 'SZ P MR '
351     C- use 'PhiVEL' as mate.
352     c diagMate = diagNum
353     c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
354     c I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
355    
356 jmc 1.6 diagName = 'UTHMASS '
357 jmc 1.26 diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
358 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
359 jmc 1.36 diagCode = 'UUr MR '
360 jmc 1.32 diagMate = diagNum + 2
361     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
362     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
363 jmc 1.6
364     diagName = 'VTHMASS '
365 jmc 1.26 diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
366 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
367 jmc 1.36 diagCode = 'VVr MR '
368 jmc 1.32 diagMate = diagNum
369     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
370     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
371 jmc 1.6
372     diagName = 'WTHMASS '
373 jmc 1.8 diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
374     diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
375 jmc 1.10 diagCode = 'WM LR '
376 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
377     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
378 jmc 1.6
379     diagName = 'USLTMASS'
380 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
381     I //sTitle, myThid )
382 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
383 jmc 1.36 diagCode = 'UUr MR '
384 jmc 1.32 diagMate = diagNum + 2
385     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
386     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
387 jmc 1.6
388     diagName = 'VSLTMASS'
389 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
390     I //sTitle, myThid )
391 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
392 jmc 1.36 diagCode = 'VVr MR '
393 jmc 1.32 diagMate = diagNum
394     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
395     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
396 jmc 1.6
397     diagName = 'WSLTMASS'
398 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
399     I //sTitle, myThid )
400 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
401 jmc 1.10 diagCode = 'WM LR '
402 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
403     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
404 jmc 1.6
405     diagName = 'UVELTH '
406 jmc 1.26 diagTitle = 'Zonal Transport of Pot Temp'
407 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
408 jmc 1.36 diagCode = 'UUR MR '
409 jmc 1.32 diagMate = diagNum + 2
410     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
411     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
412 jmc 1.6
413     diagName = 'VVELTH '
414 jmc 1.26 diagTitle = 'Meridional Transport of Pot Temp'
415 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
416 jmc 1.36 diagCode = 'VVR MR '
417 jmc 1.32 diagMate = diagNum
418     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
419     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
420 jmc 1.6
421     diagName = 'WVELTH '
422 jmc 1.26 diagTitle = 'Vertical Transport of Pot Temp'
423 jmc 1.8 diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
424 jmc 1.10 diagCode = 'WM LR '
425 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
426     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
427 jmc 1.6
428     diagName = 'UVELSLT '
429 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
430     I //sTitle, myThid )
431 jmc 1.8 diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
432 jmc 1.36 diagCode = 'UUR MR '
433 jmc 1.32 diagMate = diagNum + 2
434     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
435     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
436 jmc 1.6
437     diagName = 'VVELSLT '
438 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
439     I //sTitle, myThid )
440 jmc 1.8 diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
441 jmc 1.36 diagCode = 'VVR MR '
442 jmc 1.32 diagMate = diagNum
443     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
444     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
445 jmc 1.6
446     diagName = 'WVELSLT '
447 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
448     I //sTitle, myThid )
449 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
450 jmc 1.10 diagCode = 'WM LR '
451 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
452     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
453 jmc 1.6
454 molod 1.23 diagName = 'UVELPHI '
455 dfer 1.27 diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
456 jmc 1.26 I //pTitle//' Anomaly', myThid )
457 molod 1.23 diagUnits = 'm^3/s^3 '
458 jmc 1.36 diagCode = 'UUr MR '
459 jmc 1.32 diagMate = diagNum + 2
460     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
461     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
462 molod 1.23
463     diagName = 'VVELPHI '
464 dfer 1.27 diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
465 jmc 1.26 I //pTitle//' Anomaly', myThid )
466 molod 1.23 diagUnits = 'm^3/s^3 '
467 jmc 1.36 diagCode = 'VVr MR '
468 jmc 1.32 diagMate = diagNum
469     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
470     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
471 molod 1.23
472 jmc 1.6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
473    
474 jmc 1.7 diagName = 'RHOAnoma'
475     diagTitle = 'Density Anomaly (=Rho-rhoConst)'
476     diagUnits = 'kg/m^3 '
477 jmc 1.36 diagCode = 'SMR MR '
478 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
479     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
480 jmc 1.7
481 jmc 1.10 diagName = 'RHOANOSQ'
482     diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
483     diagUnits = 'kg^2/m^6 '
484 jmc 1.36 diagCode = 'SMRP MR '
485 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
486     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
487 jmc 1.10
488     diagName = 'URHOMASS'
489     diagTitle = 'Zonal Transport of Density'
490     diagUnits = 'kg/m^2/s '
491 jmc 1.36 diagCode = 'UUr MR '
492 jmc 1.32 diagMate = diagNum + 2
493     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
494     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
495 jmc 1.10
496     diagName = 'VRHOMASS'
497     diagTitle = 'Meridional Transport of Density'
498     diagUnits = 'kg/m^2/s '
499 jmc 1.36 diagCode = 'VVr MR '
500 jmc 1.32 diagMate = diagNum
501     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
502     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
503 jmc 1.10
504     diagName = 'WRHOMASS'
505     diagTitle = 'Vertical Transport of Potential Density'
506     diagUnits = 'kg/m^2/s '
507     diagCode = 'WM LR '
508 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
509     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
510 jmc 1.10
511 jmc 1.2 diagName = 'PHIHYD '
512 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
513     I //pTitle//' Anomaly', myThid )
514 jmc 1.8 diagUnits = 'm^2/s^2 '
515 jmc 1.36 diagCode = 'SMR MR '
516 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
517     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
518 jmc 1.2
519 molod 1.15 diagName = 'PHIHYDSQ'
520 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
521     I //pTitle//' Anomaly', myThid )
522 molod 1.15 diagUnits = 'm^4/s^4 '
523 jmc 1.36 diagCode = 'SMRP MR '
524 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
525     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
526 molod 1.15
527 jmc 1.2 diagName = 'PHIBOT '
528 jmc 1.26 c diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
529     diagTitle = DIAGS_MK_TITLE( fTitle
530     I //pTitle//' Anomaly', myThid )
531 jmc 1.8 diagUnits = 'm^2/s^2 '
532 dimitri 1.4 diagCode = 'SM M1 '
533 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
534     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
535 jmc 1.2
536     diagName = 'PHIBOTSQ'
537 jmc 1.26 c diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
538     diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
539     I //pTitle//' Anomaly', myThid )
540 jmc 1.8 diagUnits = 'm^4/s^4 '
541 jmc 1.32 diagCode = 'SM P M1 '
542     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
543     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
544 jmc 1.2
545 jmc 1.33 #ifdef ALLOW_NONHYDROSTATIC
546     diagName = 'PHI_NH '
547     diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
548     diagUnits = 'm^2/s^2 '
549 jmc 1.36 diagCode = 'SMR MR '
550 jmc 1.33 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
551     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
552     #endif /* ALLOW_NONHYDROSTATIC */
553    
554 jmc 1.31 diagName = 'MXLDEPTH'
555     diagTitle = 'Mixed-Layer Depth (>0)'
556     diagUnits = 'm '
557     diagCode = 'SM M1 '
558 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
559     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
560 jmc 1.31
561 jmc 1.1 diagName = 'DRHODR '
562 jmc 1.8 diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
563     diagUnits = 'kg/m^4 '
564     IF ( usingPCoords ) diagUnits = 's^2/m^2 '
565 jmc 1.1 diagCode = 'SM LR '
566 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
567     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
568 jmc 1.1
569 jmc 1.22 diagName = 'CONVADJ '
570     diagTitle = 'Convective Adjustment Index [0-1] '
571     diagUnits = 'fraction '
572 jmc 1.36 diagCode = 'SMR LR '
573 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
574     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
575 jmc 1.1
576 jmc 1.6 C-- surface fluxes:
577 jmc 1.28 diagName = 'oceTAUX '
578 jmc 1.26 diagTitle = 'zonal surface wind stress, >0 increases uVel'
579 jmc 1.8 diagUnits = 'N/m^2 '
580 jmc 1.32 diagCode = 'UU U1 '
581     diagMate = diagNum + 2
582     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
583     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
584 dimitri 1.5
585 jmc 1.28 diagName = 'oceTAUY '
586 jmc 1.26 diagTitle = 'meridional surf. wind stress, >0 increases vVel'
587 jmc 1.8 diagUnits = 'N/m^2 '
588 jmc 1.32 diagCode = 'VV U1 '
589     diagMate = diagNum
590     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
591     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
592 dimitri 1.5
593 jmc 1.29 diagName = 'atmPload'
594     diagTitle = 'Atmospheric pressure loading'
595     diagUnits = 'Pa '
596     diagCode = 'SM U1 '
597 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
598     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
599 jmc 1.29
600     diagName = 'sIceLoad'
601     diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
602     diagUnits = 'kg/m^2 '
603     diagCode = 'SM U1 '
604 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
605     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
606 jmc 1.29
607     diagName = 'oceFWflx'
608     diagTitle = 'net surface Fresh-Water flux into the ocean'
609     & //' (+=down), >0 decreases salinity'
610     diagUnits = 'kg/m^2/s '
611     diagCode = 'SM U1 '
612 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
613     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
614 jmc 1.29
615     diagName = 'oceSflux'
616     diagTitle = 'net surface Salt flux into the ocean (+=down),'
617     & //' >0 increases salinity'
618     diagUnits = 'g/m^2/s '
619     diagCode = 'SM U1 '
620 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
621     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
622 jmc 1.29
623     diagName = 'oceQnet '
624     diagTitle = 'net surface heat flux into the ocean (+=down),'
625     & //' >0 increases theta'
626 jmc 1.8 diagUnits = 'W/m^2 '
627 dimitri 1.5 diagCode = 'SM U1 '
628 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
629     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
630 dimitri 1.5
631 jmc 1.29 diagName = 'oceQsw '
632     diagTitle = 'net Short-Wave radiation (+=down),'
633     & //' >0 increases theta'
634 jmc 1.8 diagUnits = 'W/m^2 '
635 dimitri 1.5 diagCode = 'SM U1 '
636 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
637     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
638 dimitri 1.5
639 jmc 1.28 diagName = 'oceFreez'
640 jmc 1.29 diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
641 jmc 1.8 diagUnits = 'W/m^2 '
642 dimitri 1.5 diagCode = 'SM U1 '
643 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
644     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
645 dimitri 1.5
646 jmc 1.29 diagName = 'TRELAX '
647     diagTitle = 'surface temperature relaxation, >0 increases theta'
648     diagUnits = 'W/m^2 '
649 dimitri 1.5 diagCode = 'SM U1 '
650 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
651     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
652 dimitri 1.5
653     diagName = 'SRELAX '
654 jmc 1.26 diagTitle = 'surface salinity relaxation, >0 increases salt'
655 jmc 1.8 diagUnits = 'g/m^2/s '
656 dimitri 1.5 diagCode = 'SM U1 '
657 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
658     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
659 dimitri 1.5
660 jmc 1.29 diagName = 'surForcT'
661     diagTitle = 'model surface forcing for Temperature,'
662     & //' >0 increases theta'
663     diagUnits = 'W/m^2 '
664     diagCode = 'SM U1 '
665 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
666     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
667 jmc 1.29
668     diagName = 'surForcS'
669     diagTitle = 'model surface forcing for Salinity,'
670     & //' >0 increases salinity'
671     diagUnits = 'g/m^2/s '
672     diagCode = 'SM U1 '
673 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
674     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
675 jmc 1.29
676     diagName = 'TFLUX '
677     diagTitle = 'total heat flux (match heat-content variations),'
678     & //' >0 increases theta'
679 heimbach 1.18 diagUnits = 'W/m^2 '
680     diagCode = 'SM U1 '
681 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
682     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
683 heimbach 1.18
684 jmc 1.29 diagName = 'SFLUX '
685     diagTitle = 'total salt flux (match salt-content variations),'
686     & //' >0 increases salt'
687     diagUnits = 'g/m^2/s '
688     diagCode = 'SM U1 '
689 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
690     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
691 jmc 1.29
692 jmc 1.1 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
693 jmc 1.29
694 jmc 1.28 diagName = 'RCENTER '
695 jmc 1.26 c diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
696     diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
697 molod 1.16 diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
698     diagCode = 'SM MR '
699 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
700     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
701 jmc 1.1
702 molod 1.17 diagName = 'RSURF '
703 jmc 1.26 c diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
704     diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
705 molod 1.17 diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
706     diagCode = 'SM M1 '
707 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
708     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
709 molod 1.17
710 molod 1.24 diagName = 'TOTUTEND'
711 jmc 1.35 diagTitle = 'Tendency of Zonal Component of Velocity'
712     diagUnits = 'm/s/day '
713 jmc 1.36 diagCode = 'UUR MR '
714 jmc 1.32 diagMate = diagNum + 2
715     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
716     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
717 molod 1.24
718     diagName = 'TOTVTEND'
719 jmc 1.35 diagTitle = 'Tendency of Meridional Component of Velocity'
720     diagUnits = 'm/s/day '
721 jmc 1.36 diagCode = 'VVR MR '
722 jmc 1.32 diagMate = diagNum
723     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
724     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
725 molod 1.24
726 molod 1.25 diagName = 'TOTTTEND'
727 jmc 1.26 diagTitle = 'Tendency of Potential Temperature'
728 jmc 1.35 diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
729 jmc 1.36 diagCode = 'SMR MR '
730 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
731     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
732 molod 1.24
733     diagName = 'TOTSTEND'
734 jmc 1.26 diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
735 jmc 1.35 diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )
736 jmc 1.36 diagCode = 'SMR MR '
737 jmc 1.32 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
738     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
739 molod 1.24
740 molod 1.17
741 jmc 1.1 RETURN
742     END

  ViewVC Help
Powered by ViewVC 1.1.22