/[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.36 - (hide annotations) (download)
Mon Jan 11 19:45:11 2010 UTC (14 years, 4 months ago) by jmc
Branch: MAIN
Changes since 1.35: +44 -42 lines
set new parser diag. code 3

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

  ViewVC Help
Powered by ViewVC 1.1.22