/[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.34 by jmc, Tue May 27 00:08:48 2008 UTC revision 1.50 by jmc, Sat Jul 19 16:46:25 2014 UTC
# Line 31  C         parse(2)  &  U  &  C-grid U-Po Line 31  C         parse(2)  &  U  &  C-grid U-Po
31  C                   &  V  &  C-grid V-Point  \\  C                   &  V  &  C-grid V-Point  \\
32  C                   &  M  &  C-grid Mass Point  \\  C                   &  M  &  C-grid Mass Point  \\
33  C                   &  Z  &  C-grid Corner Point  \\  C                   &  Z  &  C-grid Corner Point  \\
34  C         parse(3)  &  --  &  --Unused--  \\  C         parse(3)  &     &  Used for Level Integrated output: cumulate levels \\
35    C                   &  r  &  same but cumulate product by model level thickness \\
36    C                   &  R  &  same but cumulate product by hFac & level thickness \\
37  C         parse(4)  &  P  &  positive definite  \\  C         parse(4)  &  P  &  positive definite  \\
38  C         parse(5 ) &  C  &  with counter array  \\  C         parse(5 ) &  C  &  with counter array  \\
39    C                   &  P  &  post-processed (not filled up) from other diags  \\
40  C                   &  D  &  disable an array for output  \\  C                   &  D  &  disable an array for output  \\
41  C         parse(6--8) & '123'  &  3-digit mate number \\  C         parse(6--8) & '123'  &  retired, formerly: 3-digit mate number \\
42  C         parse(9)  &  U  &  model-level plus 1/2  \\  C         parse(9)  &  U  &  model-level plus 1/2  \\
43  C                   &  M  &  model-level middle  \\  C                   &  M  &  model-level middle  \\
44  C                   &  L  &  model-level minus 1/2  \\  C                   &  L  &  model-level minus 1/2  \\
# Line 46  C                   &  L  &  levels = MA Line 49  C                   &  L  &  levels = MA
49  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
50  C                   &  G  &  levels = Ground_level Number \\  C                   &  G  &  levels = Ground_level Number \\
51  C                   &  I  &  levels = sea-Ice_level Number \\  C                   &  I  &  levels = sea-Ice_level Number \\
52    C                   &  X  &  free levels option (need to be set explicitly) \\
53  C       \end{tabular}  C       \end{tabular}
54  C     \end{center}  C     \end{center}
55    
# Line 72  C     sTitle     :: "salt" title Line 76  C     sTitle     :: "salt" title
76        CHARACTER*16   diagUnits        CHARACTER*16   diagUnits
77        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
78        CHARACTER*2    rUnit2c        CHARACTER*2    rUnit2c
79        CHARACTER*4    tUnit4c, sUnit4c        CHARACTER*4    tUnit4c
80          CHARACTER*5    sUnit5c
81        CHARACTER*(10) rTitle, eTitle, fTitle        CHARACTER*(10) rTitle, eTitle, fTitle
82        CHARACTER*(20) pTitle, sTitle        CHARACTER*(20) pTitle, sTitle
83    
# Line 98  C--------------------------------------- Line 103  C---------------------------------------
103        ENDIF        ENDIF
104        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
105          tUnit4c= 'K   '          tUnit4c= 'K   '
106          sUnit4c= 'g/kg'          sUnit5c= 'kg/kg'
107          sTitle = ' Specific Humidity  '          sTitle = ' Specific Humidity  '
108  #ifdef ALLOW_FIZHI          IF (useAIM) sUnit5c= 'g/kg '
109          IF (useFIZHI) sUnit4c= 'kg/kg'        ELSEIF ( eosType.EQ.'TEOS10' ) THEN
110  #endif /* ALLOW_FIZHI */          tUnit4c= 'degC'
111            sUnit5c= 'g/kg '
112    c       tTitle = 'Conservative Temp.   '
113            sTitle = ' Absolute Salinity  '
114        ELSE        ELSE
115          tUnit4c= 'degC'          tUnit4c= 'degC'
116          sUnit4c= 'psu '          sUnit5c= 'psu  '
117    c       tTitle = 'Potential Temperature'
118          sTitle = ' Salinity           '          sTitle = ' Salinity           '
119        ENDIF        ENDIF
120  C-    free-surface (eTitle) and fixed-boundary (fTitle) position:  C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
# Line 123  C-    free-surface (eTitle) and fixed-bo Line 132  C-    free-surface (eTitle) and fixed-bo
132  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
133  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
134    
   
135        diagName  = 'ETAN    '        diagName  = 'ETAN    '
136        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
137  c     IF ( fluidIsWater .AND. usingZCoords )  c     IF ( fluidIsWater .AND. usingZCoords )
# Line 152  c    &diagTitle = 'Sea Surface Elevation Line 160  c    &diagTitle = 'Sea Surface Elevation
160        diagName  = 'THETA   '        diagName  = 'THETA   '
161        diagTitle = 'Potential Temperature'        diagTitle = 'Potential Temperature'
162        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
163        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
164        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
165       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
166    
# Line 165  c    I          diagName, diagCode, diag Line 173  c    I          diagName, diagCode, diag
173    
174        diagName  = 'SALT    '        diagName  = 'SALT    '
175        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
176        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
177        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
178        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
179       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
180    
181        diagName  = 'RELHUM  '        diagName  = 'RELHUM  '
182        diagTitle = 'Relative Humidity'        diagTitle = 'Relative Humidity'
183        diagUnits = 'percent         '        diagUnits = 'percent         '
184        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
185        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
186       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
187    
188  c     diagName  = 'SSS     '  c     diagName  = 'SSS     '
189  c     diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
190  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
191  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
192  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
193  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 187  c    I          diagName, diagCode, diag Line 195  c    I          diagName, diagCode, diag
195        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
196        diagName  = 'SALTanom'        diagName  = 'SALTanom'
197        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
198        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
199        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
200        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
201       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
202        ENDIF        ENDIF
# Line 196  c    I          diagName, diagCode, diag Line 204  c    I          diagName, diagCode, diag
204        diagName  = 'UVEL    '        diagName  = 'UVEL    '
205        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
206        diagUnits = 'm/s             '        diagUnits = 'm/s             '
207        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
208        diagMate  = diagNum + 2        diagMate  = diagNum + 2
209        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
210       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 204  c    I          diagName, diagCode, diag Line 212  c    I          diagName, diagCode, diag
212        diagName  = 'VVEL    '        diagName  = 'VVEL    '
213        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
214        diagUnits = 'm/s             '        diagUnits = 'm/s             '
215        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
216        diagMate  = diagNum        diagMate  = diagNum
217        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
218       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 219  c    I          diagName, diagCode, diag Line 227  c    I          diagName, diagCode, diag
227        diagName  = 'THETASQ '        diagName  = 'THETASQ '
228        diagTitle = 'Square of Potential Temperature'        diagTitle = 'Square of Potential Temperature'
229        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
230        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
231        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
232       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
233    
234        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
235        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
236        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
237        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
238        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
239       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
240    
241        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
242        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
243        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
244        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
245        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
246        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
247       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
248        ENDIF        ENDIF
# Line 242  c    I          diagName, diagCode, diag Line 250  c    I          diagName, diagCode, diag
250        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
251        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
252        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
253        diagCode  = 'UU P    MR      '        diagCode  = 'UURP    MR      '
254        diagMate  = diagNum + 2        diagMate  = diagNum + 2
255        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
256       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 250  c    I          diagName, diagCode, diag Line 258  c    I          diagName, diagCode, diag
258        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
259        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
260        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
261        diagCode  = 'VV P    MR      '        diagCode  = 'VVRP    MR      '
262        diagMate  = diagNum        diagMate  = diagNum
263        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
264       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 265  c    I          diagName, diagCode, diag Line 273  c    I          diagName, diagCode, diag
273        diagName  = 'UE_VEL_C'        diagName  = 'UE_VEL_C'
274        diagTitle = 'Eastward Velocity (m/s) (cell center)'        diagTitle = 'Eastward Velocity (m/s) (cell center)'
275        diagUnits = 'm/s             '        diagUnits = 'm/s             '
276        diagCode  = 'UM      MR      '        diagCode  = 'UMR     MR      '
277        diagMate  = diagNum + 2        diagMate  = diagNum + 2
278        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
279       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 273  c    I          diagName, diagCode, diag Line 281  c    I          diagName, diagCode, diag
281        diagName  = 'VN_VEL_C'        diagName  = 'VN_VEL_C'
282        diagTitle = 'Northward Velocity (m/s) (cell center)'        diagTitle = 'Northward Velocity (m/s) (cell center)'
283        diagUnits = 'm/s             '        diagUnits = 'm/s             '
284        diagCode  = 'VM      MR      '        diagCode  = 'VMR     MR      '
285        diagMate  = diagNum        diagMate  = diagNum
286        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
287       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 281  c    I          diagName, diagCode, diag Line 289  c    I          diagName, diagCode, diag
289        diagName  = 'UV_VEL_C'        diagName  = 'UV_VEL_C'
290        diagTitle ='Product of horizontal Comp of velocity (cell center)'        diagTitle ='Product of horizontal Comp of velocity (cell center)'
291        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
292        diagCode  = 'UM      MR      '        diagCode  = 'UMR     MR      '
293        diagMate  = diagNum + 1        diagMate  = diagNum + 1
294        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
295       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 289  c    I          diagName, diagCode, diag Line 297  c    I          diagName, diagCode, diag
297        diagName  = 'UV_VEL_Z'        diagName  = 'UV_VEL_Z'
298        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
299        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
300        diagCode  = 'UZ      MR      '        diagCode  = 'UZR     MR      '
301        diagMate  = diagNum + 1        diagMate  = diagNum + 1
302        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
303       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 311  c    I          diagName, diagCode, diag Line 319  c    I          diagName, diagCode, diag
319        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
320        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
321        diagUnits = 'm/s             '        diagUnits = 'm/s             '
322        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
323        diagMate  = diagNum + 2        diagMate  = diagNum + 2
324        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
325       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 319  c    I          diagName, diagCode, diag Line 327  c    I          diagName, diagCode, diag
327        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
328        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
329        diagUnits = 'm/s             '        diagUnits = 'm/s             '
330        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
331        diagMate  = diagNum        diagMate  = diagNum
332        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
333       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 331  c    I          diagName, diagCode, diag Line 339  c    I          diagName, diagCode, diag
339        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
340       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
341    
342          diagName  = 'PhiVEL  '
343          diagTitle = 'Horizontal Velocity Potential (m^2/s)'
344          diagUnits = 'm^2/s           '
345          diagCode  = 'SMR P   MR      '
346    C-    use 'UVELMASS' as mate.
347          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
348         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
349    
350          diagName  = 'PsiVEL  '
351          diagTitle = 'Horizontal Velocity Stream-Function'
352          diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
353          diagCode  = 'SZ  P   MR      '
354    C-    use 'PhiVEL' as mate.
355          diagMate  = diagNum
356          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
357         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
358    
359        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
360        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
361        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
362        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
363        diagMate  = diagNum + 2        diagMate  = diagNum + 2
364        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
365       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 342  c    I          diagName, diagCode, diag Line 367  c    I          diagName, diagCode, diag
367        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
368        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
369        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
370        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
371        diagMate  = diagNum        diagMate  = diagNum
372        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
373       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 357  c    I          diagName, diagCode, diag Line 382  c    I          diagName, diagCode, diag
382        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
383        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
384       I                           //sTitle, myThid )       I                           //sTitle, myThid )
385        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
386        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
387        diagMate  = diagNum + 2        diagMate  = diagNum + 2
388        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
389       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 366  c    I          diagName, diagCode, diag Line 391  c    I          diagName, diagCode, diag
391        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
392        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
393       I                           //sTitle, myThid )       I                           //sTitle, myThid )
394        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
395        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
396        diagMate  = diagNum        diagMate  = diagNum
397        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
398       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 375  c    I          diagName, diagCode, diag Line 400  c    I          diagName, diagCode, diag
400        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
401        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
402       I                           //sTitle, myThid )       I                           //sTitle, myThid )
403        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
404        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
405        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
406       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 383  c    I          diagName, diagCode, diag Line 408  c    I          diagName, diagCode, diag
408        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
409        diagTitle = 'Zonal Transport of Pot Temp'        diagTitle = 'Zonal Transport of Pot Temp'
410        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
411        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
412        diagMate  = diagNum + 2        diagMate  = diagNum + 2
413        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
414       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 391  c    I          diagName, diagCode, diag Line 416  c    I          diagName, diagCode, diag
416        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
417        diagTitle = 'Meridional Transport of Pot Temp'        diagTitle = 'Meridional Transport of Pot Temp'
418        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
419        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
420        diagMate  = diagNum        diagMate  = diagNum
421        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
422       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 406  c    I          diagName, diagCode, diag Line 431  c    I          diagName, diagCode, diag
431        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
432        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
433       I                          //sTitle, myThid )       I                          //sTitle, myThid )
434        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
435        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
436        diagMate  = diagNum + 2        diagMate  = diagNum + 2
437        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
438       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 415  c    I          diagName, diagCode, diag Line 440  c    I          diagName, diagCode, diag
440        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
441        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
442       I                          //sTitle, myThid )       I                          //sTitle, myThid )
443        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
444        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
445        diagMate  = diagNum        diagMate  = diagNum
446        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
447       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 424  c    I          diagName, diagCode, diag Line 449  c    I          diagName, diagCode, diag
449        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
450        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
451       I                          //sTitle, myThid )       I                          //sTitle, myThid )
452        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
453        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
454        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
455       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 433  c    I          diagName, diagCode, diag Line 458  c    I          diagName, diagCode, diag
458        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
459       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
460        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
461        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
462        diagMate  = diagNum + 2        diagMate  = diagNum + 2
463        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
464       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
465    
466        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
467        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Merid. Mass-Weight Transp of '
468       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
469        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
470        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
471        diagMate  = diagNum        diagMate  = diagNum
472        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
473       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 452  C---+----1----+----2----+----3----+----4 Line 477  C---+----1----+----2----+----3----+----4
477        diagName  = 'RHOAnoma'        diagName  = 'RHOAnoma'
478        diagTitle = 'Density Anomaly (=Rho-rhoConst)'        diagTitle = 'Density Anomaly (=Rho-rhoConst)'
479        diagUnits = 'kg/m^3          '        diagUnits = 'kg/m^3          '
480        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
481        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
482       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
483    
484        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
485        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
486        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6        '
487        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
488        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
489       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
490    
491        diagName  = 'URHOMASS'        diagName  = 'URHOMASS'
492        diagTitle = 'Zonal Transport of Density'        diagTitle = 'Zonal Transport of Density'
493        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
494        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
495        diagMate  = diagNum + 2        diagMate  = diagNum + 2
496        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
497       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 474  C---+----1----+----2----+----3----+----4 Line 499  C---+----1----+----2----+----3----+----4
499        diagName  = 'VRHOMASS'        diagName  = 'VRHOMASS'
500        diagTitle = 'Meridional Transport of Density'        diagTitle = 'Meridional Transport of Density'
501        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
502        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
503        diagMate  = diagNum        diagMate  = diagNum
504        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
505       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
506    
507        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
508        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Density'
509          diagUnits = 'kg/m^2/s        '
510          diagCode  = 'WM      LR      '
511          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
512         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
513    
514          diagName  = 'WdRHO_P '
515          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
516          diagUnits = 'kg/m^2/s        '
517          diagCode  = 'WM      LR      '
518          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
519         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
520    
521          diagName  = 'WdRHOdP '
522          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
523        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
524        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
525        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 490  C---+----1----+----2----+----3----+----4 Line 529  C---+----1----+----2----+----3----+----4
529        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
530       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
531        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
532        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
533        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
534       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
535    
# Line 498  C---+----1----+----2----+----3----+----4 Line 537  C---+----1----+----2----+----3----+----4
537        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
538       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
539        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
540        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
541        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
542       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
543    
# Line 520  c     diagTitle = 'Square of ocean botto Line 559  c     diagTitle = 'Square of ocean botto
559        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
560       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
561    
562    #ifdef NONLIN_FRSURF
563          diagName  = 'PHIHYDcR'
564          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
565         I                       //pTitle//' Anomaly @ const r', myThid )
566          diagUnits = 'm^2/s^2         '
567          diagCode  = 'SMR     MR      '
568          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
569         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
570    #endif
571    
572  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
573        diagName  = 'PHI_NH  '        diagName  = 'PHI_NH  '
574        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
575        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
576        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
577        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
578       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
579  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
# Line 547  c     diagTitle = 'Square of ocean botto Line 596  c     diagTitle = 'Square of ocean botto
596        diagName  = 'CONVADJ '        diagName  = 'CONVADJ '
597        diagTitle = 'Convective Adjustment Index [0-1] '        diagTitle = 'Convective Adjustment Index [0-1] '
598        diagUnits = 'fraction        '        diagUnits = 'fraction        '
599        diagCode  = 'SM      LR      '        diagCode  = 'SMR     LR      '
600        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
601       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
602    
# Line 686  c     diagTitle = 'Free-Surface r-Positi Line 735  c     diagTitle = 'Free-Surface r-Positi
735       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
736    
737        diagName  = 'TOTUTEND'        diagName  = 'TOTUTEND'
738        diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Zonal Component of Velocity'
739        diagUnits = 'm/s^2           '        diagUnits = 'm/s/day         '
740        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
741        diagMate  = diagNum + 2        diagMate  = diagNum + 2
742        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
743       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
744    
745        diagName  = 'TOTVTEND'        diagName  = 'TOTVTEND'
746        diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Meridional Component of Velocity'
747        diagUnits = 'm/s^2           '        diagUnits = 'm/s/day         '
748        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
749        diagMate  = diagNum        diagMate  = diagNum
750        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
751       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
752    
753        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
754        diagTitle = 'Tendency of Potential Temperature'        diagTitle = 'Tendency of Potential Temperature'
755        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
756        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
757        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
758       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
759    
760        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
761        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
762        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/day', myThid )
763          diagCode  = 'SMR     MR      '
764          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
765         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
766    
767          diagName  = 'MoistCor'
768          diagTitle = 'Heating correction due to moist thermodynamics'
769          diagUnits = 'W/m^2           '
770        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
771        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
772       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
773    
774    #ifdef ALLOW_GENERIC_ADVDIFF
775          diagName  = 'gT_Forc '
776          diagTitle = 'Potential Temp. forcing tendency'
777          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
778          diagCode  = 'SMR     MR      '
779          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
780         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
781    
782          diagName  = 'gS_Forc '
783          diagTitle = DIAGS_MK_TITLE(
784         &            sTitle//'forcing tendency', myThid )
785          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
786          diagCode  = 'SMR     MR      '
787          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
788         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
789    
790          diagName  = 'AB_gT   '
791          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
792          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
793          diagCode  = 'SMR     MR      '
794          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
795         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
796    
797          diagName  = 'AB_gS   '
798          diagTitle = DIAGS_MK_TITLE(
799         &            sTitle//'tendency from Adams-Bashforth', myThid )
800          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
801          diagCode  = 'SMR     MR      '
802          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
803         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
804    
805          diagName  = 'gTinAB  '
806          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
807          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
808          diagCode  = 'SMR     MR      '
809          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
810         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
811    
812          diagName  = 'gSinAB  '
813          diagTitle = DIAGS_MK_TITLE(
814         &            sTitle//'tendency going in Adams-Bashforth', myThid )
815          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
816          diagCode  = 'SMR     MR      '
817          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
818         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
819    #endif /* ALLOW_GENERIC_ADVDIFF */
820    
821          diagName  = 'AB_gU   '
822          diagTitle = 'U momentum tendency from Adams-Bashforth'
823          diagUnits = 'm/s^2           '
824          diagCode  = 'UUR     MR      '
825          diagMate  = diagNum + 2
826          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
827         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
828    
829          diagName  = 'AB_gV   '
830          diagTitle = 'V momentum tendency from Adams-Bashforth'
831          diagUnits = 'm/s^2           '
832          diagCode  = 'VVR     MR      '
833          diagMate  = diagNum
834          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
835         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
836    
837    #ifdef ALLOW_NONHYDROSTATIC
838          diagName  = 'AB_gW   '
839          diagTitle = 'W momentum tendency from Adams-Bashforth'
840          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
841          diagCode  = 'WM      LR      '
842          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
843         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
844    #endif /* ALLOW_NONHYDROSTATIC */
845    
846    #ifdef ALLOW_EDDYPSI
847          diagName  = 'TAUXEDDY'
848          diagTitle = 'Zonal Eddy Stress'
849          diagUnits = 'N/m**2          '
850          diagCode  = 'UU P    MR      '
851          diagMate  = diagNum
852          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
853         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
854    
855          diagName  = 'TAUYEDDY'
856          diagTitle = 'Meridional Eddy Stress'
857          diagUnits = 'N/m**2          '
858          diagCode  = 'VV P    MR      '
859          diagMate  = diagNum
860          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
861         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
862    
863    # ifdef ALLOW_GMREDI
864          diagName  = 'UVELMEAN'
865          diagTitle = 'Zonal Component of Mean Velocity (m/s)'
866          diagUnits = 'm/s             '
867          diagCode  = 'UUR     MR      '
868          diagMate  = diagNum + 2
869          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
870         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
871    
872          diagName  = 'VVELMEAN'
873          diagTitle = 'Meridional Component of Mean Velocity (m/s)'
874          diagUnits = 'm/s             '
875          diagCode  = 'VVR     MR      '
876          diagMate  = diagNum
877          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
878         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
879    # endif /* ALLOW_GMREDI */
880    #endif /* ALLOW_EDDYPSI */
881    
882        RETURN        RETURN
883        END        END
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22