/[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.8 by jmc, Fri Jan 28 01:06:12 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  = 'SALT    '        diagName  = 'SALT    '
124        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
125        diagUnits = 'g/kg            '        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
126        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
127        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
128       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
129    
130        diagName  = 'UVEL    '        diagName  = 'UVEL    '
131        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
132        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
133        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
134        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
135       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
136    
137        diagName  = 'VVEL    '        diagName  = 'VVEL    '
138        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
139        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
140        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
141        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
142       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
143    
144        diagName  = 'WVEL    '        diagName  = 'WVEL    '
145        diagTitle = 'Vertical Component of Velocity (m/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
146        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
147        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
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  = 'THETASQ '        diagName  = 'THETASQ '
152        diagTitle = 'Square of Potential Temperature (K^2)'        diagTitle = 'Square of Potential Temperature (K^2)'
153        diagUnits = 'K2              '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
154        diagCode  = 'SM      MR      '        diagCode  = 'SM      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  = 'SALTSQ  '        diagName  = 'SALTSQ  '
159        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)'
160        diagUnits = 'g2/kg2          '        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
161        diagCode  = 'SM      MR      '        diagCode  = 'SM      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  = 'UVELSQ  '        diagName  = 'UVELSQ  '
166        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
167        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
168        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
169        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
170       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
171    
172        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
173        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
174        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
175        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
176        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
177       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
178    
179        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
180        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
181        diagUnits = 'm2/sec2         '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', 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  = 'UVELVVEL'        diagName  = 'UV_VEL_C'
187          diagTitle ='Product of horizontal Comp of velocity (cell center)'
188          diagUnits = 'm^2/s^2         '
189          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
190          CALL DIAGNOSTICS_ADD2LIST( diagNum,
191         I          diagName, diagCode, diagUnits, diagTitle, myThid )
192    
193          diagName  = 'UV_VEL_Z'
194        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
195        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
196        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
197        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
198       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
199    
200    
201        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
202        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
203        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
204        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
205        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
206       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
207    
208        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
209        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
210        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
211        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
212        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
213       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
214    
215        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
216        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
217        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
218        diagCode  = 'WM      MR      '        diagCode  = 'WM      MR      '
219        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
220       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
221    
222        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
223        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'
224        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
225        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
226        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
227       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
228    
229        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
230        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'
231        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
232        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
233        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
234       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
235    
236        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
237        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
238        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
239        diagCode  = 'WM      MR      '        diagCode  = 'WM      MR      '
240        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
241       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
242    
243        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
244        diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'
245        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
246        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
247        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
248       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
249    
250        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
251        diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'
252        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
253        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
254        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
255       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
256    
257        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
258        diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'
259        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
260        diagCode  = 'WM      MR      '        diagCode  = 'WM      MR      '
261        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
262       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
263    
264        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
265        diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
266        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
267        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
268        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
269       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
270    
271        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
272        diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
273        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
274        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
275        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
276       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
277    
278        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
279        diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
280        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
281        diagCode  = 'WM      MR      '        diagCode  = 'WM      MR      '
282        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
283       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
284    
285        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
286        diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'
287        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
288        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
289        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
290       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
291    
292        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
293        diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'
294        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
295        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
296        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
297       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
298    
299        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
300        diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'
301        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
302        diagCode  = 'WM      MR      '        diagCode  = 'WM      MR      '
303        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
304       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 288  C---+----1----+----2----+----3----+----4 Line 314  C---+----1----+----2----+----3----+----4
314    
315        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
316        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
317        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
318        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
319        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
320       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
321    
322        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
323        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
324        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
325        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
326        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
327       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
328    
329        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
330        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
331        diagUnits = 'm4/s4           '        diagUnits = 'm^4/s^4         '
332        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
333        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
334       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
335    
336        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
337        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
338        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
339          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
340        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
341        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
342       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 317  C---+----1----+----2----+----3----+----4 Line 344  C---+----1----+----2----+----3----+----4
344  C--   Momentum terms:  C--   Momentum terms:
345        diagName  = 'VISCA4  '        diagName  = 'VISCA4  '
346        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
347        diagUnits = 'm4/sec          '        diagUnits = 'm^4/s           '
348        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
349        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
350       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
351    
352        diagName  = 'VISCAH  '        diagName  = 'VISCAH  '
353        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
354        diagUnits = 'm2/sec          '        diagUnits = 'm^2/s           '
355        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
356        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
357       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 332  C--   Momentum terms: Line 359  C--   Momentum terms:
359  C--   surface fluxes:  C--   surface fluxes:
360        diagName  = 'TAUX    '        diagName  = 'TAUX    '
361        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
362        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
363        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
364        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
365       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367        diagName  = 'TAUY    '        diagName  = 'TAUY    '
368        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel '
369        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
370        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
371        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
372       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
373    
374        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
375        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'net surface heat flux, >0 increases theta       '
376        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
377        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
378        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
379       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
380    
381        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
382        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
383        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
384        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
385        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
386       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
387    
388        diagName  = 'TICE    '        diagName  = 'TICE    '
389        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
390        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
391        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
392        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
393       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
394    
395        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
396        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt        '
397        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
398        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
399        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
400       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
401    
402        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
403        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt  '
404        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
405        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
406        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
407       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22