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