/[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.7 by jmc, Mon Jan 3 02:32:33 2005 UTC revision 1.12 by dimitri, Fri May 6 03:46:40 2005 UTC
# Line 64  C     !LOCAL VARIABLES: Line 64  C     !LOCAL VARIABLES:
64        CHARACTER*16  diagCode        CHARACTER*16  diagCode
65        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
66        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
67          CHARACTER*2   rUnit2c
68          CHARACTER*4   tUnit4c, sUnit4c
69    
70          CHARACTER*(16) DIAGS_MK_UNITS
71          EXTERNAL DIAGS_MK_UNITS
72    
73  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
74  C     For each output variable,  C     For each output variable,
# Line 71  C     specify Name (cdiag, 8c), Descript Line 76  C     specify Name (cdiag, 8c), Descript
76  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
77  C----------------------------------------------------------------------  C----------------------------------------------------------------------
78    
79          IF ( usingPCoords ) THEN
80            rUnit2c= 'Pa'
81          ELSE
82            rUnit2c= 'm '
83          ENDIF
84          IF ( fluidIsAir ) THEN
85            tUnit4c= 'K   '
86            sUnit4c= 'g/kg'
87          ELSE
88            tUnit4c= 'degC'
89            sUnit4c= 'psu '
90          ENDIF
91    
92  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
93  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
94    
95        diagName  = 'ETAN    '        diagName  = 'ETAN    '
96        diagTitle = 'Perturbation of Surface (pressure, height) (hPa,m)'        diagTitle = 'Perturbation of Surface (pressure, height) (Pa,m)'
97        diagUnits = 'm or hPa        '        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
98        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
99        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
100       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
101    
102        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
103        diagTitle = 'Square of Perturbation of Sfc (hPa^2,m^2)'        diagTitle = 'Square of Perturbation of Sfc (Pa^2,m^2)'
104        diagUnits = 'm2 or hPa2      '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
105        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
106        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
107       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
108    
109        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
110        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
111        diagUnits = 'm2 or Pa2 /sec2 '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
112        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
113        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
114       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
115    
116        diagName  = 'THETA   '        diagName  = 'THETA   '
117        diagTitle = 'Potential Temperature (K)'        diagTitle = 'Potential Temperature (degC,K)'
118        diagUnits = 'K               '        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
119        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
120        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
121       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
122    
123          diagName  = 'SST     '
124          diagTitle = 'Sea Surface Temperature (degC,K)'
125          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
126          diagCode  = 'SM      M1      '
127          CALL DIAGNOSTICS_ADD2LIST( diagNum,
128         I          diagName, diagCode, diagUnits, diagTitle, myThid )
129    
130        diagName  = 'SALT    '        diagName  = 'SALT    '
131        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
132        diagUnits = 'g/kg            '        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
133          diagCode  = 'SM      MR      '
134          CALL DIAGNOSTICS_ADD2LIST( diagNum,
135         I          diagName, diagCode, diagUnits, diagTitle, myThid )
136    
137          diagName  = 'SSS     '
138          diagTitle = 'Sea Surface Salinity (g/kg)'
139          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
140          diagCode  = 'SM      M1      '
141          CALL DIAGNOSTICS_ADD2LIST( diagNum,
142         I          diagName, diagCode, diagUnits, diagTitle, myThid )
143    
144          diagName  = 'SALTanom'
145          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
146          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
147        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
148        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
149       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
150    
151        diagName  = 'UVEL    '        diagName  = 'UVEL    '
152        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
153        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
154        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
155        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
156       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
157    
158        diagName  = 'VVEL    '        diagName  = 'VVEL    '
159        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
160        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
161        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
162        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
163       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
164    
165          diagName  = 'UVEL_k2 '
166          diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'
167          diagUnits = 'm/s             '
168          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
169          CALL DIAGNOSTICS_ADD2LIST( diagNum,
170         I          diagName, diagCode, diagUnits, diagTitle, myThid )
171    
172          diagName  = 'VVEL_k2 '
173          diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'
174          diagUnits = 'm/s             '
175          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '
176          CALL DIAGNOSTICS_ADD2LIST( diagNum,
177         I          diagName, diagCode, diagUnits, diagTitle, myThid )
178    
179        diagName  = 'WVEL    '        diagName  = 'WVEL    '
180        diagTitle = 'Vertical Component of Velocity (m/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
181        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
182        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
183        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
184       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
185    
186        diagName  = 'THETASQ '        diagName  = 'THETASQ '
187        diagTitle = 'Square of Potential Temperature (K^2)'        diagTitle = 'Square of Potential Temperature (K^2)'
188        diagUnits = 'K2              '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
189        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
190        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
191       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
192    
193        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
194        diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'        diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
195        diagUnits = 'g2/kg2          '        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
196          diagCode  = 'SM      MR      '
197          CALL DIAGNOSTICS_ADD2LIST( diagNum,
198         I          diagName, diagCode, diagUnits, diagTitle, myThid )
199    
200          diagName  = 'SALTSQan'
201          diagTitle = 'Square of Salt anomaly (=SALTSQ-35^2 (g^2/kg^2)'
202          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
203        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
204        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
205       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
206    
207        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
208        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
209        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
210        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
211        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
212       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
213    
214        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
215        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
216        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
217        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
218        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
219       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
220    
221        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
222        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
223        diagUnits = 'm2/sec2         '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
224        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
225        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
226       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
227    
228        diagName  = 'UVELVVEL'        diagName  = 'UV_VEL_C'
229          diagTitle ='Product of horizontal Comp of velocity (cell center)'
230          diagUnits = 'm^2/s^2         '
231          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
232          CALL DIAGNOSTICS_ADD2LIST( diagNum,
233         I          diagName, diagCode, diagUnits, diagTitle, myThid )
234    
235          diagName  = 'UV_VEL_Z'
236        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
237        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
238        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
239        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
240       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
241    
242          diagName  = 'WU_VEL  '
243          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
244          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
245          diagCode  = 'WU      LR      '
246          CALL DIAGNOSTICS_ADD2LIST( diagNum,
247         I          diagName, diagCode, diagUnits, diagTitle, myThid )
248    
249          diagName  = 'WV_VEL  '
250          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
251          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
252          diagCode  = 'WV      LR      '
253          CALL DIAGNOSTICS_ADD2LIST( diagNum,
254         I          diagName, diagCode, diagUnits, diagTitle, myThid )
255    
256        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
257        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
258        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
259        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
260        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
261       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
262    
263        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
264        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
265        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
266        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
267        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
268       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
269    
270        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
271        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
272        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
273        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
274        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
275       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
276    
277        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
278        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'
279        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
280        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
281        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
282       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
283    
284        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
285        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'
286        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
287        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
288        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
289       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
290    
291        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
292        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
293        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
294        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
295        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
296       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
297    
298        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
299        diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'
300        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
301        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
302        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
303       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
304    
305        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
306        diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'
307        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
308        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
309        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
310       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
311    
312        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
313        diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'
314        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
315        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
316        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
317       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
318    
319        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
320        diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
321        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
322        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
323        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
324       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
325    
326        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
327        diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
328        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
329        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
330        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
331       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
332    
333        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
334        diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
335        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
336        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
337        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
338       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
339    
340        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
341        diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'
342        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
343        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
344        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
345       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
346    
347        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
348        diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'
349        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
350        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
351        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
352       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
353    
354        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
355        diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'
356        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
357        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
358        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
359       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
360    
# Line 286  C---+----1----+----2----+----3----+----4 Line 367  C---+----1----+----2----+----3----+----4
367        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
368       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
369    
370          diagName  = 'RHOANOSQ'
371          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
372          diagUnits = 'kg^2/m^6          '
373          diagCode  = 'SM      MR      '
374          CALL DIAGNOSTICS_ADD2LIST( diagNum,
375         I          diagName, diagCode, diagUnits, diagTitle, myThid )
376    
377          diagName  = 'URHOMASS'
378          diagTitle = 'Zonal Transport of Density'
379          diagUnits = 'kg/m^2/s        '
380          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
381          CALL DIAGNOSTICS_ADD2LIST( diagNum,
382         I          diagName, diagCode, diagUnits, diagTitle, myThid )
383    
384          diagName  = 'VRHOMASS'
385          diagTitle = 'Meridional Transport of Density'
386          diagUnits = 'kg/m^2/s        '
387          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
388          CALL DIAGNOSTICS_ADD2LIST( diagNum,
389         I          diagName, diagCode, diagUnits, diagTitle, myThid )
390    
391          diagName  = 'WRHOMASS'
392          diagTitle = 'Vertical Transport of Potential Density'
393          diagUnits = 'kg/m^2/s        '
394          diagCode  = 'WM      LR      '
395          CALL DIAGNOSTICS_ADD2LIST( diagNum,
396         I          diagName, diagCode, diagUnits, diagTitle, myThid )
397    
398        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
399        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
400        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
401        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
402        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
403       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
404    
405        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
406        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
407        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
408        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
409        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
410       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
411    
412        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
413        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
414        diagUnits = 'm4/s4           '        diagUnits = 'm^4/s^4         '
415        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
416        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
417       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
418    
419        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
420        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
421        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
422          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
423        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
424        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
425       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 317  C---+----1----+----2----+----3----+----4 Line 427  C---+----1----+----2----+----3----+----4
427  C--   Momentum terms:  C--   Momentum terms:
428        diagName  = 'VISCA4  '        diagName  = 'VISCA4  '
429        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
430        diagUnits = 'm4/sec          '        diagUnits = 'm^4/s           '
431        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
432        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
433       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
434    
435        diagName  = 'VISCAH  '        diagName  = 'VISCAH  '
436        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
437        diagUnits = 'm2/sec          '        diagUnits = 'm^2/s           '
438        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
439        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
440       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 332  C--   Momentum terms: Line 442  C--   Momentum terms:
442  C--   surface fluxes:  C--   surface fluxes:
443        diagName  = 'TAUX    '        diagName  = 'TAUX    '
444        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
445        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
446        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
447        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
448       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
449    
450        diagName  = 'TAUY    '        diagName  = 'TAUY    '
451        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel '
452        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
453        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
454        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
455       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
456    
457        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
458        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'net surface heat flux, >0 increases theta       '
459        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
460        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
461        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
462       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
463    
464        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
465        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
466        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
467        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
468        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
469       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
470    
471        diagName  = 'TICE    '        diagName  = 'TICE    '
472        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
473        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
474        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
475        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
476       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
477    
478        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
479        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt        '
480        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
481        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
482        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
483       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
484    
485        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
486        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt  '
487        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
488        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
489        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
490       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22