/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

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

revision 1.31 by jmc, Thu May 31 18:22:00 2007 UTC revision 1.34 by jmc, Tue May 27 00:08:48 2008 UTC
# Line 66  C     fTitle     :: fixed boundary title Line 66  C     fTitle     :: fixed boundary title
66  C     pTitle     :: "Phi"  title  C     pTitle     :: "Phi"  title
67  C     sTitle     :: "salt" title  C     sTitle     :: "salt" title
68        INTEGER        diagNum        INTEGER        diagNum
69          INTEGER        diagMate
70        CHARACTER*8    diagName        CHARACTER*8    diagName
71        CHARACTER*16   diagCode        CHARACTER*16   diagCode
72        CHARACTER*16   diagUnits        CHARACTER*16   diagUnits
# Line 122  C-    free-surface (eTitle) and fixed-bo Line 123  C-    free-surface (eTitle) and fixed-bo
123  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
124  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
125    
126    
127        diagName  = 'ETAN    '        diagName  = 'ETAN    '
128        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
129  c     IF ( fluidIsWater .AND. usingZCoords )  c     IF ( fluidIsWater .AND. usingZCoords )
130  c    &diagTitle = 'Sea Surface Elevation'  c    &diagTitle = 'Sea Surface Elevation'
131        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
132        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
133        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
134       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
135    
136        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
137        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
138       I                          //' Anomaly', myThid )       I                          //' Anomaly', myThid )
139        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
140        diagCode  = 'SM      M1      '        diagCode  = 'SM P    M1      '
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  = 'DETADT2 '        diagName  = 'DETADT2 '
145        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
146       I                          //' Anomaly Tendency', myThid )       I                          //' Anomaly Tendency', myThid )
147        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
148        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
149        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
150       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
151    
152        diagName  = 'THETA   '        diagName  = 'THETA   '
153        diagTitle = 'Potential Temperature'        diagTitle = 'Potential Temperature'
154        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
155        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
156        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
157       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
158    
159  c     diagName  = 'SST     '  c     diagName  = 'SST     '
160  c     diagTitle = 'Sea Surface Temperature (degC,K)'  c     diagTitle = 'Sea Surface Temperature (degC,K)'
161  c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
162  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
163  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
164  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
165    
166        diagName  = 'SALT    '        diagName  = 'SALT    '
167        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
168        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
169        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
170        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
171       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
172    
173        diagName  = 'RELHUM  '        diagName  = 'RELHUM  '
174        diagTitle = 'Relative Humidity'        diagTitle = 'Relative Humidity'
175        diagUnits = 'percent         '        diagUnits = 'percent         '
176        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
177        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
178       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
179    
180  c     diagName  = 'SSS     '  c     diagName  = 'SSS     '
181  c     diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
182  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
183  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
184  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
185  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
186    
187        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
188        diagName  = 'SALTanom'        diagName  = 'SALTanom'
189        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
190        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
191        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
192        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
193       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
194        ENDIF        ENDIF
195    
196        diagName  = 'UVEL    '        diagName  = 'UVEL    '
197        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
198        diagUnits = 'm/s             '        diagUnits = 'm/s             '
199        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
200        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
201       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
202         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
203    
204        diagName  = 'VVEL    '        diagName  = 'VVEL    '
205        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
206        diagUnits = 'm/s             '        diagUnits = 'm/s             '
207        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
208        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
209       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
210         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
211    
212        diagName  = 'WVEL    '        diagName  = 'WVEL    '
213        diagTitle = 'Vertical Component of Velocity (r_units/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
214        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
215        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
216        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
217       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
218    
219        diagName  = 'THETASQ '        diagName  = 'THETASQ '
220        diagTitle = 'Square of Potential Temperature'        diagTitle = 'Square of Potential Temperature'
221        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
222        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
223        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
224       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
225    
226        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
227        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
228        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
229        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
230        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
231       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
232    
233        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
234        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
235        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
236        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
237        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
238        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
239       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
240        ENDIF        ENDIF
241    
242        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
243        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
244        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
245        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU P    MR      '
246        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
247       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
248         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
249    
250        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
251        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
252        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
253        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV P    MR      '
254        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
255       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
256         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
257    
258        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
259        diagTitle = 'Square of Vertical Comp of Velocity'        diagTitle = 'Square of Vertical Comp of Velocity'
260        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
261        diagCode  = 'WM      LR      '        diagCode  = 'WM P    LR      '
262        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
263       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
264    
265          diagName  = 'UE_VEL_C'
266          diagTitle = 'Eastward Velocity (m/s) (cell center)'
267          diagUnits = 'm/s             '
268          diagCode  = 'UM      MR      '
269          diagMate  = diagNum + 2
270          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
271         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
272    
273          diagName  = 'VN_VEL_C'
274          diagTitle = 'Northward Velocity (m/s) (cell center)'
275          diagUnits = 'm/s             '
276          diagCode  = 'VM      MR      '
277          diagMate  = diagNum
278          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
279         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
280    
281        diagName  = 'UV_VEL_C'        diagName  = 'UV_VEL_C'
282        diagTitle ='Product of horizontal Comp of velocity (cell center)'        diagTitle ='Product of horizontal Comp of velocity (cell center)'
283        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
284        WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '        diagCode  = 'UM      MR      '
285        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 1
286       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
287         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
288    
289        diagName  = 'UV_VEL_Z'        diagName  = 'UV_VEL_Z'
290        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
291        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
292        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        diagCode  = 'UZ      MR      '
293        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 1
294       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
295         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
296    
297        diagName  = 'WU_VEL  '        diagName  = 'WU_VEL  '
298        diagTitle = 'Vertical Transport of Zonal Momentum'        diagTitle = 'Vertical Transport of Zonal Momentum'
299        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
300        diagCode  = 'WU      LR      '        diagCode  = 'WU      LR      '
301        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
302       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
303    
304        diagName  = 'WV_VEL  '        diagName  = 'WV_VEL  '
305        diagTitle ='Vertical Transport of Meridional Momentum'        diagTitle ='Vertical Transport of Meridional Momentum'
306        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
307        diagCode  = 'WV      LR      '        diagCode  = 'WV      LR      '
308        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
309       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
310    
311        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
312        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
313        diagUnits = 'm/s             '        diagUnits = 'm/s             '
314        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
315        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
316       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
317         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
318    
319        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
320        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
321        diagUnits = 'm/s             '        diagUnits = 'm/s             '
322        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
323        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
324       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
325         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
326    
327        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
328        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
329        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
330        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
331        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
332       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
333    
334        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
335        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
336        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
337        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
338        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
339       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
340         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
341    
342        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
343        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
344        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
345        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
346        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
347       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
348         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
349    
350        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
351        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
352        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
353        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
354        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
355       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
356    
357        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
358        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
359       I                           //sTitle, myThid )       I                           //sTitle, myThid )
360        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
361        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
362        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
363       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
364         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
365    
366        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
367        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
368       I                           //sTitle, myThid )       I                           //sTitle, myThid )
369        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
370        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
371        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
372       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
373         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
374    
375        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
376        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
377       I                           //sTitle, myThid )       I                           //sTitle, myThid )
378        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
379        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
380        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
381       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
382    
383        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
384        diagTitle = 'Zonal Transport of Pot Temp'        diagTitle = 'Zonal Transport of Pot Temp'
385        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
386        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
387        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
388       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
389         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
390    
391        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
392        diagTitle = 'Meridional Transport of Pot Temp'        diagTitle = 'Meridional Transport of Pot Temp'
393        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
394        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
395        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
396       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
397         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
398    
399        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
400        diagTitle = 'Vertical Transport of Pot Temp'        diagTitle = 'Vertical Transport of Pot Temp'
401        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
402        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
403        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
404       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
405    
406        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
407        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
408       I                          //sTitle, myThid )       I                          //sTitle, myThid )
409        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
410        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
411        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
412       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
413         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
414    
415        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
416        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
417       I                          //sTitle, myThid )       I                          //sTitle, myThid )
418        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
419        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
420        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
421       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
422         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
423    
424        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
425        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
426       I                          //sTitle, myThid )       I                          //sTitle, myThid )
427        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
428        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
429        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
430       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
431    
432        diagName  = 'UVELPHI '        diagName  = 'UVELPHI '
433        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
434       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
435        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
436        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
437        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
438       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
439         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
440    
441        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
442        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
443       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
444        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
445        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
446        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
447       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
448         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
449    
450  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
451    
# Line 417  C---+----1----+----2----+----3----+----4 Line 453  C---+----1----+----2----+----3----+----4
453        diagTitle = 'Density Anomaly (=Rho-rhoConst)'        diagTitle = 'Density Anomaly (=Rho-rhoConst)'
454        diagUnits = 'kg/m^3          '        diagUnits = 'kg/m^3          '
455        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
456        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
457       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
458    
459        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
460        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
461        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6          '
462        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
463        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
464       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
465    
466        diagName  = 'URHOMASS'        diagName  = 'URHOMASS'
467        diagTitle = 'Zonal Transport of Density'        diagTitle = 'Zonal Transport of Density'
468        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
469        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
470        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
471       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
472         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
473    
474        diagName  = 'VRHOMASS'        diagName  = 'VRHOMASS'
475        diagTitle = 'Meridional Transport of Density'        diagTitle = 'Meridional Transport of Density'
476        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
477        WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '        diagCode  = 'VV      MR      '
478        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
479       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
480         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
481    
482        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
483        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Potential Density'
484        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
485        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
486        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
487       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
488    
489        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
490        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
491       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
492        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
493        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
494        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
495       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
496    
497        diagName  = 'PHIHYDSQ'        diagName  = 'PHIHYDSQ'
498        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
499       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
500        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
501        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
502        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
503       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
504    
505        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
506  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
# Line 470  c     diagTitle = 'ocean bottom pressure Line 508  c     diagTitle = 'ocean bottom pressure
508       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
509        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
510        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
511        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
512       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
513    
514        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
515  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
516        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
517       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
518        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
519        diagCode  = 'SM      M1      '        diagCode  = 'SM P    M1      '
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    #ifdef ALLOW_NONHYDROSTATIC
524          diagName  = 'PHI_NH  '
525          diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
526          diagUnits = 'm^2/s^2         '
527          diagCode  = 'SM      MR      '
528          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
529         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
530    #endif /* ALLOW_NONHYDROSTATIC */
531    
532        diagName  = 'MXLDEPTH'        diagName  = 'MXLDEPTH'
533        diagTitle = 'Mixed-Layer Depth (>0)'        diagTitle = 'Mixed-Layer Depth (>0)'
534        diagUnits = 'm               '        diagUnits = 'm               '
535        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
536        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
537       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
538    
539        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
540        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
541        diagUnits = 'kg/m^4          '        diagUnits = 'kg/m^4          '
542        IF ( usingPCoords ) diagUnits = 's^2/m^2         '        IF ( usingPCoords ) diagUnits = 's^2/m^2         '
543        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
544        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
545       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
546    
547        diagName  = 'CONVADJ '        diagName  = 'CONVADJ '
548        diagTitle = 'Convective Adjustment Index [0-1] '        diagTitle = 'Convective Adjustment Index [0-1] '
549        diagUnits = 'fraction        '        diagUnits = 'fraction        '
550        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
551        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
552       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
553    
554  C--   surface fluxes:  C--   surface fluxes:
555        diagName  = 'oceTAUX '        diagName  = 'oceTAUX '
556        diagTitle = 'zonal surface wind stress, >0 increases uVel'        diagTitle = 'zonal surface wind stress, >0 increases uVel'
557        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
558        diagCode  = 'SU      U1      '        diagCode  = 'UU      U1      '
559        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
560       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
561         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
562    
563        diagName  = 'oceTAUY '        diagName  = 'oceTAUY '
564        diagTitle = 'meridional surf. wind stress, >0 increases vVel'        diagTitle = 'meridional surf. wind stress, >0 increases vVel'
565        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
566        diagCode  = 'SV      U1      '        diagCode  = 'VV      U1      '
567        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
568       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
569         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
570    
571        diagName  = 'atmPload'        diagName  = 'atmPload'
572        diagTitle = 'Atmospheric pressure loading'        diagTitle = 'Atmospheric pressure loading'
573        diagUnits = 'Pa              '        diagUnits = 'Pa              '
574        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
575        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
576       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
577    
578        diagName  = 'sIceLoad'        diagName  = 'sIceLoad'
579        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
580        diagUnits = 'kg/m^2          '        diagUnits = 'kg/m^2          '
581        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
582        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
583       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
584    
585        diagName  = 'oceFWflx'        diagName  = 'oceFWflx'
586        diagTitle = 'net surface Fresh-Water flux into the ocean'        diagTitle = 'net surface Fresh-Water flux into the ocean'
587       &          //' (+=down), >0 decreases salinity'       &          //' (+=down), >0 decreases salinity'
588        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
589        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
590        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
591       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
592    
593        diagName  = 'oceSflux'        diagName  = 'oceSflux'
594        diagTitle = 'net surface Salt flux into the ocean (+=down),'        diagTitle = 'net surface Salt flux into the ocean (+=down),'
595       &          //' >0 increases salinity'       &          //' >0 increases salinity'
596        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
597        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
598        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
599       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
600    
601        diagName  = 'oceQnet '        diagName  = 'oceQnet '
602        diagTitle = 'net surface heat flux into the ocean (+=down),'        diagTitle = 'net surface heat flux into the ocean (+=down),'
603       &          //' >0 increases theta'       &          //' >0 increases theta'
604        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
605        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
606        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
607       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
608    
609        diagName  = 'oceQsw  '        diagName  = 'oceQsw  '
610        diagTitle = 'net Short-Wave radiation (+=down),'        diagTitle = 'net Short-Wave radiation (+=down),'
611       &          //' >0 increases theta'       &          //' >0 increases theta'
612        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
613        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
614        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
615       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
616    
617        diagName  = 'oceFreez'        diagName  = 'oceFreez'
618        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
619        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
620        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
621        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
622       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
623    
624        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
625        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
626        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
627        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
628        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
629       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
630    
631        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
632        diagTitle = 'surface salinity relaxation, >0 increases salt'        diagTitle = 'surface salinity relaxation, >0 increases salt'
633        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
634        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
635        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
636       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
637    
638        diagName  = 'surForcT'        diagName  = 'surForcT'
639        diagTitle = 'model surface forcing for Temperature,'        diagTitle = 'model surface forcing for Temperature,'
640       &          //' >0 increases theta'       &          //' >0 increases theta'
641        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
642        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
643        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
644       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
645    
646        diagName  = 'surForcS'        diagName  = 'surForcS'
647        diagTitle = 'model surface forcing for Salinity,'        diagTitle = 'model surface forcing for Salinity,'
648       &          //' >0 increases salinity'       &          //' >0 increases salinity'
649        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
650        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
651        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
652       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
653    
654        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
655        diagTitle = 'total heat flux (match heat-content variations),'        diagTitle = 'total heat flux (match heat-content variations),'
656       &          //' >0 increases theta'       &          //' >0 increases theta'
657        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
658        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
659        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
660       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
661    
662        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
663        diagTitle = 'total salt flux (match salt-content variations),'        diagTitle = 'total salt flux (match salt-content variations),'
664       &          //' >0 increases salt'       &          //' >0 increases salt'
665        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
666        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
667        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
668       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
669    
670  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
671    
# Line 625  c     diagTitle = 'Cell-Center r-Positio Line 674  c     diagTitle = 'Cell-Center r-Positio
674        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
675        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
676        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
677        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
678       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
679    
680        diagName  = 'RSURF   '        diagName  = 'RSURF   '
681  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
682        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
683        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
684        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
685        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
686       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
687    
688        diagName  = 'TOTUTEND'        diagName  = 'TOTUTEND'
689        diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
690        diagUnits = 'm/s^2           '        diagUnits = 'm/s^2           '
691        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
692        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
693       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
694         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
695    
696        diagName  = 'TOTVTEND'        diagName  = 'TOTVTEND'
697        diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
698        diagUnits = 'm/s^2           '        diagUnits = 'm/s^2           '
699        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
700        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
701       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
702         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
703    
704        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
705        diagTitle = 'Tendency of Potential Temperature'        diagTitle = 'Tendency of Potential Temperature'
706        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
707        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
708        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
709       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
710    
711        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
712        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
713        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
714        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
715        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
716       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
717    
718    
719        RETURN        RETURN
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22