/[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.20 by molod, Tue Aug 16 21:29:38 2005 UTC revision 1.34 by jmc, Tue May 27 00:08:48 2008 UTC
# Line 60  C     !INPUT PARAMETERS: Line 60  C     !INPUT PARAMETERS:
60  CEOP  CEOP
61    
62  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
63        INTEGER       diagNum  C     rTitle     :: r-coordinate title
64        CHARACTER*8   diagName  C     eTitle     :: free-surface title
65        CHARACTER*16  diagCode  C     fTitle     :: fixed boundary title
66        CHARACTER*16  diagUnits  C     pTitle     :: "Phi"  title
67    C     sTitle     :: "salt" title
68          INTEGER        diagNum
69          INTEGER        diagMate
70          CHARACTER*8    diagName
71          CHARACTER*16   diagCode
72          CHARACTER*16   diagUnits
73        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
74        CHARACTER*2   rUnit2c        CHARACTER*2    rUnit2c
75        CHARACTER*4   tUnit4c, sUnit4c        CHARACTER*4    tUnit4c, sUnit4c
76          CHARACTER*(10) rTitle, eTitle, fTitle
77          CHARACTER*(20) pTitle, sTitle
78    
79        CHARACTER*(16) DIAGS_MK_UNITS        CHARACTER*(16) DIAGS_MK_UNITS
80        EXTERNAL DIAGS_MK_UNITS        EXTERNAL DIAGS_MK_UNITS
81          CHARACTER*(80) DIAGS_MK_TITLE
82          EXTERNAL DIAGS_MK_TITLE
83    
84  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
85  C     For each output variable,  C     For each output variable,
86  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
87  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
88  C----------------------------------------------------------------------  C----------------------------------------------------------------------
89    
90        IF ( usingPCoords ) THEN        IF ( usingPCoords ) THEN
91          rUnit2c= 'Pa'          rUnit2c= 'Pa'
92            rTitle = ' Pressure '
93            pTitle = ' Geopotential       '
94        ELSE        ELSE
95          rUnit2c= 'm '          rUnit2c= 'm '
96            rTitle = ' Height   '
97            pTitle = 'Pressure Pot.(p/rho)'
98        ENDIF        ENDIF
99        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
100          tUnit4c= 'K   '          tUnit4c= 'K   '
101          sUnit4c= 'g/kg'          sUnit4c= 'g/kg'
102            sTitle = ' Specific Humidity  '
103  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
104        IF(useFIZHI)THEN          IF (useFIZHI) sUnit4c= 'kg/kg'
         sUnit4c= 'kg/kg'  
       ENDIF  
105  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
106        ELSE        ELSE
107          tUnit4c= 'degC'          tUnit4c= 'degC'
108          sUnit4c= 'psu '          sUnit4c= 'psu '
109            sTitle = ' Salinity           '
110          ENDIF
111    C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
112          IF ( fluidIsAir ) THEN
113            eTitle = ' Surface  '
114            fTitle = ' Top      '
115          ELSEIF ( usingPCoords ) THEN
116            eTitle = ' Bottom   '
117            fTitle = ' Surface  '
118          ELSE
119            eTitle = ' Surface  '
120            fTitle = ' Bottom   '
121        ENDIF        ENDIF
122    
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 = 'Perturbation of Surface (pressure, height) (Pa,m)'        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
129    c     IF ( fluidIsWater .AND. usingZCoords )
130    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 = 'Square of Perturbation of Sfc (Pa^2,m^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
138         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 = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
146         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 (degC,K)'        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        diagName  = 'SST     '  c     diagName  = 'SST     '
160        diagTitle = 'Sea Surface Temperature (degC,K)'  c     diagTitle = 'Sea Surface Temperature (degC,K)'
161        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
162        diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
163        CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
164       I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
165    
166        diagName  = 'SALT    '        diagName  = 'SALT    '
167        diagTitle = 'Salt or Water Vapor mixing ratio'        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        diagName  = 'SSS     '  c     diagName  = 'SSS     '
181        diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
182        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
183        diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
184        CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
185       I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
186    
187          IF ( fluidIsWater ) THEN
188        diagName  = 'SALTanom'        diagName  = 'SALTanom'
189        diagTitle = 'Salt anomaly (=SALT-35)'        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
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 )
       diagName  = 'UVEL_k2 '  
       diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'  
       diagUnits = 'm/s             '  
       WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
   
       diagName  = 'VVEL_k2 '  
       diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'  
       diagUnits = 'm/s             '  
       WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, 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 (K^2)'        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 = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'        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
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
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 (m^2/s^2)'        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 (m^2/s^2)'        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 (m^2/s^2)'        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 (m/s)'        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 (K.m/s)'        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 (K.m/s)'        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 = 'Zonal Mass-Weight Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
359         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 = 'Meridional Mass-Weight Transp of Salt'        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
368         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 = 'Vertical Mass-Weight Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
377         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 Transp of Pot Temp (K.m/s)'        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 Transp of Pot Temp (K.m/s)'        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 Transp of Pot Temp (K.m/s)'        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 = 'Zonal Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
408         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 = 'Meridional Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
417         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 = 'Vertical Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
426         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 '
433          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
434         I                 //pTitle//' Anomaly', myThid )
435          diagUnits = 'm^3/s^3         '
436          diagCode  = 'UU      MR      '
437          diagMate  = diagNum + 2
438          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
439         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
440    
441          diagName  = 'VVELPHI '
442          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
443         I                 //pTitle//' Anomaly', myThid )
444          diagUnits = 'm^3/s^3         '
445          diagCode  = 'VV      MR      '
446          diagMate  = diagNum
447          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 377  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 = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
491         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 = 'Square of Hyd. (ocean) press / (atmos) geoPotential'        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
499         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        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
507          diagTitle = DIAGS_MK_TITLE( fTitle
508         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        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
517         I                           //pTitle//' Anomaly', myThid )
518        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
519          diagCode  = 'SM P    M1      '
520          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
521         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'
533          diagTitle = 'Mixed-Layer Depth (>0)'
534          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 )
   
 C--   Momentum terms:  
       diagName  = 'VISCA4  '  
       diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '  
       diagUnits = 'm^4/s           '  
       diagCode  = 'SM      MR      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
