/[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.6 by jmc, Sat Dec 18 23:53:04 2004 UTC revision 1.22 by jmc, Thu Sep 15 14:51:02 2005 UTC
# Line 24  C       \begin{tabular}[h]{|c|c|}\hline Line 24  C       \begin{tabular}[h]{|c|c|}\hline
24  C         \textbf{Positions}  &  \textbf{Characters}  C         \textbf{Positions}  &  \textbf{Characters}
25  C         &  \textbf{Meanings} \\\hline  C         &  \textbf{Meanings} \\\hline
26  C         parse(1)  &  S  &  scalar \\  C         parse(1)  &  S  &  scalar \\
27  C                   &  U  &  U-vector component  \\  C                   &  U  &  vector component in X direction \\
28  C                   &  V  &  V-vector component  \\  C                   &  V  &  vector component in Y direction \\
29    C                   &  W  &  vector component in vertical direction \\
30  C         parse(2)  &  U  &  C-grid U-Point  \\  C         parse(2)  &  U  &  C-grid U-Point  \\
31  C                   &  V  &  C-grid V-Point  \\  C                   &  V  &  C-grid V-Point  \\
32  C                   &  M  &  C-grid Mass Point  \\  C                   &  M  &  C-grid Mass Point  \\
# Line 43  C                   &  1  &  levels = 1 Line 44  C                   &  1  &  levels = 1
44  C                   &  R  &  levels = Nr  \\  C                   &  R  &  levels = Nr  \\
45  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\
46  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
47    C                   &  G  &  levels = Ground_level Number \\
48    C                   &  I  &  levels = sea-Ice_level Number \\
49  C       \end{tabular}  C       \end{tabular}
50  C     \end{center}  C     \end{center}
51    
# Line 50  C     !USES: Line 53  C     !USES:
53        IMPLICIT NONE        IMPLICIT NONE
54  #include "SIZE.h"  #include "SIZE.h"
55  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 c #include "EESUPPORT.h"  
56  #include "PARAMS.h"  #include "PARAMS.h"
57    
58  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 63  C     !LOCAL VARIABLES: Line 65  C     !LOCAL VARIABLES:
65        CHARACTER*16  diagCode        CHARACTER*16  diagCode
66        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
67        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
68          CHARACTER*2   rUnit2c
69          CHARACTER*4   tUnit4c, sUnit4c
70    
71          CHARACTER*(16) DIAGS_MK_UNITS
72          EXTERNAL DIAGS_MK_UNITS
73    
74  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
75  C     For each output variable,  C     For each output variable,
# Line 70  C     specify Name (cdiag, 8c), Descript Line 77  C     specify Name (cdiag, 8c), Descript
77  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
78  C----------------------------------------------------------------------  C----------------------------------------------------------------------
79    
80          IF ( usingPCoords ) THEN
81            rUnit2c= 'Pa'
82          ELSE
83            rUnit2c= 'm '
84          ENDIF
85          IF ( fluidIsAir ) THEN
86            tUnit4c= 'K   '
87            sUnit4c= 'g/kg'
88    #ifdef ALLOW_FIZHI
89          IF(useFIZHI)THEN
90            sUnit4c= 'kg/kg'
91          ENDIF
92    #endif /* ALLOW_FIZHI */
93          ELSE
94            tUnit4c= 'degC'
95            sUnit4c= 'psu '
96          ENDIF
97    
98  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
99  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
100    
101        diagName  = 'ETAN    '        diagName  = 'ETAN    '
102        diagTitle = 'Perturbation of Surface (pressure, height) (hPa,m)'        diagTitle = 'Perturbation of Surface (pressure, height) (Pa,m)'
103        diagUnits = 'm or hPa        '        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
104        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
105        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
106       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
107    
108        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
109        diagTitle = 'Square of Perturbation of Sfc (hPa^2,m^2)'        diagTitle = 'Square of Perturbation of Sfc position (r_unit^2)'
110        diagUnits = 'm2 or hPa2      '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
111        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
112        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
113       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
114    
115        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
116        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
117        diagUnits = 'm2 or Pa2 /sec2 '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
118        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
119        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
120       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
121    
122        diagName  = 'THETA   '        diagName  = 'THETA   '
123        diagTitle = 'Potential Temperature (K)'        diagTitle = 'Potential Temperature (degC,K)'
124        diagUnits = 'K               '        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
125        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
126        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
127       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
128    
129    c     diagName  = 'SST     '
130    c     diagTitle = 'Sea Surface Temperature (degC,K)'
131    c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
132    c     diagCode  = 'SM      M1      '
133    c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
134    c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
135    
136        diagName  = 'SALT    '        diagName  = 'SALT    '
137        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = 'Salt or Water Vapor mixing ratio'
138        diagUnits = 'g/kg            '        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
139          diagCode  = 'SM      MR      '
140          CALL DIAGNOSTICS_ADD2LIST( diagNum,
141         I          diagName, diagCode, diagUnits, diagTitle, myThid )
142    
143          diagName  = 'RELHUM  '
144          diagTitle = 'Relative Humidity               '
145          diagUnits = 'percent         '
146          diagCode  = 'SM      MR      '
147          CALL DIAGNOSTICS_ADD2LIST( diagNum,
148         I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150    c     diagName  = 'SSS     '
151    c     diagTitle = 'Sea Surface Salinity '
152    c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
153    c     diagCode  = 'SM      M1      '
154    c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
155    c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157          diagName  = 'SALTanom'
158          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
159          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
160        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
161        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
162       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
163    
164        diagName  = 'UVEL    '        diagName  = 'UVEL    '
165        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
166        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
167        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'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  = 'VVEL    '        diagName  = 'VVEL    '
172        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
173        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
174        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
175        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
176       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
177    
178        diagName  = 'WVEL    '        diagName  = 'WVEL    '
179        diagTitle = 'Vertical Component of Velocity (m/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
180        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
181        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
182        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
183       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185        diagName  = 'THETASQ '        diagName  = 'THETASQ '
186        diagTitle = 'Square of Potential Temperature (K^2)'        diagTitle = 'Square of Potential Temperature (K^2)'
187        diagUnits = 'K2              '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
188        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
189        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
190       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
191    
192        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
193        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)'
194        diagUnits = 'g2/kg2          '        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
195          diagCode  = 'SM      MR      '
196          CALL DIAGNOSTICS_ADD2LIST( diagNum,
197         I          diagName, diagCode, diagUnits, diagTitle, myThid )
198    
199          diagName  = 'SALTSQan'
200          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
201          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
202        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
203        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
204       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
205    
206        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
207        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
208        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
209        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
210        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
211       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
212    
213        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
214        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
215        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
216        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
217        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
218       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
219    
220        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
221        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
222        diagUnits = 'm2/sec2         '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
223        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
224        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
225       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
226    
227        diagName  = 'UVELVVEL'        diagName  = 'UV_VEL_C'
228          diagTitle ='Product of horizontal Comp of velocity (cell center)'
229          diagUnits = 'm^2/s^2         '
230          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
231          CALL DIAGNOSTICS_ADD2LIST( diagNum,
232         I          diagName, diagCode, diagUnits, diagTitle, myThid )
233    
234          diagName  = 'UV_VEL_Z'
235        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
236        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
237        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
238        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
239       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
240    
241          diagName  = 'WU_VEL  '
242          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
243          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
244          diagCode  = 'WU      LR      '
245          CALL DIAGNOSTICS_ADD2LIST( diagNum,
246         I          diagName, diagCode, diagUnits, diagTitle, myThid )
247    
248          diagName  = 'WV_VEL  '
249          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
250          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
251          diagCode  = 'WV      LR      '
252          CALL DIAGNOSTICS_ADD2LIST( diagNum,
253         I          diagName, diagCode, diagUnits, diagTitle, myThid )
254    
255        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
256        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
257        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
258        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
259        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
260       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
261    
262        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
263        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
264        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
265        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
266        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
267       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
268    
269        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
270        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
271        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
272        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
273        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
274       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
275    
276        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
277        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'
278        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
279        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
280        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
281       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
282    
283        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
284        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'
285        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
286        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
287        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
288       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
289    
290        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
291        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
292        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
293        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
294        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
295       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
296    
297        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
298        diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Mass-Weight Transp of Salt '
299        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
300        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
301        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
302       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
303    
304        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
305        diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Mass-Weight Transp of Salt'
306        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
307        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
308        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
309       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
310    
311        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
312        diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Mass-Weight Transp of Salt '
313        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
314        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
315        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
316       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
317    
318        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
319        diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
320        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
321        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
322        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
323       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
324    
325        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
326        diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
327        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
328        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
329        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
330       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
331    
332        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
333        diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
334        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
335        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
336        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
337       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
338    
339        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
340        diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Transp of Salt '
341        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
342        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
343        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
344       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
345    
346        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
347        diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Transp of Salt '
348        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
349        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
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  = 'WVELSLT '        diagName  = 'WVELSLT '
354        diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Transp of Salt '
355        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
356        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
357        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
358       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
361    
362          diagName  = 'RHOAnoma'
363          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
364          diagUnits = 'kg/m^3          '
365          diagCode  = 'SM      MR      '
366          CALL DIAGNOSTICS_ADD2LIST( diagNum,
367         I          diagName, diagCode, diagUnits, diagTitle, myThid )
368    
369          diagName  = 'RHOANOSQ'
370          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
371          diagUnits = 'kg^2/m^6          '
372          diagCode  = 'SM      MR      '
373          CALL DIAGNOSTICS_ADD2LIST( diagNum,
374         I          diagName, diagCode, diagUnits, diagTitle, myThid )
375    
376          diagName  = 'URHOMASS'
377          diagTitle = 'Zonal Transport of Density'
378          diagUnits = 'kg/m^2/s        '
379          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
380          CALL DIAGNOSTICS_ADD2LIST( diagNum,
381         I          diagName, diagCode, diagUnits, diagTitle, myThid )
382    
383          diagName  = 'VRHOMASS'
384          diagTitle = 'Meridional Transport of Density'
385          diagUnits = 'kg/m^2/s        '
386          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
387          CALL DIAGNOSTICS_ADD2LIST( diagNum,
388         I          diagName, diagCode, diagUnits, diagTitle, myThid )
389    
390          diagName  = 'WRHOMASS'
391          diagTitle = 'Vertical Transport of Potential Density'
392          diagUnits = 'kg/m^2/s        '
393          diagCode  = 'WM      LR      '
394          CALL DIAGNOSTICS_ADD2LIST( diagNum,
395         I          diagName, diagCode, diagUnits, diagTitle, myThid )
396    
397        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
398        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
399        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
400          diagCode  = 'SM      MR      '
401          CALL DIAGNOSTICS_ADD2LIST( diagNum,
402         I          diagName, diagCode, diagUnits, diagTitle, myThid )
403    
404          diagName  = 'PHIHYDSQ'
405          diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'
406          diagUnits = 'm^4/s^4         '
407        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
408        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
409       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
410    
411        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
412        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
413        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
414        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
415        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
416       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
417    
418        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
419        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
420        diagUnits = 'm4/s4           '        diagUnits = 'm^4/s^4         '
421        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
422        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
423       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
424    
425        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
426        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
427        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
428          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
429        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
430        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
431       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
432    
433  C--   Momentum terms:        diagName  = 'CONVADJ '
434        diagName  = 'VISCA4  '        diagTitle = 'Convective Adjustment Index [0-1] '
435        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '        diagUnits = 'fraction        '
436        diagUnits = 'm4/sec          '        diagCode  = 'SM      LR      '
       diagCode  = 'SM      MR      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
   
       diagName  = 'VISCAH  '  
       diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '  
       diagUnits = 'm2/sec          '  
       diagCode  = 'SM      MR      '  
437        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
438       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
439    
440  C--   surface fluxes:  C--   surface fluxes:
441        diagName  = 'TAUX    '        diagName  = 'TAUX    '
442        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
443        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
444        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
445        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
446       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
447    
448        diagName  = 'TAUY    '        diagName  = 'TAUY    '
449        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel '
450        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
451        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
452        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
453       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
454    
455        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
456        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'net surface heat flux, >0 increases theta       '
457        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
458        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
459        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
460       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
461    
462        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
463        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
464        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
465        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
466        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
467       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
468    
469        diagName  = 'TICE    '        diagName  = 'TICE    '
470        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
471        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
472        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
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  = 'SFLUX   '        diagName  = 'SFLUX   '
477        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt        '
478        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
479        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
480        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
481       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
482    
483        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
484        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt  '
485        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
486          diagCode  = 'SM      U1      '
487          CALL DIAGNOSTICS_ADD2LIST( diagNum,
488         I          diagName, diagCode, diagUnits, diagTitle, myThid )
489    
490    #ifdef SHORTWAVE_HEATING
491          diagName  = 'SWFLUX  '
492          diagTitle = 'net upward SW radiation, >0 increases theta     '
493          diagUnits = 'W/m^2           '
494        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
495        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
496       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
497    #endif
498    
499  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
500          diagName  = 'PRESSURE'
501          diagTitle = 'Atmospheric Pressure (Pa)     '
502          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
503          diagCode  = 'SM      MR      '
504          CALL DIAGNOSTICS_ADD2LIST( diagNum,
505         I          diagName, diagCode, diagUnits, diagTitle, myThid )
506    
507          diagName  = 'RSURF   '
508          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
509          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
510          diagCode  = 'SM      M1      '
511          CALL DIAGNOSTICS_ADD2LIST( diagNum,
512         I          diagName, diagCode, diagUnits, diagTitle, myThid )
513    
514    
515        RETURN        RETURN
516        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22