/[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.32 - (hide annotations) (download)
Tue Feb 5 15:34:32 2008 UTC (16 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint59p, checkpoint59o, checkpoint59n
Changes since 1.31: +203 -177 lines
call DIAGNOSTICS_ADDTOLIST (= new version) to replace DIAGNOSTICS_ADD2LIST call.

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

  ViewVC Help
Powered by ViewVC 1.1.22