/[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.25 by molod, Mon Feb 27 17:55:14 2006 UTC revision 1.26 by jmc, Sun Jun 25 23:06:24 2006 UTC
# Line 60  C     !INPUT PARAMETERS: Line 60  C     !INPUT PARAMETERS:
60  CEOP  CEOP
61    
62  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
63        INTEGER       diagNum  C     rTitle     :: r-coordinate title
64        CHARACTER*8   diagName  C     eTitle     :: free-surface title
65        CHARACTER*16  diagCode  C     fTitle     :: fixed boundary title
66        CHARACTER*16  diagUnits  C     pTitle     :: "Phi"  title
67    C     sTitle     :: "salt" title
68          INTEGER        diagNum
69          CHARACTER*8    diagName
70          CHARACTER*16   diagCode
71          CHARACTER*16   diagUnits
72        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
73        CHARACTER*2   rUnit2c        CHARACTER*2    rUnit2c
74        CHARACTER*4   tUnit4c, sUnit4c        CHARACTER*4    tUnit4c, sUnit4c
75          CHARACTER*(10) rTitle, eTitle, fTitle
76          CHARACTER*(20) pTitle, sTitle
77    
78        CHARACTER*(16) DIAGS_MK_UNITS        CHARACTER*(16) DIAGS_MK_UNITS
79        EXTERNAL DIAGS_MK_UNITS        EXTERNAL DIAGS_MK_UNITS
80          CHARACTER*(80) DIAGS_MK_TITLE
81          EXTERNAL DIAGS_MK_TITLE
82    
83  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
84  C     For each output variable,  C     For each output variable,
85  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
86  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
87  C----------------------------------------------------------------------  C----------------------------------------------------------------------
88    
89        IF ( usingPCoords ) THEN        IF ( usingPCoords ) THEN
90          rUnit2c= 'Pa'          rUnit2c= 'Pa'
91            rTitle = ' Pressure '
92            pTitle = ' Geopotential       '
93        ELSE        ELSE
94          rUnit2c= 'm '          rUnit2c= 'm '
95            rTitle = ' Height   '
96            pTitle = 'Pressure Pot.(p/rho)'
97        ENDIF        ENDIF
98        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
99          tUnit4c= 'K   '          tUnit4c= 'K   '
100          sUnit4c= 'g/kg'          sUnit4c= 'g/kg'
101            sTitle = ' Specific Humidity  '
102  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
103        IF(useFIZHI)THEN          IF (useFIZHI) sUnit4c= 'kg/kg'
         sUnit4c= 'kg/kg'  
       ENDIF  
104  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
105        ELSE        ELSE
106          tUnit4c= 'degC'          tUnit4c= 'degC'
107          sUnit4c= 'psu '          sUnit4c= 'psu '
108            sTitle = ' Salinity           '
109          ENDIF
110    C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
111          IF ( fluidIsAir ) THEN
112            eTitle = ' Surface  '
113            fTitle = ' Top      '
114          ELSEIF ( usingPCoords ) THEN
115            eTitle = ' Bottom   '
116            fTitle = ' Surface  '
117          ELSE
118            eTitle = ' Surface  '
119            fTitle = ' Bottom   '
120        ENDIF        ENDIF
121    
122  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
123  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
124    
125        diagName  = 'ETAN    '        diagName  = 'ETAN    '
126        diagTitle = 'Perturbation of Surface (pressure, height) (Pa,m)'        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
127    c     IF ( fluidIsWater .AND. usingZCoords )
128    c    &diagTitle = 'Sea Surface Elevation'
129        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
130        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
131        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
132       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
133    
134        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
135        diagTitle = 'Square of Perturbation of Sfc position (r_unit^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
136         I                          //' Anomaly', myThid )
137        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
138        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
139        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
140       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
141    
142        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
143        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
144         I                          //' Anomaly Tendency', myThid )
145        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
146        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
147        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
148       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150        diagName  = 'THETA   '        diagName  = 'THETA   '
151        diagTitle = 'Potential Temperature (degC,K)'        diagTitle = 'Potential Temperature'
152        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
153        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
154        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 134  c     CALL DIAGNOSTICS_ADD2LIST( diagNum Line 162  c     CALL DIAGNOSTICS_ADD2LIST( diagNum
162  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
163    
164        diagName  = 'SALT    '        diagName  = 'SALT    '
165        diagTitle = 'Salt or Water Vapor mixing ratio'        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
166        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
167        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
168        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
169       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
170    
171        diagName  = 'RELHUM  '        diagName  = 'RELHUM  '
172        diagTitle = 'Relative Humidity               '        diagTitle = 'Relative Humidity'
173        diagUnits = 'percent         '        diagUnits = 'percent         '
174        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
175        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 154  c     diagCode  = 'SM      M1      ' Line 182  c     diagCode  = 'SM      M1      '
182  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
183  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185          IF ( fluidIsWater ) THEN
186        diagName  = 'SALTanom'        diagName  = 'SALTanom'
187        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
188        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, 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          ENDIF
193    
194        diagName  = 'UVEL    '        diagName  = 'UVEL    '
195        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
# Line 183  c    I          diagName, diagCode, diag Line 213  c    I          diagName, diagCode, diag
213       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
214    
215        diagName  = 'THETASQ '        diagName  = 'THETASQ '
216        diagTitle = 'Square of Potential Temperature (K^2)'        diagTitle = 'Square of Potential Temperature'
217        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
218        diagCode  = 'SM      MR      '        diagCode  = 'SM      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  = 'SALTSQ  '        diagName  = 'SALTSQ  '
223        diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
224        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
225        diagCode  = 'SM      MR      '        diagCode  = 'SM      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          IF ( fluidIsWater ) THEN
230        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
231        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
232        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
234        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
235       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
236          ENDIF
237    
238        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
239        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
# Line 218  c    I          diagName, diagCode, diag Line 250  c    I          diagName, diagCode, diag
250       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
251    
252        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
253        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Vertical Comp of Velocity'
254        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
255        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
256        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 239  c    I          diagName, diagCode, diag Line 271  c    I          diagName, diagCode, diag
271       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
272    
273        diagName  = 'WU_VEL  '        diagName  = 'WU_VEL  '
274        diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Vertical Transport of Zonal Momentum'
275        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
276        diagCode  = 'WU      LR      '        diagCode  = 'WU      LR      '
277        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
278       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
279    
280        diagName  = 'WV_VEL  '        diagName  = 'WV_VEL  '
281        diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'        diagTitle ='Vertical Transport of Meridional Momentum'
282        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
283        diagCode  = 'WV      LR      '        diagCode  = 'WV      LR      '
284        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 267  c    I          diagName, diagCode, diag Line 299  c    I          diagName, diagCode, diag
299       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
300    
301        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
302        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
303        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
304        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
305        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
306       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
307    
308        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
309        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
310        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
311        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
312        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
313       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
314    
315        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
316        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
317        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
318        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
319        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 295  c    I          diagName, diagCode, diag Line 327  c    I          diagName, diagCode, diag
327       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
328    
329        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
330        diagTitle = 'Zonal Mass-Weight Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
331         I                           //sTitle, myThid )
332        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
333        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
334        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
335       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
336    
337        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
338        diagTitle = 'Meridional Mass-Weight Transp of Salt'        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
339         I                           //sTitle, myThid )
340        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
341        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
342        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
343       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
344    
345        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
346        diagTitle = 'Vertical Mass-Weight Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
347         I                           //sTitle, myThid )
348        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
349        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
350        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
351       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
354        diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'        diagTitle = 'Zonal Transport of Pot Temp'
355        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
356        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
357        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
358       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
361        diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'        diagTitle = 'Meridional Transport of Pot Temp'
362        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
363        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
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  = 'WVELTH  '        diagName  = 'WVELTH  '
368        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Transport of Pot Temp'
369        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
370        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
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  = 'UVELSLT '        diagName  = 'UVELSLT '
375        diagTitle = 'Zonal Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
376         I                          //sTitle, myThid )
377        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
378        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
379        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
380       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
381    
382        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
383        diagTitle = 'Meridional Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
384         I                          //sTitle, myThid )
385        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
386        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
387        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
388       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
389    
390        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
391        diagTitle = 'Vertical Transp of Salt '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
392         I                          //sTitle, myThid )
393        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
394        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
395        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
396       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
397    
398        diagName  = 'UVELPHI '        diagName  = 'UVELPHI '
399        diagTitle = 'Zonal Transp of Geopotential Height (m2/s2.m/s)'        diagTitle = DIAGS_MK_TITLE( 'Zonal Transp of '
400         I                 //pTitle//' Anomaly', myThid )
401        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
402        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
403        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
404       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
405    
406        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
407        diagTitle = 'Meridional Transp of Geopotential Height (m2/s2.m/s)'        diagTitle = DIAGS_MK_TITLE( 'Meridional Transp of '
408         I                 //pTitle//' Anomaly', myThid )
409        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
410        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
411        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 409  C---+----1----+----2----+----3----+----4 Line 449  C---+----1----+----2----+----3----+----4
449       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
450    
451        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
452        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
453         I                           //pTitle//' Anomaly', myThid )
454        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
455        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
456        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
457       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
458    
459        diagName  = 'PHIHYDSQ'        diagName  = 'PHIHYDSQ'
460        diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
461         I                           //pTitle//' Anomaly', myThid )
462        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
463        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
464        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
465       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
466    
467        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
468        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
469          diagTitle = DIAGS_MK_TITLE( fTitle
470         I                           //pTitle//' Anomaly', myThid )
471        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
472        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
473        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
474       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
475    
476        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
477        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
478          diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
479         I                           //pTitle//' Anomaly', myThid )
480        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
481        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
482        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 453  C---+----1----+----2----+----3----+----4 Line 499  C---+----1----+----2----+----3----+----4
499    
500  C--   surface fluxes:  C--   surface fluxes:
501        diagName  = 'TAUX    '        diagName  = 'TAUX    '
502        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel'
503        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
504        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
505        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
506       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
507    
508        diagName  = 'TAUY    '        diagName  = 'TAUY    '
509        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel'
510        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
511        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
512        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 488  C--   surface fluxes: Line 534  C--   surface fluxes:
534       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
535    
536        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
537        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt'
538        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
539        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
540        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
541       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
542    
543        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
544        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt'
545        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
546        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
547        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
548       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
549    
 #ifdef SHORTWAVE_HEATING  
550        diagName  = 'SWFLUX  '        diagName  = 'SWFLUX  '
551        diagTitle = 'net upward SW radiation, >0 increases theta     '        diagTitle = 'net upward SW radiation, >0 increases theta'
552        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
553        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
554        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
555       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
 #endif  
556    
557  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
558        diagName  = 'PRESSURE'        diagName  = 'PRESSURE'
559        diagTitle = 'Atmospheric Pressure (Pa)     '  C-    will be changed to an isomorphic name, such as:
560    c     diagName  = 'RCENTER '
561    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
562          diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
563        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
564        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
565        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
566       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
567    
568        diagName  = 'RSURF   '        diagName  = 'RSURF   '
569        diagTitle = 'Surface (Pressure, Height) (Pa,m)'  c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
570          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
571        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
572        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
573        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 540  C---+----1----+----2----+----3----+----4 Line 588  C---+----1----+----2----+----3----+----4
588       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
589    
590        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
591        diagTitle = 'Tendency of Potential Temperature (degC,K)'        diagTitle = 'Tendency of Potential Temperature'
592        diagUnits = 'deg K/sec'        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
593        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
594        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
595       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
596    
597        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
598        diagTitle = 'Tendency of Salt or Water Vapor mixing ratio'        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
599        diagUnits = 'kg/kg-sec'        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
600        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
601        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
602       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22