/[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.2 by jmc, Wed Dec 15 00:20:37 2004 UTC revision 1.38 by jmc, Tue Jun 14 00:18:37 2011 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, sUnit4c
80          CHARACTER*(10) rTitle, eTitle, fTitle
81          CHARACTER*(20) pTitle, sTitle
82    
83          CHARACTER*(16) DIAGS_MK_UNITS
84          EXTERNAL DIAGS_MK_UNITS
85          CHARACTER*(80) DIAGS_MK_TITLE
86          EXTERNAL DIAGS_MK_TITLE
87    
88  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
89  C     For each output variable,  C     For each output variable,
90  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)  C     specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
91  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
92  C----------------------------------------------------------------------  C----------------------------------------------------------------------
93    
94          IF ( usingPCoords ) THEN
95            rUnit2c= 'Pa'
96            rTitle = ' Pressure '
97            pTitle = ' Geopotential       '
98          ELSE
99            rUnit2c= 'm '
100            rTitle = ' Height   '
101            pTitle = 'Pressure Pot.(p/rho)'
102          ENDIF
103          IF ( fluidIsAir ) THEN
104            tUnit4c= 'K   '
105            sUnit4c= 'g/kg'
106            sTitle = ' Specific Humidity  '
107    #ifdef ALLOW_FIZHI
108            IF (useFIZHI) sUnit4c= 'kg/kg'
109    #endif /* ALLOW_FIZHI */
110          ELSE
111            tUnit4c= 'degC'
112            sUnit4c= 'psu '
113            sTitle = ' Salinity           '
114          ENDIF
115    C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
116          IF ( fluidIsAir ) THEN
117            eTitle = ' Surface  '
118            fTitle = ' Top      '
119          ELSEIF ( usingPCoords ) THEN
120            eTitle = ' Bottom   '
121            fTitle = ' Surface  '
122          ELSE
123            eTitle = ' Surface  '
124            fTitle = ' Bottom   '
125          ENDIF
126    
127    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
128    C-    state variables of the main code (and related quadratic var):
129    
130    
131          diagName  = 'ETAN    '
132          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
133    c     IF ( fluidIsWater .AND. usingZCoords )
134    c    &diagTitle = 'Sea Surface Elevation'
135          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
136          diagCode  = 'SM      M1      '
137          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
138         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
139    
140          diagName  = 'ETANSQ  '
141          diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
142         I                          //' Anomaly', myThid )
143          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
144          diagCode  = 'SM P    M1      '
145          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
146         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
147    
148          diagName  = 'DETADT2 '
149          diagTitle = DIAGS_MK_TITLE( 'Square of '//eTitle//rTitle
150         I                          //' Anomaly Tendency', myThid )
151          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
152          diagCode  = 'SM      M1      '
153          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
154         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
155    
156          diagName  = 'THETA   '
157          diagTitle = 'Potential Temperature'
158          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
159          diagCode  = 'SMR     MR      '
160          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
161         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
162    
163    c     diagName  = 'SST     '
164    c     diagTitle = 'Sea Surface Temperature (degC,K)'
165    c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
166    c     diagCode  = 'SM      M1      '
167    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
168    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
169    
170          diagName  = 'SALT    '
171          diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
172          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
173          diagCode  = 'SMR     MR      '
174          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
175         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
176    
177          diagName  = 'RELHUM  '
178          diagTitle = 'Relative Humidity'
179          diagUnits = 'percent         '
180          diagCode  = 'SMR     MR      '
181          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
182         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
183    
184    c     diagName  = 'SSS     '
185    c     diagTitle = 'Sea Surface Salinity '
186    c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
187    c     diagCode  = 'SM      M1      '
188    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
189    c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
190    
191          IF ( fluidIsWater ) THEN
192          diagName  = 'SALTanom'
193          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
194          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
195          diagCode  = 'SMR     MR      '
196          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
197         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
198          ENDIF
199    
200          diagName  = 'UVEL    '
201          diagTitle = 'Zonal Component of Velocity (m/s)'
202          diagUnits = 'm/s             '
203          diagCode  = 'UUR     MR      '
204          diagMate  = diagNum + 2
205          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
206         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
207    
208          diagName  = 'VVEL    '
209          diagTitle = 'Meridional Component of Velocity (m/s)'
210          diagUnits = 'm/s             '
211          diagCode  = 'VVR     MR      '
212          diagMate  = diagNum
213          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
214         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
215    
216          diagName  = 'WVEL    '
217          diagTitle = 'Vertical Component of Velocity (r_units/s)'
218          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
219          diagCode  = 'WM      LR      '
220          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
221         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
222    
223          diagName  = 'THETASQ '
224          diagTitle = 'Square of Potential Temperature'
225          diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
226          diagCode  = 'SMRP    MR      '
227          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
229    
230          diagName  = 'SALTSQ  '
231          diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
232          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233          diagCode  = 'SMRP    MR      '
234          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
235         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
236    
237          IF ( fluidIsWater ) THEN
238          diagName  = 'SALTSQan'
239          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
240          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
241          diagCode  = 'SMRP    MR      '
242          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
243         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
244          ENDIF
245    
246          diagName  = 'UVELSQ  '
247          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
248          diagUnits = 'm^2/s^2         '
249          diagCode  = 'UURP    MR      '
250          diagMate  = diagNum + 2
251          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
252         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
253    
254          diagName  = 'VVELSQ  '
255          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
256          diagUnits = 'm^2/s^2         '
257          diagCode  = 'VVRP    MR      '
258          diagMate  = diagNum
259          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
260         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
261    
262          diagName  = 'WVELSQ  '
263          diagTitle = 'Square of Vertical Comp of Velocity'
264          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
265          diagCode  = 'WM P    LR      '
266          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
267         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
268    
269          diagName  = 'UE_VEL_C'
270          diagTitle = 'Eastward Velocity (m/s) (cell center)'
271          diagUnits = 'm/s             '
272          diagCode  = 'UMR     MR      '
273          diagMate  = diagNum + 2
274          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
275         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
276    
277          diagName  = 'VN_VEL_C'
278          diagTitle = 'Northward Velocity (m/s) (cell center)'
279          diagUnits = 'm/s             '
280          diagCode  = 'VMR     MR      '
281          diagMate  = diagNum
282          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
283         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
284    
285          diagName  = 'UV_VEL_C'
286          diagTitle ='Product of horizontal Comp of velocity (cell center)'
287          diagUnits = 'm^2/s^2         '
288          diagCode  = 'UMR     MR      '
289          diagMate  = diagNum + 1
290          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
291         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
292    
293          diagName  = 'UV_VEL_Z'
294          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
295          diagUnits = 'm^2/s^2         '
296          diagCode  = 'UZR     MR      '
297          diagMate  = diagNum + 1
298          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
299         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
300    
301          diagName  = 'WU_VEL  '
302          diagTitle = 'Vertical Transport of Zonal Momentum'
303          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
304          diagCode  = 'WU      LR      '
305          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
306         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
307    
308          diagName  = 'WV_VEL  '
309          diagTitle ='Vertical Transport of Meridional Momentum'
310          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
311          diagCode  = 'WV      LR      '
312          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
313         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
314    
315          diagName  = 'UVELMASS'
316          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
317          diagUnits = 'm/s             '
318          diagCode  = 'UUr     MR      '
319          diagMate  = diagNum + 2
320          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
321         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
322    
323          diagName  = 'VVELMASS'
324          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
325          diagUnits = 'm/s             '
326          diagCode  = 'VVr     MR      '
327          diagMate  = diagNum
328          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
329         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
330    
331          diagName  = 'WVELMASS'
332          diagTitle = 'Vertical Mass-Weighted Comp of Velocity'
333          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
334          diagCode  = 'WM      LR      '
335          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
336         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
337    
338          diagName  = 'PhiVEL  '
339          diagTitle = 'Horizontal Velocity Potential (m^2/s)'
340          diagUnits = 'm^2/s           '
341          diagCode  = 'SMR P   MR      '
342    C-    use 'UVELMASS' as mate.
343          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
344         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
345    
346    c     diagName  = 'PsiVEL  '
347    c     diagTitle = 'Horizontal Velocity Stream-Function'
348    c     diagUnits = 'm^2/s           '
349    c     diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
350    c     diagCode  = 'SZ  P   MR      '
351    C-    use 'PhiVEL' as mate.
352    c     diagMate  = diagNum
353    c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
354    c    I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
355    
356          diagName  = 'UTHMASS '
357          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
358          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
359          diagCode  = 'UUr     MR      '
360          diagMate  = diagNum + 2
361          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
362         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
363    
364          diagName  = 'VTHMASS '
365          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
366          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
367          diagCode  = 'VVr     MR      '
368          diagMate  = diagNum
369          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
370         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
371    
372          diagName  = 'WTHMASS '
373          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
374          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
375          diagCode  = 'WM      LR      '
376          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
377         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
378    
379          diagName  = 'USLTMASS'
380          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
381         I                           //sTitle, myThid )
382          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
383          diagCode  = 'UUr     MR      '
384          diagMate  = diagNum + 2
385          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
386         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
387    
388          diagName  = 'VSLTMASS'
389          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
390         I                           //sTitle, myThid )
391          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
392          diagCode  = 'VVr     MR      '
393          diagMate  = diagNum
394          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
395         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
396    
397          diagName  = 'WSLTMASS'
398          diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
399         I                           //sTitle, myThid )
400          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
401          diagCode  = 'WM      LR      '
402          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
403         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
404    
405          diagName  = 'UVELTH  '
406          diagTitle = 'Zonal Transport of Pot Temp'
407          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
408          diagCode  = 'UUR     MR      '
409          diagMate  = diagNum + 2
410          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
411         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
412    
413          diagName  = 'VVELTH  '
414          diagTitle = 'Meridional Transport of Pot Temp'
415          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
416          diagCode  = 'VVR     MR      '
417          diagMate  = diagNum
418          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
419         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
420    
421          diagName  = 'WVELTH  '
422          diagTitle = 'Vertical Transport of Pot Temp'
423          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
424          diagCode  = 'WM      LR      '
425          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
426         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
427    
428          diagName  = 'UVELSLT '
429          diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
430         I                          //sTitle, myThid )
431          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
432          diagCode  = 'UUR     MR      '
433          diagMate  = diagNum + 2
434          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
435         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
436    
437          diagName  = 'VVELSLT '
438          diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
439         I                          //sTitle, myThid )
440          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
441          diagCode  = 'VVR     MR      '
442          diagMate  = diagNum
443          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
444         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
445    
446          diagName  = 'WVELSLT '
447          diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
448         I                          //sTitle, myThid )
449          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
450          diagCode  = 'WM      LR      '
451          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
452         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
453    
454          diagName  = 'UVELPHI '
455          diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
456         I                 //pTitle//' Anomaly', myThid )
457          diagUnits = 'm^3/s^3         '
458          diagCode  = 'UUr     MR      '
459          diagMate  = diagNum + 2
460          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
461         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
462    
463          diagName  = 'VVELPHI '
464          diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
465         I                 //pTitle//' Anomaly', myThid )
466          diagUnits = 'm^3/s^3         '
467          diagCode  = 'VVr     MR      '
468          diagMate  = diagNum
469          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
470         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
471    
472    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
473    
474          diagName  = 'RHOAnoma'
475          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
476          diagUnits = 'kg/m^3          '
477          diagCode  = 'SMR     MR      '
478          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
479         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
480    
481          diagName  = 'RHOANOSQ'
482          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
483          diagUnits = 'kg^2/m^6          '
484          diagCode  = 'SMRP    MR      '
485          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
486         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
487    
488          diagName  = 'URHOMASS'
489          diagTitle = 'Zonal Transport of Density'
490          diagUnits = 'kg/m^2/s        '
491          diagCode  = 'UUr     MR      '
492          diagMate  = diagNum + 2
493          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
494         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
495    
496          diagName  = 'VRHOMASS'
497          diagTitle = 'Meridional Transport of Density'
498          diagUnits = 'kg/m^2/s        '
499          diagCode  = 'VVr     MR      '
500          diagMate  = diagNum
501          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
502         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
503    
504          diagName  = 'WRHOMASS'
505          diagTitle = 'Vertical Transport of Potential Density'
506          diagUnits = 'kg/m^2/s        '
507          diagCode  = 'WM      LR      '
508          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
509         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
510    
511        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
512        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
513        diagUnits = 'm2/s2           '       I                           //pTitle//' Anomaly', myThid )
514        diagCode  = 'SM      MR      '        diagUnits = 'm^2/s^2         '
515        CALL DIAGNOSTICS_ADD2LIST(        diagCode  = 'SMR     MR      '
516       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
517         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
518    
519          diagName  = 'PHIHYDSQ'
520          diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
521         I                           //pTitle//' Anomaly', myThid )
522          diagUnits = 'm^4/s^4         '
523          diagCode  = 'SMRP    MR      '
524          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
525         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
526    
527        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
528        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'  c     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
529        diagUnits = 'm2/s2           '        diagTitle = DIAGS_MK_TITLE( fTitle
530        diagCode  = 'SM      MR      '       I                           //pTitle//' Anomaly', myThid )
531        CALL DIAGNOSTICS_ADD2LIST(        diagUnits = 'm^2/s^2         '
532       I          diagName, diagCode, diagUnits, diagTitle, myThid )        diagCode  = 'SM      M1      '
533          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
534         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
535    
536        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
537        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'  c     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
538        diagUnits = 'm4/s4           '        diagTitle = DIAGS_MK_TITLE( 'Square of '//fTitle
539         I                           //pTitle//' Anomaly', myThid )
540          diagUnits = 'm^4/s^4         '
541          diagCode  = 'SM P    M1      '
542          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
543         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
544    
545    #ifdef ALLOW_NONHYDROSTATIC
546          diagName  = 'PHI_NH  '
547          diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
548          diagUnits = 'm^2/s^2         '
549          diagCode  = 'SMR     MR      '
550          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
551         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
552    #endif /* ALLOW_NONHYDROSTATIC */
553    
554          diagName  = 'MXLDEPTH'
555          diagTitle = 'Mixed-Layer Depth (>0)'
556          diagUnits = 'm               '
557        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
558        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
559       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
560    
561        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
562        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
563        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
564          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
565        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
566        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
567       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
568    
569        diagName  = 'DETADT2 '        diagName  = 'CONVADJ '
570        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = 'Convective Adjustment Index [0-1] '
571        diagUnits = 'm2 or Pa2 /sec2 '        diagUnits = 'fraction        '
572        diagCode  = 'SM      M1      '        diagCode  = 'SMR     LR      '
573        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
574       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
575    
576    C--   surface fluxes:
577          diagName  = 'oceTAUX '
578          diagTitle = 'zonal surface wind stress, >0 increases uVel'
579          diagUnits = 'N/m^2           '
580          diagCode  = 'UU      U1      '
581          diagMate  = diagNum + 2
582          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
583         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
584    
585          diagName  = 'oceTAUY '
586          diagTitle = 'meridional surf. wind stress, >0 increases vVel'
587          diagUnits = 'N/m^2           '
588          diagCode  = 'VV      U1      '
589          diagMate  = diagNum
590          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
591         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
592    
593          diagName  = 'atmPload'
594          diagTitle = 'Atmospheric pressure loading'
595          diagUnits = 'Pa              '
596          diagCode  = 'SM      U1      '
597          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
598         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
599    
600          diagName  = 'sIceLoad'
601          diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
602          diagUnits = 'kg/m^2          '
603          diagCode  = 'SM      U1      '
604          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
605         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
606    
607          diagName  = 'oceFWflx'
608          diagTitle = 'net surface Fresh-Water flux into the ocean'
609         &          //' (+=down), >0 decreases salinity'
610          diagUnits = 'kg/m^2/s        '
611          diagCode  = 'SM      U1      '
612          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
613         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
614    
615          diagName  = 'oceSflux'
616          diagTitle = 'net surface Salt flux into the ocean (+=down),'
617         &          //' >0 increases salinity'
618          diagUnits = 'g/m^2/s         '
619          diagCode  = 'SM      U1      '
620          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
621         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
622    
623          diagName  = 'oceQnet '
624          diagTitle = 'net surface heat flux into the ocean (+=down),'
625         &          //' >0 increases theta'
626          diagUnits = 'W/m^2           '
627          diagCode  = 'SM      U1      '
628          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
629         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
630    
631          diagName  = 'oceQsw  '
632          diagTitle = 'net Short-Wave radiation (+=down),'
633         &          //' >0 increases theta'
634          diagUnits = 'W/m^2           '
635          diagCode  = 'SM      U1      '
636          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
637         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
638    
639          diagName  = 'oceFreez'
640          diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
641          diagUnits = 'W/m^2           '
642          diagCode  = 'SM      U1      '
643          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
644         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
645    
646          diagName  = 'TRELAX  '
647          diagTitle = 'surface temperature relaxation, >0 increases theta'
648          diagUnits = 'W/m^2           '
649          diagCode  = 'SM      U1      '
650          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
651         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
652    
653          diagName  = 'SRELAX  '
654          diagTitle = 'surface salinity relaxation, >0 increases salt'
655          diagUnits = 'g/m^2/s         '
656          diagCode  = 'SM      U1      '
657          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
658         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
659    
660          diagName  = 'surForcT'
661          diagTitle = 'model surface forcing for Temperature,'
662         &          //' >0 increases theta'
663          diagUnits = 'W/m^2           '
664          diagCode  = 'SM      U1      '
665          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
666         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
667    
668          diagName  = 'surForcS'
669          diagTitle = 'model surface forcing for Salinity,'
670         &          //' >0 increases salinity'
671          diagUnits = 'g/m^2/s         '
672          diagCode  = 'SM      U1      '
673          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
674         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
675    
676          diagName  = 'TFLUX   '
677          diagTitle = 'total heat flux (match heat-content variations),'
678         &          //' >0 increases theta'
679          diagUnits = 'W/m^2           '
680          diagCode  = 'SM      U1      '
681          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
682         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
683    
684          diagName  = 'SFLUX   '
685          diagTitle = 'total salt flux (match salt-content variations),'
686         &          //' >0 increases salt'
687          diagUnits = 'g/m^2/s         '
688          diagCode  = 'SM      U1      '
689          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
690         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
691    
692  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
693    
694          diagName  = 'RCENTER '
695    c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
696          diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
697          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
698          diagCode  = 'SM      MR      '
699          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
700         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
701    
702          diagName  = 'RSURF   '
703    c     diagTitle = 'Free-Surface r-Position (Pressure, Height) (Pa,m)'
704          diagTitle = DIAGS_MK_TITLE( eTitle//rTitle, myThid )
705          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
706          diagCode  = 'SM      M1      '
707          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
708         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
709    
710          diagName  = 'TOTUTEND'
711          diagTitle = 'Tendency of Zonal Component of Velocity'
712          diagUnits = 'm/s/day         '
713          diagCode  = 'UUR     MR      '
714          diagMate  = diagNum + 2
715          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
716         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
717    
718          diagName  = 'TOTVTEND'
719          diagTitle = 'Tendency of Meridional Component of Velocity'
720          diagUnits = 'm/s/day         '
721          diagCode  = 'VVR     MR      '
722          diagMate  = diagNum
723          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
724         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
725    
726          diagName  = 'TOTTTEND'
727          diagTitle = 'Tendency of Potential Temperature'
728          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
729          diagCode  = 'SMR     MR      '
730          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
731         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
732    
733          diagName  = 'TOTSTEND'
734          diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
735          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )
736          diagCode  = 'SMR     MR      '
737          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
738         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
739    
740    
741        RETURN        RETURN
742        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22