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