/[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.32 by jmc, Tue Feb 5 15:34:32 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  = 'UV_VEL_C'        diagName  = 'UV_VEL_C'
266        diagTitle ='Product of horizontal Comp of velocity (cell center)'        diagTitle ='Product of horizontal Comp of velocity (cell center)'
267        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
268        WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '        diagCode  = 'UM      MR      '
269        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 1
270       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
271         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
272    
273        diagName  = 'UV_VEL_Z'        diagName  = 'UV_VEL_Z'
274        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
275        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
276        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        diagCode  = 'UZ      MR      '
277        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 1
278       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
279         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
280    
281        diagName  = 'WU_VEL  '        diagName  = 'WU_VEL  '
282        diagTitle = 'Vertical Transport of Zonal Momentum'        diagTitle = 'Vertical Transport of Zonal Momentum'
283        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
284        diagCode  = 'WU      LR      '        diagCode  = 'WU      LR      '
285        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
286       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
287    
288        diagName  = 'WV_VEL  '        diagName  = 'WV_VEL  '
289        diagTitle ='Vertical Transport of Meridional Momentum'        diagTitle ='Vertical Transport of Meridional Momentum'
290        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
291        diagCode  = 'WV      LR      '        diagCode  = 'WV      LR      '
292        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
293       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
294    
295        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
296        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
297        diagUnits = 'm/s             '        diagUnits = 'm/s             '
298        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
299        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
300       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
301         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
302    
303        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
304        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
305        diagUnits = 'm/s             '        diagUnits = 'm/s             '
306        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
307        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
308       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
309         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
310    
311        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
312        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
313        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
314        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
315        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
316       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
317    
318        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
319        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
320        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
321        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
322        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
323       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
324         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
325    
326        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
327        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
328        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
329        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
330        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
331       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
332         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
333    
334        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
335        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
336        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
337        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
338        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
339       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
340    
341        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
342        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
343       I                           //sTitle, myThid )       I                           //sTitle, myThid )
344        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
345        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
346        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
347       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
348         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
349    
350        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
351        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
352       I                           //sTitle, myThid )       I                           //sTitle, myThid )
353        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
354        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
355        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
356       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
357         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
358    
359        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
360        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
361       I                           //sTitle, myThid )       I                           //sTitle, myThid )
362        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
363        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
364        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
365       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
366    
367        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
368        diagTitle = 'Zonal Transport of Pot Temp'        diagTitle = 'Zonal Transport of Pot Temp'
369        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
370        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
371        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
372       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
373         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
374    
375        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
376        diagTitle = 'Meridional Transport of Pot Temp'        diagTitle = 'Meridional Transport of Pot Temp'
377        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
378        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
379        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
380       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
381         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
382    
383        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
384        diagTitle = 'Vertical Transport of Pot Temp'        diagTitle = 'Vertical Transport of Pot Temp'
385        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
386        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
387        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
388       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
389    
390        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
391        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
392       I                          //sTitle, myThid )       I                          //sTitle, myThid )
393        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
394        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
395        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
396       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
397         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
398    
399        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
400        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
401       I                          //sTitle, myThid )       I                          //sTitle, myThid )
402        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
403        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
404        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
405       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
406         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
407    
408        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
409        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
410       I                          //sTitle, myThid )       I                          //sTitle, myThid )
411        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
412        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
413        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
414       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
415    
416        diagName  = 'UVELPHI '        diagName  = 'UVELPHI '
417        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
418       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
419        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
420        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
421        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
422       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
423         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
424    
425        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
426        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
427       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
428        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
429        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
430        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
431       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
432         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
433    
434  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
435    
# Line 417  C---+----1----+----2----+----3----+----4 Line 437  C---+----1----+----2----+----3----+----4
437        diagTitle = 'Density Anomaly (=Rho-rhoConst)'        diagTitle = 'Density Anomaly (=Rho-rhoConst)'
438        diagUnits = 'kg/m^3          '        diagUnits = 'kg/m^3          '
439        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
440        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
441       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
442    
443        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
444        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
445        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6          '
446        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
447        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
448       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
449    
450        diagName  = 'URHOMASS'        diagName  = 'URHOMASS'
451        diagTitle = 'Zonal Transport of Density'        diagTitle = 'Zonal Transport of Density'
452        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
453        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
454        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
455       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
456         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
457    
458        diagName  = 'VRHOMASS'        diagName  = 'VRHOMASS'
459        diagTitle = 'Meridional Transport of Density'        diagTitle = 'Meridional Transport of Density'
460        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
461        WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '        diagCode  = 'VV      MR      '
462        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
463       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
464         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
465    
466        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
467        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Potential Density'
468        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
469        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
470        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
471       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
472    
473        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
474        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
475       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
476        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
477        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
478        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
479       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
480    
481        diagName  = 'PHIHYDSQ'        diagName  = 'PHIHYDSQ'
482        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
483       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
484        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
485        diagCode  = 'SM      MR      '        diagCode  = 'SM P    MR      '
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  = 'PHIBOT  '        diagName  = 'PHIBOT  '
490  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 492  c     diagTitle = 'ocean bottom pressure
492       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
493        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
494        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
495        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
496       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
497    
498        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
499  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
500        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
501       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
502        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
503        diagCode  = 'SM      M1      '        diagCode  = 'SM P    M1      '
504        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
505       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
506    
507        diagName  = 'MXLDEPTH'        diagName  = 'MXLDEPTH'
508        diagTitle = 'Mixed-Layer Depth (>0)'        diagTitle = 'Mixed-Layer Depth (>0)'
509        diagUnits = 'm               '        diagUnits = 'm               '
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  = 'DRHODR  '        diagName  = 'DRHODR  '
515        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
516        diagUnits = 'kg/m^4          '        diagUnits = 'kg/m^4          '
517        IF ( usingPCoords ) diagUnits = 's^2/m^2         '        IF ( usingPCoords ) diagUnits = 's^2/m^2         '
518        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
519        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
520       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
521    
522        diagName  = 'CONVADJ '        diagName  = 'CONVADJ '
523        diagTitle = 'Convective Adjustment Index [0-1] '        diagTitle = 'Convective Adjustment Index [0-1] '
524        diagUnits = 'fraction        '        diagUnits = 'fraction        '
525        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
526        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
527       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
528    
529  C--   surface fluxes:  C--   surface fluxes:
530        diagName  = 'oceTAUX '        diagName  = 'oceTAUX '
531        diagTitle = 'zonal surface wind stress, >0 increases uVel'        diagTitle = 'zonal surface wind stress, >0 increases uVel'
532        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
533        diagCode  = 'SU      U1      '        diagCode  = 'UU      U1      '
534        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
535       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
536         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
537    
538        diagName  = 'oceTAUY '        diagName  = 'oceTAUY '
539        diagTitle = 'meridional surf. wind stress, >0 increases vVel'        diagTitle = 'meridional surf. wind stress, >0 increases vVel'
540        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
541        diagCode  = 'SV      U1      '        diagCode  = 'VV      U1      '
542        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
543       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
544         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
545    
546        diagName  = 'atmPload'        diagName  = 'atmPload'
547        diagTitle = 'Atmospheric pressure loading'        diagTitle = 'Atmospheric pressure loading'
548        diagUnits = 'Pa              '        diagUnits = 'Pa              '
549        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
550        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
551       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
552    
553        diagName  = 'sIceLoad'        diagName  = 'sIceLoad'
554        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
555        diagUnits = 'kg/m^2          '        diagUnits = 'kg/m^2          '
556        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
557        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
558       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
559    
560        diagName  = 'oceFWflx'        diagName  = 'oceFWflx'
561        diagTitle = 'net surface Fresh-Water flux into the ocean'        diagTitle = 'net surface Fresh-Water flux into the ocean'
562       &          //' (+=down), >0 decreases salinity'       &          //' (+=down), >0 decreases salinity'
563        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
564        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
565        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
566       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
567    
568        diagName  = 'oceSflux'        diagName  = 'oceSflux'
569        diagTitle = 'net surface Salt flux into the ocean (+=down),'        diagTitle = 'net surface Salt flux into the ocean (+=down),'
570       &          //' >0 increases salinity'       &          //' >0 increases salinity'
571        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
572        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
573        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
574       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
575    
576        diagName  = 'oceQnet '        diagName  = 'oceQnet '
577        diagTitle = 'net surface heat flux into the ocean (+=down),'        diagTitle = 'net surface heat flux into the ocean (+=down),'
578       &          //' >0 increases theta'       &          //' >0 increases theta'
579        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
580        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
581        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
582       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
583    
584        diagName  = 'oceQsw  '        diagName  = 'oceQsw  '
585        diagTitle = 'net Short-Wave radiation (+=down),'        diagTitle = 'net Short-Wave radiation (+=down),'
586       &          //' >0 increases theta'       &          //' >0 increases theta'
587        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
588        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
589        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
590       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
591    
592        diagName  = 'oceFreez'        diagName  = 'oceFreez'
593        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
594        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
595        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
596        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
597       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
598    
599        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
600        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
601        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
602        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
603        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
604       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
605    
606        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
607        diagTitle = 'surface salinity relaxation, >0 increases salt'        diagTitle = 'surface salinity relaxation, >0 increases salt'
608        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
609        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
610        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
611       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
612    
613        diagName  = 'surForcT'        diagName  = 'surForcT'
614        diagTitle = 'model surface forcing for Temperature,'        diagTitle = 'model surface forcing for Temperature,'
615       &          //' >0 increases theta'       &          //' >0 increases theta'
616        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
617        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
618        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
619       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
620    
621        diagName  = 'surForcS'        diagName  = 'surForcS'
622        diagTitle = 'model surface forcing for Salinity,'        diagTitle = 'model surface forcing for Salinity,'
623       &          //' >0 increases salinity'       &          //' >0 increases salinity'
624        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
625        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
626        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
627       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
628    
629        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
630        diagTitle = 'total heat flux (match heat-content variations),'        diagTitle = 'total heat flux (match heat-content variations),'
631       &          //' >0 increases theta'       &          //' >0 increases theta'
632        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
633        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
634        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
635       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
636    
637        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
638        diagTitle = 'total salt flux (match salt-content variations),'        diagTitle = 'total salt flux (match salt-content variations),'
639       &          //' >0 increases salt'       &          //' >0 increases salt'
640        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
641        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
642        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
643       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
644    
645  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
646    
# Line 625  c     diagTitle = 'Cell-Center r-Positio Line 649  c     diagTitle = 'Cell-Center r-Positio
649        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
650        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
651        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
652        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
653       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
654    
655        diagName  = 'RSURF   '        diagName  = 'RSURF   '
656  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
657        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
658        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
659        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
660        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
661       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
662    
663        diagName  = 'TOTUTEND'        diagName  = 'TOTUTEND'
664        diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
665        diagUnits = 'm/s^2           '        diagUnits = 'm/s^2           '
666        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
667        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
668       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
669         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
670    
671        diagName  = 'TOTVTEND'        diagName  = 'TOTVTEND'
672        diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
673        diagUnits = 'm/s^2           '        diagUnits = 'm/s^2           '
674        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
675        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
676       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
677         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
678    
679        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
680        diagTitle = 'Tendency of Potential Temperature'        diagTitle = 'Tendency of Potential Temperature'
681        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
682        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
683        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
684       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
685    
686        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
687        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
688        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
689        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
690        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
691       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
692    
693    
694        RETURN        RETURN
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22