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