/[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.1 by jmc, Mon Dec 13 21:43:54 2004 UTC revision 1.50 by jmc, Sat Jul 19 16:46:25 2014 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  \\
33  C                   &  Z  &  C-grid Corner Point  \\  C                   &  Z  &  C-grid Corner Point  \\
34  C         parse(3)  &  --  &  --Unused--  \\  C         parse(3)  &     &  Used for Level Integrated output: cumulate levels \\
35    C                   &  r  &  same but cumulate product by model level thickness \\
36    C                   &  R  &  same but cumulate product by hFac & level thickness \\
37  C         parse(4)  &  P  &  positive definite  \\  C         parse(4)  &  P  &  positive definite  \\
38  C         parse(5 ) &  C  &  with counter array  \\  C         parse(5 ) &  C  &  with counter array  \\
39    C                   &  P  &  post-processed (not filled up) from other diags  \\
40  C                   &  D  &  disable an array for output  \\  C                   &  D  &  disable an array for output  \\
41  C         parse(6--8) & '123'  &  3-digit mate number \\  C         parse(6--8) & '123'  &  retired, formerly: 3-digit mate number \\
42  C         parse(9)  &  U  &  model-level plus 1/2  \\  C         parse(9)  &  U  &  model-level plus 1/2  \\
43  C                   &  M  &  model-level middle  \\  C                   &  M  &  model-level middle  \\
44  C                   &  L  &  model-level minus 1/2  \\  C                   &  L  &  model-level minus 1/2  \\
# Line 43  C                   &  1  &  levels = 1 Line 47  C                   &  1  &  levels = 1
47  C                   &  R  &  levels = Nr  \\  C                   &  R  &  levels = Nr  \\
48  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\
49  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
50    C                   &  G  &  levels = Ground_level Number \\
51    C                   &  I  &  levels = sea-Ice_level Number \\
52    C                   &  X  &  free levels option (need to be set explicitly) \\
53  C       \end{tabular}  C       \end{tabular}
54  C     \end{center}  C     \end{center}
55    
# Line 50  C     !USES: Line 57  C     !USES:
57        IMPLICIT NONE        IMPLICIT NONE
58  #include "SIZE.h"  #include "SIZE.h"
59  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 c #include "EESUPPORT.h"  
60  #include "PARAMS.h"  #include "PARAMS.h"
61    
62  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 58  C     !INPUT PARAMETERS: Line 64  C     !INPUT PARAMETERS:
64  CEOP  CEOP
65    
66  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
67        CHARACTER*8   diagName  C     rTitle     :: r-coordinate title
68        CHARACTER*16  diagCode  C     eTitle     :: free-surface title
69        CHARACTER*16  diagUnits  C     fTitle     :: fixed boundary title
70    C     pTitle     :: "Phi"  title
71    C     sTitle     :: "salt" title
72          INTEGER        diagNum
73          INTEGER        diagMate
74          CHARACTER*8    diagName
75          CHARACTER*16   diagCode
76          CHARACTER*16   diagUnits
77        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
78          CHARACTER*2    rUnit2c
79          CHARACTER*4    tUnit4c
80          CHARACTER*5    sUnit5c
81          CHARACTER*(10) rTitle, eTitle, fTitle
82          CHARACTER*(20) pTitle, sTitle
83    
84          CHARACTER*(16) DIAGS_MK_UNITS
85          EXTERNAL DIAGS_MK_UNITS
86          CHARACTER*(80) DIAGS_MK_TITLE
87          EXTERNAL DIAGS_MK_TITLE
88    
89  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
90  C     For each output variable,  C     For each output variable,
91  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
92  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
93  C----------------------------------------------------------------------  C----------------------------------------------------------------------
94    
95        diagName  = 'DRHODR  '        IF ( usingPCoords ) THEN
96        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '          rUnit2c= 'Pa'
97        diagUnits = 'kg/m4           '          rTitle = ' Pressure '
98        diagCode  = 'SM      LR      '          pTitle = ' Geopotential       '
99        CALL DIAGNOSTICS_ADD2LIST(        ELSE
100       I          diagName, diagCode, diagUnits, diagTitle, myThid )          rUnit2c= 'm '
101            rTitle = ' Height   '
102            pTitle = 'Pressure Pot.(p/rho)'
103          ENDIF
104          IF ( fluidIsAir ) THEN
105            tUnit4c= 'K   '
106            sUnit5c= 'kg/kg'
107            sTitle = ' Specific Humidity  '
108            IF (useAIM) sUnit5c= 'g/kg '
109          ELSEIF ( eosType.EQ.'TEOS10' ) THEN
110            tUnit4c= 'degC'
111            sUnit5c= 'g/kg '
112    c       tTitle = 'Conservative Temp.   '
113            sTitle = ' Absolute Salinity  '
114          ELSE
115            tUnit4c= 'degC'
116            sUnit5c= 'psu  '
117    c       tTitle = 'Potential Temperature'
118            sTitle = ' Salinity           '
119          ENDIF
120    C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
121          IF ( fluidIsAir ) THEN
122            eTitle = ' Surface  '
123            fTitle = ' Top      '
124          ELSEIF ( usingPCoords ) THEN
125            eTitle = ' Bottom   '
126            fTitle = ' Surface  '
127          ELSE
128            eTitle = ' Surface  '
129            fTitle = ' Bottom   '
130          ENDIF
131    
132    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
133    C-    state variables of the main code (and related quadratic var):
134    
135          diagName  = 'ETAN    '
136          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
137    c     IF ( fluidIsWater .AND. usingZCoords )
138    c    &diagTitle = 'Sea Surface Elevation'
139          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
140          diagCode  = 'SM      M1      '
141          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
142         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
143    
144          diagName  = 'ETANSQ  '
145          diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
146         I                          //' Anomaly', myThid )
147          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
148          diagCode  = 'SM P    M1      '
149          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
150         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
151    
152        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
153        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
154        diagUnits = 'm2 or Pa2 /sec2 '       I                          //' Anomaly Tendency', myThid )
155          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
156          diagCode  = 'SM      M1      '
157          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
158         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
159    
160          diagName  = 'THETA   '
161          diagTitle = 'Potential Temperature'
162          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
163          diagCode  = 'SMR     MR      '
164          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
165         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
166    
167    c     diagName  = 'SST     '
168    c     diagTitle = 'Sea Surface Temperature (degC,K)'
169    c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
170    c     diagCode  = 'SM      M1      '
171    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
172    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
173    
174          diagName  = 'SALT    '
175          diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
176          diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
177          diagCode  = 'SMR     MR      '
178          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
179         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
180    
181          diagName  = 'RELHUM  '
182          diagTitle = 'Relative Humidity'
183          diagUnits = 'percent         '
184          diagCode  = 'SMR     MR      '
185          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
186         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
187    
188    c     diagName  = 'SSS     '
189    c     diagTitle = 'Sea Surface Salinity '
190    c     diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
191    c     diagCode  = 'SM      M1      '
192    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
193    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
194    
195          IF ( fluidIsWater ) THEN
196          diagName  = 'SALTanom'
197          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
198          diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
199          diagCode  = 'SMR     MR      '
200          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
201         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
202          ENDIF
203    
204          diagName  = 'UVEL    '
205          diagTitle = 'Zonal Component of Velocity (m/s)'
206          diagUnits = 'm/s             '
207          diagCode  = 'UUR     MR      '
208          diagMate  = diagNum + 2
209          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
210         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
211    
212          diagName  = 'VVEL    '
213          diagTitle = 'Meridional Component of Velocity (m/s)'
214          diagUnits = 'm/s             '
215          diagCode  = 'VVR     MR      '
216          diagMate  = diagNum
217          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
218         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
219    
220          diagName  = 'WVEL    '
221          diagTitle = 'Vertical Component of Velocity (r_units/s)'
222          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
223          diagCode  = 'WM      LR      '
224          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
225         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
226    
227          diagName  = 'THETASQ '
228          diagTitle = 'Square of Potential Temperature'
229          diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
230          diagCode  = 'SMRP    MR      '
231          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
232         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
233    
234          diagName  = 'SALTSQ  '
235          diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
236          diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
237          diagCode  = 'SMRP    MR      '
238          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
239         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
240    
241          IF ( fluidIsWater ) THEN
242          diagName  = 'SALTSQan'
243          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
244          diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
245          diagCode  = 'SMRP    MR      '
246          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
247         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
248          ENDIF
249    
250          diagName  = 'UVELSQ  '
251          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
252          diagUnits = 'm^2/s^2         '
253          diagCode  = 'UURP    MR      '
254          diagMate  = diagNum + 2
255          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
256         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
257    
258          diagName  = 'VVELSQ  '
259          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
260          diagUnits = 'm^2/s^2         '
261          diagCode  = 'VVRP    MR      '
262          diagMate  = diagNum
263          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
264         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
265    
266          diagName  = 'WVELSQ  '
267          diagTitle = 'Square of Vertical Comp of Velocity'
268          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
269          diagCode  = 'WM P    LR      '
270          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
271         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
272    
273          diagName  = 'UE_VEL_C'
274          diagTitle = 'Eastward Velocity (m/s) (cell center)'
275          diagUnits = 'm/s             '
276          diagCode  = 'UMR     MR      '
277          diagMate  = diagNum + 2
278          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
279         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
280    
281          diagName  = 'VN_VEL_C'
282          diagTitle = 'Northward Velocity (m/s) (cell center)'
283          diagUnits = 'm/s             '
284          diagCode  = 'VMR     MR      '
285          diagMate  = diagNum
286          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
287         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
288    
289          diagName  = 'UV_VEL_C'
290          diagTitle ='Product of horizontal Comp of velocity (cell center)'
291          diagUnits = 'm^2/s^2         '
292          diagCode  = 'UMR     MR      '
293          diagMate  = diagNum + 1
294          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
295         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
296    
297          diagName  = 'UV_VEL_Z'
298          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
299          diagUnits = 'm^2/s^2         '
300          diagCode  = 'UZR     MR      '
301          diagMate  = diagNum + 1
302          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
303         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
304    
305          diagName  = 'WU_VEL  '
306          diagTitle = 'Vertical Transport of Zonal Momentum'
307          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
308          diagCode  = 'WU      LR      '
309          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
310         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
311    
312          diagName  = 'WV_VEL  '
313          diagTitle ='Vertical Transport of Meridional Momentum'
314          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
315          diagCode  = 'WV      LR      '
316          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
317         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
318    
319          diagName  = 'UVELMASS'
320          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
321          diagUnits = 'm/s             '
322          diagCode  = 'UUr     MR      '
323          diagMate  = diagNum + 2
324          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
325         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
326    
327          diagName  = 'VVELMASS'
328          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
329          diagUnits = 'm/s             '
330          diagCode  = 'VVr     MR      '
331          diagMate  = diagNum
332          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
333         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
334    
335          diagName  = 'WVELMASS'
336          diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
337          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
338          diagCode  = 'WM      LR      '
339          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
340         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
341    
342          diagName  = 'PhiVEL  '
343          diagTitle = 'Horizontal Velocity Potential (m^2/s)'
344          diagUnits = 'm^2/s           '
345          diagCode  = 'SMR P   MR      '
346    C-    use 'UVELMASS' as mate.
347          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
348         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
349    
350          diagName  = 'PsiVEL  '
351          diagTitle = 'Horizontal Velocity Stream-Function'
352          diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
353          diagCode  = 'SZ  P   MR      '
354    C-    use 'PhiVEL' as mate.
355          diagMate  = diagNum
356          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
357         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
358    
359          diagName  = 'UTHMASS '
360          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
361          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
362          diagCode  = 'UUr     MR      '
363          diagMate  = diagNum + 2
364          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
365         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
366    
367          diagName  = 'VTHMASS '
368          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
369          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
370          diagCode  = 'VVr     MR      '
371          diagMate  = diagNum
372          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
373         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
374    
375          diagName  = 'WTHMASS '
376          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
377          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
378          diagCode  = 'WM      LR      '
379          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
380         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
381    
382          diagName  = 'USLTMASS'
383          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
384         I                           //sTitle, myThid )
385          diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
386          diagCode  = 'UUr     MR      '
387          diagMate  = diagNum + 2
388          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
389         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
390    
391          diagName  = 'VSLTMASS'
392          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
393         I                           //sTitle, myThid )
394          diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
395          diagCode  = 'VVr     MR      '
396          diagMate  = diagNum
397          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
398         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
399    
400          diagName  = 'WSLTMASS'
401          diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
402         I                           //sTitle, myThid )
403          diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
404          diagCode  = 'WM      LR      '
405          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
406         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
407    
408          diagName  = 'UVELTH  '
409          diagTitle = 'Zonal Transport of Pot Temp'
410          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
411          diagCode  = 'UUR     MR      '
412          diagMate  = diagNum + 2
413          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
414         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
415    
416          diagName  = 'VVELTH  '
417          diagTitle = 'Meridional Transport of Pot Temp'
418          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
419          diagCode  = 'VVR     MR      '
420          diagMate  = diagNum
421          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
422         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
423    
424          diagName  = 'WVELTH  '
425          diagTitle = 'Vertical Transport of Pot Temp'
426          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
427          diagCode  = 'WM      LR      '
428          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
429         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
430    
431          diagName  = 'UVELSLT '
432          diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
433         I                          //sTitle, myThid )
434          diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
435          diagCode  = 'UUR     MR      '
436          diagMate  = diagNum + 2
437          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
438         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
439    
440          diagName  = 'VVELSLT '
441          diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
442         I                          //sTitle, myThid )
443          diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
444          diagCode  = 'VVR     MR      '
445          diagMate  = diagNum
446          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
447         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
448    
449          diagName  = 'WVELSLT '
450          diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
451         I                          //sTitle, myThid )
452          diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
453          diagCode  = 'WM      LR      '
454          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
455         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
456    
457          diagName  = 'UVELPHI '
458          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
459         I                 //pTitle//' Anomaly', myThid )
460          diagUnits = 'm^3/s^3         '
461          diagCode  = 'UUr     MR      '
462          diagMate  = diagNum + 2
463          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
464         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
465    
466          diagName  = 'VVELPHI '
467          diagTitle = DIAGS_MK_TITLE( 'Merid. Mass-Weight Transp of '
468         I                 //pTitle//' Anomaly', myThid )
469          diagUnits = 'm^3/s^3         '
470          diagCode  = 'VVr     MR      '
471          diagMate  = diagNum
472          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
473         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
474    
475    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
476    
477          diagName  = 'RHOAnoma'
478          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
479          diagUnits = 'kg/m^3          '
480          diagCode  = 'SMR     MR      '
481          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
482         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
483    
484          diagName  = 'RHOANOSQ'
485          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
486          diagUnits = 'kg^2/m^6        '
487          diagCode  = 'SMRP    MR      '
488          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
489         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
490    
491          diagName  = 'URHOMASS'
492          diagTitle = 'Zonal Transport of Density'
493          diagUnits = 'kg/m^2/s        '
494          diagCode  = 'UUr     MR      '
495          diagMate  = diagNum + 2
496          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
497         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
498    
499          diagName  = 'VRHOMASS'
500          diagTitle = 'Meridional Transport of Density'
501          diagUnits = 'kg/m^2/s        '
502          diagCode  = 'VVr     MR      '
503          diagMate  = diagNum
504          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
505         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
506    
507          diagName  = 'WRHOMASS'
508          diagTitle = 'Vertical Transport of Density'
509          diagUnits = 'kg/m^2/s        '
510          diagCode  = 'WM      LR      '
511          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
512         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
513    
514          diagName  = 'WdRHO_P '
515          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
516          diagUnits = 'kg/m^2/s        '
517          diagCode  = 'WM      LR      '
518          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
519         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
520    
521          diagName  = 'WdRHOdP '
522          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
523          diagUnits = 'kg/m^2/s        '
524          diagCode  = 'WM      LR      '
525          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
526         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
527    
528          diagName  = 'PHIHYD  '
529          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
530         I                           //pTitle//' Anomaly', myThid )
531          diagUnits = 'm^2/s^2         '
532          diagCode  = 'SMR     MR      '
533          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
534         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
535    
536          diagName  = 'PHIHYDSQ'
537          diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
538         I                           //pTitle//' Anomaly', myThid )
539          diagUnits = 'm^4/s^4         '
540          diagCode  = 'SMRP    MR      '
541          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
542         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
543    
544          diagName  = 'PHIBOT  '
545    c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
546          diagTitle = DIAGS_MK_TITLE( fTitle
547         I                           //pTitle//' Anomaly', myThid )
548          diagUnits = 'm^2/s^2         '
549        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
550        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
551       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
552    
553          diagName  = 'PHIBOTSQ'
554    c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
555          diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
556         I                           //pTitle//' Anomaly', myThid )
557          diagUnits = 'm^4/s^4         '
558          diagCode  = 'SM P    M1      '
559          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
560         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
561    
562    #ifdef NONLIN_FRSURF
563          diagName  = 'PHIHYDcR'
564          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
565         I                       //pTitle//' Anomaly @ const r', myThid )
566          diagUnits = 'm^2/s^2         '
567          diagCode  = 'SMR     MR      '
568          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
569         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
570    #endif
571    
572    #ifdef ALLOW_NONHYDROSTATIC
573          diagName  = 'PHI_NH  '
574          diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
575          diagUnits = 'm^2/s^2         '
576          diagCode  = 'SMR     MR      '
577          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
578         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
579    #endif /* ALLOW_NONHYDROSTATIC */
580    
581          diagName  = 'MXLDEPTH'
582          diagTitle = 'Mixed-Layer Depth (>0)'
583          diagUnits = 'm               '
584          diagCode  = 'SM      M1      '
585          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
586         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
587    
588          diagName  = 'DRHODR  '
589          diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
590          diagUnits = 'kg/m^4          '
591          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
592          diagCode  = 'SM      LR      '
593          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
594         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
595    
596          diagName  = 'CONVADJ '
597          diagTitle = 'Convective Adjustment Index [0-1] '
598          diagUnits = 'fraction        '
599          diagCode  = 'SMR     LR      '
600          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
601         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
602    
603    C--   surface fluxes:
604          diagName  = 'oceTAUX '
605          diagTitle = 'zonal surface wind stress, >0 increases uVel'
606          diagUnits = 'N/m^2           '
607          diagCode  = 'UU      U1      '
608          diagMate  = diagNum + 2
609          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
610         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
611    
612          diagName  = 'oceTAUY '
613          diagTitle = 'meridional surf. wind stress, >0 increases vVel'
614          diagUnits = 'N/m^2           '
615          diagCode  = 'VV      U1      '
616          diagMate  = diagNum
617          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
618         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
619    
620          diagName  = 'atmPload'
621          diagTitle = 'Atmospheric pressure loading'
622          diagUnits = 'Pa              '
623          diagCode  = 'SM      U1      '
624          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
625         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
626    
627          diagName  = 'sIceLoad'
628          diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
629          diagUnits = 'kg/m^2          '
630          diagCode  = 'SM      U1      '
631          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
632         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
633    
634          diagName  = 'oceFWflx'
635          diagTitle = 'net surface Fresh-Water flux into the ocean'
636         &          //' (+=down), >0 decreases salinity'
637          diagUnits = 'kg/m^2/s        '
638          diagCode  = 'SM      U1      '
639          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
640         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
641    
642          diagName  = 'oceSflux'
643          diagTitle = 'net surface Salt flux into the ocean (+=down),'
644         &          //' >0 increases salinity'
645          diagUnits = 'g/m^2/s         '
646          diagCode  = 'SM      U1      '
647          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
648         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
649    
650          diagName  = 'oceQnet '
651          diagTitle = 'net surface heat flux into the ocean (+=down),'
652         &          //' >0 increases theta'
653          diagUnits = 'W/m^2           '
654          diagCode  = 'SM      U1      '
655          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
656         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
657    
658          diagName  = 'oceQsw  '
659          diagTitle = 'net Short-Wave radiation (+=down),'
660         &          //' >0 increases theta'
661          diagUnits = 'W/m^2           '
662          diagCode  = 'SM      U1      '
663          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
664         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
665    
666          diagName  = 'oceFreez'
667          diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
668          diagUnits = 'W/m^2           '
669          diagCode  = 'SM      U1      '
670          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
671         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
672    
673          diagName  = 'TRELAX  '
674          diagTitle = 'surface temperature relaxation, >0 increases theta'
675          diagUnits = 'W/m^2           '
676          diagCode  = 'SM      U1      '
677          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
678         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
679    
680          diagName  = 'SRELAX  '
681          diagTitle = 'surface salinity relaxation, >0 increases salt'
682          diagUnits = 'g/m^2/s         '
683          diagCode  = 'SM      U1      '
684          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
685         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
686    
687          diagName  = 'surForcT'
688          diagTitle = 'model surface forcing for Temperature,'
689         &          //' >0 increases theta'
690          diagUnits = 'W/m^2           '
691          diagCode  = 'SM      U1      '
692          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
693         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
694    
695          diagName  = 'surForcS'
696          diagTitle = 'model surface forcing for Salinity,'
697         &          //' >0 increases salinity'
698          diagUnits = 'g/m^2/s         '
699          diagCode  = 'SM      U1      '
700          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
701         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
702    
703          diagName  = 'TFLUX   '
704          diagTitle = 'total heat flux (match heat-content variations),'
705         &          //' >0 increases theta'
706          diagUnits = 'W/m^2           '
707          diagCode  = 'SM      U1      '
708          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
709         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
710    
711          diagName  = 'SFLUX   '
712          diagTitle = 'total salt flux (match salt-content variations),'
713         &          //' >0 increases salt'
714          diagUnits = 'g/m^2/s         '
715          diagCode  = 'SM      U1      '
716          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
717         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
718    
719  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
720    
721          diagName  = 'RCENTER '
722    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
723          diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
724          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
725          diagCode  = 'SM      MR      '
726          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
727         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
728    
729          diagName  = 'RSURF   '
730    c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
731          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
732          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
733          diagCode  = 'SM      M1      '
734          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
735         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
736    
737          diagName  = 'TOTUTEND'
738          diagTitle = 'Tendency of Zonal Component of Velocity'
739          diagUnits = 'm/s/day         '
740          diagCode  = 'UUR     MR      '
741          diagMate  = diagNum + 2
742          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
743         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
744    
745          diagName  = 'TOTVTEND'
746          diagTitle = 'Tendency of Meridional Component of Velocity'
747          diagUnits = 'm/s/day         '
748          diagCode  = 'VVR     MR      '
749          diagMate  = diagNum
750          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
751         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
752    
753          diagName  = 'TOTTTEND'
754          diagTitle = 'Tendency of Potential Temperature'
755          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
756          diagCode  = 'SMR     MR      '
757          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
758         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
759    
760          diagName  = 'TOTSTEND'
761          diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
762          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/day', myThid )
763          diagCode  = 'SMR     MR      '
764          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
765         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
766    
767          diagName  = 'MoistCor'
768          diagTitle = 'Heating correction due to moist thermodynamics'
769          diagUnits = 'W/m^2           '
770          diagCode  = 'SM      MR      '
771          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
772         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
773    
774    #ifdef ALLOW_GENERIC_ADVDIFF
775          diagName  = 'gT_Forc '
776          diagTitle = 'Potential Temp. forcing tendency'
777          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
778          diagCode  = 'SMR     MR      '
779          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
780         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
781    
782          diagName  = 'gS_Forc '
783          diagTitle = DIAGS_MK_TITLE(
784         &            sTitle//'forcing tendency', myThid )
785          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
786          diagCode  = 'SMR     MR      '
787          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
788         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
789    
790          diagName  = 'AB_gT   '
791          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
792          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
793          diagCode  = 'SMR     MR      '
794          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
795         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
796    
797          diagName  = 'AB_gS   '
798          diagTitle = DIAGS_MK_TITLE(
799         &            sTitle//'tendency from Adams-Bashforth', myThid )
800          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
801          diagCode  = 'SMR     MR      '
802          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
803         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
804    
805          diagName  = 'gTinAB  '
806          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
807          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
808          diagCode  = 'SMR     MR      '
809          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
810         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
811    
812          diagName  = 'gSinAB  '
813          diagTitle = DIAGS_MK_TITLE(
814         &            sTitle//'tendency going in Adams-Bashforth', myThid )
815          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
816          diagCode  = 'SMR     MR      '
817          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
818         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
819    #endif /* ALLOW_GENERIC_ADVDIFF */
820    
821          diagName  = 'AB_gU   '
822          diagTitle = 'U momentum tendency from Adams-Bashforth'
823          diagUnits = 'm/s^2           '
824          diagCode  = 'UUR     MR      '
825          diagMate  = diagNum + 2
826          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
827         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
828    
829          diagName  = 'AB_gV   '
830          diagTitle = 'V momentum tendency from Adams-Bashforth'
831          diagUnits = 'm/s^2           '
832          diagCode  = 'VVR     MR      '
833          diagMate  = diagNum
834          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
835         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
836    
837    #ifdef ALLOW_NONHYDROSTATIC
838          diagName  = 'AB_gW   '
839          diagTitle = 'W momentum tendency from Adams-Bashforth'
840          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
841          diagCode  = 'WM      LR      '
842          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
843         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
844    #endif /* ALLOW_NONHYDROSTATIC */
845    
846    #ifdef ALLOW_EDDYPSI
847          diagName  = 'TAUXEDDY'
848          diagTitle = 'Zonal Eddy Stress'
849          diagUnits = 'N/m**2          '
850          diagCode  = 'UU P    MR      '
851          diagMate  = diagNum
852          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
853         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
854    
855          diagName  = 'TAUYEDDY'
856          diagTitle = 'Meridional Eddy Stress'
857          diagUnits = 'N/m**2          '
858          diagCode  = 'VV P    MR      '
859          diagMate  = diagNum
860          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
861         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
862    
863    # ifdef ALLOW_GMREDI
864          diagName  = 'UVELMEAN'
865          diagTitle = 'Zonal Component of Mean Velocity (m/s)'
866          diagUnits = 'm/s             '
867          diagCode  = 'UUR     MR      '
868          diagMate  = diagNum + 2
869          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
870         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
871    
872          diagName  = 'VVELMEAN'
873          diagTitle = 'Meridional Component of Mean Velocity (m/s)'
874          diagUnits = 'm/s             '
875          diagCode  = 'VVR     MR      '
876          diagMate  = diagNum
877          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
878         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
879    # endif /* ALLOW_GMREDI */
880    #endif /* ALLOW_EDDYPSI */
881    
882        RETURN        RETURN
883        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22