/[MITgcm]/MITgcm/pkg/mom_common/mom_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/mom_common/mom_diagnostics_init.F

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


Revision 1.10 - (hide annotations) (download)
Tue Apr 22 22:03:17 2008 UTC (16 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint60, checkpoint61, checkpoint62, checkpoint59q, checkpoint59r, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.9: +321 -285 lines
call DIAGNOSTICS_ADDTOLIST (= new version) to replace DIAGNOSTICS_ADD2LIST
 call and add diagnostics for vertical momentum tendencies.

1 jmc 1.10 C $Header: /u/gcmpack/MITgcm/pkg/mom_common/mom_diagnostics_init.F,v 1.9 2006/07/13 21:30:12 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "MOM_COMMON_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: MOM_DIAGNOSTICS_INIT
8     C !INTERFACE:
9     SUBROUTINE MOM_DIAGNOSTICS_INIT( myThid )
10     C !DESCRIPTION:
11     C Routine to initialize Momentum diagnostics
12     C (from mom_common, mom_fluxform & mom_vecinv packages)
13    
14     C !USES:
15     IMPLICIT NONE
16     C === Global variables ===
17     #include "SIZE.h"
18     #include "EEPARAMS.h"
19     #include "PARAMS.h"
20    
21     C !INPUT/OUTPUT PARAMETERS:
22     C === Routine arguments ===
23     C myThid - Number of this instance of GAD_INIT
24     INTEGER myThid
25     CEOP
26    
27     #ifdef ALLOW_DIAGNOSTICS
28     #ifdef ALLOW_MOM_COMMON
29     C !LOCAL VARIABLES:
30     C === Local variables ===
31     C msgBuf - Informational/error meesage buffer
32     c CHARACTER*(MAX_LEN_MBUF) msgBuf
33    
34     INTEGER diagNum
35 jmc 1.10 INTEGER diagMate
36 jmc 1.1 CHARACTER*8 diagName
37     CHARACTER*16 diagCode
38     CHARACTER*16 diagUnits
39     CHARACTER*(80) diagTitle
40    
41     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
42    
43 jmc 1.3 C--- Add diagnostics to the (long) list of available diagnostics:
44    
45     C-- add diagnostics for viscosity coefficients:
46 baylor 1.5 diagName = 'VISCAHZ '
47     diagTitle = 'Harmonic Visc Coefficient (m2/s) (Zeta Pt)'
48     diagUnits = 'm^2/s '
49     diagCode = 'SZ MR '
50 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
51     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
52 baylor 1.5
53     diagName = 'VISCA4Z '
54     diagTitle = 'Biharmonic Visc Coefficient (m4/s) (Zeta Pt)'
55     diagUnits = 'm^4/s '
56     diagCode = 'SZ MR '
57 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
58     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
59 baylor 1.5
60 baylor 1.4 diagName = 'VISCAHD '
61     diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (Div Pt)'
62 jmc 1.3 diagUnits = 'm^2/s '
63     diagCode = 'SM MR '
64 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
65     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
66 jmc 1.3
67 baylor 1.4 diagName = 'VISCA4D '
68     diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (Div Pt)'
69 jmc 1.3 diagUnits = 'm^4/s '
70     diagCode = 'SM MR '
71 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
72     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
73 baylor 1.4
74 baylor 1.8 diagName = 'VISCAHW '
75     diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (W Pt)'
76     diagUnits = 'm^2/s '
77     diagCode = 'WM LR '
78 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
79     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
80 baylor 1.8
81     diagName = 'VISCA4W '
82     diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (W Pt)'
83     diagUnits = 'm^4/s '
84     diagCode = 'WM LR '
85 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
86     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
87 baylor 1.8
88 baylor 1.5 diagName = 'VAHZMAX '
89     diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Zeta Pt)'
90     diagUnits = 'm^2/s '
91     diagCode = 'SZ MR '
92 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
93     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
94 baylor 1.5
95     diagName = 'VA4ZMAX '
96     diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Zeta Pt)'
97     diagUnits = 'm^4/s '
98     diagCode = 'SZ MR '
99 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
100     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
101 baylor 1.5
102     diagName = 'VAHDMAX '
103     diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Div Pt)'
104     diagUnits = 'm^2/s '
105     diagCode = 'SM MR '
106 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
107     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
108 baylor 1.5
109     diagName = 'VA4DMAX '
110     diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Div Pt)'
111     diagUnits = 'm^4/s '
112     diagCode = 'SM MR '
113 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
114     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
115 baylor 1.5
116     diagName = 'VAHZMIN '
117     diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Zeta Pt)'
118     diagUnits = 'm^2/s '
119     diagCode = 'SZ MR '
120 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
121     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
122 baylor 1.5
123     diagName = 'VA4ZMIN '
124     diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Zeta Pt)'
125     diagUnits = 'm^4/s '
126     diagCode = 'SZ MR '
127 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
128     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
129 baylor 1.5
130     diagName = 'VAHDMIN '
131     diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Div Pt)'
132     diagUnits = 'm^2/s '
133     diagCode = 'SM MR '
134 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
135     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
136 baylor 1.5
137     diagName = 'VA4DMIN '
138     diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Div Pt)'
139     diagUnits = 'm^4/s '
140     diagCode = 'SM MR '
141 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
142     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
143 baylor 1.5
144     diagName = 'VAHZLTH '
145     diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Zeta Pt)'
146     diagUnits = 'm^2/s '
147     diagCode = 'SZ MR '
148 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
149     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
150 baylor 1.5
151     diagName = 'VA4ZLTH '
152     diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Zeta Pt)'
153     diagUnits = 'm^4/s '
154     diagCode = 'SZ MR '
155 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
156     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
157 baylor 1.5
158     diagName = 'VAHDLTH '
159     diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Div Pt)'
160     diagUnits = 'm^2/s '
161     diagCode = 'SM MR '
162 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
163     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
164 baylor 1.5
165     diagName = 'VA4DLTH '
166     diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Div Pt)'
167     diagUnits = 'm^4/s '
168     diagCode = 'SM MR '
169 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
170     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
171 baylor 1.5
172     diagName = 'VAHZLTHD'
173     diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Zeta Pt)'
174     diagUnits = 'm^2/s '
175     diagCode = 'SZ MR '
176 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
177     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
178 baylor 1.5
179     diagName = 'VA4ZLTHD'
180     diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Zeta Pt)'
181     diagUnits = 'm^4/s '
182     diagCode = 'SZ MR '
183 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
184     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
185 baylor 1.5
186     diagName = 'VAHDLTHD'
187     diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Div Pt)'
188     diagUnits = 'm^2/s '
189     diagCode = 'SM MR '
190 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
191     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
192 baylor 1.5
193     diagName = 'VA4DLTHD'
194     diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Div Pt)'
195     diagUnits = 'm^4/s '
196     diagCode = 'SM MR '
197 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
198     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
199 baylor 1.5
200     diagName = 'VAHZSMAG'
201     diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Zeta Pt)'
202 baylor 1.4 diagUnits = 'm^2/s '
203     diagCode = 'SZ MR '
204 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
205     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
206 baylor 1.4
207 baylor 1.5 diagName = 'VA4ZSMAG'
208     diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Zeta Pt)'
209 baylor 1.4 diagUnits = 'm^4/s '
210     diagCode = 'SZ MR '
211 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
212     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
213 jmc 1.1
214 baylor 1.5 diagName = 'VAHDSMAG'
215     diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Div Pt)'
216     diagUnits = 'm^2/s '
217     diagCode = 'SM MR '
218 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
219     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
220 baylor 1.5
221     diagName = 'VA4DSMAG'
222     diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Div Pt)'
223     diagUnits = 'm^4/s '
224     diagCode = 'SM MR '
225 jmc 1.10 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
226     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
227 baylor 1.5
228    
229 jmc 1.10 C-- add diagnostics of simple derivative quantities:
230     diagName = 'momKE '
231     diagTitle = 'Kinetic Energy (in momentum Eq.)'
232     diagCode = 'SM MR '
233     diagUnits = 'm^2/s^2 '
234     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
235     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
236    
237     #ifdef ALLOW_MOM_VECINV
238     diagName = 'momHDiv '
239     diagTitle = 'Horizontal Divergence (in momentum Eq.)'
240     diagCode = 'SM MR '
241     diagUnits = 's^-1 '
242     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
243     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
244 baylor 1.5
245 jmc 1.10 diagName = 'momVort3'
246     diagTitle = '3rd component (vertical) of Vorticity'
247     diagCode = 'SZ MR '
248     diagUnits = 's^-1 '
249     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
250     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
251    
252     diagName = 'Strain '
253     diagTitle = 'Horizontal Strain of Horizontal Velocities'
254     diagCode = 'SZ MR '
255     diagUnits = 's^-1 '
256     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
257     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
258 jmc 1.1
259 jmc 1.10 diagName = 'Tension '
260     diagTitle = 'Horizontal Tension of Horizontal Velocities'
261     diagCode = 'SM MR '
262     diagUnits = 's^-1 '
263     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
264     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
265 jmc 1.1 #endif /* ALLOW_MOM_VECINV */
266    
267 jmc 1.10 C-- add diagnostics of tendencies from each terms
268     diagUnits = 'm/s^2 '
269 jmc 1.1
270 jmc 1.10 diagName = 'UBotDrag'
271     diagTitle = 'U momentum tendency from Bottom Drag'
272     diagCode = 'UU MR '
273     diagMate = diagNum + 2
274     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
275     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
276     diagName = 'VBotDrag'
277     diagTitle = 'V momentum tendency from Bottom Drag'
278     diagCode = 'VV MR '
279     diagMate = diagNum
280     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
281     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
282     diagName = 'USidDrag'
283     diagTitle = 'U momentum tendency from Side Drag'
284     diagCode = 'UU MR '
285     diagMate = diagNum + 2
286     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
287     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
288     diagName = 'VSidDrag'
289     diagTitle = 'V momentum tendency from Side Drag'
290     diagCode = 'VV MR '
291     diagMate = diagNum
292     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
293     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
294    
295     diagName = 'Um_Diss '
296     diagTitle = 'U momentum tendency from Dissipation'
297     diagCode = 'UU MR '
298     diagMate = diagNum + 2
299     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
300     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
301     diagName = 'Vm_Diss '
302     diagTitle = 'V momentum tendency from Dissipation'
303     diagCode = 'VV MR '
304     diagMate = diagNum
305     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
306     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
307    
308     diagName = 'Um_Advec'
309     diagTitle = 'U momentum tendency from Advection terms'
310     diagCode = 'UU MR '
311     diagMate = diagNum + 2
312     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
313     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
314     diagName = 'Vm_Advec'
315     diagTitle = 'V momentum tendency from Advection terms'
316     diagCode = 'VV MR '
317     diagMate = diagNum
318     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
319     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
320    
321     diagName = 'Um_Cori '
322     diagTitle = 'U momentum tendency from Coriolis term'
323     diagCode = 'UU MR '
324     diagMate = diagNum + 2
325     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
326     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
327     diagName = 'Vm_Cori '
328     diagTitle = 'V momentum tendency from Coriolis term'
329     diagCode = 'VV MR '
330     diagMate = diagNum
331     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
332     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
333    
334     diagName = 'Um_Ext '
335     diagTitle = 'U momentum tendency from external forcing'
336     diagCode = 'UU MR '
337     diagMate = diagNum + 2
338     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
339     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
340     diagName = 'Vm_Ext '
341     diagTitle = 'V momentum tendency from external forcing'
342     diagCode = 'VV MR '
343     diagMate = diagNum
344     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
345     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
346 baylor 1.6
347 jmc 1.1 #ifdef ALLOW_MOM_VECINV
348 jmc 1.10 diagName = 'Um_AdvZ3'
349     diagTitle = 'U momentum tendency from Vorticity Advection'
350     diagCode = 'UU MR '
351     diagMate = diagNum + 2
352     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
353     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
354     diagName = 'Vm_AdvZ3'
355     diagTitle = 'V momentum tendency from Vorticity Advection'
356     diagCode = 'VV MR '
357     diagMate = diagNum
358     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
359     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
360    
361     diagName = 'Um_AdvRe'
362     diagTitle = 'U momentum tendency from vertical Advection'
363     & //' (Explicit part)'
364     diagCode = 'UU MR '
365     diagMate = diagNum + 2
366     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
367     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
368     diagName = 'Vm_AdvRe'
369     diagTitle = 'V momentum tendency from vertical Advection'
370     & //' (Explicit part)'
371     diagCode = 'VV MR '
372     diagMate = diagNum
373     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
374     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
375 baylor 1.6
376 jmc 1.1 #endif /* ALLOW_MOM_VECINV */
377    
378 jmc 1.9 #ifdef ALLOW_NONHYDROSTATIC
379 jmc 1.10 C- vertical momentum tendencies
380     IF ( usingPCoords ) THEN
381     diagUnits = 'Pa/s^2 '
382     ELSE
383     diagUnits = 'm/s^2 '
384     ENDIF
385     diagName = 'Wm_Diss '
386     diagTitle = 'W momentum tendency from Dissipation'
387     diagCode = 'WM LR '
388     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
389     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
390    
391     diagName = 'Wm_Advec'
392     diagTitle = 'W momentum tendency from Advection terms'
393     diagCode = 'WM LR '
394     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
395     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
396    
397     diagName = 'WSidDrag'
398     diagTitle = 'Vertical momentum tendency from Side Drag'
399     diagCode = 'WM LR '
400     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
401     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
402 jmc 1.9 #endif /* ALLOW_NONHYDROSTATIC */
403    
404 jmc 1.10 C-- add diagnostics of advective & viscous flux :
405     IF ( usingPCoords ) THEN
406     diagUnits = 'Pa.m^3/s^2 '
407     ELSE
408     diagUnits = 'm^4/s^2 '
409     ENDIF
410 jmc 1.1
411     C- Advective flux of Zonal momentum:
412     #ifdef ALLOW_MOM_FLUXFORM
413 jmc 1.10 diagName = 'ADVx_Um '
414     diagTitle = 'Zonal Advective Flux of U momentum'
415     diagCode = 'UM MR '
416     diagMate = diagNum + 2
417     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
418     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
419     diagName = 'ADVy_Um '
420     diagTitle = 'Meridional Advective Flux of U momentum'
421     diagCode = 'VZ MR '
422     diagMate = diagNum
423     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
424     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
425     diagName = 'ADVrE_Um'
426     diagTitle = 'Vertical Advective Flux of U momentum'
427     & //' (Explicit part)'
428     diagCode = 'WU LR '
429     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
430     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
431 jmc 1.1 #endif /* ALLOW_MOM_FLUXFORM */
432 jmc 1.10 c diagName = 'ADVrI_Um'
433     c diagTitle = 'Vertical Advective Flux of U momentum'
434     c & //' (Implicit part)'
435     c diagCode = 'WU LR '
436     c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
437     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
438 jmc 1.1
439 jmc 1.10 C- Advective flux of Meridional momentum:
440 jmc 1.1 #ifdef ALLOW_MOM_FLUXFORM
441 jmc 1.10 diagName = 'ADVx_Vm '
442     diagTitle = 'Zonal Advective Flux of V momentum'
443     diagCode = 'UZ MR '
444     diagMate = diagNum + 2
445     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
446     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
447     diagName = 'ADVy_Vm '
448     diagTitle = 'Meridional Advective Flux of V momentum'
449     diagCode = 'VM MR '
450     diagMate = diagNum
451     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
452     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
453     diagName = 'ADVrE_Vm'
454     diagTitle = 'Vertical Advective Flux of V momentum'
455     & //' (Explicit part)'
456     diagCode = 'WV LR '
457     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
458     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
459 jmc 1.1 #endif /* ALLOW_MOM_FLUXFORM */
460 jmc 1.10 c diagName = 'ADVrI_Vm'
461     c diagTitle = 'Vertical Advective Flux of V momentum'
462     c & //' (Implicit part)'
463     c diagCode = 'WV LR '
464     c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
465     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
466 jmc 1.1
467 jmc 1.10 C- Viscous flux of Zonal momentum :
468 jmc 1.1 #ifdef ALLOW_MOM_FLUXFORM
469 jmc 1.10 diagName = 'VISCx_Um'
470     diagTitle = 'Zonal Viscous Flux of U momentum'
471     diagCode = 'UM MR '
472     diagMate = diagNum + 2
473     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
474     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
475     diagName = 'VISCy_Um'
476     diagTitle = 'Meridional Viscous Flux of U momentum'
477     diagCode = 'VZ MR '
478     diagMate = diagNum
479     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
480     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
481     diagName = 'VISrE_Um'
482     diagTitle = 'Vertical Viscous Flux of U momentum'
483     & //' (Explicit part)'
484     diagCode = 'WU LR '
485     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
486     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
487 jmc 1.1 #endif /* ALLOW_MOM_FLUXFORM */
488 jmc 1.10 diagName = 'VISrI_Um'
489     diagTitle = 'Vertical Viscous Flux of U momentum'
490     & //' (Implicit part)'
491     diagCode = 'WU LR '
492     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
493     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
494 jmc 1.1
495 jmc 1.10 C- Viscous flux of Meridional momentum :
496 jmc 1.1 #ifdef ALLOW_MOM_FLUXFORM
497 jmc 1.10 diagName = 'VISCx_Vm'
498     diagTitle = 'Zonal Viscous Flux of V momentum'
499     diagCode = 'UZ MR '
500     diagMate = diagNum + 2
501     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
502     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
503     diagName = 'VISCy_Vm'
504     diagTitle = 'Meridional Viscous Flux of V momentum'
505     diagCode = 'VM MR '
506     diagMate = diagNum
507     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
508     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
509     diagName = 'VISrE_Vm'
510     diagTitle = 'Vertical Viscous Flux of V momentum'
511     & //' (Explicit part)'
512     diagCode = 'WV LR '
513     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
514     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
515 jmc 1.1 #endif /* ALLOW_MOM_FLUXFORM */
516 jmc 1.10 diagName = 'VISrI_Vm'
517     diagTitle = 'Vertical Viscous Flux of V momentum'
518     & //' (Implicit part)'
519     diagCode = 'WV LR '
520     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
521     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
522 jmc 1.1
523     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
524    
525     #endif /* ALLOW_MOM_COMMON */
526     #endif /* ALLOW_DIAGNOSTICS */
527    
528     RETURN
529     END

  ViewVC Help
Powered by ViewVC 1.1.22