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