/[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.19 by molod, Tue Aug 16 20:36:10 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 (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 )
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          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  = 'SSS     '
144          diagTitle = 'Sea Surface Salinity '
145          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
146          diagCode  = 'SM      M1      '
147          CALL DIAGNOSTICS_ADD2LIST( diagNum,
148         I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150          diagName  = 'SALTanom'
151          diagTitle = 'Salt anomaly (=SALT-35)'
152          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
153          diagCode  = 'SM      MR      '
154          CALL DIAGNOSTICS_ADD2LIST( diagNum,
155         I          diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157          diagName  = 'UVEL    '
158          diagTitle = 'Zonal Component of Velocity (m/s)'
159          diagUnits = 'm/s             '
160          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
161          CALL DIAGNOSTICS_ADD2LIST( diagNum,
162         I          diagName, diagCode, diagUnits, diagTitle, myThid )
163    
164          diagName  = 'VVEL    '
165          diagTitle = 'Meridional Component of Velocity (m/s)'
166          diagUnits = 'm/s             '
167          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
168          CALL DIAGNOSTICS_ADD2LIST( diagNum,
169         I          diagName, diagCode, diagUnits, diagTitle, myThid )
170    
171          diagName  = 'UVEL_k2 '
172          diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'
173          diagUnits = 'm/s             '
174          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
175          CALL DIAGNOSTICS_ADD2LIST( diagNum,
176         I          diagName, diagCode, diagUnits, diagTitle, myThid )
177    
178          diagName  = 'VVEL_k2 '
179          diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'
180          diagUnits = 'm/s             '
181          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '
182          CALL DIAGNOSTICS_ADD2LIST( diagNum,
183         I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185          diagName  = 'WVEL    '
186          diagTitle = 'Vertical Component of Velocity (r_units/s)'
187          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
188          diagCode  = 'WM      LR      '
189          CALL DIAGNOSTICS_ADD2LIST( diagNum,
190         I          diagName, diagCode, diagUnits, diagTitle, myThid )
191    
192          diagName  = 'THETASQ '
193          diagTitle = 'Square of Potential Temperature (K^2)'
194          diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
195          diagCode  = 'SM      MR      '
196          CALL DIAGNOSTICS_ADD2LIST( diagNum,
197         I          diagName, diagCode, diagUnits, diagTitle, myThid )
198    
199          diagName  = 'SALTSQ  '
200          diagTitle = 'Square of Salt or Water Vapor mix rat (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  = 'SALTSQan'
207          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (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  = 'UVELSQ  '
214          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
215          diagUnits = 'm^2/s^2         '
216          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
217          CALL DIAGNOSTICS_ADD2LIST( diagNum,
218         I          diagName, diagCode, diagUnits, diagTitle, myThid )
219    
220          diagName  = 'VVELSQ  '
221          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
222          diagUnits = 'm^2/s^2         '
223          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
224          CALL DIAGNOSTICS_ADD2LIST( diagNum,
225         I          diagName, diagCode, diagUnits, diagTitle, myThid )
226    
227          diagName  = 'WVELSQ  '
228          diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
229          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
230          diagCode  = 'WM      LR      '
231          CALL DIAGNOSTICS_ADD2LIST( diagNum,
232         I          diagName, diagCode, diagUnits, diagTitle, myThid )
233    
234          diagName  = 'UV_VEL_C'
235          diagTitle ='Product of horizontal Comp of velocity (cell center)'
236          diagUnits = 'm^2/s^2         '
237          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
238          CALL DIAGNOSTICS_ADD2LIST( diagNum,
239         I          diagName, diagCode, diagUnits, diagTitle, myThid )
240    
241          diagName  = 'UV_VEL_Z'
242          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
243          diagUnits = 'm^2/s^2         '
244          WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
245          CALL DIAGNOSTICS_ADD2LIST( diagNum,
246         I          diagName, diagCode, diagUnits, diagTitle, myThid )
247    
248          diagName  = 'WU_VEL  '
249          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
250          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
251          diagCode  = 'WU      LR      '
252          CALL DIAGNOSTICS_ADD2LIST( diagNum,
253         I          diagName, diagCode, diagUnits, diagTitle, myThid )
254    
255          diagName  = 'WV_VEL  '
256          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
257          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
258          diagCode  = 'WV      LR      '
259          CALL DIAGNOSTICS_ADD2LIST( diagNum,
260         I          diagName, diagCode, diagUnits, diagTitle, myThid )
261    
262          diagName  = 'UVELMASS'
263          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
264          diagUnits = 'm/s             '
265          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
266          CALL DIAGNOSTICS_ADD2LIST( diagNum,
267         I          diagName, diagCode, diagUnits, diagTitle, myThid )
268    
269          diagName  = 'VVELMASS'
270          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
271          diagUnits = 'm/s             '
272          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
273          CALL DIAGNOSTICS_ADD2LIST( diagNum,
274         I          diagName, diagCode, diagUnits, diagTitle, myThid )
275    
276          diagName  = 'WVELMASS'
277          diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
278          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
279          diagCode  = 'WM      LR      '
280          CALL DIAGNOSTICS_ADD2LIST( diagNum,
281         I          diagName, diagCode, diagUnits, diagTitle, myThid )
282    
283          diagName  = 'UTHMASS '
284          diagTitle = 'Zonal 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)') 'UU   ',diagNum+2,'MR      '
287          CALL DIAGNOSTICS_ADD2LIST( diagNum,
288         I          diagName, diagCode, diagUnits, diagTitle, myThid )
289    
290          diagName  = 'VTHMASS '
291          diagTitle = 'Meridional 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)') 'VV   ', diagNum ,'MR      '
294          CALL DIAGNOSTICS_ADD2LIST( diagNum,
295         I          diagName, diagCode, diagUnits, diagTitle, myThid )
296    
297          diagName  = 'WTHMASS '
298          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
299          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
300          diagCode  = 'WM      LR      '
301          CALL DIAGNOSTICS_ADD2LIST( diagNum,
302         I          diagName, diagCode, diagUnits, diagTitle, myThid )
303    
304          diagName  = 'USLTMASS'
305          diagTitle = 'Zonal Mass-Weight Transp of Salt '
306          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
307          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
308          CALL DIAGNOSTICS_ADD2LIST( diagNum,
309         I          diagName, diagCode, diagUnits, diagTitle, myThid )
310    
311          diagName  = 'VSLTMASS'
312          diagTitle = 'Meridional Mass-Weight Transp of Salt'
313          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
314          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
315          CALL DIAGNOSTICS_ADD2LIST( diagNum,
316         I          diagName, diagCode, diagUnits, diagTitle, myThid )
317    
318          diagName  = 'WSLTMASS'
319          diagTitle = 'Vertical Mass-Weight Transp of Salt '
320          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
321          diagCode  = 'WM      LR      '
322          CALL DIAGNOSTICS_ADD2LIST( diagNum,
323         I          diagName, diagCode, diagUnits, diagTitle, myThid )
324    
325          diagName  = 'UVELTH  '
326          diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
327          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
328          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
329          CALL DIAGNOSTICS_ADD2LIST( diagNum,
330         I          diagName, diagCode, diagUnits, diagTitle, myThid )
331    
332          diagName  = 'VVELTH  '
333          diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
334          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
335          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
336          CALL DIAGNOSTICS_ADD2LIST( diagNum,
337         I          diagName, diagCode, diagUnits, diagTitle, myThid )
338    
339          diagName  = 'WVELTH  '
340          diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
341          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
342          diagCode  = 'WM      LR      '
343          CALL DIAGNOSTICS_ADD2LIST( diagNum,
344         I          diagName, diagCode, diagUnits, diagTitle, myThid )
345    
346          diagName  = 'UVELSLT '
347          diagTitle = 'Zonal Transp of Salt '
348          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
349          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
350          CALL DIAGNOSTICS_ADD2LIST( diagNum,
351         I          diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353          diagName  = 'VVELSLT '
354          diagTitle = 'Meridional Transp of Salt '
355          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
356          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
357          CALL DIAGNOSTICS_ADD2LIST( diagNum,
358         I          diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360          diagName  = 'WVELSLT '
361          diagTitle = 'Vertical Transp of Salt '
362          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
363          diagCode  = 'WM      LR      '
364          CALL DIAGNOSTICS_ADD2LIST( diagNum,
365         I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
368    
369          diagName  = 'RHOAnoma'
370          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
371          diagUnits = 'kg/m^3          '
372          diagCode  = 'SM      MR      '
373          CALL DIAGNOSTICS_ADD2LIST( diagNum,
374         I          diagName, diagCode, diagUnits, diagTitle, myThid )
375    
376          diagName  = 'RHOANOSQ'
377          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
378          diagUnits = 'kg^2/m^6          '
379          diagCode  = 'SM      MR      '
380          CALL DIAGNOSTICS_ADD2LIST( diagNum,
381         I          diagName, diagCode, diagUnits, diagTitle, myThid )
382    
383          diagName  = 'URHOMASS'
384          diagTitle = 'Zonal Transport of Density'
385          diagUnits = 'kg/m^2/s        '
386          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
387          CALL DIAGNOSTICS_ADD2LIST( diagNum,
388         I          diagName, diagCode, diagUnits, diagTitle, myThid )
389    
390          diagName  = 'VRHOMASS'
391          diagTitle = 'Meridional Transport of Density'
392          diagUnits = 'kg/m^2/s        '
393          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
394          CALL DIAGNOSTICS_ADD2LIST( diagNum,
395         I          diagName, diagCode, diagUnits, diagTitle, myThid )
396    
397          diagName  = 'WRHOMASS'
398          diagTitle = 'Vertical Transport of Potential Density'
399          diagUnits = 'kg/m^2/s        '
400          diagCode  = 'WM      LR      '
401          CALL DIAGNOSTICS_ADD2LIST( diagNum,
402         I          diagName, diagCode, diagUnits, diagTitle, myThid )
403    
404        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
405        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
406        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
407          diagCode  = 'SM      MR      '
408          CALL DIAGNOSTICS_ADD2LIST( diagNum,
409         I          diagName, diagCode, diagUnits, diagTitle, myThid )
410    
411          diagName  = 'PHIHYDSQ'
412          diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'
413          diagUnits = 'm^4/s^4         '
414        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
415        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
416       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
417    
418        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
419        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
420        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
421        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
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  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
426        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
427        diagUnits = 'm4/s4           '        diagUnits = 'm^4/s^4         '
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  = 'DRHODR  '        diagName  = 'DRHODR  '
433        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
434        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
435          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
436        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
437        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
438       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
439    
440        diagName  = 'DETADT2 '  C--   Momentum terms:
441        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagName  = 'VISCA4  '
442        diagUnits = 'm2 or Pa2 /sec2 '        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
443        diagCode  = 'SM      M1      '        diagUnits = 'm^4/s           '
444          diagCode  = 'SM      MR      '
445        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
446       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
447    
448          diagName  = 'VISCAH  '
449          diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
450          diagUnits = 'm^2/s           '
451          diagCode  = 'SM      MR      '
452          CALL DIAGNOSTICS_ADD2LIST( diagNum,
453         I          diagName, diagCode, diagUnits, diagTitle, myThid )
454    
455    C--   surface fluxes:
456        diagName  = 'TAUX    '        diagName  = 'TAUX    '
457        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
458        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
459        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
460        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
461       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
462    
463        diagName  = 'TAUY    '        diagName  = 'TAUY    '
464        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel '
465        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
466        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
467        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
468       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
469    
470        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
471        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'net surface heat flux, >0 increases theta       '
472        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
473        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
474        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
475       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
476    
477        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
478        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
479        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
480        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
481        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
482       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
483    
484        diagName  = 'TICE    '        diagName  = 'TICE    '
485        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
486        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
487        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
488        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
489       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
490    
491        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
492        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt        '
493        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
494        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
495        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
496       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
497    
498        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
499        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt  '
500        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
501        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
502        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
503       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
504    
505    #ifdef SHORTWAVE_HEATING
506          diagName  = 'SWFLUX  '
507          diagTitle = 'net upward SW radiation, >0 increases theta     '
508          diagUnits = 'W/m^2           '
509          diagCode  = 'SM      U1      '
510          CALL DIAGNOSTICS_ADD2LIST( diagNum,
511         I          diagName, diagCode, diagUnits, diagTitle, myThid )
512    #endif
513    
514  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
515          diagName  = 'PRESSURE'
516          diagTitle = 'Atmospheric Pressure (Pa)     '
517          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
518          diagCode  = 'SM      MR      '
519          CALL DIAGNOSTICS_ADD2LIST( diagNum,
520         I          diagName, diagCode, diagUnits, diagTitle, myThid )
521    
522          diagName  = 'RSURF   '
523          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
524          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
525          diagCode  = 'SM      M1      '
526          CALL DIAGNOSTICS_ADD2LIST( diagNum,
527         I          diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529    
530        RETURN        RETURN
531        END        END
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22