/[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.6 by jmc, Sat Dec 18 23:53:04 2004 UTC revision 1.16 by molod, Wed Jun 15 12:59: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          ELSE
89            tUnit4c= 'degC'
90            sUnit4c= 'psu '
91          ENDIF
92    
93  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
94  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
95    
96        diagName  = 'ETAN    '        diagName  = 'ETAN    '
97        diagTitle = 'Perturbation of Surface (pressure, height) (hPa,m)'        diagTitle = 'Perturbation of Surface (pressure, height) (Pa,m)'
98        diagUnits = 'm or hPa        '        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
99        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
100        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
101       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
102    
103        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
104        diagTitle = 'Square of Perturbation of Sfc (hPa^2,m^2)'        diagTitle = 'Square of Perturbation of Sfc (Pa^2,m^2)'
105        diagUnits = 'm2 or hPa2      '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
106        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
107        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
108       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
109    
110        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
111        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
112        diagUnits = 'm2 or Pa2 /sec2 '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
113        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
114        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
115       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
116    
117        diagName  = 'THETA   '        diagName  = 'THETA   '
118        diagTitle = 'Potential Temperature (K)'        diagTitle = 'Potential Temperature (degC,K)'
119        diagUnits = 'K               '        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
120        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
121        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
122       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
123    
124          diagName  = 'SST     '
125          diagTitle = 'Sea Surface Temperature (degC,K)'
126          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
127          diagCode  = 'SM      M1      '
128          CALL DIAGNOSTICS_ADD2LIST( diagNum,
129         I          diagName, diagCode, diagUnits, diagTitle, myThid )
130    
131        diagName  = 'SALT    '        diagName  = 'SALT    '
132        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
133        diagUnits = 'g/kg            '        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
134          diagCode  = 'SM      MR      '
135          CALL DIAGNOSTICS_ADD2LIST( diagNum,
136         I          diagName, diagCode, diagUnits, diagTitle, myThid )
137    
138          diagName  = 'SSS     '
139          diagTitle = 'Sea Surface Salinity (g/kg)'
140          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
141          diagCode  = 'SM      M1      '
142          CALL DIAGNOSTICS_ADD2LIST( diagNum,
143         I          diagName, diagCode, diagUnits, diagTitle, myThid )
144    
145          diagName  = 'SALTanom'
146          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
147          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
148        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
149        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
150       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
151    
152        diagName  = 'UVEL    '        diagName  = 'UVEL    '
153        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
154        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
155        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
156        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
157       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
158    
159        diagName  = 'VVEL    '        diagName  = 'VVEL    '
160        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
161        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
162        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
163        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
164       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
165    
166          diagName  = 'UVEL_k2 '
167          diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'
168          diagUnits = 'm/s             '
169          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
170          CALL DIAGNOSTICS_ADD2LIST( diagNum,
171         I          diagName, diagCode, diagUnits, diagTitle, myThid )
172    
173          diagName  = 'VVEL_k2 '
174          diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'
175          diagUnits = 'm/s             '
176          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '
177          CALL DIAGNOSTICS_ADD2LIST( diagNum,
178         I          diagName, diagCode, diagUnits, diagTitle, myThid )
179    
180        diagName  = 'WVEL    '        diagName  = 'WVEL    '
181        diagTitle = 'Vertical Component of Velocity (m/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
182        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
183        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
184        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
185       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
186    
187        diagName  = 'THETASQ '        diagName  = 'THETASQ '
188        diagTitle = 'Square of Potential Temperature (K^2)'        diagTitle = 'Square of Potential Temperature (K^2)'
189        diagUnits = 'K2              '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
190        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
191        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
192       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
193    
194        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
195        diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'        diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
196        diagUnits = 'g2/kg2          '        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
197          diagCode  = 'SM      MR      '
198          CALL DIAGNOSTICS_ADD2LIST( diagNum,
199         I          diagName, diagCode, diagUnits, diagTitle, myThid )
200    
201          diagName  = 'SALTSQan'
202          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
203          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
204        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
205        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
206       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
207    
208        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
209        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
210        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
211        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
212        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
213       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
214    
215        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
216        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
217        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
218        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
219        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
220       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
221    
222        diagName  = 'WVELSQ  '        diagName  = 'WVELSQ  '
223        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
224        diagUnits = 'm2/sec2         '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
225        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
226        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
227       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
228    
229        diagName  = 'UVELVVEL'        diagName  = 'UV_VEL_C'
230          diagTitle ='Product of horizontal Comp of velocity (cell center)'
231          diagUnits = 'm^2/s^2         '
232          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
233          CALL DIAGNOSTICS_ADD2LIST( diagNum,
234         I          diagName, diagCode, diagUnits, diagTitle, myThid )
235    
236          diagName  = 'UV_VEL_Z'
237        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
238        diagUnits = 'm2/sec2         '        diagUnits = 'm^2/s^2         '
239        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
240        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
241       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
242    
243          diagName  = 'WU_VEL  '
244          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
245          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
246          diagCode  = 'WU      LR      '
247          CALL DIAGNOSTICS_ADD2LIST( diagNum,
248         I          diagName, diagCode, diagUnits, diagTitle, myThid )
249    
250          diagName  = 'WV_VEL  '
251          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
252          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
253          diagCode  = 'WV      LR      '
254          CALL DIAGNOSTICS_ADD2LIST( diagNum,
255         I          diagName, diagCode, diagUnits, diagTitle, myThid )
256    
257        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
258        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
259        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
260        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
261        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
262       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
263    
264        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
265        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
266        diagUnits = 'm/sec           '        diagUnits = 'm/s             '
267        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
268        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
269       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
270    
271        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
272        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
273        diagUnits = 'm/sec           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
274        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
275        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
276       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
277    
278        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
279        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'
280        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
281        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
282        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
283       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
284    
285        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
286        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'
287        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
288        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
289        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
290       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
291    
292        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
293        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
294        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
295        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
296        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
297       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
298    
299        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
300        diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'
301        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
302        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
303        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
304       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
305    
306        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
307        diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'
308        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
309        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
310        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
311       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
312    
313        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
314        diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'
315        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
316        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
317        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
318       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
319    
320        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
321        diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'        diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
322        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
323        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
324        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
325       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
326    
327        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
328        diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'        diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
329        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
330        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
331        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
332       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
333    
334        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
335        diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
336        diagUnits = 'm-K/sec         '        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
337        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
338        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
339       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
340    
341        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
342        diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'        diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'
343        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
344        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
345        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
346       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
347    
348        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
349        diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'        diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'
350        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
351        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
352        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
353       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
354    
355        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
356        diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'
357        diagUnits = 'm-kg/s-kg       '        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
358        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
359        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
360       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
361    
362  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
363    
364          diagName  = 'RHOAnoma'
365          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
366          diagUnits = 'kg/m^3          '
367          diagCode  = 'SM      MR      '
368          CALL DIAGNOSTICS_ADD2LIST( diagNum,
369         I          diagName, diagCode, diagUnits, diagTitle, myThid )
370    
371          diagName  = 'RHOANOSQ'
372          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
373          diagUnits = 'kg^2/m^6          '
374          diagCode  = 'SM      MR      '
375          CALL DIAGNOSTICS_ADD2LIST( diagNum,
376         I          diagName, diagCode, diagUnits, diagTitle, myThid )
377    
378          diagName  = 'URHOMASS'
379          diagTitle = 'Zonal Transport of Density'
380          diagUnits = 'kg/m^2/s        '
381          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
382          CALL DIAGNOSTICS_ADD2LIST( diagNum,
383         I          diagName, diagCode, diagUnits, diagTitle, myThid )
384    
385          diagName  = 'VRHOMASS'
386          diagTitle = 'Meridional Transport of Density'
387          diagUnits = 'kg/m^2/s        '
388          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
389          CALL DIAGNOSTICS_ADD2LIST( diagNum,
390         I          diagName, diagCode, diagUnits, diagTitle, myThid )
391    
392          diagName  = 'WRHOMASS'
393          diagTitle = 'Vertical Transport of Potential Density'
394          diagUnits = 'kg/m^2/s        '
395          diagCode  = 'WM      LR      '
396          CALL DIAGNOSTICS_ADD2LIST( diagNum,
397         I          diagName, diagCode, diagUnits, diagTitle, myThid )
398    
399        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
400        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
401        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
402          diagCode  = 'SM      MR      '
403          CALL DIAGNOSTICS_ADD2LIST( diagNum,
404         I          diagName, diagCode, diagUnits, diagTitle, myThid )
405    
406          diagName  = 'PHIHYDSQ'
407          diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'
408          diagUnits = 'm^4/s^4         '
409        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
410        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
411       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
412    
413        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
414        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
415        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
416        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
417        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
418       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
419    
420        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
421        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
422        diagUnits = 'm4/s4           '        diagUnits = 'm^4/s^4         '
423        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
424        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
425       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
426    
427        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
428        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
429        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
430          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
431        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
432        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
433       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 309  C---+----1----+----2----+----3----+----4 Line 435  C---+----1----+----2----+----3----+----4
435  C--   Momentum terms:  C--   Momentum terms:
436        diagName  = 'VISCA4  '        diagName  = 'VISCA4  '
437        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
438        diagUnits = 'm4/sec          '        diagUnits = 'm^4/s           '
439        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
440        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
441       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
442    
443        diagName  = 'VISCAH  '        diagName  = 'VISCAH  '
444        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '        diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
445        diagUnits = 'm2/sec          '        diagUnits = 'm^2/s           '
446        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
447        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
448       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 324  C--   Momentum terms: Line 450  C--   Momentum terms:
450  C--   surface fluxes:  C--   surface fluxes:
451        diagName  = 'TAUX    '        diagName  = 'TAUX    '
452        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
453        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
454        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
455        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
456       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
457    
458        diagName  = 'TAUY    '        diagName  = 'TAUY    '
459        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel '
460        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
461        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
462        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
463       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
464    
465        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
466        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'net surface heat flux, >0 increases theta       '
467        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
468        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
469        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
470       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
471    
472        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
473        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
474        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
475        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
476        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
477       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
478    
479        diagName  = 'TICE    '        diagName  = 'TICE    '
480        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
481        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
482        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
483        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
484       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
485    
486        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
487        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt        '
488        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
489        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
490        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
491       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
492    
493        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
494        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt  '
495        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
496        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
497        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
498       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
499    
500  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
501          diagName  = 'PRESSURE'
502          diagTitle = 'Atmospheric Pressure (Pa)     '
503          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
504          diagCode  = 'SM      MR      '
505          CALL DIAGNOSTICS_ADD2LIST( diagNum,
506         I          diagName, diagCode, diagUnits, diagTitle, myThid )
507    
508        RETURN        RETURN
509        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22