546    
547        diagName  = 'VISCAH  '        diagName  = 'CONVADJ '
548        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '        diagTitle = 'Convective Adjustment Index [0-1] '
549        diagUnits = 'm^2/s           '        diagUnits = 'fraction        '
550        diagCode  = 'SM      MR      '        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  = 'TAUX    '        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  = 'TAUY    '        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'
572          diagTitle = 'Atmospheric pressure loading'
573          diagUnits = 'Pa              '
574          diagCode  = 'SM      U1      '
575          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
576         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
577    
578        diagName  = 'TFLUX   '        diagName  = 'sIceLoad'
579        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
580          diagUnits = 'kg/m^2          '
581          diagCode  = 'SM      U1      '
582          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
583         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
584    
585          diagName  = 'oceFWflx'
586          diagTitle = 'net surface Fresh-Water flux into the ocean'
587         &          //' (+=down), >0 decreases salinity'
588          diagUnits = 'kg/m^2/s        '
589          diagCode  = 'SM      U1      '
590          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
591         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
592    
593          diagName  = 'oceSflux'
594          diagTitle = 'net surface Salt flux into the ocean (+=down),'
595         &          //' >0 increases salinity'
596          diagUnits = 'g/m^2/s         '
597          diagCode  = 'SM      U1      '
598          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
599         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
600    
601          diagName  = 'oceQnet '
602          diagTitle = 'net surface heat flux into the ocean (+=down),'
603         &          //' >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  = 'TRELAX  '        diagName  = 'oceQsw  '
610        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'net Short-Wave radiation (+=down),'
611         &          //' >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  = 'TICE    '        diagName  = 'oceFreez'
618        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        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  = 'SFLUX   '        diagName  = 'TRELAX  '
625        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'surface temperature relaxation, >0 increases theta'
626        diagUnits = 'g/m^2/s         '        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         '
634          diagCode  = 'SM      U1      '
635          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
636         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
637    
638          diagName  = 'surForcT'
639          diagTitle = 'model surface forcing for Temperature,'
640         &          //' >0 increases theta'
641          diagUnits = 'W/m^2           '
642          diagCode  = 'SM      U1      '
643          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
644         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
645    
646          diagName  = 'surForcS'
647          diagTitle = 'model surface forcing for Salinity,'
648         &          //' >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  #ifdef SHORTWAVE_HEATING        diagName  = 'TFLUX   '
655        diagName  = 'SWFLUX  '        diagTitle = 'total heat flux (match heat-content variations),'
656        diagTitle = 'net upward SW radiation, >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  #endif  
662          diagName  = 'SFLUX   '
663          diagTitle = 'total salt flux (match salt-content variations),'
664         &          //' >0 increases salt'
665          diagUnits = 'g/m^2/s         '
666          diagCode  = 'SM      U1      '
667          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
668         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        diagName  = 'PRESSURE'  
672        diagTitle = 'Atmospheric Pressure (Pa)     '        diagName  = 'RCENTER '
673    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
674          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        diagTitle = 'Surface (Pressure, Height) (Pa,m)'  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
682          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'
689          diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
690          diagUnits = 'm/s^2           '
691          diagCode  = 'UU      MR      '
692          diagMate  = diagNum + 2
693          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
694         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
695    
696          diagName  = 'TOTVTEND'
697          diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
698          diagUnits = 'm/s^2           '
699          diagCode  = 'VV      MR      '
700          diagMate  = diagNum
701          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
702         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
703    
704          diagName  = 'TOTTTEND'
705          diagTitle = 'Tendency of Potential Temperature'
706          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
707          diagCode  = 'SM      MR      '
708          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
709         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
710    
711          diagName  = 'TOTSTEND'
712          diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
713          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
714          diagCode  = 'SM      MR      '
715          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
716         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
717    
718    
719        RETURN        RETURN
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22