/[MITgcm]/MITgcm/pkg/gmredi/gmredi_diagnostics_init.F
ViewVC logotype

Contents of /MITgcm/pkg/gmredi/gmredi_diagnostics_init.F

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


Revision 1.22 - (show annotations) (download)
Thu Sep 24 21:31:35 2015 UTC (8 years, 7 months ago) by dfer
Branch: MAIN
CVS Tags: checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, HEAD
Changes since 1.21: +38 -38 lines
Small fixes

1 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_diagnostics_init.F,v 1.21 2014/11/24 10:12:15 dfer Exp $
2 C $Name: $
3
4 #include "GMREDI_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: GMREDI_DIAGNOSTICS_INIT
8 C !INTERFACE:
9 SUBROUTINE GMREDI_DIAGNOSTICS_INIT( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE GMREDI_DIAGNOSTICS_INIT
14 C | o Routine to initialize list of all available diagnostics
15 C | for GM/Redi package
16 C *==========================================================*
17 C \ev
18 C !USES:
19 IMPLICIT NONE
20
21 C === Global variables ===
22 #include "EEPARAMS.h"
23 c #include "SIZE.h"
24 c #include "PARAMS.h"
25 c #include "GMREDI.h"
26
27 C !INPUT/OUTPUT PARAMETERS:
28 C === Routine arguments ===
29 C myThid :: my Thread Id number
30 INTEGER myThid
31 CEOP
32
33 #ifdef ALLOW_DIAGNOSTICS
34 C !LOCAL VARIABLES:
35 C === Local variables ===
36 C diagNum :: diagnostics number in the (long) list of available diag.
37 C diagMate :: diag. mate number in the (long) list of available diag.
38 C diagName :: local short name (8c) of a diagnostics
39 C diagCode :: local parser field with characteristics of the diagnostics
40 C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
41 C diagUnits :: local string (16c): physical units of a diagnostic field
42 C diagTitle :: local string (80c): description of field in diagnostic
43 INTEGER diagNum
44 INTEGER diagMate
45 CHARACTER*8 diagName
46 CHARACTER*16 diagCode
47 CHARACTER*16 diagUnits
48 CHARACTER*(80) diagTitle
49
50 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
51
52 c IF ( useDiagnotics ) THEN
53
54 diagName = 'GM_VisbK'
55 diagTitle =
56 & 'Mixing coefficient from Visbeck etal parameterization'
57 diagUnits = 'm^2/s '
58 diagCode = 'SM P M1 '
59 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
60 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
61
62 diagName = 'GM_hTrsL'
63 diagTitle = 'Base depth (>0) of the Transition Layer'
64 diagUnits = 'm '
65 diagCode = 'SM P M1 '
66 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
67 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
68
69 diagName = 'GM_baseS'
70 diagTitle = 'Slope at the base of the Transition Layer'
71 diagUnits = '1 '
72 diagCode = 'SM P M1 '
73 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
74 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
75
76 diagName = 'GM_rLamb'
77 diagTitle =
78 & 'Slope vertical gradient at Trans. Layer Base (=recip.Lambda)'
79 diagUnits = '1/m '
80 diagCode = 'SM P M1 '
81 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
82 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
83
84 #ifndef GM_EXCLUDE_SUBMESO
85 diagName = 'SubMesLf'
86 diagTitle = 'Sub-Meso horiz. Length Scale (Lf)'
87 diagUnits = 'm '
88 diagCode = 'SM P M1 '
89 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
90 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
91
92 diagName = 'SubMpsiX'
93 diagTitle =
94 & 'Sub-Meso transp.stream-funct. magnitude (Psi0): U component'
95 diagUnits = 'm^2/s '
96 diagCode = 'UU M1 '
97 diagMate = diagNum + 2
98 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
99 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
100
101 diagName = 'SubMpsiY'
102 diagTitle =
103 & 'Sub-Meso transp.stream-funct. magnitude (Psi0): V component'
104 diagUnits = 'm^2/s '
105 diagCode = 'VV M1 '
106 diagMate = diagNum
107 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
108 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
109 #endif
110
111 diagName = 'GM_Kux '
112 diagTitle = 'K_11 element (U.point, X.dir) of GM-Redi tensor'
113 diagUnits = 'm^2/s '
114 diagCode = 'UU P MR '
115 diagMate = diagNum + 2
116 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
117 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
118
119 diagName = 'GM_Kvy '
120 diagTitle = 'K_22 element (V.point, Y.dir) of GM-Redi tensor'
121 diagUnits = 'm^2/s '
122 diagCode = 'VV P MR '
123 diagMate = diagNum
124 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
125 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
126
127 diagName = 'GM_Kuz '
128 diagTitle = 'K_13 element (U.point, Z.dir) of GM-Redi tensor'
129 diagUnits = 'm^2/s '
130 diagCode = 'UU MR '
131 diagMate = diagNum + 2
132 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
133 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
134
135 diagName = 'GM_Kvz '
136 diagTitle = 'K_23 element (V.point, Z.dir) of GM-Redi tensor'
137 diagUnits = 'm^2/s '
138 diagCode = 'VV MR '
139 diagMate = diagNum
140 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
141 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
142
143 diagName = 'GM_Kwx '
144 diagTitle = 'K_31 element (W.point, X.dir) of GM-Redi tensor'
145 diagUnits = 'm^2/s '
146 diagCode = 'UM LR '
147 diagMate = diagNum + 2
148 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
149 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
150
151 diagName = 'GM_Kwy '
152 diagTitle = 'K_32 element (W.point, Y.dir) of GM-Redi tensor'
153 diagUnits = 'm^2/s '
154 diagCode = 'VM LR '
155 diagMate = diagNum
156 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
157 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
158
159 diagName = 'GM_Kwz '
160 diagTitle = 'K_33 element (W.point, Z.dir) of GM-Redi tensor'
161 diagUnits = 'm^2/s '
162 diagCode = 'WM P LR '
163 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
164 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
165
166 diagName = 'GM_PsiX '
167 diagTitle = 'GM Bolus transport stream-function : U component'
168 diagUnits = 'm^2/s '
169 diagCode = 'UU LR '
170 diagMate = diagNum + 2
171 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
172 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
173
174 diagName = 'GM_PsiY '
175 diagTitle = 'GM Bolus transport stream-function : V component'
176 diagUnits = 'm^2/s '
177 diagCode = 'VV LR '
178 diagMate = diagNum
179 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
180 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
181
182 diagName = 'GM_KuzTz'
183 diagTitle = 'Redi Off-diagonal Temperature flux: X component'
184 diagUnits = 'degC.m^3/s '
185 diagCode = 'UU MR '
186 diagMate = diagNum + 2
187 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
188 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
189
190 diagName = 'GM_KvzTz'
191 diagTitle = 'Redi Off-diagonal Temperature flux: Y component'
192 diagUnits = 'degC.m^3/s '
193 diagCode = 'VV MR '
194 diagMate = diagNum
195 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
196 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
197
198 diagName = 'GM_KwzTz'
199 diagTitle = 'Redi main-diagonal vertical Temperature flux'
200 diagUnits = 'degC.m^3/s '
201 diagCode = 'WM LR '
202 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
203 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
204
205 diagName = 'GM_ubT '
206 diagTitle = 'Zonal Mass-Weight Bolus Transp of Pot Temp'
207 diagUnits = 'degC.m^3/s '
208 diagCode = 'UUr MR '
209 diagMate = diagNum + 2
210 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
211 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
212
213 diagName = 'GM_vbT '
214 diagTitle = 'Meridional Mass-Weight Bolus Transp of Pot Temp'
215 diagUnits = 'degC.m^3/s '
216 diagCode = 'VVr MR '
217 diagMate = diagNum
218 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
219 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
220
221 #ifdef ALLOW_EDDYPSI
222 diagName = 'GMEdTauX'
223 diagTitle = 'eddy-induced stress X-comp.'
224 diagUnits = 'N/m^2 '
225 diagCode = 'UM LR '
226 diagMate = diagNum + 2
227 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
229
230 diagName = 'GMEdTauY'
231 diagTitle = 'eddy-induced stress Y-comp.'
232 diagUnits = 'N/m^2 '
233 diagCode = 'VM LR '
234 diagMate = diagNum
235 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
236 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
237 #endif
238
239 #ifdef GM_BOLUS_BVP
240 diagName = 'GM_BVPcW'
241 diagTitle = 'WKB wave speed (at Western edge location)'
242 diagUnits = 'm/s '
243 diagCode = 'SU P M1 '
244 diagMate = diagNum + 2
245 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
246 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
247
248 diagName = 'GM_BVPcS'
249 diagTitle = 'WKB wave speed (at Southern edge location)'
250 diagUnits = 'm/s '
251 diagCode = 'SV P M1 '
252 diagMate = diagNum
253 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
254 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
255 #endif
256
257 #ifdef GM_K3D
258 diagName = 'GM_K3D '
259 diagTitle = '3D diffusivity'
260 diagUnits = 'm**2/s '
261 diagCode = 'SM P MR '
262 diagMate = diagNum
263 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
264 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
265
266 diagName = 'GM_KPV '
267 diagTitle = '3D eddy transport diffusivity'
268 diagUnits = 'm**2/s '
269 diagCode = 'SM P MR '
270 diagMate = diagNum
271 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
272 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
273
274 diagName = 'GM_K3D_T'
275 diagTitle = '3D Redi diffusivity'
276 diagUnits = 'm**2/s '
277 diagCode = 'SM P MR '
278 diagMate = diagNum
279 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
280 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
281
282 #ifdef HAVE_LAPACK
283 diagName = 'GM_A3D'
284 diagTitle = '3D lower diagona'
285 diagUnits = '1/m**2 '
286 diagCode = 'SM MR '
287 diagMate = diagNum
288 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
289 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
290
291 diagName = 'GM_B3D'
292 diagTitle = '3D lower diagona'
293 diagUnits = '1/m**2 '
294 diagCode = 'SM MR '
295 diagMate = diagNum
296 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
297 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
298
299 diagName = 'GM_C3D'
300 diagTitle = '3D lower diagona'
301 diagUnits = '1/m**2 '
302 diagCode = 'SM MR '
303 diagMate = diagNum
304 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
305 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
306 #endif
307
308 diagName = 'GM_MODEC'
309 diagTitle = 'Eigen vector for the first barcolinic mode'
310 diagUnits = 'dimensionless '
311 diagCode = 'SM MR '
312 diagMate = diagNum
313 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
314 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
315
316 diagName = 'GM_URMS'
317 diagTitle = 'rms Eddy Velocity'
318 diagUnits = 'm/s '
319 diagCode = 'SM P MR '
320 diagMate = diagNum
321 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
322 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
323
324 diagName = 'GM_UMC '
325 diagTitle = 'ubar-c'
326 diagUnits = 'm/s '
327 diagCode = 'SM MR '
328 diagMate = diagNum
329 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
330 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
331
332 diagName = 'GM_SFLYR'
333 diagTitle = 'mixed layer depth'
334 diagUnits = 'm '
335 diagCode = 'SM P MR '
336 diagMate = diagNum
337 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
338 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
339
340 diagName = 'GM_USTAR'
341 diagTitle = 'u^*'
342 diagUnits = 'm/s '
343 diagCode = 'UU MR '
344 diagMate = diagNum
345 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
346 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
347
348 diagName = 'GM_VSTAR'
349 diagTitle = 'v^*'
350 diagUnits = 'm/s '
351 diagCode = 'VV MR '
352 diagMate = diagNum
353 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
354 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
355
356 diagName = 'GM_RDEF'
357 diagTitle = 'Deformation Radius'
358 diagUnits = 'm '
359 diagCode = 'SM P M1 '
360 diagMate = diagNum
361 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
362 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
363
364 diagName = 'GM_RURMS'
365 diagTitle =
366 I 'Mixing length for calculating u_rms: MIN(Rdef,GM_K3D_Rmax)'
367 diagUnits = 'm '
368 diagCode = 'SM P M1 '
369 diagMate = diagNum
370 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
371 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
372
373 diagName = 'GM_RRHNS'
374 diagTitle = 'Rhines scale: sqrt(urms/beta)'
375 diagUnits = 'm '
376 diagCode = 'SM P M1 '
377 diagMate = diagNum
378 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
379 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
380
381 diagName = 'GM_RMIX'
382 diagTitle = 'Unmodulated Mixing Length: MIN(Rdef, RRhines)'
383 diagUnits = 'm '
384 diagCode = 'SM P M1 '
385 diagMate = diagNum
386 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
387 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
388
389 diagName = 'GM_SUPP'
390 diagTitle = 'Suppression Factor for K3D'
391 diagUnits = 'none '
392 diagCode = 'SM P MR '
393 diagMate = diagNum
394 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
395 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
396
397 diagName = 'GM_dqdx '
398 diagTitle = 'dq/dx'
399 diagUnits = '1/(m*s) '
400 diagCode = 'UU MR '
401 diagMate = diagNum
402 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
403 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
404
405 diagName = 'GM_dqdy '
406 diagTitle = 'dq/dy'
407 diagUnits = '1/(m*s) '
408 diagCode = 'VV MR '
409 diagMate = diagNum
410 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
411 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
412
413 diagName = 'GM_Kdqdx'
414 diagTitle = 'K*dq/dx'
415 diagUnits = 'm/s*2 '
416 diagCode = 'UU MR '
417 diagMate = diagNum
418 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
419 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
420
421 diagName = 'GM_Kdqdy'
422 diagTitle = 'K*dq/dy'
423 diagUnits = 'm/s**2 '
424 diagCode = 'VV MR '
425 diagMate = diagNum
426 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
427 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
428
429 diagName = 'GM_Xix '
430 diagTitle = '-k dq/dx expansion'
431 diagUnits = 'm/s**2 '
432 diagCode = 'UU MR '
433 diagMate = diagNum
434 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
435 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
436
437 diagName = 'GM_Xiy '
438 diagTitle = '-k dq/dy expansion'
439 diagUnits = 'm/s**2 '
440 diagCode = 'VV MR '
441 diagMate = diagNum
442 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
443 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
444
445 diagName = 'GM_Sx '
446 diagTitle = 'Zonal isopycnal slope'
447 diagUnits = 'none '
448 diagCode = 'UU LR '
449 diagMate = diagNum
450 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
451 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
452
453 diagName = 'GM_Sy '
454 diagTitle = 'Meridional isopycnal slope'
455 diagUnits = 'none '
456 diagCode = 'VV LR '
457 diagMate = diagNum
458 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
459 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
460
461 diagName = 'GM_TFLXX'
462 diagTitle = 'Zonal thickness flux'
463 diagUnits = '1/(m*s) '
464 diagCode = 'UU MR '
465 diagMate = diagNum
466 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
467 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
468
469 diagName = 'GM_TFLXY'
470 diagTitle = 'meridional thickness flux'
471 diagUnits = '1/(m*s) '
472 diagCode = 'VV MR '
473 diagMate = diagNum
474 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
475 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
476
477 diagName = 'GM_C'
478 diagTitle = 'Doppler shifted long Rossby wave speed'
479 diagUnits = 'm/s '
480 diagCode = 'SM M1 '
481 diagMate = diagNum
482 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
483 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
484
485 diagName = 'GM_UBARO'
486 diagTitle = 'Barotropic velocity'
487 diagUnits = 'm/s '
488 diagCode = 'SM M1 '
489 diagMate = diagNum
490 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
491 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
492
493 diagName = 'GM_EADY '
494 diagTitle = 'Eady Growth rate'
495 diagUnits = '1/s '
496 diagCode = 'SM M1 '
497 diagMate = diagNum
498 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
499 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
500
501 diagName = 'GM_UBAR '
502 diagTitle = 'Mean zonal velocity'
503 diagUnits = 'm/s '
504 diagCode = 'SM MR '
505 diagMate = diagNum
506 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
507 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
508
509 diagName = 'GM_M4 '
510 diagTitle = 'M**4'
511 diagUnits = '1/s**4 '
512 diagCode = 'SM MR '
513 diagMate = diagNum
514 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
515 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
516
517 diagName = 'GM_N2 '
518 diagTitle = 'N**2'
519 diagUnits = '1/s**2 '
520 diagCode = 'SM MR '
521 diagMate = diagNum
522 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
523 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
524
525 diagName = 'GM_M4_N2'
526 diagTitle = 'M**4/N**2'
527 diagUnits = '1/s**2 '
528 diagCode = 'SM MR '
529 diagMate = diagNum
530 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
531 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
532
533 diagName = 'GM_SLOPE'
534 diagTitle = 'SLOPE'
535 diagUnits = 'm/m '
536 diagCode = 'SM MR '
537 diagMate = diagNum
538 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
539 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
540
541 diagName = 'GM_RENRM'
542 diagTitle = 'renormalisation factor'
543 diagUnits = 'none '
544 diagCode = 'SM M1 '
545 diagMate = diagNum
546 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
547 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
548
549 #ifdef GM_K3D_PASSIVE
550 diagName = 'GM_PSTAR'
551 diagTitle = 'Diagnosed bolus streamfunction: V component'
552 diagUnits = 'm**2/s '
553 diagCode = 'VV LR '
554 diagMate = diagNum
555 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
556 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
557 #endif
558 #endif
559
560 c ENDIF
561
562 #endif /* ALLOW_DIAGNOSTICS */
563
564 RETURN
565 END

  ViewVC Help
Powered by ViewVC 1.1.22