/[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.43 by jmc, Wed Nov 30 20:58:35 2011 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  \\
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 43  C                   &  1  &  levels = 1 Line 47  C                   &  1  &  levels = 1
47  C                   &  R  &  levels = Nr  \\  C                   &  R  &  levels = Nr  \\
48  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\
49  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
50    C                   &  G  &  levels = Ground_level Number \\
51    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 50  C     !USES: Line 57  C     !USES:
57        IMPLICIT NONE        IMPLICIT NONE
58  #include "SIZE.h"  #include "SIZE.h"
59  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 c #include "EESUPPORT.h"  
60  #include "PARAMS.h"  #include "PARAMS.h"
61    
62  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 58  C     !INPUT PARAMETERS: Line 64  C     !INPUT PARAMETERS:
64  CEOP  CEOP
65    
66  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
67        INTEGER       diagNum  C     rTitle     :: r-coordinate title
68        CHARACTER*8   diagName  C     eTitle     :: free-surface title
69        CHARACTER*16  diagCode  C     fTitle     :: fixed boundary title
70        CHARACTER*16  diagUnits  C     pTitle     :: "Phi"  title
71    C     sTitle     :: "salt" title
72          INTEGER        diagNum
73          INTEGER        diagMate
74          CHARACTER*8    diagName
75          CHARACTER*16   diagCode
76          CHARACTER*16   diagUnits
77        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
78          CHARACTER*2    rUnit2c
79          CHARACTER*4    tUnit4c, sUnit4c
80          CHARACTER*(10) rTitle, eTitle, fTitle
81          CHARACTER*(20) pTitle, sTitle
82    
83          CHARACTER*(16) DIAGS_MK_UNITS
84          EXTERNAL DIAGS_MK_UNITS
85          CHARACTER*(80) DIAGS_MK_TITLE
86          EXTERNAL DIAGS_MK_TITLE
87    
88  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
89  C     For each output variable,  C     For each output variable,
90  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
91  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
92  C----------------------------------------------------------------------  C----------------------------------------------------------------------
93    
94          IF ( usingPCoords ) THEN
95            rUnit2c= 'Pa'
96            rTitle = ' Pressure '
97            pTitle = ' Geopotential       '
98          ELSE
99            rUnit2c= 'm '
100            rTitle = ' Height   '
101            pTitle = 'Pressure Pot.(p/rho)'
102          ENDIF
103          IF ( fluidIsAir ) THEN
104            tUnit4c= 'K   '
105            sUnit4c= 'g/kg'
106            sTitle = ' Specific Humidity  '
107    #ifdef ALLOW_FIZHI
108            IF (useFIZHI) sUnit4c= 'kg/kg'
109    #endif /* ALLOW_FIZHI */
110          ELSE
111            tUnit4c= 'degC'
112            sUnit4c= 'psu '
113            sTitle = ' Salinity           '
114          ENDIF
115    C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
116          IF ( fluidIsAir ) THEN
117            eTitle = ' Surface  '
118            fTitle = ' Top      '
119          ELSEIF ( usingPCoords ) THEN
120            eTitle = ' Bottom   '
121            fTitle = ' Surface  '
122          ELSE
123            eTitle = ' Surface  '
124            fTitle = ' Bottom   '
125          ENDIF
126    
127  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
128  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
129    
130    
131        diagName  = 'ETAN    '        diagName  = 'ETAN    '
132        diagTitle = 'Perturbation of Surface (pressure, height) (hPa,m)'        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
133        diagUnits = 'm or hPa        '  c     IF ( fluidIsWater .AND. usingZCoords )
134    c    &diagTitle = 'Sea Surface Elevation'
135          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
136        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
137        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
138       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
139    
140        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
141        diagTitle = 'Square of Perturbation of Sfc (hPa^2,m^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
142        diagUnits = 'm2 or hPa2      '       I                          //' Anomaly', myThid )
143        diagCode  = 'SM      M1      '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
144        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'SM P    M1      '
145       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
146         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
147    
148        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
149        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
150        diagUnits = 'm2 or Pa2 /sec2 '       I                          //' Anomaly Tendency', myThid )
151          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
152        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
153        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
154       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
155    
156        diagName  = 'THETA   '        diagName  = 'THETA   '
157        diagTitle = 'Potential Temperature (K)'        diagTitle = 'Potential Temperature'
158        diagUnits = 'K               '        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
159        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
160        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
161       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
162    
163    c     diagName  = 'SST     '
164    c     diagTitle = 'Sea Surface Temperature (degC,K)'
165    c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
166    c     diagCode  = 'SM      M1      '
167    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
168    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
169    
170        diagName  = 'SALT    '        diagName  = 'SALT    '
171        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
172        diagUnits = 'g/kg            '        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
173        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
174        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
175       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
176    
177          diagName  = 'RELHUM  '
178          diagTitle = 'Relative Humidity'
179          diagUnits = 'percent         '
180          diagCode  = 'SMR     MR      '
181          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
182         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
183    
184    c     diagName  = 'SSS     '
185    c     diagTitle = 'Sea Surface Salinity '
186    c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
187    c     diagCode  = 'SM      M1      '
188    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
189    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
190    
191          IF ( fluidIsWater ) THEN
192          diagName  = 'SALTanom'
193          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
194          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
195          diagCode  = 'SMR     MR      '
196          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
197         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
198          ENDIF
199    
200        diagName  = 'UVEL    '        diagName  = 'UVEL    '
201        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
202        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
203        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UUR     MR      '
204        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
205       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
206         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
207    
208        diagName  = 'VVEL    '        diagName  = 'VVEL    '
209        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
210        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
211        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVR     MR      '
212        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
213       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
214         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
215    
216        diagName  = 'WVEL    '        diagName  = 'WVEL    '
217        diagTitle = 'Vertical Component of Velocity (m/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
218        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
219        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
220        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
221       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
222    
223        diagName  = 'THETASQ '        diagName  = 'THETASQ '
224        diagTitle = 'Square of Potential Temperature (K^2)'        diagTitle = 'Square of Potential Temperature'
225        diagUnits = 'K2              '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
226        diagCode  = 'SM      MR      '        diagCode  = 'SMRP    MR      '
227        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
229    
230        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
231        diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
232        diagUnits = 'g2/kg2          '        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233        diagCode  = 'SM      MR      '        diagCode  = 'SMRP    MR      '
234        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
235       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
236    
237          IF ( fluidIsWater ) THEN
238          diagName  = 'SALTSQan'
239          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
240          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
241          diagCode  = 'SMRP    MR      '
242          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
243         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
244          ENDIF
245    
246        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
247        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
248        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
249        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UURP    MR      '
250        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
251       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
252         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
253    
254        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
255        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
256        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
257        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVRP    MR      '
258        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
259       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
260         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
261    
262        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
263        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Vertical Comp of Velocity'
264        diagUnits = 'm2/sec2         '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
265        diagCode  = 'WM      LR      '        diagCode  = 'WM P    LR      '
266        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
267       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
268    
269          diagName  = 'UE_VEL_C'
270          diagTitle = 'Eastward Velocity (m/s) (cell center)'
271          diagUnits = 'm/s             '
272          diagCode  = 'UMR     MR      '
273          diagMate  = diagNum + 2
274          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
275         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
276    
277          diagName  = 'VN_VEL_C'
278          diagTitle = 'Northward Velocity (m/s) (cell center)'
279          diagUnits = 'm/s             '
280          diagCode  = 'VMR     MR      '
281          diagMate  = diagNum
282          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
283         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
284    
285          diagName  = 'UV_VEL_C'
286          diagTitle ='Product of horizontal Comp of velocity (cell center)'
287          diagUnits = 'm^2/s^2         '
288          diagCode  = 'UMR     MR      '
289          diagMate  = diagNum + 1
290          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
291         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
292    
293        diagName  = 'UVELVVEL'        diagName  = 'UV_VEL_Z'
294        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
295        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
296        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        diagCode  = 'UZR     MR      '
297        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 1
298       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
299         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
300    
301          diagName  = 'WU_VEL  '
302          diagTitle = 'Vertical Transport of Zonal Momentum'
303          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
304          diagCode  = 'WU      LR      '
305          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
306         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
307    
308          diagName  = 'WV_VEL  '
309          diagTitle ='Vertical Transport of Meridional Momentum'
310          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
311          diagCode  = 'WV      LR      '
312          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
313         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
314    
315        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
316        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
317        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
318        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UUr     MR      '
319        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
320       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
321         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
322    
323        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
324        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
325        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
326        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVr     MR      '
327        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
328       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
329         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
330    
331        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
332        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
333        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
334        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
335        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
336       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
337    
338          diagName  = 'PhiVEL  '
339          diagTitle = 'Horizontal Velocity Potential (m^2/s)'
340          diagUnits = 'm^2/s           '
341          diagCode  = 'SMR P   MR      '
342    C-    use 'UVELMASS' as mate.
343          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
344         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
345    
346          diagName  = 'PsiVEL  '
347          diagTitle = 'Horizontal Velocity Stream-Function'
348          diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
349          diagCode  = 'SZ  P   MR      '
350    C-    use 'PhiVEL' as mate.
351          diagMate  = diagNum
352          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
353         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
354    
355        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
356        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
357        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
358        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UUr     MR      '
359        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
360       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
361         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
362    
363        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
364        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
365        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
366        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVr     MR      '
367        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
368       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
369         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
370    
371        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
372        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
373        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
374        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
375        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
376       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
377    
378        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
379        diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
380        diagUnits = 'm-kg/s-kg       '       I                           //sTitle, myThid )
381        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
382        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'UUr     MR      '
383       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum + 2
384          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
385         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
386    
387        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
388        diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
389        diagUnits = 'm-kg/s-kg       '       I                           //sTitle, myThid )
390        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
391        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VVr     MR      '
392       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
393          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
394         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
395    
396        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
397        diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
398        diagUnits = 'm-kg/s-kg       '       I                           //sTitle, myThid )
399        diagCode  = 'WM      MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
400        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'WM      LR      '
401       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
402         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
403    
404        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
405        diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Transport of Pot Temp'
406        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
407        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UUR     MR      '
408        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
409       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
410         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
411    
412        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
413        diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Transport of Pot Temp'
414        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
415        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVR     MR      '
416        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
417       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
418         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
419    
420        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
421        diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Transport of Pot Temp'
422        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
423        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
424        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
425       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
426    
427        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
428        diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
429        diagUnits = 'm-kg/s-kg       '       I                          //sTitle, myThid )
430        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
431        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'UUR     MR      '
432       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum + 2
433          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
434         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
435    
436        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
437        diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
438        diagUnits = 'm-kg/s-kg       '       I                          //sTitle, myThid )
439        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
440        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VVR     MR      '
441       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
442          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
443         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
444    
445        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
446        diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
447        diagUnits = 'm-kg/s-kg       '       I                          //sTitle, myThid )
448        diagCode  = 'WM      MR      '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
449        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'WM      LR      '
450       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
451         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
452    
453          diagName  = 'UVELPHI '
454          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
455         I                 //pTitle//' Anomaly', myThid )
456          diagUnits = 'm^3/s^3         '
457          diagCode  = 'UUr     MR      '
458          diagMate  = diagNum + 2
459          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
460         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
461    
462          diagName  = 'VVELPHI '
463          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
464         I                 //pTitle//' Anomaly', myThid )
465          diagUnits = 'm^3/s^3         '
466          diagCode  = 'VVr     MR      '
467          diagMate  = diagNum
468          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
469         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
470    
471  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
472    
473          diagName  = 'RHOAnoma'
474          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
475          diagUnits = 'kg/m^3          '
476          diagCode  = 'SMR     MR      '
477          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
478         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
479    
480          diagName  = 'RHOANOSQ'
481          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
482          diagUnits = 'kg^2/m^6          '
483          diagCode  = 'SMRP    MR      '
484          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
485         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
486    
487          diagName  = 'URHOMASS'
488          diagTitle = 'Zonal Transport of Density'
489          diagUnits = 'kg/m^2/s        '
490          diagCode  = 'UUr     MR      '
491          diagMate  = diagNum + 2
492          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
493         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
494    
495          diagName  = 'VRHOMASS'
496          diagTitle = 'Meridional Transport of Density'
497          diagUnits = 'kg/m^2/s        '
498          diagCode  = 'VVr     MR      '
499          diagMate  = diagNum
500          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
501         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
502    
503          diagName  = 'WRHOMASS'
504          diagTitle = 'Vertical Transport of Density'
505          diagUnits = 'kg/m^2/s        '
506          diagCode  = 'WM      LR      '
507          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
508         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
509    
510          diagName  = 'WdRHO_P '
511          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
512          diagUnits = 'kg/m^2/s        '
513          diagCode  = 'WM      LR      '
514          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
515         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
516    
517          diagName  = 'WdRHOdP '
518          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
519          diagUnits = 'kg/m^2/s        '
520          diagCode  = 'WM      LR      '
521          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
522         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
523    
524        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
525        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
526        diagUnits = 'm2/s2           '       I                           //pTitle//' Anomaly', myThid )
527        diagCode  = 'SM      MR      '        diagUnits = 'm^2/s^2         '
528        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'SMR     MR      '
529       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
530         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
531    
532          diagName  = 'PHIHYDSQ'
533          diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
534         I                           //pTitle//' Anomaly', myThid )
535          diagUnits = 'm^4/s^4         '
536          diagCode  = 'SMRP    MR      '
537          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
538         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
539    
540        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
541        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
542        diagUnits = 'm2/s2           '        diagTitle = DIAGS_MK_TITLE( fTitle
543         I                           //pTitle//' Anomaly', myThid )
544          diagUnits = 'm^2/s^2         '
545        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
546        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
547       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
548    
549        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
550        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
551        diagUnits = 'm4/s4           '        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
552         I                           //pTitle//' Anomaly', myThid )
553          diagUnits = 'm^4/s^4         '
554          diagCode  = 'SM P    M1      '
555          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
556         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
557    
558    #ifdef NONLIN_FRSURF
559          diagName  = 'PHIHYDcR'
560          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
561         I                       //pTitle//' Anomaly @ const r', myThid )
562          diagUnits = 'm^2/s^2         '
563          diagCode  = 'SMR     MR      '
564          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
565         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
566    #endif
567    
568    #ifdef ALLOW_NONHYDROSTATIC
569          diagName  = 'PHI_NH  '
570          diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
571          diagUnits = 'm^2/s^2         '
572          diagCode  = 'SMR     MR      '
573          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
574         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
575    #endif /* ALLOW_NONHYDROSTATIC */
576    
577          diagName  = 'MXLDEPTH'
578          diagTitle = 'Mixed-Layer Depth (>0)'
579          diagUnits = 'm               '
580        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
581        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
582       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
583    
584        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
585        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
586        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
587          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
588        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
589        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
590       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 )  
591    
592        diagName  = 'VISCAH  '        diagName  = 'CONVADJ '
593        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '        diagTitle = 'Convective Adjustment Index [0-1] '
594        diagUnits = 'm2/sec          '        diagUnits = 'fraction        '
595        diagCode  = 'SM      MR      '        diagCode  = 'SMR     LR      '
596        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
597       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
598    
599  C--   surface fluxes:  C--   surface fluxes:
600        diagName  = 'TAUX    '        diagName  = 'oceTAUX '
601        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel'
602        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
603        diagCode  = 'SU      U1      '        diagCode  = 'UU      U1      '
604        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
605       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
606         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
607        diagName  = 'TAUY    '  
608        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagName  = 'oceTAUY '
609        diagUnits = 'N/m2            '        diagTitle = 'meridional surf. wind stress, >0 increases vVel'
610        diagCode  = 'SV      U1      '        diagUnits = 'N/m^2           '
611        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagCode  = 'VV      U1      '
612       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagMate  = diagNum
613          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
614         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
615    
616          diagName  = 'atmPload'
617          diagTitle = 'Atmospheric pressure loading'
618          diagUnits = 'Pa              '
619          diagCode  = 'SM      U1      '
620          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
621         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
622    
623        diagName  = 'TFLUX   '        diagName  = 'sIceLoad'
624        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
625        diagUnits = 'W/m2            '        diagUnits = 'kg/m^2          '
626          diagCode  = 'SM      U1      '
627          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
628         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
629    
630          diagName  = 'oceFWflx'
631          diagTitle = 'net surface Fresh-Water flux into the ocean'
632         &          //' (+=down), >0 decreases salinity'
633          diagUnits = 'kg/m^2/s        '
634          diagCode  = 'SM      U1      '
635          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
636         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
637    
638          diagName  = 'oceSflux'
639          diagTitle = 'net surface Salt flux into the ocean (+=down),'
640         &          //' >0 increases salinity'
641          diagUnits = 'g/m^2/s         '
642          diagCode  = 'SM      U1      '
643          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
644         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
645    
646          diagName  = 'oceQnet '
647          diagTitle = 'net surface heat flux into the ocean (+=down),'
648         &          //' >0 increases theta'
649          diagUnits = 'W/m^2           '
650        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
651        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
652       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
653    
654          diagName  = 'oceQsw  '
655          diagTitle = 'net Short-Wave radiation (+=down),'
656         &          //' >0 increases theta'
657          diagUnits = 'W/m^2           '
658          diagCode  = 'SM      U1      '
659          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
660         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
661    
662          diagName  = 'oceFreez'
663          diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
664          diagUnits = 'W/m^2           '
665          diagCode  = 'SM      U1      '
666          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
667         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
668    
669        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
670        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
671        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
672        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
673        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
674       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
675    
676        diagName  = 'TICE    '        diagName  = 'SRELAX  '
677        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'surface salinity relaxation, >0 increases salt'
678        diagUnits = 'W/m2            '        diagUnits = 'g/m^2/s         '
679        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
680        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
681       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
682    
683        diagName  = 'SFLUX   '        diagName  = 'surForcT'
684        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'model surface forcing for Temperature,'
685        diagUnits = 'g/m2/s          '       &          //' >0 increases theta'
686          diagUnits = 'W/m^2           '
687        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
688        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
689       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
690    
691        diagName  = 'SRELAX  '        diagName  = 'surForcS'
692        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'model surface forcing for Salinity,'
693        diagUnits = 'g/m2/s          '       &          //' >0 increases salinity'
694          diagUnits = 'g/m^2/s         '
695          diagCode  = 'SM      U1      '
696          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
697         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
698    
699          diagName  = 'TFLUX   '
700          diagTitle = 'total heat flux (match heat-content variations),'
701         &          //' >0 increases theta'
702          diagUnits = 'W/m^2           '
703          diagCode  = 'SM      U1      '
704          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
705         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
706    
707          diagName  = 'SFLUX   '
708          diagTitle = 'total salt flux (match salt-content variations),'
709         &          //' >0 increases salt'
710          diagUnits = 'g/m^2/s         '
711        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
712        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
713       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
714    
715  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
716    
717          diagName  = 'RCENTER '
718    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
719          diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
720          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
721          diagCode  = 'SM      MR      '
722          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
723         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
724    
725          diagName  = 'RSURF   '
726    c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
727          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
728          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
729          diagCode  = 'SM      M1      '
730          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
731         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
732    
733          diagName  = 'TOTUTEND'
734          diagTitle = 'Tendency of Zonal Component of Velocity'
735          diagUnits = 'm/s/day         '
736          diagCode  = 'UUR     MR      '
737          diagMate  = diagNum + 2
738          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
739         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
740    
741          diagName  = 'TOTVTEND'
742          diagTitle = 'Tendency of Meridional Component of Velocity'
743          diagUnits = 'm/s/day         '
744          diagCode  = 'VVR     MR      '
745          diagMate  = diagNum
746          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
747         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
748    
749          diagName  = 'TOTTTEND'
750          diagTitle = 'Tendency of Potential Temperature'
751          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
752          diagCode  = 'SMR     MR      '
753          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
754         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
755    
756          diagName  = 'TOTSTEND'
757          diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
758          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )
759          diagCode  = 'SMR     MR      '
760          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
761         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
762    
763          diagName  = 'AB_gU   '
764          diagTitle = 'U momentum tendency from Adams-Bashforth'
765          diagUnits = 'm/s^2           '
766          diagCode  = 'UUR     MR      '
767          diagMate  = diagNum + 2
768          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
769         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
770    
771          diagName  = 'AB_gV   '
772          diagTitle = 'V momentum tendency from Adams-Bashforth'
773          diagUnits = 'm/s^2           '
774          diagCode  = 'VVR     MR      '
775          diagMate  = diagNum
776          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
777         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
778    
779    #ifdef ALLOW_NONHYDROSTATIC
780          diagName  = 'AB_gW   '
781          diagTitle = 'W momentum tendency from Adams-Bashforth'
782          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
783          diagCode  = 'WM      LR      '
784          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
785         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
786    #endif /* ALLOW_NONHYDROSTATIC */
787    
788    #ifdef ALLOW_GENERIC_ADVDIFF
789          diagName  = 'AB_gT   '
790          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
791          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
792          diagCode  = 'SMR     MR      '
793          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
794         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
795    
796          diagName  = 'AB_gS   '
797          diagTitle = DIAGS_MK_TITLE(
798         &            sTitle//'tendency from Adams-Bashforth', myThid )
799          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
800          diagCode  = 'SMR     MR      '
801          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
802         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
803    
804          diagName  = 'gTinAB  '
805          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
806          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
807          diagCode  = 'SMR     MR      '
808          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
809         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
810    
811          diagName  = 'gSinAB  '
812          diagTitle = DIAGS_MK_TITLE(
813         &            sTitle//'tendency going in Adams-Bashforth', myThid )
814          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
815          diagCode  = 'SMR     MR      '
816          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
817         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
818    #endif /* ALLOW_GENERIC_ADVDIFF */
819    
820        RETURN        RETURN
821        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.43
  ViewVC Help
Powered by ViewVC 1.1.22