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

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

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

revision 1.1 by jmc, Sun Sep 4 19:22:48 2005 UTC revision 1.9 by jmc, Thu Jul 13 21:30:12 2006 UTC
# Line 39  c     CHARACTER*(MAX_LEN_MBUF) msgBuf Line 39  c     CHARACTER*(MAX_LEN_MBUF) msgBuf
39    
40  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
41    
42  C--   Add diagnostics to the (long) list  C---  Add diagnostics to the (long) list of available diagnostics:
43    
44    C--   add diagnostics for viscosity coefficients:
45          diagName  = 'VISCAHZ '
46          diagTitle = 'Harmonic Visc Coefficient (m2/s) (Zeta Pt)'
47          diagUnits = 'm^2/s           '
48          diagCode  = 'SZ      MR      '
49          CALL DIAGNOSTICS_ADD2LIST( diagNum,
50         I          diagName, diagCode, diagUnits, diagTitle, myThid )
51    
52          diagName  = 'VISCA4Z '
53          diagTitle = 'Biharmonic Visc Coefficient (m4/s) (Zeta Pt)'
54          diagUnits = 'm^4/s           '
55          diagCode  = 'SZ      MR      '
56          CALL DIAGNOSTICS_ADD2LIST( diagNum,
57         I          diagName, diagCode, diagUnits, diagTitle, myThid )
58    
59          diagName  = 'VISCAHD '
60          diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (Div Pt)'
61          diagUnits = 'm^2/s           '
62          diagCode  = 'SM      MR      '
63          CALL DIAGNOSTICS_ADD2LIST( diagNum,
64         I          diagName, diagCode, diagUnits, diagTitle, myThid )
65    
66          diagName  = 'VISCA4D '
67          diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (Div Pt)'
68          diagUnits = 'm^4/s           '
69          diagCode  = 'SM      MR      '
70          CALL DIAGNOSTICS_ADD2LIST( diagNum,
71         I          diagName, diagCode, diagUnits, diagTitle, myThid )
72    
73          diagName  = 'VISCAHW '
74          diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (W Pt)'
75          diagUnits = 'm^2/s           '
76          diagCode  = 'WM      LR      '
77          CALL DIAGNOSTICS_ADD2LIST( diagNum,
78         I          diagName, diagCode, diagUnits, diagTitle, myThid )
79    
80          diagName  = 'VISCA4W '
81          diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (W Pt)'
82          diagUnits = 'm^4/s           '
83          diagCode  = 'WM      LR      '
84          CALL DIAGNOSTICS_ADD2LIST( diagNum,
85         I          diagName, diagCode, diagUnits, diagTitle, myThid )
86    
87          diagName  = 'VAHZMAX '
88          diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Zeta Pt)'
89          diagUnits = 'm^2/s           '
90          diagCode  = 'SZ      MR      '
91          CALL DIAGNOSTICS_ADD2LIST( diagNum,
92         I          diagName, diagCode, diagUnits, diagTitle, myThid )
93    
94          diagName  = 'VA4ZMAX '
95          diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Zeta Pt)'
96          diagUnits = 'm^4/s           '
97          diagCode  = 'SZ      MR      '
98          CALL DIAGNOSTICS_ADD2LIST( diagNum,
99         I          diagName, diagCode, diagUnits, diagTitle, myThid )
100    
101          diagName  = 'VAHDMAX '
102          diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Div Pt)'
103          diagUnits = 'm^2/s           '
104          diagCode  = 'SM      MR      '
105          CALL DIAGNOSTICS_ADD2LIST( diagNum,
106         I          diagName, diagCode, diagUnits, diagTitle, myThid )
107    
108          diagName  = 'VA4DMAX '
109          diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Div Pt)'
110          diagUnits = 'm^4/s           '
111          diagCode  = 'SM      MR      '
112          CALL DIAGNOSTICS_ADD2LIST( diagNum,
113         I          diagName, diagCode, diagUnits, diagTitle, myThid )
114    
115          diagName  = 'VAHZMIN '
116          diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Zeta Pt)'
117          diagUnits = 'm^2/s           '
118          diagCode  = 'SZ      MR      '
119          CALL DIAGNOSTICS_ADD2LIST( diagNum,
120         I          diagName, diagCode, diagUnits, diagTitle, myThid )
121    
122          diagName  = 'VA4ZMIN '
123          diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Zeta Pt)'
124          diagUnits = 'm^4/s           '
125          diagCode  = 'SZ      MR      '
126          CALL DIAGNOSTICS_ADD2LIST( diagNum,
127         I          diagName, diagCode, diagUnits, diagTitle, myThid )
128    
129          diagName  = 'VAHDMIN '
130          diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Div Pt)'
131          diagUnits = 'm^2/s           '
132          diagCode  = 'SM      MR      '
133          CALL DIAGNOSTICS_ADD2LIST( diagNum,
134         I          diagName, diagCode, diagUnits, diagTitle, myThid )
135    
136          diagName  = 'VA4DMIN '
137          diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Div Pt)'
138          diagUnits = 'm^4/s           '
139          diagCode  = 'SM      MR      '
140          CALL DIAGNOSTICS_ADD2LIST( diagNum,
141         I          diagName, diagCode, diagUnits, diagTitle, myThid )
142    
143          diagName  = 'VAHZLTH '
144          diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Zeta Pt)'
145          diagUnits = 'm^2/s           '
146          diagCode  = 'SZ      MR      '
147          CALL DIAGNOSTICS_ADD2LIST( diagNum,
148         I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150          diagName  = 'VA4ZLTH '
151          diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Zeta Pt)'
152          diagUnits = 'm^4/s           '
153          diagCode  = 'SZ      MR      '
154          CALL DIAGNOSTICS_ADD2LIST( diagNum,
155         I          diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157          diagName  = 'VAHDLTH '
158          diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Div Pt)'
159          diagUnits = 'm^2/s           '
160          diagCode  = 'SM      MR      '
161          CALL DIAGNOSTICS_ADD2LIST( diagNum,
162         I          diagName, diagCode, diagUnits, diagTitle, myThid )
163    
164          diagName  = 'VA4DLTH '
165          diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Div Pt)'
166          diagUnits = 'm^4/s           '
167          diagCode  = 'SM      MR      '
168          CALL DIAGNOSTICS_ADD2LIST( diagNum,
169         I          diagName, diagCode, diagUnits, diagTitle, myThid )
170    
171          diagName  = 'VAHZLTHD'
172          diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Zeta Pt)'
173          diagUnits = 'm^2/s           '
174          diagCode  = 'SZ      MR      '
175          CALL DIAGNOSTICS_ADD2LIST( diagNum,
176         I          diagName, diagCode, diagUnits, diagTitle, myThid )
177    
178          diagName  = 'VA4ZLTHD'
179          diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Zeta Pt)'
180          diagUnits = 'm^4/s           '
181          diagCode  = 'SZ      MR      '
182          CALL DIAGNOSTICS_ADD2LIST( diagNum,
183         I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185          diagName  = 'VAHDLTHD'
186          diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Div Pt)'
187          diagUnits = 'm^2/s           '
188          diagCode  = 'SM      MR      '
189          CALL DIAGNOSTICS_ADD2LIST( diagNum,
190         I          diagName, diagCode, diagUnits, diagTitle, myThid )
191    
192          diagName  = 'VA4DLTHD'
193          diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Div Pt)'
194          diagUnits = 'm^4/s           '
195          diagCode  = 'SM      MR      '
196          CALL DIAGNOSTICS_ADD2LIST( diagNum,
197         I          diagName, diagCode, diagUnits, diagTitle, myThid )
198    
199          diagName  = 'VAHZSMAG'
200          diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Zeta Pt)'
201          diagUnits = 'm^2/s           '
202          diagCode  = 'SZ      MR      '
203          CALL DIAGNOSTICS_ADD2LIST( diagNum,
204         I          diagName, diagCode, diagUnits, diagTitle, myThid )
205    
206          diagName  = 'VA4ZSMAG'
207          diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Zeta Pt)'
208          diagUnits = 'm^4/s           '
209          diagCode  = 'SZ      MR      '
210          CALL DIAGNOSTICS_ADD2LIST( diagNum,
211         I          diagName, diagCode, diagUnits, diagTitle, myThid )
212    
213          diagName  = 'VAHDSMAG'
214          diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Div Pt)'
215          diagUnits = 'm^2/s           '
216          diagCode  = 'SM      MR      '
217          CALL DIAGNOSTICS_ADD2LIST( diagNum,
218         I          diagName, diagCode, diagUnits, diagTitle, myThid )
219    
220          diagName  = 'VA4DSMAG'
221          diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Div Pt)'
222          diagUnits = 'm^4/s           '
223          diagCode  = 'SM      MR      '
224          CALL DIAGNOSTICS_ADD2LIST( diagNum,
225         I          diagName, diagCode, diagUnits, diagTitle, myThid )
226    
227    
228    
229  C--    add diagnostics of simple derivative quantities:  C--    add diagnostics of simple derivative quantities:
230         diagName  = 'momKE   '         diagName  = 'momKE   '
# Line 47  C--    add diagnostics of simple derivat Line 232  C--    add diagnostics of simple derivat
232         diagCode  = 'SM      MR      '         diagCode  = 'SM      MR      '
233         diagUnits = 'm^2/s^2         '         diagUnits = 'm^2/s^2         '
234         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
235       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
236    
237  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
238         diagName  = 'momHDiv '         diagName  = 'momHDiv '
# Line 55  C--    add diagnostics of simple derivat Line 240  C--    add diagnostics of simple derivat
240         diagCode  = 'SM      MR      '         diagCode  = 'SM      MR      '
241         diagUnits = 's^-1            '         diagUnits = 's^-1            '
242         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
243       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
244    
245         diagName  = 'momVort3'         diagName  = 'momVort3'
246         diagTitle = '3rd component (vertical) of Vorticity'         diagTitle = '3rd component (vertical) of Vorticity'
247         diagCode  = 'SZ      MR      '         diagCode  = 'SZ      MR      '
248         diagUnits = 's^-1            '         diagUnits = 's^-1            '
249         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
250       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
251    
252           diagName  = 'Strain  '
253           diagTitle = 'Horizontal Strain of Horizontal Velocities'
254           diagCode  = 'SZ      MR      '
255           diagUnits = 's^-1            '
256           CALL DIAGNOSTICS_ADD2LIST( diagNum,
257         I              diagName, diagCode, diagUnits, diagTitle, myThid )
258    
259           diagName  = 'Tension '
260           diagTitle = 'Horizontal Tension of Horizontal Velocities'
261           diagCode  = 'SM      MR      '
262           diagUnits = 's^-1            '
263           CALL DIAGNOSTICS_ADD2LIST( diagNum,
264         I              diagName, diagCode, diagUnits, diagTitle, myThid )
265  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
266    
267  C--    add diagnostics of tendencies from each terms  C--    add diagnostics of tendencies from each terms
268         diagUnits = 'm/s^2           '         diagUnits = 'm/s^2           '
269    
270           diagName  = 'UBotDrag'
271           diagTitle = 'U momentum tendency from Bottom Drag'
272           WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
273           CALL DIAGNOSTICS_ADD2LIST( diagNum,
274         I              diagName, diagCode, diagUnits, diagTitle, myThid )
275           diagName  = 'VBotDrag'
276           diagTitle = 'V momentum tendency from Bottom Drag'
277           WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
278           CALL DIAGNOSTICS_ADD2LIST( diagNum,
279         I              diagName, diagCode, diagUnits, diagTitle, myThid )
280           diagName  = 'USidDrag'
281           diagTitle = 'U momentum tendency from Side Drag'
282           WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
283           CALL DIAGNOSTICS_ADD2LIST( diagNum,
284         I              diagName, diagCode, diagUnits, diagTitle, myThid )
285           diagName  = 'VSidDrag'
286           diagTitle = 'V momentum tendency from Side Drag'
287           WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
288           CALL DIAGNOSTICS_ADD2LIST( diagNum,
289         I              diagName, diagCode, diagUnits, diagTitle, myThid )
290    
291         diagName  = 'Um_Diss '         diagName  = 'Um_Diss '
292         diagTitle = 'U momentum tendency from Dissipation'         diagTitle = 'U momentum tendency from Dissipation'
293         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
294         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
295       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
296         diagName  = 'Vm_Diss '         diagName  = 'Vm_Diss '
297         diagTitle = 'V momentum tendency from Dissipation'         diagTitle = 'V momentum tendency from Dissipation'
298         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
299         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
300       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
301    
302         diagName  = 'Um_Advec'         diagName  = 'Um_Advec'
303         diagTitle = 'U momentum tendency from Advection terms'         diagTitle = 'U momentum tendency from Advection terms'
304         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
305         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
306       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
307         diagName  = 'Vm_Advec'         diagName  = 'Vm_Advec'
308         diagTitle = 'V momentum tendency from Advection terms'         diagTitle = 'V momentum tendency from Advection terms'
309         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
310         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
311       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
312    
313         diagName  = 'Um_Cori '         diagName  = 'Um_Cori '
314         diagTitle = 'U momentum tendency from Coriolis term'         diagTitle = 'U momentum tendency from Coriolis term'
315         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
316         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
317       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
318         diagName  = 'Vm_Cori '         diagName  = 'Vm_Cori '
319         diagTitle = 'V momentum tendency from Coriolis term'         diagTitle = 'V momentum tendency from Coriolis term'
320         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
321         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
322       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
323    
324           diagName  = 'Um_Ext  '
325           diagTitle = 'U momentum tendency from external forcing'
326           WRITE(diagCode,'(A,I3.3,A)') 'UU   ', diagNum, 'MR      '
327           CALL DIAGNOSTICS_ADD2LIST( diagNum,
328         I              diagName, diagCode, diagUnits, diagTitle, myThid )
329           diagName  = 'Vm_Ext  '
330           diagTitle = 'V momentum tendency from external forcing'
331           WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
332           CALL DIAGNOSTICS_ADD2LIST( diagNum,
333         I              diagName, diagCode, diagUnits, diagTitle, myThid )
334    
335  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
336         diagName  = 'Um_AdvZ3'         diagName  = 'Um_AdvZ3'
337         diagTitle = 'U momentum tendency from Vorticity Advection'         diagTitle = 'U momentum tendency from Vorticity Advection'
338         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
339         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
340       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
341         diagName  = 'Vm_AdvZ3'         diagName  = 'Vm_AdvZ3'
342         diagTitle = 'V momentum tendency from Vorticity Advection'         diagTitle = 'V momentum tendency from Vorticity Advection'
343         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
344         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
345       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
346    
347         diagName  = 'Um_AdvRe'         diagName  = 'Um_AdvRe'
348         diagTitle = 'U momentum tendency from vertical Advection'         diagTitle = 'U momentum tendency from vertical Advection'
349       &           //' (Explicit part)'       &           //' (Explicit part)'
350         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
351         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
352       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
353         diagName  = 'Vm_AdvRe'         diagName  = 'Vm_AdvRe'
354         diagTitle = 'V momentum tendency from vertical Advection'         diagTitle = 'V momentum tendency from vertical Advection'
355       &           //' (Explicit part)'       &           //' (Explicit part)'
356         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
357         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
358       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
361    
362    #ifdef ALLOW_NONHYDROSTATIC
363    C-     vertical momentum tendencies
364           IF ( usingPCoords ) THEN
365             diagUnits = 'Pa/s^2          '
366           ELSE
367             diagUnits = 'm/s^2           '
368           ENDIF
369           diagName  = 'WSidDrag'
370           diagTitle = 'Vertical momentum tendency from Side Drag'
371           diagCode  = 'WM      LR      '
372           CALL DIAGNOSTICS_ADD2LIST( diagNum,
373         I              diagName, diagCode, diagUnits, diagTitle, myThid )
374    #endif /* ALLOW_NONHYDROSTATIC */
375    
376  C--    add diagnostics of advective & viscous flux :  C--    add diagnostics of advective & viscous flux :
377         IF ( usingPCoords ) THEN         IF ( usingPCoords ) THEN
378           diagUnits = 'Pa.m^3/s^2      '           diagUnits = 'Pa.m^3/s^2      '
# Line 140  C-     Advective flux of Zonal momentum: Line 386  C-     Advective flux of Zonal momentum:
386         diagTitle = 'Zonal      Advective Flux of U momentum'         diagTitle = 'Zonal      Advective Flux of U momentum'
387         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '
388         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
389       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
390         diagName  = 'ADVy_Um '         diagName  = 'ADVy_Um '
391         diagTitle = 'Meridional Advective Flux of U momentum'         diagTitle = 'Meridional Advective Flux of U momentum'
392         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '
393         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
394       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
395         diagName  = 'ADVrE_Um'         diagName  = 'ADVrE_Um'
396         diagTitle = 'Vertical   Advective Flux of U momentum'         diagTitle = 'Vertical   Advective Flux of U momentum'
397       &           //' (Explicit part)'       &           //' (Explicit part)'
398         diagCode  = 'WU      LR      '         diagCode  = 'WU      LR      '
399         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
400       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
401  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
402  c      diagName  = 'ADVrI_Um'  c      diagName  = 'ADVrI_Um'
403  c      diagTitle = 'Vertical   Advective Flux of U momentum'  c      diagTitle = 'Vertical   Advective Flux of U momentum'
404  c    &           //' (Implicit part)'  c    &           //' (Implicit part)'
405  c      diagCode  = 'WU      LR      '  c      diagCode  = 'WU      LR      '
406  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,
407  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )
408    
409  C-     Advective flux of Meridional momentum:  C-     Advective flux of Meridional momentum:
410  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
# Line 166  C-     Advective flux of Meridional mome Line 412  C-     Advective flux of Meridional mome
412         diagTitle = 'Zonal      Advective Flux of V momentum'         diagTitle = 'Zonal      Advective Flux of V momentum'
413         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '
414         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
415       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
416         diagName  = 'ADVy_Vm '         diagName  = 'ADVy_Vm '
417         diagTitle = 'Meridional Advective Flux of V momentum'         diagTitle = 'Meridional Advective Flux of V momentum'
418         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '
419         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
420       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
421         diagName  = 'ADVrE_Vm'         diagName  = 'ADVrE_Vm'
422         diagTitle = 'Vertical   Advective Flux of V momentum'         diagTitle = 'Vertical   Advective Flux of V momentum'
423       &           //' (Explicit part)'       &           //' (Explicit part)'
424         diagCode  = 'WV      LR      '         diagCode  = 'WV      LR      '
425         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
426       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
427  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
428  c      diagName  = 'ADVrI_Vm'  c      diagName  = 'ADVrI_Vm'
429  c      diagTitle = 'Vertical   Advective Flux of V momentum'  c      diagTitle = 'Vertical   Advective Flux of V momentum'
430  c    &           //' (Implicit part)'  c    &           //' (Implicit part)'
431  c      diagCode  = 'WV      LR      '  c      diagCode  = 'WV      LR      '
432  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,
433  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )
434    
435  c-     Viscous flux of Zonal momentum :  c-     Viscous flux of Zonal momentum :
436  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
# Line 192  c-     Viscous flux of Zonal momentum : Line 438  c-     Viscous flux of Zonal momentum :
438         diagTitle = 'Zonal      Viscous Flux of U momentum'         diagTitle = 'Zonal      Viscous Flux of U momentum'
439         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '
440         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
441       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
442         diagName  = 'VISCy_Um'         diagName  = 'VISCy_Um'
443         diagTitle = 'Meridional Viscous Flux of U momentum'         diagTitle = 'Meridional Viscous Flux of U momentum'
444         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '
445         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
446       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
447         diagName  = 'VISrE_Um'         diagName  = 'VISrE_Um'
448         diagTitle = 'Vertical   Viscous Flux of U momentum'         diagTitle = 'Vertical   Viscous Flux of U momentum'
449       &           //' (Explicit part)'       &           //' (Explicit part)'
450         diagCode  = 'WU      LR      '         diagCode  = 'WU      LR      '
451         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
452       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
453  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
454         diagName  = 'VISrI_Um'         diagName  = 'VISrI_Um'
455         diagTitle = 'Vertical   Viscous Flux of U momentum'         diagTitle = 'Vertical   Viscous Flux of U momentum'
456       &           //' (Implicit part)'       &           //' (Implicit part)'
457         diagCode  = 'WU      LR      '         diagCode  = 'WU      LR      '
458         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
459       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
460    
461  c-     Viscous flux of Meridional momentum :  c-     Viscous flux of Meridional momentum :
462  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
# Line 218  c-     Viscous flux of Meridional moment Line 464  c-     Viscous flux of Meridional moment
464         diagTitle = 'Zonal      Viscous Flux of V momentum'         diagTitle = 'Zonal      Viscous Flux of V momentum'
465         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '
466         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
467       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
468         diagName  = 'VISCy_Vm'         diagName  = 'VISCy_Vm'
469         diagTitle = 'Meridional Viscous Flux of V momentum'         diagTitle = 'Meridional Viscous Flux of V momentum'
470         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '
471         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
472       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
473         diagName  = 'VISrE_Vm'         diagName  = 'VISrE_Vm'
474         diagTitle = 'Vertical   Viscous Flux of V momentum'         diagTitle = 'Vertical   Viscous Flux of V momentum'
475       &           //' (Explicit part)'       &           //' (Explicit part)'
476         diagCode  = 'WV      LR      '         diagCode  = 'WV      LR      '
477         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
478       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
479  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
480         diagName  = 'VISrI_Vm'         diagName  = 'VISrI_Vm'
481         diagTitle = 'Vertical   Viscous Flux of V momentum'         diagTitle = 'Vertical   Viscous Flux of V momentum'
482       &           //' (Implicit part)'       &           //' (Implicit part)'
483         diagCode  = 'WV      LR      '         diagCode  = 'WV      LR      '
484         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
485       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
486    
487  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
488    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22