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

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

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


Revision 1.19 - (hide annotations) (download)
Fri Mar 28 04:22:11 2014 UTC (10 years, 2 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64w, checkpoint64v
Changes since 1.18: +9 -1 lines
K3D: Added upper and lower bounds to renormalisation factor as namelist variables. Fixed small bug with interpolation to U & V grid. Also added extra documentation and a diagnostic.

1 m_bates 1.19 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_diagnostics_init.F,v 1.18 2014/01/01 23:20:48 m_bates Exp $
2 jmc 1.1 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 jmc 1.7 C diagMate :: diag. mate number in the (long) list of available diag.
38 jmc 1.1 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 jmc 1.7 INTEGER diagMate
45 jmc 1.1 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 jmc 1.7 diagTitle =
56 jmc 1.1 & 'Mixing coefficient from Visbeck etal parameterization'
57     diagUnits = 'm^2/s '
58     diagCode = 'SM P M1 '
59 jmc 1.7 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
60     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
61 jmc 1.1
62 jmc 1.5 diagName = 'GM_hTrsL'
63     diagTitle = 'Base depth (>0) of the Transition Layer'
64     diagUnits = 'm '
65     diagCode = 'SM P M1 '
66 jmc 1.7 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
67     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
68 jmc 1.5
69     diagName = 'GM_baseS'
70     diagTitle = 'Slope at the base of the Transition Layer'
71     diagUnits = '1 '
72     diagCode = 'SM P M1 '
73 jmc 1.7 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
74     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
75 jmc 1.5
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 jmc 1.7 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
82     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
83 jmc 1.5
84 jmc 1.9 #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 jmc 1.1 diagName = 'GM_Kux '
112     diagTitle = 'K_11 element (U.point, X.dir) of GM-Redi tensor'
113     diagUnits = 'm^2/s '
114 jmc 1.7 diagCode = 'UU P MR '
115     diagMate = diagNum + 2
116     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
117     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
118 jmc 1.1
119     diagName = 'GM_Kvy '
120     diagTitle = 'K_22 element (V.point, Y.dir) of GM-Redi tensor'
121     diagUnits = 'm^2/s '
122 jmc 1.7 diagCode = 'VV P MR '
123     diagMate = diagNum
124     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
125     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
126 jmc 1.1
127     diagName = 'GM_Kuz '
128     diagTitle = 'K_13 element (U.point, Z.dir) of GM-Redi tensor'
129     diagUnits = 'm^2/s '
130 jmc 1.7 diagCode = 'UU MR '
131     diagMate = diagNum + 2
132     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
133     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
134 jmc 1.1
135     diagName = 'GM_Kvz '
136     diagTitle = 'K_23 element (V.point, Z.dir) of GM-Redi tensor'
137     diagUnits = 'm^2/s '
138 jmc 1.7 diagCode = 'VV MR '
139     diagMate = diagNum
140     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
141     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
142 jmc 1.1
143     diagName = 'GM_Kwx '
144     diagTitle = 'K_31 element (W.point, X.dir) of GM-Redi tensor'
145     diagUnits = 'm^2/s '
146 jmc 1.7 diagCode = 'UM LR '
147     diagMate = diagNum + 2
148     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
149     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
150 jmc 1.1
151     diagName = 'GM_Kwy '
152     diagTitle = 'K_32 element (W.point, Y.dir) of GM-Redi tensor'
153     diagUnits = 'm^2/s '
154 jmc 1.7 diagCode = 'VM LR '
155     diagMate = diagNum
156     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
157     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
158 jmc 1.1
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 jmc 1.7 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
164     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
165 jmc 1.1
166     diagName = 'GM_PsiX '
167 jmc 1.9 diagTitle = 'GM Bolus transport stream-function : U component'
168 jmc 1.1 diagUnits = 'm^2/s '
169 jmc 1.7 diagCode = 'UU LR '
170     diagMate = diagNum + 2
171     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
172     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
173 jmc 1.1
174     diagName = 'GM_PsiY '
175 jmc 1.9 diagTitle = 'GM Bolus transport stream-function : V component'
176 jmc 1.1 diagUnits = 'm^2/s '
177 jmc 1.7 diagCode = 'VV LR '
178     diagMate = diagNum
179     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
180     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
181 jmc 1.1
182 jmc 1.2 diagName = 'GM_KuzTz'
183 dfer 1.4 diagTitle = 'Redi Off-diagonal Temperature flux: X component'
184 jmc 1.2 diagUnits = 'degC.m^3/s '
185 jmc 1.7 diagCode = 'UU MR '
186     diagMate = diagNum + 2
187     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
188     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
189 jmc 1.2
190     diagName = 'GM_KvzTz'
191 dfer 1.4 diagTitle = 'Redi Off-diagonal Temperature flux: Y component'
192 jmc 1.2 diagUnits = 'degC.m^3/s '
193 jmc 1.7 diagCode = 'VV MR '
194     diagMate = diagNum
195     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
196     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
197 jmc 1.2
198 jmc 1.10 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 dfer 1.4 diagName = 'GM_ubT '
206     diagTitle = 'Zonal Mass-Weight Bolus Transp of Pot Temp'
207     diagUnits = 'degC.m^3/s '
208 jmc 1.7 diagCode = 'UUr MR '
209     diagMate = diagNum + 2
210     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
211     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
212 dfer 1.4
213     diagName = 'GM_vbT '
214     diagTitle = 'Meridional Mass-Weight Bolus Transp of Pot Temp'
215     diagUnits = 'degC.m^3/s '
216 jmc 1.7 diagCode = 'VVr MR '
217     diagMate = diagNum
218     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
219     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
220 dfer 1.4
221 gforget 1.6 #ifdef ALLOW_EDDYPSI
222 heimbach 1.3 diagName = 'GMEdTauX'
223 m_bates 1.16 diagTitle = 'eddy-induced stress X-comp.'
224 heimbach 1.3 diagUnits = 'N/m^2 '
225 jmc 1.7 diagCode = 'UM LR '
226     diagMate = diagNum + 2
227     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
229 heimbach 1.3
230     diagName = 'GMEdTauY'
231 m_bates 1.16 diagTitle = 'eddy-induced stress Y-comp.'
232 heimbach 1.3 diagUnits = 'N/m^2 '
233 jmc 1.7 diagCode = 'VM LR '
234     diagMate = diagNum
235     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
236     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
237 heimbach 1.3 #endif
238    
239 jmc 1.8 #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 m_bates 1.11 #ifdef GM_K3D
258 m_bates 1.18 diagName = 'GM_K3D '
259 m_bates 1.11 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 m_bates 1.18 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 m_bates 1.12 #ifdef use_lapack
283 m_bates 1.11 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 m_bates 1.12 #endif
307 m_bates 1.11
308 m_bates 1.12 diagName = 'GM_MODEC'
309 m_bates 1.11 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 m_bates 1.14 diagName = 'GM_RURMS'
365     diagTitle = 'Mixing length for calculating u_rms'
366     diagUnits = 'm '
367     diagCode = 'SM P M1 '
368     diagMate = diagNum
369     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
370     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
371    
372     diagName = 'GM_RRHNS'
373 m_bates 1.11 diagTitle = 'Rhines Radius'
374     diagUnits = 'm '
375     diagCode = 'SM P M1 '
376     diagMate = diagNum
377     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
378     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
379    
380     diagName = 'GM_RMIX'
381     diagTitle = 'Unmodulated Mixing Length'
382     diagUnits = 'm '
383     diagCode = 'SM P M1 '
384     diagMate = diagNum
385     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
386     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
387    
388     diagName = 'GM_SUPP'
389     diagTitle = 'Suppression Factor for K3D'
390     diagUnits = 'none '
391     diagCode = 'SM P MR '
392     diagMate = diagNum
393     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
394     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
395    
396     diagName = 'GM_dqdx '
397     diagTitle = 'dq/dx'
398     diagUnits = '1/(m*s) '
399     diagCode = 'UU MR '
400     diagMate = diagNum
401     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
402     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
403    
404     diagName = 'GM_dqdy '
405     diagTitle = 'dq/dy'
406     diagUnits = '1/(m*s) '
407     diagCode = 'VV MR '
408     diagMate = diagNum
409     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
410     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
411    
412 m_bates 1.17 diagName = 'GM_Kdqdx'
413     diagTitle = 'K*dq/dx'
414     diagUnits = 'm/s*2 '
415     diagCode = 'UU MR '
416     diagMate = diagNum
417     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
418     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
419    
420     diagName = 'GM_Kdqdy'
421     diagTitle = 'K*dq/dy'
422     diagUnits = 'm/s**2 '
423     diagCode = 'VV MR '
424     diagMate = diagNum
425     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
426     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
427    
428 m_bates 1.11 diagName = 'GM_Xix '
429     diagTitle = '-k dq/dx expansion'
430 m_bates 1.13 diagUnits = 'm/s**2 '
431 m_bates 1.11 diagCode = 'UU MR '
432     diagMate = diagNum
433     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
434     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
435    
436     diagName = 'GM_Xiy '
437     diagTitle = '-k dq/dy expansion'
438 m_bates 1.13 diagUnits = 'm/s**2 '
439 m_bates 1.11 diagCode = 'VV MR '
440     diagMate = diagNum
441     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
442     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
443    
444     diagName = 'GM_Sx '
445     diagTitle = 'Zonal isopycnal slope'
446     diagUnits = 'none '
447     diagCode = 'UU LR '
448     diagMate = diagNum
449     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
450     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
451    
452     diagName = 'GM_Sy '
453     diagTitle = 'Meridional isopycnal slope'
454     diagUnits = 'none '
455     diagCode = 'VV LR '
456     diagMate = diagNum
457     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
458     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
459    
460     diagName = 'GM_TFLXX'
461     diagTitle = 'Zonal thickness flux'
462     diagUnits = '1/(m*s) '
463     diagCode = 'UU MR '
464     diagMate = diagNum
465     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
466     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
467    
468     diagName = 'GM_TFLXY'
469     diagTitle = 'meridional thickness flux'
470     diagUnits = '1/(m*s) '
471     diagCode = 'VV MR '
472     diagMate = diagNum
473     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
474     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
475    
476     diagName = 'GM_C'
477     diagTitle = 'Doppler shifted long Rossby wave speed'
478     diagUnits = 'm/s '
479     diagCode = 'SM M1 '
480     diagMate = diagNum
481     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
482     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
483    
484     diagName = 'GM_UBARO'
485     diagTitle = 'Barotropic velocity'
486     diagUnits = 'm/s '
487     diagCode = 'SM M1 '
488     diagMate = diagNum
489     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
490     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
491    
492     diagName = 'GM_EADY '
493     diagTitle = 'Eady Growth rate'
494     diagUnits = '1/s '
495     diagCode = 'SM M1 '
496     diagMate = diagNum
497     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
498     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
499    
500     diagName = 'GM_UBAR '
501     diagTitle = 'Mean zonal velocity'
502     diagUnits = 'm/s '
503     diagCode = 'SM MR '
504     diagMate = diagNum
505     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
506     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
507    
508 m_bates 1.15 diagName = 'GM_M4 '
509     diagTitle = 'M**4'
510     diagUnits = '1/s**4 '
511     diagCode = 'SM MR '
512     diagMate = diagNum
513     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
514     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
515    
516     diagName = 'GM_N2 '
517     diagTitle = 'N**2'
518     diagUnits = '1/s**2 '
519     diagCode = 'SM MR '
520     diagMate = diagNum
521     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
522     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
523    
524 m_bates 1.17 diagName = 'GM_M4_N2'
525     diagTitle = 'M**4/N**2'
526     diagUnits = '1/s**2 '
527     diagCode = 'SM MR '
528     diagMate = diagNum
529     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
530     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
531    
532 m_bates 1.15 diagName = 'GM_SLOPE'
533     diagTitle = 'SLOPE'
534     diagUnits = 'm/m '
535     diagCode = 'SM MR '
536     diagMate = diagNum
537     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
538     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
539    
540 m_bates 1.11 diagName = 'GM_KDEF '
541     diagTitle = 'Rossby wavenumber'
542     diagUnits = '1/m '
543     diagCode = 'SM M1 '
544     diagMate = diagNum
545     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
546     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
547    
548 m_bates 1.19 diagName = 'GM_RENRM'
549     diagTitle = 'renormalisation factor'
550     diagUnits = 'none '
551     diagCode = 'SM M1 '
552     diagMate = diagNum
553     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
554     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
555    
556 m_bates 1.11 #endif
557    
558 jmc 1.1 c ENDIF
559    
560     #endif /* ALLOW_DIAGNOSTICS */
561    
562     RETURN
563     END

  ViewVC Help
Powered by ViewVC 1.1.22