/[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.9 by jmc, Thu Jul 13 21:30:12 2006 UTC revision 1.10 by jmc, Tue Apr 22 22:03:17 2008 UTC
# Line 32  C     msgBuf      - Informational/error Line 32  C     msgBuf      - Informational/error
32  c     CHARACTER*(MAX_LEN_MBUF) msgBuf  c     CHARACTER*(MAX_LEN_MBUF) msgBuf
33    
34        INTEGER       diagNum        INTEGER       diagNum
35          INTEGER       diagMate
36        CHARACTER*8   diagName        CHARACTER*8   diagName
37        CHARACTER*16  diagCode        CHARACTER*16  diagCode
38        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
# Line 46  C--   add diagnostics for viscosity coef Line 47  C--   add diagnostics for viscosity coef
47        diagTitle = 'Harmonic Visc Coefficient (m2/s) (Zeta Pt)'        diagTitle = 'Harmonic Visc Coefficient (m2/s) (Zeta Pt)'
48        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
49        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
50        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
51       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
52    
53        diagName  = 'VISCA4Z '        diagName  = 'VISCA4Z '
54        diagTitle = 'Biharmonic Visc Coefficient (m4/s) (Zeta Pt)'        diagTitle = 'Biharmonic Visc Coefficient (m4/s) (Zeta Pt)'
55        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
56        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
57        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
58       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
59    
60        diagName  = 'VISCAHD '        diagName  = 'VISCAHD '
61        diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (Div Pt)'        diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (Div Pt)'
62        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
63        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
64        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
65       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
66    
67        diagName  = 'VISCA4D '        diagName  = 'VISCA4D '
68        diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (Div Pt)'        diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (Div Pt)'
69        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
70        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
71        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
72       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
73    
74        diagName  = 'VISCAHW '        diagName  = 'VISCAHW '
75        diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (W Pt)'        diagTitle = 'Harmonic Viscosity Coefficient (m2/s) (W Pt)'
76        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
77        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
78        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
79       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
80    
81        diagName  = 'VISCA4W '        diagName  = 'VISCA4W '
82        diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (W Pt)'        diagTitle = 'Biharmonic Viscosity Coefficient (m4/s) (W Pt)'
83        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
84        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
85        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
86       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
87    
88        diagName  = 'VAHZMAX '        diagName  = 'VAHZMAX '
89        diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Zeta Pt)'        diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Zeta Pt)'
90        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
91        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
92        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
93       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
94    
95        diagName  = 'VA4ZMAX '        diagName  = 'VA4ZMAX '
96        diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Zeta Pt)'        diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Zeta Pt)'
97        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
98        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
99        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
100       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
101    
102        diagName  = 'VAHDMAX '        diagName  = 'VAHDMAX '
103        diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Div Pt)'        diagTitle = 'CFL-MAX Harm Visc Coefficient (m2/s) (Div Pt)'
104        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
105        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
106        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
107       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
108    
109        diagName  = 'VA4DMAX '        diagName  = 'VA4DMAX '
110        diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Div Pt)'        diagTitle = 'CFL-MAX Biharm Visc Coefficient (m4/s) (Div Pt)'
111        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
112        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
113        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
114       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
115    
116        diagName  = 'VAHZMIN '        diagName  = 'VAHZMIN '
117        diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Zeta Pt)'        diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Zeta Pt)'
118        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
119        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
120        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
121       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
122    
123        diagName  = 'VA4ZMIN '        diagName  = 'VA4ZMIN '
124        diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Zeta Pt)'        diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Zeta Pt)'
125        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
126        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
127        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
128       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
129    
130        diagName  = 'VAHDMIN '        diagName  = 'VAHDMIN '
131        diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Div Pt)'        diagTitle = 'RE-MIN Harm Visc Coefficient (m2/s) (Div Pt)'
132        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
133        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
134        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
135       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
136    
137        diagName  = 'VA4DMIN '        diagName  = 'VA4DMIN '
138        diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Div Pt)'        diagTitle = 'RE-MIN Biharm Visc Coefficient (m4/s) (Div Pt)'
139        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
140        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
141        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
142       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
143    
144        diagName  = 'VAHZLTH '        diagName  = 'VAHZLTH '
145        diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Zeta Pt)'        diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Zeta Pt)'
146        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
147        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
148        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
149       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
150    
151        diagName  = 'VA4ZLTH '        diagName  = 'VA4ZLTH '
152        diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Zeta Pt)'        diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Zeta Pt)'
153        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
154        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
155        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
156       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
157    
158        diagName  = 'VAHDLTH '        diagName  = 'VAHDLTH '
159        diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Div Pt)'        diagTitle = 'Leith Harm Visc Coefficient (m2/s) (Div Pt)'
160        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
161        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
162        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
163       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
164    
165        diagName  = 'VA4DLTH '        diagName  = 'VA4DLTH '
166        diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Div Pt)'        diagTitle = 'Leith Biharm Visc Coefficient (m4/s) (Div Pt)'
167        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
168        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
169        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
170       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
171    
172        diagName  = 'VAHZLTHD'        diagName  = 'VAHZLTHD'
173        diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Zeta Pt)'        diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Zeta Pt)'
174        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
175        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
176        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
177       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
178    
179        diagName  = 'VA4ZLTHD'        diagName  = 'VA4ZLTHD'
180        diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Zeta Pt)'        diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Zeta Pt)'
181        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
182        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
183        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
184       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
185    
186        diagName  = 'VAHDLTHD'        diagName  = 'VAHDLTHD'
187        diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Div Pt)'        diagTitle = 'LeithD Harm Visc Coefficient (m2/s) (Div Pt)'
188        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
189        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
190        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
191       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
192    
193        diagName  = 'VA4DLTHD'        diagName  = 'VA4DLTHD'
194        diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Div Pt)'        diagTitle = 'LeithD Biharm Visc Coefficient (m4/s) (Div Pt)'
195        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
196        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
197        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
198       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
199    
200        diagName  = 'VAHZSMAG'        diagName  = 'VAHZSMAG'
201        diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Zeta Pt)'        diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Zeta Pt)'
202        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
203        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
204        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
205       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
206    
207        diagName  = 'VA4ZSMAG'        diagName  = 'VA4ZSMAG'
208        diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Zeta Pt)'        diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Zeta Pt)'
209        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
210        diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
211        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
212       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
213    
214        diagName  = 'VAHDSMAG'        diagName  = 'VAHDSMAG'
215        diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Div Pt)'        diagTitle = 'Smagorinsky Harm Visc Coefficient (m2/s) (Div Pt)'
216        diagUnits = 'm^2/s           '        diagUnits = 'm^2/s           '
217        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
218        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
219       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
220    
221        diagName  = 'VA4DSMAG'        diagName  = 'VA4DSMAG'
222        diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Div Pt)'        diagTitle = 'Smagorinsky Biharm Visc Coeff. (m4/s) (Div Pt)'
223        diagUnits = 'm^4/s           '        diagUnits = 'm^4/s           '
224        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
225        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
226       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
227    
228    
229    C--   add diagnostics of simple derivative quantities:
230  C--    add diagnostics of simple derivative quantities:        diagName  = 'momKE   '
231         diagName  = 'momKE   '        diagTitle = 'Kinetic Energy (in momentum Eq.)'
232         diagTitle = 'Kinetic Energy (in momentum Eq.)'        diagCode  = 'SM      MR      '
233         diagCode  = 'SM      MR      '        diagUnits = 'm^2/s^2         '
234         diagUnits = 'm^2/s^2         '        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
235         CALL DIAGNOSTICS_ADD2LIST( diagNum,       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
      I              diagName, diagCode, diagUnits, diagTitle, myThid )  
