/[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.28 - (hide annotations) (download)
Sun Nov 19 22:03:13 2006 UTC (17 years, 5 months ago) by jmc
Branch: MAIN
Changes since 1.27: +6 -6 lines
rename some diagnostics:
s/'PRESSURE'/'RCENTER '/g
s/'TICE    '/'oceFreez'/g
s/'TAUX    '/'oceTAUX '/g
s/'TAUY    '/'oceTAUY '/g
s/'SWFLUX  '/'oceQsw  '/g
s/DIFx/DFxE/g
s/DIFy/DFyE/g

1 jmc 1.28 C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_main_init.F,v 1.27 2006/08/29 17:00:30 dfer 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     CHARACTER*8 diagName
70     CHARACTER*16 diagCode
71     CHARACTER*16 diagUnits
72 jmc 1.1 CHARACTER*(80) diagTitle
73 jmc 1.26 CHARACTER*2 rUnit2c
74     CHARACTER*4 tUnit4c, sUnit4c
75     CHARACTER*(10) rTitle, eTitle, fTitle
76     CHARACTER*(20) pTitle, sTitle
77 jmc 1.8
78     CHARACTER*(16) DIAGS_MK_UNITS
79     EXTERNAL DIAGS_MK_UNITS
80 jmc 1.26 CHARACTER*(80) DIAGS_MK_TITLE
81     EXTERNAL DIAGS_MK_TITLE
82 jmc 1.1
83     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
84     C For each output variable,
85 jmc 1.26 C specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
86 jmc 1.1 C and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
87     C----------------------------------------------------------------------
88    
89 jmc 1.8 IF ( usingPCoords ) THEN
90     rUnit2c= 'Pa'
91 jmc 1.26 rTitle = ' Pressure '
92     pTitle = ' Geopotential '
93 jmc 1.8 ELSE
94     rUnit2c= 'm '
95 jmc 1.26 rTitle = ' Height '
96     pTitle = 'Pressure Pot.(p/rho)'
97 jmc 1.8 ENDIF
98     IF ( fluidIsAir ) THEN
99     tUnit4c= 'K '
100     sUnit4c= 'g/kg'
101 jmc 1.26 sTitle = ' Specific Humidity '
102 molod 1.19 #ifdef ALLOW_FIZHI
103 jmc 1.26 IF (useFIZHI) sUnit4c= 'kg/kg'
104 molod 1.19 #endif /* ALLOW_FIZHI */
105 jmc 1.8 ELSE
106     tUnit4c= 'degC'
107     sUnit4c= 'psu '
108 jmc 1.26 sTitle = ' Salinity '
109     ENDIF
110     C- free-surface (eTitle) and fixed-boundary (fTitle) position:
111     IF ( fluidIsAir ) THEN
112     eTitle = ' Surface '
113     fTitle = ' Top '
114     ELSEIF ( usingPCoords ) THEN
115     eTitle = ' Bottom '
116     fTitle = ' Surface '
117     ELSE
118     eTitle = ' Surface '
119     fTitle = ' Bottom '
120 jmc 1.8 ENDIF
121    
122 jmc 1.6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
123     C- state variables of the main code (and related quadratic var):
124    
125     diagName = 'ETAN '
126 jmc 1.26 diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
127     c IF ( fluidIsWater .AND. usingZCoords )
128     c &diagTitle = 'Sea Surface Elevation'
129 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
130 jmc 1.6 diagCode = 'SM M1 '
131     CALL DIAGNOSTICS_ADD2LIST( diagNum,
132     I diagName, diagCode, diagUnits, diagTitle, myThid )
133    
134     diagName = 'ETANSQ '
135 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
136     I //' Anomaly', myThid )
137 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
138 jmc 1.6 diagCode = 'SM M1 '
139     CALL DIAGNOSTICS_ADD2LIST( diagNum,
140     I diagName, diagCode, diagUnits, diagTitle, myThid )
141    
142     diagName = 'DETADT2 '
143 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
144     I //' Anomaly Tendency', myThid )
145 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
146 jmc 1.6 diagCode = 'SM M1 '
147     CALL DIAGNOSTICS_ADD2LIST( diagNum,
148     I diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150     diagName = 'THETA '
151 jmc 1.26 diagTitle = 'Potential Temperature'
152 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
153 jmc 1.6 diagCode = 'SM MR '
154     CALL DIAGNOSTICS_ADD2LIST( diagNum,
155     I diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157 jmc 1.21 c diagName = 'SST '
158     c diagTitle = 'Sea Surface Temperature (degC,K)'
159     c diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
160     c diagCode = 'SM M1 '
161     c CALL DIAGNOSTICS_ADD2LIST( diagNum,
162     c I diagName, diagCode, diagUnits, diagTitle, myThid )
163 dimitri 1.12
164 jmc 1.6 diagName = 'SALT '
165 jmc 1.26 diagTitle = DIAGS_MK_TITLE( sTitle, myThid )
166 jmc 1.8 diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
167 jmc 1.6 diagCode = 'SM MR '
168     CALL DIAGNOSTICS_ADD2LIST( diagNum,
169     I diagName, diagCode, diagUnits, diagTitle, myThid )
170    
171 molod 1.20 diagName = 'RELHUM '
172 jmc 1.26 diagTitle = 'Relative Humidity'
173 molod 1.20 diagUnits = 'percent '
174     diagCode = 'SM MR '
175     CALL DIAGNOSTICS_ADD2LIST( diagNum,
176     I diagName, diagCode, diagUnits, diagTitle, myThid )
177    
178 jmc 1.21 c diagName = 'SSS '
179     c diagTitle = 'Sea Surface Salinity '
180     c diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
181     c diagCode = 'SM M1 '
182     c CALL DIAGNOSTICS_ADD2LIST( diagNum,
183     c I diagName, diagCode, diagUnits, diagTitle, myThid )
184 dimitri 1.12
185 jmc 1.26 IF ( fluidIsWater ) THEN
186 dimitri 1.11 diagName = 'SALTanom'
187 jmc 1.21 diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
188 dimitri 1.11 diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
189     diagCode = 'SM MR '
190     CALL DIAGNOSTICS_ADD2LIST( diagNum,
191     I diagName, diagCode, diagUnits, diagTitle, myThid )
192 jmc 1.26 ENDIF
193 dimitri 1.11
194 jmc 1.6 diagName = 'UVEL '
195     diagTitle = 'Zonal Component of Velocity (m/s)'
196 jmc 1.8 diagUnits = 'm/s '
197 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
198     CALL DIAGNOSTICS_ADD2LIST( diagNum,
199     I diagName, diagCode, diagUnits, diagTitle, myThid )
200    
201     diagName = 'VVEL '
202     diagTitle = 'Meridional Component of Velocity (m/s)'
203 jmc 1.8 diagUnits = 'm/s '
204 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
205     CALL DIAGNOSTICS_ADD2LIST( diagNum,
206     I diagName, diagCode, diagUnits, diagTitle, myThid )
207    
208     diagName = 'WVEL '
209 jmc 1.8 diagTitle = 'Vertical Component of Velocity (r_units/s)'
210     diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
211 jmc 1.6 diagCode = 'WM LR '
212     CALL DIAGNOSTICS_ADD2LIST( diagNum,
213     I diagName, diagCode, diagUnits, diagTitle, myThid )
214    
215     diagName = 'THETASQ '
216 jmc 1.26 diagTitle = 'Square of Potential Temperature'
217 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
218 jmc 1.6 diagCode = 'SM MR '
219     CALL DIAGNOSTICS_ADD2LIST( diagNum,
220     I diagName, diagCode, diagUnits, diagTitle, myThid )
221    
222     diagName = 'SALTSQ '
223 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
224 jmc 1.8 diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
225 jmc 1.6 diagCode = 'SM MR '
226     CALL DIAGNOSTICS_ADD2LIST( diagNum,
227     I diagName, diagCode, diagUnits, diagTitle, myThid )
228    
229 jmc 1.26 IF ( fluidIsWater ) THEN
230 dimitri 1.11 diagName = 'SALTSQan'
231 dimitri 1.13 diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
232 dimitri 1.11 diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233     diagCode = 'SM MR '
234     CALL DIAGNOSTICS_ADD2LIST( diagNum,
235     I diagName, diagCode, diagUnits, diagTitle, myThid )
236 jmc 1.26 ENDIF
237 dimitri 1.11
238 jmc 1.6 diagName = 'UVELSQ '
239     diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
240 jmc 1.8 diagUnits = 'm^2/s^2 '
241 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
242     CALL DIAGNOSTICS_ADD2LIST( diagNum,
243     I diagName, diagCode, diagUnits, diagTitle, myThid )
244    
245     diagName = 'VVELSQ '
246     diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
247 jmc 1.8 diagUnits = 'm^2/s^2 '
248 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
249     CALL DIAGNOSTICS_ADD2LIST( diagNum,
250     I diagName, diagCode, diagUnits, diagTitle, myThid )
251    
252     diagName = 'WVELSQ '
253 jmc 1.26 diagTitle = 'Square of Vertical Comp of Velocity'
254 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
255 jmc 1.6 diagCode = 'WM LR '
256     CALL DIAGNOSTICS_ADD2LIST( diagNum,
257     I diagName, diagCode, diagUnits, diagTitle, myThid )
258    
259 jmc 1.8 diagName = 'UV_VEL_C'
260     diagTitle ='Product of horizontal Comp of velocity (cell center)'
261     diagUnits = 'm^2/s^2 '
262     WRITE(diagCode,'(A,I3.3,A)') 'UM ', diagNum+1 ,'MR '
263     CALL DIAGNOSTICS_ADD2LIST( diagNum,
264     I diagName, diagCode, diagUnits, diagTitle, myThid )
265    
266     diagName = 'UV_VEL_Z'
267 jmc 1.6 diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
268 jmc 1.8 diagUnits = 'm^2/s^2 '
269 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UZ ', diagNum+1 ,'MR '
270     CALL DIAGNOSTICS_ADD2LIST( diagNum,
271     I diagName, diagCode, diagUnits, diagTitle, myThid )
272    
273 jmc 1.9 diagName = 'WU_VEL '
274 jmc 1.26 diagTitle = 'Vertical Transport of Zonal Momentum'
275 jmc 1.9 diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
276     diagCode = 'WU LR '
277     CALL DIAGNOSTICS_ADD2LIST( diagNum,
278     I diagName, diagCode, diagUnits, diagTitle, myThid )
279    
280     diagName = 'WV_VEL '
281 jmc 1.26 diagTitle ='Vertical Transport of Meridional Momentum'
282 jmc 1.9 diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
283     diagCode = 'WV LR '
284     CALL DIAGNOSTICS_ADD2LIST( diagNum,
285     I diagName, diagCode, diagUnits, diagTitle, myThid )
286 jmc 1.8
287 jmc 1.6 diagName = 'UVELMASS'
288     diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
289 jmc 1.8 diagUnits = 'm/s '
290 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
291     CALL DIAGNOSTICS_ADD2LIST( diagNum,
292     I diagName, diagCode, diagUnits, diagTitle, myThid )
293    
294     diagName = 'VVELMASS'
295     diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
296 jmc 1.8 diagUnits = 'm/s '
297 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
298     CALL DIAGNOSTICS_ADD2LIST( diagNum,
299     I diagName, diagCode, diagUnits, diagTitle, myThid )
300    
301     diagName = 'WVELMASS'
302 jmc 1.26 diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
303 jmc 1.8 diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
304 jmc 1.10 diagCode = 'WM LR '
305 jmc 1.6 CALL DIAGNOSTICS_ADD2LIST( diagNum,
306     I diagName, diagCode, diagUnits, diagTitle, myThid )
307    
308     diagName = 'UTHMASS '
309 jmc 1.26 diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
310 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
311 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
312     CALL DIAGNOSTICS_ADD2LIST( diagNum,
313     I diagName, diagCode, diagUnits, diagTitle, myThid )
314    
315     diagName = 'VTHMASS '
316 jmc 1.26 diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
317 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
318 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
319     CALL DIAGNOSTICS_ADD2LIST( diagNum,
320     I diagName, diagCode, diagUnits, diagTitle, myThid )
321    
322     diagName = 'WTHMASS '
323 jmc 1.8 diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
324     diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
325 jmc 1.10 diagCode = 'WM LR '
326 jmc 1.6 CALL DIAGNOSTICS_ADD2LIST( diagNum,
327     I diagName, diagCode, diagUnits, diagTitle, myThid )
328    
329     diagName = 'USLTMASS'
330 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
331     I //sTitle, myThid )
332 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
333 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
334     CALL DIAGNOSTICS_ADD2LIST( diagNum,
335     I diagName, diagCode, diagUnits, diagTitle, myThid )
336    
337     diagName = 'VSLTMASS'
338 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
339     I //sTitle, myThid )
340 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
341 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
342     CALL DIAGNOSTICS_ADD2LIST( diagNum,
343     I diagName, diagCode, diagUnits, diagTitle, myThid )
344    
345     diagName = 'WSLTMASS'
346 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
347     I //sTitle, myThid )
348 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
349 jmc 1.10 diagCode = 'WM LR '
350 jmc 1.6 CALL DIAGNOSTICS_ADD2LIST( diagNum,
351     I diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353     diagName = 'UVELTH '
354 jmc 1.26 diagTitle = 'Zonal Transport of Pot Temp'
355 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
356 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
357     CALL DIAGNOSTICS_ADD2LIST( diagNum,
358     I diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360     diagName = 'VVELTH '
361 jmc 1.26 diagTitle = 'Meridional Transport of Pot Temp'
362 jmc 1.8 diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
363 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
364     CALL DIAGNOSTICS_ADD2LIST( diagNum,
365     I diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367     diagName = 'WVELTH '
368 jmc 1.26 diagTitle = 'Vertical Transport of Pot Temp'
369 jmc 1.8 diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
370 jmc 1.10 diagCode = 'WM LR '
371 jmc 1.6 CALL DIAGNOSTICS_ADD2LIST( diagNum,
372     I diagName, diagCode, diagUnits, diagTitle, myThid )
373    
374     diagName = 'UVELSLT '
375 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
376     I //sTitle, myThid )
377 jmc 1.8 diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
378 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
379     CALL DIAGNOSTICS_ADD2LIST( diagNum,
380     I diagName, diagCode, diagUnits, diagTitle, myThid )
381    
382     diagName = 'VVELSLT '
383 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
384     I //sTitle, myThid )
385 jmc 1.8 diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
386 jmc 1.6 WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
387     CALL DIAGNOSTICS_ADD2LIST( diagNum,
388     I diagName, diagCode, diagUnits, diagTitle, myThid )
389    
390     diagName = 'WVELSLT '
391 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
392     I //sTitle, myThid )
393 jmc 1.8 diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
394 jmc 1.10 diagCode = 'WM LR '
395 jmc 1.6 CALL DIAGNOSTICS_ADD2LIST( diagNum,
396     I diagName, diagCode, diagUnits, diagTitle, myThid )
397    
398 molod 1.23 diagName = 'UVELPHI '
399 dfer 1.27 diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
400 jmc 1.26 I //pTitle//' Anomaly', myThid )
401 molod 1.23 diagUnits = 'm^3/s^3 '
402     WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
403     CALL DIAGNOSTICS_ADD2LIST( diagNum,
404     I diagName, diagCode, diagUnits, diagTitle, myThid )
405    
406     diagName = 'VVELPHI '
407 dfer 1.27 diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
408 jmc 1.26 I //pTitle//' Anomaly', myThid )
409 molod 1.23 diagUnits = 'm^3/s^3 '
410     WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
411     CALL DIAGNOSTICS_ADD2LIST( diagNum,
412     I diagName, diagCode, diagUnits, diagTitle, myThid )
413    
414 jmc 1.6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
415    
416 jmc 1.7 diagName = 'RHOAnoma'
417     diagTitle = 'Density Anomaly (=Rho-rhoConst)'
418     diagUnits = 'kg/m^3 '
419     diagCode = 'SM MR '
420     CALL DIAGNOSTICS_ADD2LIST( diagNum,
421     I diagName, diagCode, diagUnits, diagTitle, myThid )
422    
423 jmc 1.10 diagName = 'RHOANOSQ'
424     diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
425     diagUnits = 'kg^2/m^6 '
426     diagCode = 'SM MR '
427     CALL DIAGNOSTICS_ADD2LIST( diagNum,
428     I diagName, diagCode, diagUnits, diagTitle, myThid )
429    
430     diagName = 'URHOMASS'
431     diagTitle = 'Zonal Transport of Density'
432     diagUnits = 'kg/m^2/s '
433     WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
434     CALL DIAGNOSTICS_ADD2LIST( diagNum,
435     I diagName, diagCode, diagUnits, diagTitle, myThid )
436    
437     diagName = 'VRHOMASS'
438     diagTitle = 'Meridional Transport of Density'
439     diagUnits = 'kg/m^2/s '
440     WRITE(diagCode,'(A,I3.3,A)') 'VV ',diagNum,'MR '
441     CALL DIAGNOSTICS_ADD2LIST( diagNum,
442     I diagName, diagCode, diagUnits, diagTitle, myThid )
443    
444     diagName = 'WRHOMASS'
445     diagTitle = 'Vertical Transport of Potential Density'
446     diagUnits = 'kg/m^2/s '
447     diagCode = 'WM LR '
448     CALL DIAGNOSTICS_ADD2LIST( diagNum,
449     I diagName, diagCode, diagUnits, diagTitle, myThid )
450    
451 jmc 1.2 diagName = 'PHIHYD '
452 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
453     I //pTitle//' Anomaly', myThid )
454 jmc 1.8 diagUnits = 'm^2/s^2 '
455 jmc 1.2 diagCode = 'SM MR '
456 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
457 jmc 1.2 I diagName, diagCode, diagUnits, diagTitle, myThid )
458    
459 molod 1.15 diagName = 'PHIHYDSQ'
460 jmc 1.26 diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
461     I //pTitle//' Anomaly', myThid )
462 molod 1.15 diagUnits = 'm^4/s^4 '
463     diagCode = 'SM MR '
464     CALL DIAGNOSTICS_ADD2LIST( diagNum,
465     I diagName, diagCode, diagUnits, diagTitle, myThid )
466    
467 jmc 1.2 diagName = 'PHIBOT '
468 jmc 1.26 c diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
469     diagTitle = DIAGS_MK_TITLE( fTitle
470     I //pTitle//' Anomaly', myThid )
471 jmc 1.8 diagUnits = 'm^2/s^2 '
472 dimitri 1.4 diagCode = 'SM M1 '
473 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
474 jmc 1.2 I diagName, diagCode, diagUnits, diagTitle, myThid )
475    
476     diagName = 'PHIBOTSQ'
477 jmc 1.26 c diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
478     diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
479     I //pTitle//' Anomaly', myThid )
480 jmc 1.8 diagUnits = 'm^4/s^4 '
481 jmc 1.2 diagCode = 'SM M1 '
482 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
483 jmc 1.2 I diagName, diagCode, diagUnits, diagTitle, myThid )
484    
485 jmc 1.1 diagName = 'DRHODR '
486 jmc 1.8 diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
487     diagUnits = 'kg/m^4 '
488     IF ( usingPCoords ) diagUnits = 's^2/m^2 '
489 jmc 1.1 diagCode = 'SM LR '
490 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
491 jmc 1.1 I diagName, diagCode, diagUnits, diagTitle, myThid )
492    
493 jmc 1.22 diagName = 'CONVADJ '
494     diagTitle = 'Convective Adjustment Index [0-1] '
495     diagUnits = 'fraction '
496     diagCode = 'SM LR '
497 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
498 jmc 1.1 I diagName, diagCode, diagUnits, diagTitle, myThid )
499    
500 jmc 1.6 C-- surface fluxes:
501 jmc 1.28 diagName = 'oceTAUX '
502 jmc 1.26 diagTitle = 'zonal surface wind stress, >0 increases uVel'
503 jmc 1.8 diagUnits = 'N/m^2 '
504 dimitri 1.5 diagCode = 'SU U1 '
505     CALL DIAGNOSTICS_ADD2LIST( diagNum,
506     I diagName, diagCode, diagUnits, diagTitle, myThid )
507    
508 jmc 1.28 diagName = 'oceTAUY '
509 jmc 1.26 diagTitle = 'meridional surf. wind stress, >0 increases vVel'
510 jmc 1.8 diagUnits = 'N/m^2 '
511 dimitri 1.5 diagCode = 'SV U1 '
512     CALL DIAGNOSTICS_ADD2LIST( diagNum,
513     I diagName, diagCode, diagUnits, diagTitle, myThid )
514    
515     diagName = 'TFLUX '
516     diagTitle = 'net surface heat flux, >0 increases theta '
517 jmc 1.8 diagUnits = 'W/m^2 '
518 dimitri 1.5 diagCode = 'SM U1 '
519     CALL DIAGNOSTICS_ADD2LIST( diagNum,
520     I diagName, diagCode, diagUnits, diagTitle, myThid )
521    
522     diagName = 'TRELAX '
523     diagTitle = 'surface temperature relaxation, >0 increases theta'
524 jmc 1.8 diagUnits = 'W/m^2 '
525 dimitri 1.5 diagCode = 'SM U1 '
526     CALL DIAGNOSTICS_ADD2LIST( diagNum,
527     I diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529 jmc 1.28 diagName = 'oceFreez'
530 dimitri 1.5 diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
531 jmc 1.8 diagUnits = 'W/m^2 '
532 dimitri 1.5 diagCode = 'SM U1 '
533     CALL DIAGNOSTICS_ADD2LIST( diagNum,
534     I diagName, diagCode, diagUnits, diagTitle, myThid )
535    
536     diagName = 'SFLUX '
537 jmc 1.26 diagTitle = 'net surface salt flux, >0 increases salt'
538 jmc 1.8 diagUnits = 'g/m^2/s '
539 dimitri 1.5 diagCode = 'SM U1 '
540     CALL DIAGNOSTICS_ADD2LIST( diagNum,
541     I diagName, diagCode, diagUnits, diagTitle, myThid )
542    
543     diagName = 'SRELAX '
544 jmc 1.26 diagTitle = 'surface salinity relaxation, >0 increases salt'
545 jmc 1.8 diagUnits = 'g/m^2/s '
546 dimitri 1.5 diagCode = 'SM U1 '
547     CALL DIAGNOSTICS_ADD2LIST( diagNum,
548     I diagName, diagCode, diagUnits, diagTitle, myThid )
549    
550 jmc 1.28 diagName = 'oceQsw '
551 jmc 1.26 diagTitle = 'net upward SW radiation, >0 increases theta'
552 heimbach 1.18 diagUnits = 'W/m^2 '
553     diagCode = 'SM U1 '
554     CALL DIAGNOSTICS_ADD2LIST( diagNum,
555     I diagName, diagCode, diagUnits, diagTitle, myThid )
556    
557 jmc 1.1 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
558 jmc 1.28 diagName = 'RCENTER '
559 jmc 1.26 C- will be changed to an isomorphic name, such as:
560     c diagName = 'RCENTER '
561     c diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
562     diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
563 molod 1.16 diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
564     diagCode = 'SM MR '
565     CALL DIAGNOSTICS_ADD2LIST( diagNum,
566     I diagName, diagCode, diagUnits, diagTitle, myThid )
567 jmc 1.1
568 molod 1.17 diagName = 'RSURF '
569 jmc 1.26 c diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
570     diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
571 molod 1.17 diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
572     diagCode = 'SM M1 '
573     CALL DIAGNOSTICS_ADD2LIST( diagNum,
574     I diagName, diagCode, diagUnits, diagTitle, myThid )
575    
576 molod 1.24 diagName = 'TOTUTEND'
577     diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
578     diagUnits = 'm/s^2 '
579     WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR '
580     CALL DIAGNOSTICS_ADD2LIST( diagNum,
581     I diagName, diagCode, diagUnits, diagTitle, myThid )
582    
583     diagName = 'TOTVTEND'
584     diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
585     diagUnits = 'm/s^2 '
586     WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR '
587     CALL DIAGNOSTICS_ADD2LIST( diagNum,
588     I diagName, diagCode, diagUnits, diagTitle, myThid )
589    
590 molod 1.25 diagName = 'TOTTTEND'
591 jmc 1.26 diagTitle = 'Tendency of Potential Temperature'
592     diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
593 molod 1.24 diagCode = 'SM MR '
594     CALL DIAGNOSTICS_ADD2LIST( diagNum,
595     I diagName, diagCode, diagUnits, diagTitle, myThid )
596    
597     diagName = 'TOTSTEND'
598 jmc 1.26 diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
599     diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
600 molod 1.24 diagCode = 'SM MR '
601     CALL DIAGNOSTICS_ADD2LIST( diagNum,
602     I diagName, diagCode, diagUnits, diagTitle, myThid )
603    
604 molod 1.17
605 jmc 1.1 RETURN
606     END

  ViewVC Help
Powered by ViewVC 1.1.22