/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Mon Dec 13 21:43:54 2004 UTC revision 1.20 by molod, Tue Aug 16 21:29:38 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 58  C     !INPUT PARAMETERS: Line 60  C     !INPUT PARAMETERS:
60  CEOP  CEOP
61    
62  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
63          INTEGER       diagNum
64        CHARACTER*8   diagName        CHARACTER*8   diagName
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 69  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        diagName  = 'DRHODR  '        IF ( usingPCoords ) THEN
81        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '          rUnit2c= 'Pa'
82        diagUnits = 'kg/m4           '        ELSE
83        diagCode  = 'SM      LR      '          rUnit2c= 'm '
84        CALL DIAGNOSTICS_ADD2LIST(        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 (Pa^2,m^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 )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
114    
115        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
116        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
117        diagUnits = 'm2 or Pa2 /sec2 '        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          diagName  = 'SST     '
130          diagTitle = 'Sea Surface Temperature (degC,K)'
131          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
132          diagCode  = 'SM      M1      '
133          CALL DIAGNOSTICS_ADD2LIST( diagNum,
134         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          diagName  = 'SSS     '
151          diagTitle = 'Sea Surface Salinity '
152          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
153        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
154        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADD2LIST( diagNum,
155         I          diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157          diagName  = 'SALTanom'
158          diagTitle = 'Salt anomaly (=SALT-35)'
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  = 'UVEL_k2 '
179          diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'
180          diagUnits = 'm/s             '
181          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
182          CALL DIAGNOSTICS_ADD2LIST( diagNum,
183         I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185          diagName  = 'VVEL_k2 '
186          diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'
187          diagUnits = 'm/s             '
188          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '
189          CALL DIAGNOSTICS_ADD2LIST( diagNum,
190         I          diagName, diagCode, diagUnits, diagTitle, myThid )
191    
192          diagName  = 'WVEL    '
193          diagTitle = 'Vertical Component of Velocity (r_units/s)'
194          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
195          diagCode  = 'WM      LR      '
196          CALL DIAGNOSTICS_ADD2LIST( diagNum,
197         I          diagName, diagCode, diagUnits, diagTitle, myThid )
198    
199          diagName  = 'THETASQ '
200          diagTitle = 'Square of Potential Temperature (K^2)'
201          diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
202          diagCode  = 'SM      MR      '
203          CALL DIAGNOSTICS_ADD2LIST( diagNum,
204         I          diagName, diagCode, diagUnits, diagTitle, myThid )
205    
206          diagName  = 'SALTSQ  '
207          diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
208          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
209          diagCode  = 'SM      MR      '
210          CALL DIAGNOSTICS_ADD2LIST( diagNum,
211         I          diagName, diagCode, diagUnits, diagTitle, myThid )
212    
213          diagName  = 'SALTSQan'
214          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
215          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
216          diagCode  = 'SM      MR      '
217          CALL DIAGNOSTICS_ADD2LIST( diagNum,
218         I          diagName, diagCode, diagUnits, diagTitle, myThid )
219    
220          diagName  = 'UVELSQ  '
221          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
222          diagUnits = 'm^2/s^2         '
223          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
224          CALL DIAGNOSTICS_ADD2LIST( diagNum,
225         I          diagName, diagCode, diagUnits, diagTitle, myThid )
226    
227          diagName  = 'VVELSQ  '
228          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
229          diagUnits = 'm^2/s^2         '
230          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
231          CALL DIAGNOSTICS_ADD2LIST( diagNum,
232         I          diagName, diagCode, diagUnits, diagTitle, myThid )
233    
234          diagName  = 'WVELSQ  '
235          diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
236          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
237          diagCode  = 'WM      LR      '
238          CALL DIAGNOSTICS_ADD2LIST( diagNum,
239         I          diagName, diagCode, diagUnits, diagTitle, myThid )
240    
241          diagName  = 'UV_VEL_C'
242          diagTitle ='Product of horizontal Comp of velocity (cell center)'
243          diagUnits = 'm^2/s^2         '
244          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
245          CALL DIAGNOSTICS_ADD2LIST( diagNum,
246         I          diagName, diagCode, diagUnits, diagTitle, myThid )
247    
248          diagName  = 'UV_VEL_Z'
249          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
250          diagUnits = 'm^2/s^2         '
251          WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
252          CALL DIAGNOSTICS_ADD2LIST( diagNum,
253         I          diagName, diagCode, diagUnits, diagTitle, myThid )
254    
255          diagName  = 'WU_VEL  '
256          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
257          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
258          diagCode  = 'WU      LR      '
259          CALL DIAGNOSTICS_ADD2LIST( diagNum,
260         I          diagName, diagCode, diagUnits, diagTitle, myThid )
261    
262          diagName  = 'WV_VEL  '
263          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
264          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
265          diagCode  = 'WV      LR      '
266          CALL DIAGNOSTICS_ADD2LIST( diagNum,
267         I          diagName, diagCode, diagUnits, diagTitle, myThid )
268    
269          diagName  = 'UVELMASS'
270          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
271          diagUnits = 'm/s             '
272          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
273          CALL DIAGNOSTICS_ADD2LIST( diagNum,
274         I          diagName, diagCode, diagUnits, diagTitle, myThid )
275    
276          diagName  = 'VVELMASS'
277          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
278          diagUnits = 'm/s             '
279          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
280          CALL DIAGNOSTICS_ADD2LIST( diagNum,
281         I          diagName, diagCode, diagUnits, diagTitle, myThid )
282    
283          diagName  = 'WVELMASS'
284          diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
285          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
286          diagCode  = 'WM      LR      '
287          CALL DIAGNOSTICS_ADD2LIST( diagNum,
288         I          diagName, diagCode, diagUnits, diagTitle, myThid )
289    
290          diagName  = 'UTHMASS '
291          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'
292          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
293          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
294          CALL DIAGNOSTICS_ADD2LIST( diagNum,
295         I          diagName, diagCode, diagUnits, diagTitle, myThid )
296    
297          diagName  = 'VTHMASS '
298          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'
299          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
300          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
301          CALL DIAGNOSTICS_ADD2LIST( diagNum,
302         I          diagName, diagCode, diagUnits, diagTitle, myThid )
303    
304          diagName  = 'WTHMASS '
305          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
306          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
307          diagCode  = 'WM      LR      '
308          CALL DIAGNOSTICS_ADD2LIST( diagNum,
309         I          diagName, diagCode, diagUnits, diagTitle, myThid )
310    
311          diagName  = 'USLTMASS'
312          diagTitle = 'Zonal Mass-Weight Transp of Salt '
313          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
314          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
315          CALL DIAGNOSTICS_ADD2LIST( diagNum,
316         I          diagName, diagCode, diagUnits, diagTitle, myThid )
317    
318          diagName  = 'VSLTMASS'
319          diagTitle = 'Meridional Mass-Weight Transp of Salt'
320          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
321          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
322          CALL DIAGNOSTICS_ADD2LIST( diagNum,
323         I          diagName, diagCode, diagUnits, diagTitle, myThid )
324    
325          diagName  = 'WSLTMASS'
326          diagTitle = 'Vertical Mass-Weight Transp of Salt '
327          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
328          diagCode  = 'WM      LR      '
329          CALL DIAGNOSTICS_ADD2LIST( diagNum,
330         I          diagName, diagCode, diagUnits, diagTitle, myThid )
331    
332          diagName  = 'UVELTH  '
333          diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
334          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
335          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
336          CALL DIAGNOSTICS_ADD2LIST( diagNum,
337         I          diagName, diagCode, diagUnits, diagTitle, myThid )
338    
339          diagName  = 'VVELTH  '
340          diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
341          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
342          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
343          CALL DIAGNOSTICS_ADD2LIST( diagNum,
344         I          diagName, diagCode, diagUnits, diagTitle, myThid )
345    
346          diagName  = 'WVELTH  '
347          diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
348          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
349          diagCode  = 'WM      LR      '
350          CALL DIAGNOSTICS_ADD2LIST( diagNum,
351         I          diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353          diagName  = 'UVELSLT '
354          diagTitle = 'Zonal Transp of Salt '
355          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
356          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
357          CALL DIAGNOSTICS_ADD2LIST( diagNum,
358         I          diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360          diagName  = 'VVELSLT '
361          diagTitle = 'Meridional Transp of Salt '
362          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
363          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
364          CALL DIAGNOSTICS_ADD2LIST( diagNum,
365         I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367          diagName  = 'WVELSLT '
368          diagTitle = 'Vertical Transp of Salt '
369          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
370          diagCode  = 'WM      LR      '
371          CALL DIAGNOSTICS_ADD2LIST( diagNum,
372       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
373    
374  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  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  '
412          diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
413          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      '
422          CALL DIAGNOSTICS_ADD2LIST( diagNum,
423         I          diagName, diagCode, diagUnits, diagTitle, myThid )
424    
425          diagName  = 'PHIBOT  '
426          diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
427          diagUnits = 'm^2/s^2         '
428          diagCode  = 'SM      M1      '
429          CALL DIAGNOSTICS_ADD2LIST( diagNum,
430         I          diagName, diagCode, diagUnits, diagTitle, myThid )
431    
432          diagName  = 'PHIBOTSQ'
433          diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
434          diagUnits = 'm^4/s^4         '
435          diagCode  = 'SM      M1      '
436          CALL DIAGNOSTICS_ADD2LIST( diagNum,
437         I          diagName, diagCode, diagUnits, diagTitle, myThid )
438    
439          diagName  = 'DRHODR  '
440          diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
441          diagUnits = 'kg/m^4          '
442          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
443          diagCode  = 'SM      LR      '
444          CALL DIAGNOSTICS_ADD2LIST( diagNum,
445         I          diagName, diagCode, diagUnits, diagTitle, myThid )
446    
447    C--   Momentum terms:
448          diagName  = 'VISCA4  '
449          diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
450          diagUnits = 'm^4/s           '
451          diagCode  = 'SM      MR      '
452          CALL DIAGNOSTICS_ADD2LIST( diagNum,
453         I          diagName, diagCode, diagUnits, diagTitle, myThid )
454    
455          diagName  = 'VISCAH  '
456          diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
457          diagUnits = 'm^2/s           '
458          diagCode  = 'SM      MR      '
459          CALL DIAGNOSTICS_ADD2LIST( diagNum,
460         I          diagName, diagCode, diagUnits, diagTitle, myThid )
461    
462    C--   surface fluxes:
463          diagName  = 'TAUX    '
464          diagTitle = 'zonal surface wind stress, >0 increases uVel    '
465          diagUnits = 'N/m^2           '
466          diagCode  = 'SU      U1      '
467          CALL DIAGNOSTICS_ADD2LIST( diagNum,
468         I          diagName, diagCode, diagUnits, diagTitle, myThid )
469    
470          diagName  = 'TAUY    '
471          diagTitle = 'meridional surf. wind stress, >0 increases vVel '
472          diagUnits = 'N/m^2           '
473          diagCode  = 'SV      U1      '
474          CALL DIAGNOSTICS_ADD2LIST( diagNum,
475         I          diagName, diagCode, diagUnits, diagTitle, myThid )
476    
477          diagName  = 'TFLUX   '
478          diagTitle = 'net surface heat flux, >0 increases theta       '
479          diagUnits = 'W/m^2           '
480          diagCode  = 'SM      U1      '
481          CALL DIAGNOSTICS_ADD2LIST( diagNum,
482         I          diagName, diagCode, diagUnits, diagTitle, myThid )
483    
484          diagName  = 'TRELAX  '
485          diagTitle = 'surface temperature relaxation, >0 increases theta'
486          diagUnits = 'W/m^2           '
487          diagCode  = 'SM      U1      '
488          CALL DIAGNOSTICS_ADD2LIST( diagNum,
489         I          diagName, diagCode, diagUnits, diagTitle, myThid )
490    
491          diagName  = 'TICE    '
492          diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
493          diagUnits = 'W/m^2           '
494          diagCode  = 'SM      U1      '
495          CALL DIAGNOSTICS_ADD2LIST( diagNum,
496         I          diagName, diagCode, diagUnits, diagTitle, myThid )
497    
498          diagName  = 'SFLUX   '
499          diagTitle = 'net surface salt flux, >0 increases salt        '
500          diagUnits = 'g/m^2/s         '
501          diagCode  = 'SM      U1      '
502          CALL DIAGNOSTICS_ADD2LIST( diagNum,
503         I          diagName, diagCode, diagUnits, diagTitle, myThid )
504    
505          diagName  = 'SRELAX  '
506          diagTitle = 'surface salinity relaxation, >0 increases salt  '
507          diagUnits = 'g/m^2/s         '
508          diagCode  = 'SM      U1      '
509          CALL DIAGNOSTICS_ADD2LIST( diagNum,
510         I          diagName, diagCode, diagUnits, diagTitle, myThid )
511    
512    #ifdef SHORTWAVE_HEATING
513          diagName  = 'SWFLUX  '
514          diagTitle = 'net upward SW radiation, >0 increases theta     '
515          diagUnits = 'W/m^2           '
516          diagCode  = 'SM      U1      '
517          CALL DIAGNOSTICS_ADD2LIST( diagNum,
518         I          diagName, diagCode, diagUnits, diagTitle, myThid )
519    #endif
520    
521    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
522          diagName  = 'PRESSURE'
523          diagTitle = 'Atmospheric Pressure (Pa)     '
524          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
525          diagCode  = 'SM      MR      '
526          CALL DIAGNOSTICS_ADD2LIST( diagNum,
527         I          diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529          diagName  = 'RSURF   '
530          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
531          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
532          diagCode  = 'SM      M1      '
533          CALL DIAGNOSTICS_ADD2LIST( diagNum,
534         I          diagName, diagCode, diagUnits, diagTitle, myThid )
535    
536    
537        RETURN        RETURN
538        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22