/[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.37 - (hide annotations) (download)
Fri Jan 15 00:24:37 2010 UTC (14 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62y, checkpoint62x
Changes since 1.36: +2 -1 lines
add internal flag to check if adding diag to the list from the right place.

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

  ViewVC Help
Powered by ViewVC 1.1.22