/[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.35 by jmc, Sat Nov 15 19:52:32 2008 UTC revision 1.47 by m_bates, Wed Sep 18 22:04:27 2013 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 152  c    &diagTitle = 'Sea Surface Elevation Line 161  c    &diagTitle = 'Sea Surface Elevation
161        diagName  = 'THETA   '        diagName  = 'THETA   '
162        diagTitle = 'Potential Temperature'        diagTitle = 'Potential Temperature'
163        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
164        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
165        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
166       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
167    
# Line 165  c    I          diagName, diagCode, diag Line 174  c    I          diagName, diagCode, diag
174    
175        diagName  = 'SALT    '        diagName  = 'SALT    '
176        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
177        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
178        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
179        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
180       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
181    
182        diagName  = 'RELHUM  '        diagName  = 'RELHUM  '
183        diagTitle = 'Relative Humidity'        diagTitle = 'Relative Humidity'
184        diagUnits = 'percent         '        diagUnits = 'percent         '
185        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
186        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
187       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
188    
189  c     diagName  = 'SSS     '  c     diagName  = 'SSS     '
190  c     diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
191  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
192  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
193  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
194  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 187  c    I          diagName, diagCode, diag Line 196  c    I          diagName, diagCode, diag
196        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
197        diagName  = 'SALTanom'        diagName  = 'SALTanom'
198        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
199        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
200        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
201        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
202       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
203        ENDIF        ENDIF
# Line 196  c    I          diagName, diagCode, diag Line 205  c    I          diagName, diagCode, diag
205        diagName  = 'UVEL    '        diagName  = 'UVEL    '
206        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
207        diagUnits = 'm/s             '        diagUnits = 'm/s             '
208        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
209        diagMate  = diagNum + 2        diagMate  = diagNum + 2
210        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
211       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 204  c    I          diagName, diagCode, diag Line 213  c    I          diagName, diagCode, diag
213        diagName  = 'VVEL    '        diagName  = 'VVEL    '
214        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
215        diagUnits = 'm/s             '        diagUnits = 'm/s             '
216        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
217        diagMate  = diagNum        diagMate  = diagNum
218        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
219       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 219  c    I          diagName, diagCode, diag Line 228  c    I          diagName, diagCode, diag
228        diagName  = 'THETASQ '        diagName  = 'THETASQ '
229        diagTitle = 'Square of Potential Temperature'        diagTitle = 'Square of Potential Temperature'
230        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
231        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
232        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
233       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
234    
235        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
236        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
237        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
238        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
239        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
240       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
241    
242        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
243        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
244        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
245        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
246        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
247        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
248       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
249        ENDIF        ENDIF
# Line 242  c    I          diagName, diagCode, diag Line 251  c    I          diagName, diagCode, diag
251        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
252        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
253        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
254        diagCode  = 'UU P    MR      '        diagCode  = 'UURP    MR      '
255        diagMate  = diagNum + 2        diagMate  = diagNum + 2
256        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
257       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 250  c    I          diagName, diagCode, diag Line 259  c    I          diagName, diagCode, diag
259        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
260        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
261        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
262        diagCode  = 'VV P    MR      '        diagCode  = 'VVRP    MR      '
263        diagMate  = diagNum        diagMate  = diagNum
264        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
265       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 265  c    I          diagName, diagCode, diag Line 274  c    I          diagName, diagCode, diag
274        diagName  = 'UE_VEL_C'        diagName  = 'UE_VEL_C'
275        diagTitle = 'Eastward Velocity (m/s) (cell center)'        diagTitle = 'Eastward Velocity (m/s) (cell center)'
276        diagUnits = 'm/s             '        diagUnits = 'm/s             '
277        diagCode  = 'UM      MR      '        diagCode  = 'UMR     MR      '
278        diagMate  = diagNum + 2        diagMate  = diagNum + 2
279        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
280       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 273  c    I          diagName, diagCode, diag Line 282  c    I          diagName, diagCode, diag
282        diagName  = 'VN_VEL_C'        diagName  = 'VN_VEL_C'
283        diagTitle = 'Northward Velocity (m/s) (cell center)'        diagTitle = 'Northward Velocity (m/s) (cell center)'
284        diagUnits = 'm/s             '        diagUnits = 'm/s             '
285        diagCode  = 'VM      MR      '        diagCode  = 'VMR     MR      '
286        diagMate  = diagNum        diagMate  = diagNum
287        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
288       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 281  c    I          diagName, diagCode, diag Line 290  c    I          diagName, diagCode, diag
290        diagName  = 'UV_VEL_C'        diagName  = 'UV_VEL_C'
291        diagTitle ='Product of horizontal Comp of velocity (cell center)'        diagTitle ='Product of horizontal Comp of velocity (cell center)'
292        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
293        diagCode  = 'UM      MR      '        diagCode  = 'UMR     MR      '
294        diagMate  = diagNum + 1        diagMate  = diagNum + 1
295        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
296       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 289  c    I          diagName, diagCode, diag Line 298  c    I          diagName, diagCode, diag
298        diagName  = 'UV_VEL_Z'        diagName  = 'UV_VEL_Z'
299        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
300        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
301        diagCode  = 'UZ      MR      '        diagCode  = 'UZR     MR      '
302        diagMate  = diagNum + 1        diagMate  = diagNum + 1
303        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
304       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 311  c    I          diagName, diagCode, diag Line 320  c    I          diagName, diagCode, diag
320        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
321        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
322        diagUnits = 'm/s             '        diagUnits = 'm/s             '
323        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
324        diagMate  = diagNum + 2        diagMate  = diagNum + 2
325        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
326       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 319  c    I          diagName, diagCode, diag Line 328  c    I          diagName, diagCode, diag
328        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
329        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
330        diagUnits = 'm/s             '        diagUnits = 'm/s             '
331        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
332        diagMate  = diagNum        diagMate  = diagNum
333        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
334       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 331  c    I          diagName, diagCode, diag Line 340  c    I          diagName, diagCode, diag
340        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
341       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
342    
343          diagName  = 'PhiVEL  '
344          diagTitle = 'Horizontal Velocity Potential (m^2/s)'
345          diagUnits = 'm^2/s           '
346          diagCode  = 'SMR P   MR      '
347    C-    use 'UVELMASS' as mate.
348          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
349         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
350    
351          diagName  = 'PsiVEL  '
352          diagTitle = 'Horizontal Velocity Stream-Function'
353          diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
354          diagCode  = 'SZ  P   MR      '
355    C-    use 'PhiVEL' as mate.
356          diagMate  = diagNum
357          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
358         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
359    
360        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
361        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
362        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
363        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
364        diagMate  = diagNum + 2        diagMate  = diagNum + 2
365        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
366       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 342  c    I          diagName, diagCode, diag Line 368  c    I          diagName, diagCode, diag
368        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
369        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
370        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
371        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
372        diagMate  = diagNum        diagMate  = diagNum
373        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
374       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 357  c    I          diagName, diagCode, diag Line 383  c    I          diagName, diagCode, diag
383        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
384        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
385       I                           //sTitle, myThid )       I                           //sTitle, myThid )
386        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
387        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
388        diagMate  = diagNum + 2        diagMate  = diagNum + 2
389        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
390       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 366  c    I          diagName, diagCode, diag Line 392  c    I          diagName, diagCode, diag
392        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
393        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
394       I                           //sTitle, myThid )       I                           //sTitle, myThid )
395        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
396        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
397        diagMate  = diagNum        diagMate  = diagNum
398        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
399       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 375  c    I          diagName, diagCode, diag Line 401  c    I          diagName, diagCode, diag
401        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
402        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
403       I                           //sTitle, myThid )       I                           //sTitle, myThid )
404        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
405        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
406        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
407       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 383  c    I          diagName, diagCode, diag Line 409  c    I          diagName, diagCode, diag
409        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
410        diagTitle = 'Zonal Transport of Pot Temp'        diagTitle = 'Zonal Transport of Pot Temp'
411        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
412        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
413        diagMate  = diagNum + 2        diagMate  = diagNum + 2
414        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
415       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 391  c    I          diagName, diagCode, diag Line 417  c    I          diagName, diagCode, diag
417        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
418        diagTitle = 'Meridional Transport of Pot Temp'        diagTitle = 'Meridional Transport of Pot Temp'
419        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
420        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
421        diagMate  = diagNum        diagMate  = diagNum
422        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
423       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 406  c    I          diagName, diagCode, diag Line 432  c    I          diagName, diagCode, diag
432        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
433        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
434       I                          //sTitle, myThid )       I                          //sTitle, myThid )
435        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
436        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
437        diagMate  = diagNum + 2        diagMate  = diagNum + 2
438        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
439       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 415  c    I          diagName, diagCode, diag Line 441  c    I          diagName, diagCode, diag
441        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
442        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
443       I                          //sTitle, myThid )       I                          //sTitle, myThid )
444        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
445        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
446        diagMate  = diagNum        diagMate  = diagNum
447        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
448       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 424  c    I          diagName, diagCode, diag Line 450  c    I          diagName, diagCode, diag
450        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
451        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
452       I                          //sTitle, myThid )       I                          //sTitle, myThid )
453        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
454        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
455        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
456       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 433  c    I          diagName, diagCode, diag Line 459  c    I          diagName, diagCode, diag
459        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
460       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
461        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
462        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
463        diagMate  = diagNum + 2        diagMate  = diagNum + 2
464        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
465       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 442  c    I          diagName, diagCode, diag Line 468  c    I          diagName, diagCode, diag
468        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
469       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
470        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
471        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
472        diagMate  = diagNum        diagMate  = diagNum
473        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
474       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 452  C---+----1----+----2----+----3----+----4 Line 478  C---+----1----+----2----+----3----+----4
478        diagName  = 'RHOAnoma'        diagName  = 'RHOAnoma'
479        diagTitle = 'Density Anomaly (=Rho-rhoConst)'        diagTitle = 'Density Anomaly (=Rho-rhoConst)'
480        diagUnits = 'kg/m^3          '        diagUnits = 'kg/m^3          '
481        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
482        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
483       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
484    
485        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
486        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
487        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6          '
488        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
489        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
490       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
491    
492        diagName  = 'URHOMASS'        diagName  = 'URHOMASS'
493        diagTitle = 'Zonal Transport of Density'        diagTitle = 'Zonal Transport of Density'
494        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
495        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
496        diagMate  = diagNum + 2        diagMate  = diagNum + 2
497        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
498       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 474  C---+----1----+----2----+----3----+----4 Line 500  C---+----1----+----2----+----3----+----4
500        diagName  = 'VRHOMASS'        diagName  = 'VRHOMASS'
501        diagTitle = 'Meridional Transport of Density'        diagTitle = 'Meridional Transport of Density'
502        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
503        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
504        diagMate  = diagNum        diagMate  = diagNum
505        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
506       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
507    
508        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
509        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Density'
510          diagUnits = 'kg/m^2/s        '
511          diagCode  = 'WM      LR      '
512          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
513         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
514    
515          diagName  = 'WdRHO_P '
516          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
517          diagUnits = 'kg/m^2/s        '
518          diagCode  = 'WM      LR      '
519          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
520         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
521    
522          diagName  = 'WdRHOdP '
523          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
524        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
525        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
526        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 490  C---+----1----+----2----+----3----+----4 Line 530  C---+----1----+----2----+----3----+----4
530        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
531       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
532        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
533        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
534        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
535       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
536    
# Line 498  C---+----1----+----2----+----3----+----4 Line 538  C---+----1----+----2----+----3----+----4
538        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
539       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
540        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
541        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
542        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
543       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
544    
# Line 520  c     diagTitle = 'Square of ocean botto Line 560  c     diagTitle = 'Square of ocean botto
560        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
561       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
562    
563    #ifdef NONLIN_FRSURF
564          diagName  = 'PHIHYDcR'
565          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
566         I                       //pTitle//' Anomaly @ const r', myThid )
567          diagUnits = 'm^2/s^2         '
568          diagCode  = 'SMR     MR      '
569          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
570         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
571    #endif
572    
573  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
574        diagName  = 'PHI_NH  '        diagName  = 'PHI_NH  '
575        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
576        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
577        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
578        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
579       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
580  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
# Line 547  c     diagTitle = 'Square of ocean botto Line 597  c     diagTitle = 'Square of ocean botto
597        diagName  = 'CONVADJ '        diagName  = 'CONVADJ '
598        diagTitle = 'Convective Adjustment Index [0-1] '        diagTitle = 'Convective Adjustment Index [0-1] '
599        diagUnits = 'fraction        '        diagUnits = 'fraction        '
600        diagCode  = 'SM      LR      '        diagCode  = 'SMR     LR      '
601        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
602       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
603    
# Line 688  c     diagTitle = 'Free-Surface r-Positi Line 738  c     diagTitle = 'Free-Surface r-Positi
738        diagName  = 'TOTUTEND'        diagName  = 'TOTUTEND'
739        diagTitle = 'Tendency of Zonal Component of Velocity'        diagTitle = 'Tendency of Zonal Component of Velocity'
740        diagUnits = 'm/s/day         '        diagUnits = 'm/s/day         '
741        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
742        diagMate  = diagNum + 2        diagMate  = diagNum + 2
743        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
744       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 696  c     diagTitle = 'Free-Surface r-Positi Line 746  c     diagTitle = 'Free-Surface r-Positi
746        diagName  = 'TOTVTEND'        diagName  = 'TOTVTEND'
747        diagTitle = 'Tendency of Meridional Component of Velocity'        diagTitle = 'Tendency of Meridional Component of Velocity'
748        diagUnits = 'm/s/day         '        diagUnits = 'm/s/day         '
749        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
750        diagMate  = diagNum        diagMate  = diagNum
751        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
752       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 704  c     diagTitle = 'Free-Surface r-Positi Line 754  c     diagTitle = 'Free-Surface r-Positi
754        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
755        diagTitle = 'Tendency of Potential Temperature'        diagTitle = 'Tendency of Potential Temperature'
756        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
757        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
758        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
759       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
760    
761        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
762        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
763        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/day', myThid )
764          diagCode  = 'SMR     MR      '
765          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
766         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
767    
768          diagName  = 'MoistCor'
769          diagTitle = 'Heating correction due to moist thermodynamics'
770          diagUnits = 'W/s^2           '
771        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
772        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
773       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
774    
775          diagName  = 'AB_gU   '
776          diagTitle = 'U momentum tendency from Adams-Bashforth'
777          diagUnits = 'm/s^2           '
778          diagCode  = 'UUR     MR      '
779          diagMate  = diagNum + 2
780          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
781         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
782    
783          diagName  = 'AB_gV   '
784          diagTitle = 'V momentum tendency from Adams-Bashforth'
785          diagUnits = 'm/s^2           '
786          diagCode  = 'VVR     MR      '
787          diagMate  = diagNum
788          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
789         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
790    
791    #ifdef ALLOW_NONHYDROSTATIC
792          diagName  = 'AB_gW   '
793          diagTitle = 'W momentum tendency from Adams-Bashforth'
794          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
795          diagCode  = 'WM      LR      '
796          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
797         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
798    #endif /* ALLOW_NONHYDROSTATIC */
799    
800    #ifdef ALLOW_GENERIC_ADVDIFF
801          diagName  = 'AB_gT   '
802          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
803          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
804          diagCode  = 'SMR     MR      '
805          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
806         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
807    
808          diagName  = 'AB_gS   '
809          diagTitle = DIAGS_MK_TITLE(
810         &            sTitle//'tendency from Adams-Bashforth', myThid )
811          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
812          diagCode  = 'SMR     MR      '
813          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
814         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
815    
816          diagName  = 'gTinAB  '
817          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
818          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
819          diagCode  = 'SMR     MR      '
820          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
821         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
822    
823          diagName  = 'gSinAB  '
824          diagTitle = DIAGS_MK_TITLE(
825         &            sTitle//'tendency going in Adams-Bashforth', myThid )
826          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
827          diagCode  = 'SMR     MR      '
828          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
829         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
830    #endif /* ALLOW_GENERIC_ADVDIFF */
831    
832    #ifdef ALLOW_EDDYPSI
833          diagName  = 'TAUXEDDY'
834          diagTitle = 'Zonal Eddy Stress'
835          diagUnits = 'N/m**2          '
836          diagCode  = 'UU P    MR      '
837          diagMate  = diagNum
838          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
839         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
840            
841          diagName  = 'TAUYEDDY'
842          diagTitle = 'Meridional Eddy Stress'
843          diagUnits = 'N/m**2          '
844          diagCode  = 'VV P    MR      '
845          diagMate  = diagNum
846          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
847         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
848    
849    #ifdef ALLOW_GMREDI
850          diagName  = 'UVELMEAN'
851          diagTitle = 'Zonal Component of Mean Velocity (m/s)'
852          diagUnits = 'm/s             '
853          diagCode  = 'UUR     MR      '
854          diagMate  = diagNum + 2
855          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
856         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
857            
858          diagName  = 'VVELMEAN'
859          diagTitle = 'Meridional Component of Mean Velocity (m/s)'
860          diagUnits = 'm/s             '
861          diagCode  = 'VVR     MR      '
862          diagMate  = diagNum
863          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
864         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
865    
866    #endif /* ALLOW_GMREDI */
867    
868    #endif /* ALLOW_EDDYPSI */
869    
870        RETURN        RETURN
871        END        END
Legend:
Removed from v.1.35  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.22