/[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.29 by jmc, Mon Nov 20 04:25:19 2006 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 66  C     fTitle     :: fixed boundary title Line 68  C     fTitle     :: fixed boundary title
68  C     pTitle     :: "Phi"  title  C     pTitle     :: "Phi"  title
69  C     sTitle     :: "salt" title  C     sTitle     :: "salt" title
70        INTEGER        diagNum        INTEGER        diagNum
71          INTEGER        diagMate
72        CHARACTER*8    diagName        CHARACTER*8    diagName
73        CHARACTER*16   diagCode        CHARACTER*16   diagCode
74        CHARACTER*16   diagUnits        CHARACTER*16   diagUnits
# Line 122  C-    free-surface (eTitle) and fixed-bo Line 125  C-    free-surface (eTitle) and fixed-bo
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 = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
131  c     IF ( fluidIsWater .AND. usingZCoords )  c     IF ( fluidIsWater .AND. usingZCoords )
132  c    &diagTitle = 'Sea Surface Elevation'  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 = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
140       I                          //' Anomaly', myThid )       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 = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
148       I                          //' Anomaly Tendency', myThid )       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'        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  c     diagName  = 'SST     '  c     diagName  = 'SST     '
162  c     diagTitle = 'Sea Surface Temperature (degC,K)'  c     diagTitle = 'Sea Surface Temperature (degC,K)'
163  c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
164  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
165  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
166  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
167    
168        diagName  = 'SALT    '        diagName  = 'SALT    '
169        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        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  = 'RELHUM  '        diagName  = 'RELHUM  '
176        diagTitle = 'Relative Humidity'        diagTitle = 'Relative Humidity'
177        diagUnits = 'percent         '        diagUnits = 'percent         '
178        diagCode  = 'SM      MR      '        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     '  c     diagName  = 'SSS     '
183  c     diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
184  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
185  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
186  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
187  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
188    
189        IF ( fluidIsWater ) THEN        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        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 )
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'        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 = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        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        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        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'        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'        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'        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'        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'        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'        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 = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
361       I                           //sTitle, myThid )       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 = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
370       I                           //sTitle, myThid )       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 = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
379       I                           //sTitle, myThid )       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 Transport of Pot Temp'        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 Transport of Pot Temp'        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 Transport of Pot Temp'        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 = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
410       I                          //sTitle, myThid )       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 = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
419       I                          //sTitle, myThid )       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 = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
428       I                          //sTitle, myThid )       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 '        diagName  = 'UVELPHI '
435        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
436       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
437        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
438        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UUr     MR      '
439        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
440       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
441         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
442    
443        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
444        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
445       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
446        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
447        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVr     MR      '
448        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
449       I          diagName, diagCode, diagUnits, diagTitle, myThid )        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 = DIAGS_MK_TITLE( 'Hydrostatic '        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
493       I                           //pTitle//' Anomaly', myThid )       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 = DIAGS_MK_TITLE( 'Square of Hyd. '        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
501       I                           //pTitle//' Anomaly', myThid )       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  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
# Line 470  c     diagTitle = 'ocean bottom pressure Line 510  c     diagTitle = 'ocean bottom pressure
510       I                           //pTitle//' Anomaly', myThid )       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  c     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        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
519       I                           //pTitle//' Anomaly', myThid )       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        diagName  = 'CONVADJ '        diagName  = 'CONVADJ '
550        diagTitle = 'Convective Adjustment Index [0-1] '        diagTitle = 'Convective Adjustment Index [0-1] '
551        diagUnits = 'fraction        '        diagUnits = 'fraction        '
552        diagCode  = 'SM      LR      '        diagCode  = 'SMR     LR      '
553        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
554       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
555    
556  C--   surface fluxes:  C--   surface fluxes:
557        diagName  = 'oceTAUX '        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  = 'oceTAUY '        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    
 #ifdef ATMOSPHERIC_LOADING  
573        diagName  = 'atmPload'        diagName  = 'atmPload'
574        diagTitle = 'Atmospheric pressure loading'        diagTitle = 'Atmospheric pressure loading'
575        diagUnits = 'Pa              '        diagUnits = 'Pa              '
576        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
577        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
578       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
579    
580        diagName  = 'sIceLoad'        diagName  = 'sIceLoad'
581        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'        diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
582        diagUnits = 'kg/m^2          '        diagUnits = 'kg/m^2          '
583        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
584        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
585       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
 #endif /* ATMOSPHERIC_LOADING */  
586    
587        diagName  = 'oceFWflx'        diagName  = 'oceFWflx'
588        diagTitle = 'net surface Fresh-Water flux into the ocean'        diagTitle = 'net surface Fresh-Water flux into the ocean'
589       &          //' (+=down), >0 decreases salinity'       &          //' (+=down), >0 decreases salinity'
590        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
591        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
592        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
593       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
594    
595        diagName  = 'oceSflux'        diagName  = 'oceSflux'
596        diagTitle = 'net surface Salt flux into the ocean (+=down),'        diagTitle = 'net surface Salt flux into the ocean (+=down),'
597       &          //' >0 increases salinity'       &          //' >0 increases salinity'
598        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
599        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
600        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
601       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
602    
603        diagName  = 'oceQnet '        diagName  = 'oceQnet '
604        diagTitle = 'net surface heat flux into the ocean (+=down),'        diagTitle = 'net surface heat flux into the ocean (+=down),'
605       &          //' >0 increases theta'       &          //' >0 increases theta'
606        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
607        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
608        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
609       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
610    
611        diagName  = 'oceQsw  '        diagName  = 'oceQsw  '
612        diagTitle = 'net Short-Wave radiation (+=down),'        diagTitle = 'net Short-Wave radiation (+=down),'
613       &          //' >0 increases theta'       &          //' >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'        diagName  = 'oceFreez'
620        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
621        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
622        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
623        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
624       I          diagName, diagCode, diagUnits, diagTitle, myThid )       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  '        diagName  = 'SRELAX  '
634        diagTitle = 'surface salinity relaxation, >0 increases salt'        diagTitle = 'surface salinity relaxation, >0 increases salt'
635        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
636        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
637        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
638       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
639    
640        diagName  = 'surForcT'        diagName  = 'surForcT'
641        diagTitle = 'model surface forcing for Temperature,'        diagTitle = 'model surface forcing for Temperature,'
642       &          //' >0 increases theta'       &          //' >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  = 'surForcS'        diagName  = 'surForcS'
649        diagTitle = 'model surface forcing for Salinity,'        diagTitle = 'model surface forcing for Salinity,'
650       &          //' >0 increases salinity'       &          //' >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  = 'TFLUX   '        diagName  = 'TFLUX   '
657        diagTitle = 'total heat flux (match heat-content variations),'        diagTitle = 'total heat flux (match heat-content variations),'
658       &          //' >0 increases theta'       &          //' >0 increases theta'
659        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
660        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
661        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
662       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
663    
664        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
665        diagTitle = 'total salt flux (match salt-content variations),'        diagTitle = 'total salt flux (match salt-content variations),'
666       &          //' >0 increases salt'       &          //' >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    
# Line 620  c     diagTitle = 'Cell-Center r-Positio Line 676  c     diagTitle = 'Cell-Center r-Positio
676        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )        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  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
684        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )        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'        diagName  = 'TOTUTEND'
691        diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Zonal Component of Velocity'
692        diagUnits = 'm/s^2           '        diagUnits = 'm/s/day         '
693        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        diagCode  = 'UUR     MR      '
694        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
695       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
696         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
697    
698        diagName  = 'TOTVTEND'        diagName  = 'TOTVTEND'
699        diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'        diagTitle = 'Tendency of Meridional Component of Velocity'
700        diagUnits = 'm/s^2           '        diagUnits = 'm/s/day         '
701        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        diagCode  = 'VVR     MR      '
702        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
703       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
704         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
705    
706        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
707        diagTitle = 'Tendency of Potential Temperature'        diagTitle = 'Tendency of Potential Temperature'
708        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
709        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
710        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
711       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
712    
713        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
714        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
715        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )
716        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
717        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
718       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
719    
720    
721        RETURN        RETURN
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22