236    
237  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
238         diagName  = 'momHDiv '        diagName  = 'momHDiv '
239         diagTitle = 'Horizontal Divergence (in momentum Eq.)'        diagTitle = 'Horizontal Divergence (in momentum Eq.)'
240         diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
241         diagUnits = 's^-1            '        diagUnits = 's^-1            '
242         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
243       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, 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_ADDTOLIST( diagNum,
250       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
251    
252         diagName  = 'Strain  '        diagName  = 'Strain  '
253         diagTitle = 'Horizontal Strain of Horizontal Velocities'        diagTitle = 'Horizontal Strain of Horizontal Velocities'
254         diagCode  = 'SZ      MR      '        diagCode  = 'SZ      MR      '
255         diagUnits = 's^-1            '        diagUnits = 's^-1            '
256         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
257       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
258    
259         diagName  = 'Tension '        diagName  = 'Tension '
260         diagTitle = 'Horizontal Tension of Horizontal Velocities'        diagTitle = 'Horizontal Tension of Horizontal Velocities'
261         diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
262         diagUnits = 's^-1            '        diagUnits = 's^-1            '
263         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
264       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, 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'        diagName  = 'UBotDrag'
271         diagTitle = 'U momentum tendency from Bottom Drag'        diagTitle = 'U momentum tendency from Bottom Drag'
272         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
273         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
274       I              diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
275         diagName  = 'VBotDrag'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
276         diagTitle = 'V momentum tendency from Bottom Drag'        diagName  = 'VBotDrag'
277         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '        diagTitle = 'V momentum tendency from Bottom Drag'
278         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      MR      '
279       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
280         diagName  = 'USidDrag'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
281         diagTitle = 'U momentum tendency from Side Drag'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
282         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagName  = 'USidDrag'
283         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagTitle = 'U momentum tendency from Side Drag'
284       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'UU      MR      '
285         diagName  = 'VSidDrag'        diagMate  = diagNum + 2
286         diagTitle = 'V momentum tendency from Side Drag'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
287         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
288         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagName  = 'VSidDrag'
289       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagTitle = 'V momentum tendency from Side Drag'
290          diagCode  = 'VV      MR      '
291         diagName  = 'Um_Diss '        diagMate  = diagNum
292         diagTitle = 'U momentum tendency from Dissipation'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
293         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
294         CALL DIAGNOSTICS_ADD2LIST( diagNum,  
295       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagName  = 'Um_Diss '
296         diagName  = 'Vm_Diss '        diagTitle = 'U momentum tendency from Dissipation'
297         diagTitle = 'V momentum tendency from Dissipation'        diagCode  = 'UU      MR      '
298         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '        diagMate  = diagNum + 2
299         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
300       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
301          diagName  = 'Vm_Diss '
302         diagName  = 'Um_Advec'        diagTitle = 'V momentum tendency from Dissipation'
303         diagTitle = 'U momentum tendency from Advection terms'        diagCode  = 'VV      MR      '
304         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagMate  = diagNum
305         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
306       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
307         diagName  = 'Vm_Advec'  
308         diagTitle = 'V momentum tendency from Advection terms'        diagName  = 'Um_Advec'
309         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '        diagTitle = 'U momentum tendency from Advection terms'
310         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'UU      MR      '
311       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum + 2
312          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
313         diagName  = 'Um_Cori '       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
314         diagTitle = 'U momentum tendency from Coriolis term'        diagName  = 'Vm_Advec'
315         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagTitle = 'V momentum tendency from Advection terms'
316         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      MR      '
317       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
318         diagName  = 'Vm_Cori '        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
319         diagTitle = 'V momentum tendency from Coriolis term'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
320         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '  
321         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagName  = 'Um_Cori '
322       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagTitle = 'U momentum tendency from Coriolis term'
323          diagCode  = 'UU      MR      '
324         diagName  = 'Um_Ext  '        diagMate  = diagNum + 2
325         diagTitle = 'U momentum tendency from external forcing'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
326         WRITE(diagCode,'(A,I3.3,A)') 'UU   ', diagNum, 'MR      '       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
327         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagName  = 'Vm_Cori '
328       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagTitle = 'V momentum tendency from Coriolis term'
329         diagName  = 'Vm_Ext  '        diagCode  = 'VV      MR      '
330         diagTitle = 'V momentum tendency from external forcing'        diagMate  = diagNum
331         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
332         CALL DIAGNOSTICS_ADD2LIST( diagNum,       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
333       I              diagName, diagCode, diagUnits, diagTitle, myThid )  
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    
347  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
348         diagName  = 'Um_AdvZ3'        diagName  = 'Um_AdvZ3'
349         diagTitle = 'U momentum tendency from Vorticity Advection'        diagTitle = 'U momentum tendency from Vorticity Advection'
350         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
351         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
352       I              diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
353         diagName  = 'Vm_AdvZ3'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
354         diagTitle = 'V momentum tendency from Vorticity Advection'        diagName  = 'Vm_AdvZ3'
355         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '        diagTitle = 'V momentum tendency from Vorticity Advection'
356         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      MR      '
357       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
358          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
359         diagName  = 'Um_AdvRe'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
360         diagTitle = 'U momentum tendency from vertical Advection'  
361       &           //' (Explicit part)'        diagName  = 'Um_AdvRe'
362         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagTitle = 'U momentum tendency from vertical Advection'
363         CALL DIAGNOSTICS_ADD2LIST( diagNum,       &          //' (Explicit part)'
364       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'UU      MR      '
365         diagName  = 'Vm_AdvRe'        diagMate  = diagNum + 2
366         diagTitle = 'V momentum tendency from vertical Advection'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
367       &           //' (Explicit part)'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
368         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '        diagName  = 'Vm_AdvRe'
369         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagTitle = 'V momentum tendency from vertical Advection'
370       I              diagName, diagCode, diagUnits, diagTitle, myThid )       &          //' (Explicit part)'
371          diagCode  = 'VV      MR      '
372          diagMate  = diagNum
373          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
374         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
375    
376  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
377    
378  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
379  C-     vertical momentum tendencies  C-    vertical momentum tendencies
380         IF ( usingPCoords ) THEN        IF ( usingPCoords ) THEN
381           diagUnits = 'Pa/s^2          '          diagUnits = 'Pa/s^2          '
382         ELSE        ELSE
383           diagUnits = 'm/s^2           '          diagUnits = 'm/s^2           '
384         ENDIF        ENDIF
385         diagName  = 'WSidDrag'        diagName  = 'Wm_Diss '
386         diagTitle = 'Vertical momentum tendency from Side Drag'        diagTitle = 'W momentum tendency from Dissipation'
387         diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
388         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
389       I              diagName, diagCode, diagUnits, diagTitle, myThid )       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  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
403    
404  C--    add diagnostics of advective & viscous flux :  C--   add diagnostics of advective & viscous flux :
405         IF ( usingPCoords ) THEN        IF ( usingPCoords ) THEN
406           diagUnits = 'Pa.m^3/s^2      '          diagUnits = 'Pa.m^3/s^2      '
407         ELSE        ELSE
408           diagUnits = 'm^4/s^2         '          diagUnits = 'm^4/s^2         '
409         ENDIF        ENDIF
410    
411  C-     Advective flux of Zonal momentum:  C-     Advective flux of Zonal momentum:
412  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
413         diagName  = 'ADVx_Um '        diagName  = 'ADVx_Um '
414         diagTitle = 'Zonal      Advective Flux of U momentum'        diagTitle = 'Zonal      Advective Flux of U momentum'
415         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '        diagCode  = 'UM      MR      '
416         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
417       I              diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
418         diagName  = 'ADVy_Um '       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
419         diagTitle = 'Meridional Advective Flux of U momentum'        diagName  = 'ADVy_Um '
420         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '        diagTitle = 'Meridional Advective Flux of U momentum'
421         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VZ      MR      '
422       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
423         diagName  = 'ADVrE_Um'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
424         diagTitle = 'Vertical   Advective Flux of U momentum'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
425       &           //' (Explicit part)'        diagName  = 'ADVrE_Um'
426         diagCode  = 'WU      LR      '        diagTitle = 'Vertical   Advective Flux of U momentum'
427         CALL DIAGNOSTICS_ADD2LIST( diagNum,       &          //' (Explicit part)'
428       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'WU      LR      '
429          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
430         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
431  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
432  c      diagName  = 'ADVrI_Um'  c     diagName  = 'ADVrI_Um'
433  c      diagTitle = 'Vertical   Advective Flux of U momentum'  c     diagTitle = 'Vertical   Advective Flux of U momentum'
434  c    &           //' (Implicit part)'  c    &          //' (Implicit part)'
435  c      diagCode  = 'WU      LR      '  c     diagCode  = 'WU      LR      '
436  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
437  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
438    
439  C-     Advective flux of Meridional momentum:  C-    Advective flux of Meridional momentum:
440  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
441         diagName  = 'ADVx_Vm '        diagName  = 'ADVx_Vm '
442         diagTitle = 'Zonal      Advective Flux of V momentum'        diagTitle = 'Zonal      Advective Flux of V momentum'
443         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '        diagCode  = 'UZ      MR      '
444         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
445       I              diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
446         diagName  = 'ADVy_Vm '       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
447         diagTitle = 'Meridional Advective Flux of V momentum'        diagName  = 'ADVy_Vm '
448         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '        diagTitle = 'Meridional Advective Flux of V momentum'
449         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VM      MR      '
450       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
451         diagName  = 'ADVrE_Vm'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
452         diagTitle = 'Vertical   Advective Flux of V momentum'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
453       &           //' (Explicit part)'        diagName  = 'ADVrE_Vm'
454         diagCode  = 'WV      LR      '        diagTitle = 'Vertical   Advective Flux of V momentum'
455         CALL DIAGNOSTICS_ADD2LIST( diagNum,       &          //' (Explicit part)'
456       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'WV      LR      '
457          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
458         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
459  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
460  c      diagName  = 'ADVrI_Vm'  c     diagName  = 'ADVrI_Vm'
461  c      diagTitle = 'Vertical   Advective Flux of V momentum'  c     diagTitle = 'Vertical   Advective Flux of V momentum'
462  c    &           //' (Implicit part)'  c    &          //' (Implicit part)'
463  c      diagCode  = 'WV      LR      '  c     diagCode  = 'WV      LR      '
464  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
465  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
466    
467  c-     Viscous flux of Zonal momentum :  C-    Viscous flux of Zonal momentum :
468  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
469         diagName  = 'VISCx_Um'        diagName  = 'VISCx_Um'
470         diagTitle = 'Zonal      Viscous Flux of U momentum'        diagTitle = 'Zonal      Viscous Flux of U momentum'
471         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '        diagCode  = 'UM      MR      '
472         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
473       I              diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
474         diagName  = 'VISCy_Um'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
475         diagTitle = 'Meridional Viscous Flux of U momentum'        diagName  = 'VISCy_Um'
476         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '        diagTitle = 'Meridional Viscous Flux of U momentum'
477         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VZ      MR      '
478       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
479         diagName  = 'VISrE_Um'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
480         diagTitle = 'Vertical   Viscous Flux of U momentum'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
481       &           //' (Explicit part)'        diagName  = 'VISrE_Um'
482         diagCode  = 'WU      LR      '        diagTitle = 'Vertical   Viscous Flux of U momentum'
483         CALL DIAGNOSTICS_ADD2LIST( diagNum,       &          //' (Explicit part)'
484       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'WU      LR      '
485          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
486         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
487  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
488         diagName  = 'VISrI_Um'        diagName  = 'VISrI_Um'
489         diagTitle = 'Vertical   Viscous Flux of U momentum'        diagTitle = 'Vertical   Viscous Flux of U momentum'
490       &           //' (Implicit part)'       &          //' (Implicit part)'
491         diagCode  = 'WU      LR      '        diagCode  = 'WU      LR      '
492         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
493       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
494    
495  c-     Viscous flux of Meridional momentum :  C-    Viscous flux of Meridional momentum :
496  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
497         diagName  = 'VISCx_Vm'        diagName  = 'VISCx_Vm'
498         diagTitle = 'Zonal      Viscous Flux of V momentum'        diagTitle = 'Zonal      Viscous Flux of V momentum'
499         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '        diagCode  = 'UZ      MR      '
500         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
501       I              diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
502         diagName  = 'VISCy_Vm'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
503         diagTitle = 'Meridional Viscous Flux of V momentum'        diagName  = 'VISCy_Vm'
504         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '        diagTitle = 'Meridional Viscous Flux of V momentum'
505         CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VM      MR      '
506       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
507         diagName  = 'VISrE_Vm'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
508         diagTitle = 'Vertical   Viscous Flux of V momentum'       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
509       &           //' (Explicit part)'        diagName  = 'VISrE_Vm'
510         diagCode  = 'WV      LR      '        diagTitle = 'Vertical   Viscous Flux of V momentum'
511         CALL DIAGNOSTICS_ADD2LIST( diagNum,       &          //' (Explicit part)'
512       I              diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'WV      LR      '
513          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
514         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
515  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
516         diagName  = 'VISrI_Vm'        diagName  = 'VISrI_Vm'
517         diagTitle = 'Vertical   Viscous Flux of V momentum'        diagTitle = 'Vertical   Viscous Flux of V momentum'
518       &           //' (Implicit part)'       &          //' (Implicit part)'
519         diagCode  = 'WV      LR      '        diagCode  = 'WV      LR      '
520         CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
521       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
522    
523  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
524    
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22