/[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.6 by jmc, Sat Dec 18 23:53:04 2004 UTC revision 1.33 by jmc, Tue Apr 22 22:01:01 2008 UTC
# Line 24  C       \begin{tabular}[h]{|c|c|}\hline Line 24  C       \begin{tabular}[h]{|c|c|}\hline
24  C         \textbf{Positions}  &  \textbf{Characters}  C         \textbf{Positions}  &  \textbf{Characters}
25  C         &  \textbf{Meanings} \\\hline  C         &  \textbf{Meanings} \\\hline
26  C         parse(1)  &  S  &  scalar \\  C         parse(1)  &  S  &  scalar \\
27  C                   &  U  &  U-vector component  \\  C                   &  U  &  vector component in X direction \\
28  C                   &  V  &  V-vector component  \\  C                   &  V  &  vector component in Y direction \\
29    C                   &  W  &  vector component in vertical direction \\
30  C         parse(2)  &  U  &  C-grid U-Point  \\  C         parse(2)  &  U  &  C-grid U-Point  \\
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  \\
# Line 43  C                   &  1  &  levels = 1 Line 44  C                   &  1  &  levels = 1
44  C                   &  R  &  levels = Nr  \\  C                   &  R  &  levels = Nr  \\
45  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\
46  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
47    C                   &  G  &  levels = Ground_level Number \\
48    C                   &  I  &  levels = sea-Ice_level Number \\
49  C       \end{tabular}  C       \end{tabular}
50  C     \end{center}  C     \end{center}
51    
# Line 50  C     !USES: Line 53  C     !USES:
53        IMPLICIT NONE        IMPLICIT NONE
54  #include "SIZE.h"  #include "SIZE.h"
55  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 c #include "EESUPPORT.h"  
56  #include "PARAMS.h"  #include "PARAMS.h"
57    
58  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 58  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
75          CHARACTER*4    tUnit4c, sUnit4c
76          CHARACTER*(10) rTitle, eTitle, fTitle
77          CHARACTER*(20) pTitle, sTitle
78    
79          CHARACTER*(16) DIAGS_MK_UNITS
80          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
91            rUnit2c= 'Pa'
92            rTitle = ' Pressure '
93            pTitle = ' Geopotential       '
94          ELSE
95            rUnit2c= 'm '
96            rTitle = ' Height   '
97            pTitle = 'Pressure Pot.(p/rho)'
98          ENDIF
99          IF ( fluidIsAir ) THEN
100            tUnit4c= 'K   '
101            sUnit4c= 'g/kg'
102            sTitle = ' Specific Humidity  '
103    #ifdef ALLOW_FIZHI
104            IF (useFIZHI) sUnit4c= 'kg/kg'
105    #endif /* ALLOW_FIZHI */
106          ELSE
107            tUnit4c= 'degC'
108            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
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) (hPa,m)'        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
129        diagUnits = 'm or hPa        '  c     IF ( fluidIsWater .AND. usingZCoords )
130    c    &diagTitle = 'Sea Surface Elevation'
131          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 (hPa^2,m^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
138        diagUnits = 'm2 or hPa2      '       I                          //' Anomaly', myThid )
139        diagCode  = 'SM      M1      '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
140        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'SM P    M1      '
141       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
142         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        diagUnits = 'm2 or Pa2 /sec2 '       I                          //' Anomaly Tendency', myThid )
147          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 (K)'        diagTitle = 'Potential Temperature'
154        diagUnits = 'K               '        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
155        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
156        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
157       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
158    
159    c     diagName  = 'SST     '
160    c     diagTitle = 'Sea Surface Temperature (degC,K)'
161    c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
162    c     diagCode  = 'SM      M1      '
163    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
164    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
165    
166        diagName  = 'SALT    '        diagName  = 'SALT    '
167        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
168        diagUnits = 'g/kg            '        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
169          diagCode  = 'SM      MR      '
170          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
171         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
172    
173          diagName  = 'RELHUM  '
174          diagTitle = 'Relative Humidity'
175          diagUnits = 'percent         '
176          diagCode  = 'SM      MR      '
177          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
178         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
179    
180    c     diagName  = 'SSS     '
181    c     diagTitle = 'Sea Surface Salinity '
182    c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
183    c     diagCode  = 'SM      M1      '
184    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
185    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
186    
187          IF ( fluidIsWater ) THEN
188          diagName  = 'SALTanom'
189          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
190          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/sec           '        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/sec           '        diagUnits = 'm/s             '
207        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
208        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
209       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
210         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
211    
212        diagName  = 'WVEL    '        diagName  = 'WVEL    '
213        diagTitle = 'Vertical Component of Velocity (m/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
214        diagUnits = 'm/sec           '        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 = 'K2              '        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 = 'g2/kg2          '        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'
235          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
236          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
237          diagCode  = 'SM P    MR      '
238          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
239         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 = 'm2/sec2         '        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 = 'm2/sec2         '        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 = 'm2/sec2         '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
261        diagCode  = 'WM      LR      '        diagCode  = 'WM P    LR      '
262        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
263       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
264    
265          diagName  = 'UV_VEL_C'
266          diagTitle ='Product of horizontal Comp of velocity (cell center)'
267          diagUnits = 'm^2/s^2         '
268          diagCode  = 'UM      MR      '
269          diagMate  = diagNum + 1
270          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
271         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
272    
273        diagName  = 'UVELVVEL'        diagName  = 'UV_VEL_Z'
274        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
275        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
276        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        diagCode  = 'UZ      MR      '
277        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 1
278       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
279         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
280    
281          diagName  = 'WU_VEL  '
282          diagTitle = 'Vertical Transport of Zonal Momentum'
283          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
284          diagCode  = 'WU      LR      '
285          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
286         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
287    
288          diagName  = 'WV_VEL  '
289          diagTitle ='Vertical Transport of Meridional Momentum'
290          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
291          diagCode  = 'WV      LR      '
292          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
293         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
294    
295        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
296        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
297        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
298        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
299        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
300       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
301         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
302    
303        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
304        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
305        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
306        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
307        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
308       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
309         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
310    
311        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
312        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
313        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
314        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
315        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
316       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
317    
318        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
319        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
320        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
321        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
322        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
323       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
324         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
325    
326        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
327        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
328        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
329        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
330        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
331       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
332         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
333    
334        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
335        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
336        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
337        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
338        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
339       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
340    
341        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
342        diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
343        diagUnits = 'm-kg/s-kg       '       I                           //sTitle, myThid )
344        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
345        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'UU      MR      '
346       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum + 2
347          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
348         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
349    
350        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
351        diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
352        diagUnits = 'm-kg/s-kg       '       I                           //sTitle, myThid )
353        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
354        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      MR      '
355       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
356          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
357         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
358    
359        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
360        diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
361        diagUnits = 'm-kg/s-kg       '       I                           //sTitle, myThid )
362        diagCode  = 'WM      MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
363        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'WM      LR      '
364       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
365         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
366    
367        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
368        diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Transport of Pot Temp'
369        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
370        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UU      MR      '
371        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
372       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
373         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
374    
375        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
376        diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Transport of Pot Temp'
377        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
378        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VV      MR      '
379        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
380       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
381         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
382    
383        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
384        diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Transport of Pot Temp'
385        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
386        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
387        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
388       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
389    
390        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
391        diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
392        diagUnits = 'm-kg/s-kg       '       I                          //sTitle, myThid )
393        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
394        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'UU      MR      '
395       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum + 2
396          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
397         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
398    
399        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
400        diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
401        diagUnits = 'm-kg/s-kg       '       I                          //sTitle, myThid )
402        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
403        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      MR      '
404       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
405          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
406         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
407    
408        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
409        diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
410        diagUnits = 'm-kg/s-kg       '       I                          //sTitle, myThid )
411        diagCode  = 'WM      MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
412        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'WM      LR      '
413       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
414         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
415    
416          diagName  = 'UVELPHI '
417          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
418         I                 //pTitle//' Anomaly', myThid )
419          diagUnits = 'm^3/s^3         '
420          diagCode  = 'UU      MR      '
421          diagMate  = diagNum + 2
422          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
423         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
424    
425          diagName  = 'VVELPHI '
426          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
427         I                 //pTitle//' Anomaly', myThid )
428          diagUnits = 'm^3/s^3         '
429          diagCode  = 'VV      MR      '
430          diagMate  = diagNum
431          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
432         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
433    
434  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
435    
436          diagName  = 'RHOAnoma'
437          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
438          diagUnits = 'kg/m^3          '
439          diagCode  = 'SM      MR      '
440          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
441         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
442    
443          diagName  = 'RHOANOSQ'
444          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
445          diagUnits = 'kg^2/m^6          '
446          diagCode  = 'SM P    MR      '
447          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
448         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
449    
450          diagName  = 'URHOMASS'
451          diagTitle = 'Zonal Transport of Density'
452          diagUnits = 'kg/m^2/s        '
453          diagCode  = 'UU      MR      '
454          diagMate  = diagNum + 2
455          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
456         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
457    
458          diagName  = 'VRHOMASS'
459          diagTitle = 'Meridional Transport of Density'
460          diagUnits = 'kg/m^2/s        '
461          diagCode  = 'VV      MR      '
462          diagMate  = diagNum
463          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
464         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
465    
466          diagName  = 'WRHOMASS'
467          diagTitle = 'Vertical Transport of Potential Density'
468          diagUnits = 'kg/m^2/s        '
469          diagCode  = 'WM      LR      '
470          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
471         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
472    
473        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
474        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
475        diagUnits = 'm2/s2           '       I                           //pTitle//' Anomaly', myThid )
476          diagUnits = 'm^2/s^2         '
477        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
478        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
479       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
480    
481          diagName  = 'PHIHYDSQ'
482          diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
483         I                           //pTitle//' Anomaly', myThid )
484          diagUnits = 'm^4/s^4         '
485          diagCode  = 'SM P    MR      '
486          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
487         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
488    
489        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
490        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
491        diagUnits = 'm2/s2           '        diagTitle = DIAGS_MK_TITLE( fTitle
492         I                           //pTitle//' Anomaly', myThid )
493          diagUnits = 'm^2/s^2         '
494        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
495        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
496       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
497    
498        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
499        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
500        diagUnits = 'm4/s4           '        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
501         I                           //pTitle//' Anomaly', myThid )
502          diagUnits = 'm^4/s^4         '
503          diagCode  = 'SM P    M1      '
504          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
505         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
506    
507    #ifdef ALLOW_NONHYDROSTATIC
508          diagName  = 'PHI_NH  '
509          diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
510          diagUnits = 'm^2/s^2         '
511          diagCode  = 'SM      MR      '
512          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
513         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
514    #endif /* ALLOW_NONHYDROSTATIC */
515    
516          diagName  = 'MXLDEPTH'
517          diagTitle = 'Mixed-Layer Depth (>0)'
518          diagUnits = 'm               '
519        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
520        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
521       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
522    
523        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
524        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
525        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
526          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
527        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
528        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
529       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 = 'm4/sec          '  
       diagCode  = 'SM      MR      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
530    
531        diagName  = 'VISCAH  '        diagName  = 'CONVADJ '
532        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '        diagTitle = 'Convective Adjustment Index [0-1] '
533        diagUnits = 'm2/sec          '        diagUnits = 'fraction        '
534        diagCode  = 'SM      MR      '        diagCode  = 'SM      LR      '
535        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
536       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
537    
538  C--   surface fluxes:  C--   surface fluxes:
539        diagName  = 'TAUX    '        diagName  = 'oceTAUX '
540        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel'
541        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
542        diagCode  = 'SU      U1      '        diagCode  = 'UU      U1      '
543        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
544       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
545         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
546        diagName  = 'TAUY    '  
547        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagName  = 'oceTAUY '
548        diagUnits = 'N/m2            '        diagTitle = 'meridional surf. wind stress, >0 increases vVel'
549        diagCode  = 'SV      U1      '        diagUnits = 'N/m^2           '
550        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      U1      '
551       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
552          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
553         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
554    
555          diagName  = 'atmPload'
556          diagTitle = 'Atmospheric pressure loading'
557          diagUnits = 'Pa              '
558          diagCode  = 'SM      U1      '
559          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
560         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
561    
562        diagName  = 'TFLUX   '        diagName  = 'sIceLoad'
563        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
564        diagUnits = 'W/m2            '        diagUnits = 'kg/m^2          '
565        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
566        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
567       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
568    
569          diagName  = 'oceFWflx'
570          diagTitle = 'net surface Fresh-Water flux into the ocean'
571         &          //' (+=down), >0 decreases salinity'
572          diagUnits = 'kg/m^2/s        '
573          diagCode  = 'SM      U1      '
574          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
575         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
576    
577          diagName  = 'oceSflux'
578          diagTitle = 'net surface Salt flux into the ocean (+=down),'
579         &          //' >0 increases salinity'
580          diagUnits = 'g/m^2/s         '
581          diagCode  = 'SM      U1      '
582          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
583         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
584    
585          diagName  = 'oceQnet '
586          diagTitle = 'net surface heat flux into the ocean (+=down),'
587         &          //' >0 increases theta'
588          diagUnits = 'W/m^2           '
589          diagCode  = 'SM      U1      '
590          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
591         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
592    
593          diagName  = 'oceQsw  '
594          diagTitle = 'net Short-Wave radiation (+=down),'
595         &          //' >0 increases theta'
596          diagUnits = 'W/m^2           '
597          diagCode  = 'SM      U1      '
598          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
599         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
600    
601          diagName  = 'oceFreez'
602          diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
603          diagUnits = 'W/m^2           '
604          diagCode  = 'SM      U1      '
605          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
606         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
607    
608        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
609        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
610        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
611        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
612        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
613       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
614    
615        diagName  = 'TICE    '        diagName  = 'SRELAX  '
616        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'surface salinity relaxation, >0 increases salt'
617        diagUnits = 'W/m2            '        diagUnits = 'g/m^2/s         '
618        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
619        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
620       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
621    
622        diagName  = 'SFLUX   '        diagName  = 'surForcT'
623        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'model surface forcing for Temperature,'
624        diagUnits = 'g/m2/s          '       &          //' >0 increases theta'
625          diagUnits = 'W/m^2           '
626        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
627        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
628       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
629    
630        diagName  = 'SRELAX  '        diagName  = 'surForcS'
631        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'model surface forcing for Salinity,'
632        diagUnits = 'g/m2/s          '       &          //' >0 increases salinity'
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  = 'TFLUX   '
639          diagTitle = 'total heat flux (match heat-content variations),'
640         &          //' >0 increases theta'
641          diagUnits = 'W/m^2           '
642        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
643        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
644       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
645    
646          diagName  = 'SFLUX   '
647          diagTitle = 'total salt flux (match salt-content variations),'
648         &          //' >0 increases salt'
649          diagUnits = 'g/m^2/s         '
650          diagCode  = 'SM      U1      '
651          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
652         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
653    
654  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
655    
656          diagName  = 'RCENTER '
657    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
658          diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
659          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
660          diagCode  = 'SM      MR      '
661          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
662         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
663    
664          diagName  = 'RSURF   '
665    c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
666          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
667          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
668          diagCode  = 'SM      M1      '
669          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
670         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
671    
672          diagName  = 'TOTUTEND'
673          diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
674          diagUnits = 'm/s^2           '
675          diagCode  = 'UU      MR      '
676          diagMate  = diagNum + 2
677          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
678         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
679    
680          diagName  = 'TOTVTEND'
681          diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
682          diagUnits = 'm/s^2           '
683          diagCode  = 'VV      MR      '
684          diagMate  = diagNum
685          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
686         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
687    
688          diagName  = 'TOTTTEND'
689          diagTitle = 'Tendency of Potential Temperature'
690          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
691          diagCode  = 'SM      MR      '
692          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
693         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
694    
695          diagName  = 'TOTSTEND'
696          diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
697          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
698          diagCode  = 'SM      MR      '
699          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
700         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
701    
702    
703        RETURN        RETURN
704        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22