/[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.31 by jmc, Thu May 31 18:22:00 2007 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 58  C     !INPUT PARAMETERS: Line 60  C     !INPUT PARAMETERS:
60  CEOP  CEOP
61    
62  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
63        CHARACTER*8   diagName  C     rTitle     :: r-coordinate title
64        CHARACTER*16  diagCode  C     eTitle     :: free-surface title
65        CHARACTER*16  diagUnits  C     fTitle     :: fixed boundary title
66    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
74          CHARACTER*4    tUnit4c, sUnit4c
75          CHARACTER*(10) rTitle, eTitle, fTitle
76          CHARACTER*(20) pTitle, sTitle
77    
78          CHARACTER*(16) DIAGS_MK_UNITS
79          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        diagName  = 'DRHODR  '        IF ( usingPCoords ) THEN
90        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '          rUnit2c= 'Pa'
91        diagUnits = 'kg/m4           '          rTitle = ' Pressure '
92        diagCode  = 'SM      LR      '          pTitle = ' Geopotential       '
93        CALL DIAGNOSTICS_ADD2LIST(        ELSE
94            rUnit2c= 'm '
95            rTitle = ' Height   '
96            pTitle = 'Pressure Pot.(p/rho)'
97          ENDIF
98          IF ( fluidIsAir ) THEN
99            tUnit4c= 'K   '
100            sUnit4c= 'g/kg'
101            sTitle = ' Specific Humidity  '
102    #ifdef ALLOW_FIZHI
103            IF (useFIZHI) sUnit4c= 'kg/kg'
104    #endif /* ALLOW_FIZHI */
105          ELSE
106            tUnit4c= 'degC'
107            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
121    
122    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
123    C-    state variables of the main code (and related quadratic var):
124    
125          diagName  = 'ETAN    '
126          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 )
130          diagCode  = 'SM      M1      '
131          CALL DIAGNOSTICS_ADD2LIST( diagNum,
132         I          diagName, diagCode, diagUnits, diagTitle, myThid )
133    
134          diagName  = 'ETANSQ  '
135          diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
136         I                          //' Anomaly', myThid )
137          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
138          diagCode  = 'SM      M1      '
139          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        diagUnits = 'm2 or Pa2 /sec2 '       I                          //' Anomaly Tendency', myThid )
145          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
146          diagCode  = 'SM      M1      '
147          CALL DIAGNOSTICS_ADD2LIST( diagNum,
148         I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150          diagName  = 'THETA   '
151          diagTitle = 'Potential Temperature'
152          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
153          diagCode  = 'SM      MR      '
154          CALL DIAGNOSTICS_ADD2LIST( diagNum,
155         I          diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157    c     diagName  = 'SST     '
158    c     diagTitle = 'Sea Surface Temperature (degC,K)'
159    c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
160    c     diagCode  = 'SM      M1      '
161    c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
162    c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
163    
164          diagName  = 'SALT    '
165          diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
166          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
167          diagCode  = 'SM      MR      '
168          CALL DIAGNOSTICS_ADD2LIST( diagNum,
169         I          diagName, diagCode, diagUnits, diagTitle, myThid )
170    
171          diagName  = 'RELHUM  '
172          diagTitle = 'Relative Humidity'
173          diagUnits = 'percent         '
174          diagCode  = 'SM      MR      '
175          CALL DIAGNOSTICS_ADD2LIST( diagNum,
176         I          diagName, diagCode, diagUnits, diagTitle, myThid )
177    
178    c     diagName  = 'SSS     '
179    c     diagTitle = 'Sea Surface Salinity '
180    c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
181    c     diagCode  = 'SM      M1      '
182    c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
183    c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185          IF ( fluidIsWater ) THEN
186          diagName  = 'SALTanom'
187          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
188          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
189          diagCode  = 'SM      MR      '
190          CALL DIAGNOSTICS_ADD2LIST( diagNum,
191         I          diagName, diagCode, diagUnits, diagTitle, myThid )
192          ENDIF
193    
194          diagName  = 'UVEL    '
195          diagTitle = 'Zonal Component of Velocity (m/s)'
196          diagUnits = 'm/s             '
197          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
198          CALL DIAGNOSTICS_ADD2LIST( diagNum,
199         I          diagName, diagCode, diagUnits, diagTitle, myThid )
200    
201          diagName  = 'VVEL    '
202          diagTitle = 'Meridional Component of Velocity (m/s)'
203          diagUnits = 'm/s             '
204          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
205          CALL DIAGNOSTICS_ADD2LIST( diagNum,
206         I          diagName, diagCode, diagUnits, diagTitle, myThid )
207    
208          diagName  = 'WVEL    '
209          diagTitle = 'Vertical Component of Velocity (r_units/s)'
210          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
211          diagCode  = 'WM      LR      '
212          CALL DIAGNOSTICS_ADD2LIST( diagNum,
213         I          diagName, diagCode, diagUnits, diagTitle, myThid )
214    
215          diagName  = 'THETASQ '
216          diagTitle = 'Square of Potential Temperature'
217          diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
218          diagCode  = 'SM      MR      '
219          CALL DIAGNOSTICS_ADD2LIST( diagNum,
220         I          diagName, diagCode, diagUnits, diagTitle, myThid )
221    
222          diagName  = 'SALTSQ  '
223          diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
224          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
225          diagCode  = 'SM      MR      '
226          CALL DIAGNOSTICS_ADD2LIST( diagNum,
227         I          diagName, diagCode, diagUnits, diagTitle, myThid )
228    
229          IF ( fluidIsWater ) THEN
230          diagName  = 'SALTSQan'
231          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
232          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233          diagCode  = 'SM      MR      '
234          CALL DIAGNOSTICS_ADD2LIST( diagNum,
235         I          diagName, diagCode, diagUnits, diagTitle, myThid )
236          ENDIF
237    
238          diagName  = 'UVELSQ  '
239          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
240          diagUnits = 'm^2/s^2         '
241          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
242          CALL DIAGNOSTICS_ADD2LIST( diagNum,
243         I          diagName, diagCode, diagUnits, diagTitle, myThid )
244    
245          diagName  = 'VVELSQ  '
246          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
247          diagUnits = 'm^2/s^2         '
248          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
249          CALL DIAGNOSTICS_ADD2LIST( diagNum,
250         I          diagName, diagCode, diagUnits, diagTitle, myThid )
251    
252          diagName  = 'WVELSQ  '
253          diagTitle = 'Square of Vertical Comp of Velocity'
254          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
255          diagCode  = 'WM      LR      '
256          CALL DIAGNOSTICS_ADD2LIST( diagNum,
257         I          diagName, diagCode, diagUnits, diagTitle, myThid )
258    
259          diagName  = 'UV_VEL_C'
260          diagTitle ='Product of horizontal Comp of velocity (cell center)'
261          diagUnits = 'm^2/s^2         '
262          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
263          CALL DIAGNOSTICS_ADD2LIST( diagNum,
264         I          diagName, diagCode, diagUnits, diagTitle, myThid )
265    
266          diagName  = 'UV_VEL_Z'
267          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
268          diagUnits = 'm^2/s^2         '
269          WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
270          CALL DIAGNOSTICS_ADD2LIST( diagNum,
271         I          diagName, diagCode, diagUnits, diagTitle, myThid )
272    
273          diagName  = 'WU_VEL  '
274          diagTitle = 'Vertical Transport of Zonal Momentum'
275          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
276          diagCode  = 'WU      LR      '
277          CALL DIAGNOSTICS_ADD2LIST( diagNum,
278         I          diagName, diagCode, diagUnits, diagTitle, myThid )
279    
280          diagName  = 'WV_VEL  '
281          diagTitle ='Vertical Transport of Meridional Momentum'
282          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
283          diagCode  = 'WV      LR      '
284          CALL DIAGNOSTICS_ADD2LIST( diagNum,
285         I          diagName, diagCode, diagUnits, diagTitle, myThid )
286    
287          diagName  = 'UVELMASS'
288          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
289          diagUnits = 'm/s             '
290          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
291          CALL DIAGNOSTICS_ADD2LIST( diagNum,
292         I          diagName, diagCode, diagUnits, diagTitle, myThid )
293    
294          diagName  = 'VVELMASS'
295          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
296          diagUnits = 'm/s             '
297          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
298          CALL DIAGNOSTICS_ADD2LIST( diagNum,
299         I          diagName, diagCode, diagUnits, diagTitle, myThid )
300    
301          diagName  = 'WVELMASS'
302          diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
303          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
304          diagCode  = 'WM      LR      '
305          CALL DIAGNOSTICS_ADD2LIST( diagNum,
306         I          diagName, diagCode, diagUnits, diagTitle, myThid )
307    
308          diagName  = 'UTHMASS '
309          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
310          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
311          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
312          CALL DIAGNOSTICS_ADD2LIST( diagNum,
313         I          diagName, diagCode, diagUnits, diagTitle, myThid )
314    
315          diagName  = 'VTHMASS '
316          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
317          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
318          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
319          CALL DIAGNOSTICS_ADD2LIST( diagNum,
320         I          diagName, diagCode, diagUnits, diagTitle, myThid )
321    
322          diagName  = 'WTHMASS '
323          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
324          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
325          diagCode  = 'WM      LR      '
326          CALL DIAGNOSTICS_ADD2LIST( diagNum,
327         I          diagName, diagCode, diagUnits, diagTitle, myThid )
328    
329          diagName  = 'USLTMASS'
330          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
331         I                           //sTitle, myThid )
332          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
333          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
334          CALL DIAGNOSTICS_ADD2LIST( diagNum,
335         I          diagName, diagCode, diagUnits, diagTitle, myThid )
336    
337          diagName  = 'VSLTMASS'
338          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
339         I                           //sTitle, myThid )
340          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
341          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
342          CALL DIAGNOSTICS_ADD2LIST( diagNum,
343         I          diagName, diagCode, diagUnits, diagTitle, myThid )
344    
345          diagName  = 'WSLTMASS'
346          diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
347         I                           //sTitle, myThid )
348          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
349          diagCode  = 'WM      LR      '
350          CALL DIAGNOSTICS_ADD2LIST( diagNum,
351         I          diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353          diagName  = 'UVELTH  '
354          diagTitle = 'Zonal Transport of Pot Temp'
355          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
356          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
357          CALL DIAGNOSTICS_ADD2LIST( diagNum,
358         I          diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360          diagName  = 'VVELTH  '
361          diagTitle = 'Meridional Transport of Pot Temp'
362          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
363          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
364          CALL DIAGNOSTICS_ADD2LIST( diagNum,
365         I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367          diagName  = 'WVELTH  '
368          diagTitle = 'Vertical Transport of Pot Temp'
369          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
370          diagCode  = 'WM      LR      '
371          CALL DIAGNOSTICS_ADD2LIST( diagNum,
372         I          diagName, diagCode, diagUnits, diagTitle, myThid )
373    
374          diagName  = 'UVELSLT '
375          diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
376         I                          //sTitle, myThid )
377          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
378          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
379          CALL DIAGNOSTICS_ADD2LIST( diagNum,
380         I          diagName, diagCode, diagUnits, diagTitle, myThid )
381    
382          diagName  = 'VVELSLT '
383          diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
384         I                          //sTitle, myThid )
385          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
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  = 'WVELSLT '
391          diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
392         I                          //sTitle, myThid )
393          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
394          diagCode  = 'WM      LR      '
395          CALL DIAGNOSTICS_ADD2LIST( diagNum,
396         I          diagName, diagCode, diagUnits, diagTitle, myThid )
397    
398          diagName  = 'UVELPHI '
399          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
400         I                 //pTitle//' Anomaly', myThid )
401          diagUnits = 'm^3/s^3         '
402          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
403          CALL DIAGNOSTICS_ADD2LIST( diagNum,
404         I          diagName, diagCode, diagUnits, diagTitle, myThid )
405    
406          diagName  = 'VVELPHI '
407          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
408         I                 //pTitle//' Anomaly', myThid )
409          diagUnits = 'm^3/s^3         '
410          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
411          CALL DIAGNOSTICS_ADD2LIST( diagNum,
412         I          diagName, diagCode, diagUnits, diagTitle, myThid )
413    
414    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
415    
416          diagName  = 'RHOAnoma'
417          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
418          diagUnits = 'kg/m^3          '
419          diagCode  = 'SM      MR      '
420          CALL DIAGNOSTICS_ADD2LIST( diagNum,
421         I          diagName, diagCode, diagUnits, diagTitle, myThid )
422    
423          diagName  = 'RHOANOSQ'
424          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
425          diagUnits = 'kg^2/m^6          '
426          diagCode  = 'SM      MR      '
427          CALL DIAGNOSTICS_ADD2LIST( diagNum,
428         I          diagName, diagCode, diagUnits, diagTitle, myThid )
429    
430          diagName  = 'URHOMASS'
431          diagTitle = 'Zonal Transport of Density'
432          diagUnits = 'kg/m^2/s        '
433          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
434          CALL DIAGNOSTICS_ADD2LIST( diagNum,
435         I          diagName, diagCode, diagUnits, diagTitle, myThid )
436    
437          diagName  = 'VRHOMASS'
438          diagTitle = 'Meridional Transport of Density'
439          diagUnits = 'kg/m^2/s        '
440          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
441          CALL DIAGNOSTICS_ADD2LIST( diagNum,
442         I          diagName, diagCode, diagUnits, diagTitle, myThid )
443    
444          diagName  = 'WRHOMASS'
445          diagTitle = 'Vertical Transport of Potential Density'
446          diagUnits = 'kg/m^2/s        '
447          diagCode  = 'WM      LR      '
448          CALL DIAGNOSTICS_ADD2LIST( diagNum,
449         I          diagName, diagCode, diagUnits, diagTitle, myThid )
450    
451          diagName  = 'PHIHYD  '
452          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
453         I                           //pTitle//' Anomaly', myThid )
454          diagUnits = 'm^2/s^2         '
455          diagCode  = 'SM      MR      '
456          CALL DIAGNOSTICS_ADD2LIST( diagNum,
457         I          diagName, diagCode, diagUnits, diagTitle, myThid )
458    
459          diagName  = 'PHIHYDSQ'
460          diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
461         I                           //pTitle//' Anomaly', myThid )
462          diagUnits = 'm^4/s^4         '
463          diagCode  = 'SM      MR      '
464          CALL DIAGNOSTICS_ADD2LIST( diagNum,
465         I          diagName, diagCode, diagUnits, diagTitle, myThid )
466    
467          diagName  = 'PHIBOT  '
468    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         '
472          diagCode  = 'SM      M1      '
473          CALL DIAGNOSTICS_ADD2LIST( diagNum,
474         I          diagName, diagCode, diagUnits, diagTitle, myThid )
475    
476          diagName  = 'PHIBOTSQ'
477    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         '
481          diagCode  = 'SM      M1      '
482          CALL DIAGNOSTICS_ADD2LIST( diagNum,
483         I          diagName, diagCode, diagUnits, diagTitle, myThid )
484    
485          diagName  = 'MXLDEPTH'
486          diagTitle = 'Mixed-Layer Depth (>0)'
487          diagUnits = 'm               '
488        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
489        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADD2LIST( diagNum,
490         I          diagName, diagCode, diagUnits, diagTitle, myThid )
491    
492          diagName  = 'DRHODR  '
493          diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
494          diagUnits = 'kg/m^4          '
495          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
496          diagCode  = 'SM      LR      '
497          CALL DIAGNOSTICS_ADD2LIST( diagNum,
498         I          diagName, diagCode, diagUnits, diagTitle, myThid )
499    
500          diagName  = 'CONVADJ '
501          diagTitle = 'Convective Adjustment Index [0-1] '
502          diagUnits = 'fraction        '
503          diagCode  = 'SM      LR      '
504          CALL DIAGNOSTICS_ADD2LIST( diagNum,
505         I          diagName, diagCode, diagUnits, diagTitle, myThid )
506    
507    C--   surface fluxes:
508          diagName  = 'oceTAUX '
509          diagTitle = 'zonal surface wind stress, >0 increases uVel'
510          diagUnits = 'N/m^2           '
511          diagCode  = 'SU      U1      '
512          CALL DIAGNOSTICS_ADD2LIST( diagNum,
513         I          diagName, diagCode, diagUnits, diagTitle, myThid )
514    
515          diagName  = 'oceTAUY '
516          diagTitle = 'meridional surf. wind stress, >0 increases vVel'
517          diagUnits = 'N/m^2           '
518          diagCode  = 'SV      U1      '
519          CALL DIAGNOSTICS_ADD2LIST( diagNum,
520         I          diagName, diagCode, diagUnits, diagTitle, myThid )
521    
522          diagName  = 'atmPload'
523          diagTitle = 'Atmospheric pressure loading'
524          diagUnits = 'Pa              '
525          diagCode  = 'SM      U1      '
526          CALL DIAGNOSTICS_ADD2LIST( diagNum,
527         I          diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529          diagName  = 'sIceLoad'
530          diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
531          diagUnits = 'kg/m^2          '
532          diagCode  = 'SM      U1      '
533          CALL DIAGNOSTICS_ADD2LIST( diagNum,
534         I          diagName, diagCode, diagUnits, diagTitle, myThid )
535    
536          diagName  = 'oceFWflx'
537          diagTitle = 'net surface Fresh-Water flux into the ocean'
538         &          //' (+=down), >0 decreases salinity'
539          diagUnits = 'kg/m^2/s        '
540          diagCode  = 'SM      U1      '
541          CALL DIAGNOSTICS_ADD2LIST( diagNum,
542         I          diagName, diagCode, diagUnits, diagTitle, myThid )
543    
544          diagName  = 'oceSflux'
545          diagTitle = 'net surface Salt flux into the ocean (+=down),'
546         &          //' >0 increases salinity'
547          diagUnits = 'g/m^2/s         '
548          diagCode  = 'SM      U1      '
549          CALL DIAGNOSTICS_ADD2LIST( diagNum,
550         I          diagName, diagCode, diagUnits, diagTitle, myThid )
551    
552          diagName  = 'oceQnet '
553          diagTitle = 'net surface heat flux into the ocean (+=down),'
554         &          //' >0 increases theta'
555          diagUnits = 'W/m^2           '
556          diagCode  = 'SM      U1      '
557          CALL DIAGNOSTICS_ADD2LIST( diagNum,
558         I          diagName, diagCode, diagUnits, diagTitle, myThid )
559    
560          diagName  = 'oceQsw  '
561          diagTitle = 'net Short-Wave radiation (+=down),'
562         &          //' >0 increases theta'
563          diagUnits = 'W/m^2           '
564          diagCode  = 'SM      U1      '
565          CALL DIAGNOSTICS_ADD2LIST( diagNum,
566         I          diagName, diagCode, diagUnits, diagTitle, myThid )
567    
568          diagName  = 'oceFreez'
569          diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
570          diagUnits = 'W/m^2           '
571          diagCode  = 'SM      U1      '
572          CALL DIAGNOSTICS_ADD2LIST( diagNum,
573         I          diagName, diagCode, diagUnits, diagTitle, myThid )
574    
575          diagName  = 'TRELAX  '
576          diagTitle = 'surface temperature relaxation, >0 increases theta'
577          diagUnits = 'W/m^2           '
578          diagCode  = 'SM      U1      '
579          CALL DIAGNOSTICS_ADD2LIST( diagNum,
580         I          diagName, diagCode, diagUnits, diagTitle, myThid )
581    
582          diagName  = 'SRELAX  '
583          diagTitle = 'surface salinity relaxation, >0 increases salt'
584          diagUnits = 'g/m^2/s         '
585          diagCode  = 'SM      U1      '
586          CALL DIAGNOSTICS_ADD2LIST( diagNum,
587         I          diagName, diagCode, diagUnits, diagTitle, myThid )
588    
589          diagName  = 'surForcT'
590          diagTitle = 'model surface forcing for Temperature,'
591         &          //' >0 increases theta'
592          diagUnits = 'W/m^2           '
593          diagCode  = 'SM      U1      '
594          CALL DIAGNOSTICS_ADD2LIST( diagNum,
595         I          diagName, diagCode, diagUnits, diagTitle, myThid )
596    
597          diagName  = 'surForcS'
598          diagTitle = 'model surface forcing for Salinity,'
599         &          //' >0 increases salinity'
600          diagUnits = 'g/m^2/s         '
601          diagCode  = 'SM      U1      '
602          CALL DIAGNOSTICS_ADD2LIST( diagNum,
603         I          diagName, diagCode, diagUnits, diagTitle, myThid )
604    
605          diagName  = 'TFLUX   '
606          diagTitle = 'total heat flux (match heat-content variations),'
607         &          //' >0 increases theta'
608          diagUnits = 'W/m^2           '
609          diagCode  = 'SM      U1      '
610          CALL DIAGNOSTICS_ADD2LIST( diagNum,
611         I          diagName, diagCode, diagUnits, diagTitle, myThid )
612    
613          diagName  = 'SFLUX   '
614          diagTitle = 'total salt flux (match salt-content variations),'
615         &          //' >0 increases salt'
616          diagUnits = 'g/m^2/s         '
617          diagCode  = 'SM      U1      '
618          CALL DIAGNOSTICS_ADD2LIST( diagNum,
619       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
620    
621  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
622    
623          diagName  = 'RCENTER '
624    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
625          diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
626          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
627          diagCode  = 'SM      MR      '
628          CALL DIAGNOSTICS_ADD2LIST( diagNum,
629         I          diagName, diagCode, diagUnits, diagTitle, myThid )
630    
631          diagName  = 'RSURF   '
632    c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
633          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
634          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
635          diagCode  = 'SM      M1      '
636          CALL DIAGNOSTICS_ADD2LIST( diagNum,
637         I          diagName, diagCode, diagUnits, diagTitle, myThid )
638    
639          diagName  = 'TOTUTEND'
640          diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
641          diagUnits = 'm/s^2           '
642          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
643          CALL DIAGNOSTICS_ADD2LIST( diagNum,
644         I          diagName, diagCode, diagUnits, diagTitle, myThid )
645    
646          diagName  = 'TOTVTEND'
647          diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
648          diagUnits = 'm/s^2           '
649          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
650          CALL DIAGNOSTICS_ADD2LIST( diagNum,
651         I          diagName, diagCode, diagUnits, diagTitle, myThid )
652    
653          diagName  = 'TOTTTEND'
654          diagTitle = 'Tendency of Potential Temperature'
655          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
656          diagCode  = 'SM      MR      '
657          CALL DIAGNOSTICS_ADD2LIST( diagNum,
658         I          diagName, diagCode, diagUnits, diagTitle, myThid )
659    
660          diagName  = 'TOTSTEND'
661          diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
662          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
663          diagCode  = 'SM      MR      '
664          CALL DIAGNOSTICS_ADD2LIST( diagNum,
665         I          diagName, diagCode, diagUnits, diagTitle, myThid )
666    
667    
668        RETURN        RETURN
669        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22