/[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.7 by jmc, Mon Jan 3 02:32:33 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 58  C     !INPUT PARAMETERS: Line 59  C     !INPUT PARAMETERS:
59  CEOP  CEOP
60    
61  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
62          INTEGER       diagNum
63        CHARACTER*8   diagName        CHARACTER*8   diagName
64        CHARACTER*16  diagCode        CHARACTER*16  diagCode
65        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
# Line 69  C     specify Name (cdiag, 8c), Descript Line 71  C     specify Name (cdiag, 8c), Descript
71  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
72  C----------------------------------------------------------------------  C----------------------------------------------------------------------
73    
74        diagName  = 'DRHODR  '  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
75        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '  C-    state variables of the main code (and related quadratic var):
76        diagUnits = 'kg/m4           '  
77        diagCode  = 'SM      LR      '        diagName  = 'ETAN    '
78        CALL DIAGNOSTICS_ADD2LIST(        diagTitle = 'Perturbation of Surface (pressure, height) (hPa,m)'
79          diagUnits = 'm or hPa        '
80          diagCode  = 'SM      M1      '
81          CALL DIAGNOSTICS_ADD2LIST( diagNum,
82         I          diagName, diagCode, diagUnits, diagTitle, myThid )
83    
84          diagName  = 'ETANSQ  '
85          diagTitle = 'Square of Perturbation of Sfc (hPa^2,m^2)'
86          diagUnits = 'm2 or hPa2      '
87          diagCode  = 'SM      M1      '
88          CALL DIAGNOSTICS_ADD2LIST( diagNum,
89       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
90    
91        diagName  = 'DETADT2 '        diagName  = 'DETADT2 '
92        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
93        diagUnits = 'm2 or Pa2 /sec2 '        diagUnits = 'm2 or Pa2 /sec2 '
94        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
95        CALL DIAGNOSTICS_ADD2LIST(        CALL DIAGNOSTICS_ADD2LIST( diagNum,
96         I          diagName, diagCode, diagUnits, diagTitle, myThid )
97    
98          diagName  = 'THETA   '
99          diagTitle = 'Potential Temperature (K)'
100          diagUnits = 'K               '
101          diagCode  = 'SM      MR      '
102          CALL DIAGNOSTICS_ADD2LIST( diagNum,
103         I          diagName, diagCode, diagUnits, diagTitle, myThid )
104    
105          diagName  = 'SALT    '
106          diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
107          diagUnits = 'g/kg            '
108          diagCode  = 'SM      MR      '
109          CALL DIAGNOSTICS_ADD2LIST( diagNum,
110         I          diagName, diagCode, diagUnits, diagTitle, myThid )
111    
112          diagName  = 'UVEL    '
113          diagTitle = 'Zonal Component of Velocity (m/s)'
114          diagUnits = 'm/sec           '
115          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
116          CALL DIAGNOSTICS_ADD2LIST( diagNum,
117         I          diagName, diagCode, diagUnits, diagTitle, myThid )
118    
119          diagName  = 'VVEL    '
120          diagTitle = 'Meridional Component of Velocity (m/s)'
121          diagUnits = 'm/sec           '
122          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
123          CALL DIAGNOSTICS_ADD2LIST( diagNum,
124         I          diagName, diagCode, diagUnits, diagTitle, myThid )
125    
126          diagName  = 'WVEL    '
127          diagTitle = 'Vertical Component of Velocity (m/s)'
128          diagUnits = 'm/sec           '
129          diagCode  = 'WM      LR      '
130          CALL DIAGNOSTICS_ADD2LIST( diagNum,
131         I          diagName, diagCode, diagUnits, diagTitle, myThid )
132    
133          diagName  = 'THETASQ '
134          diagTitle = 'Square of Potential Temperature (K^2)'
135          diagUnits = 'K2              '
136          diagCode  = 'SM      MR      '
137          CALL DIAGNOSTICS_ADD2LIST( diagNum,
138         I          diagName, diagCode, diagUnits, diagTitle, myThid )
139    
140          diagName  = 'SALTSQ  '
141          diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
142          diagUnits = 'g2/kg2          '
143          diagCode  = 'SM      MR      '
144          CALL DIAGNOSTICS_ADD2LIST( diagNum,
145         I          diagName, diagCode, diagUnits, diagTitle, myThid )
146    
147          diagName  = 'UVELSQ  '
148          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
149          diagUnits = 'm2/sec2         '
150          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
151          CALL DIAGNOSTICS_ADD2LIST( diagNum,
152         I          diagName, diagCode, diagUnits, diagTitle, myThid )
153    
154          diagName  = 'VVELSQ  '
155          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
156          diagUnits = 'm2/sec2         '
157          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
158          CALL DIAGNOSTICS_ADD2LIST( diagNum,
159         I          diagName, diagCode, diagUnits, diagTitle, myThid )
160    
161          diagName  = 'WVELSQ  '
162          diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
163          diagUnits = 'm2/sec2         '
164          diagCode  = 'WM      LR      '
165          CALL DIAGNOSTICS_ADD2LIST( diagNum,
166         I          diagName, diagCode, diagUnits, diagTitle, myThid )
167    
168          diagName  = 'UVELVVEL'
169          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
170          diagUnits = 'm2/sec2         '
171          WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
172          CALL DIAGNOSTICS_ADD2LIST( diagNum,
173         I          diagName, diagCode, diagUnits, diagTitle, myThid )
174    
175          diagName  = 'UVELMASS'
176          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
177          diagUnits = 'm/sec           '
178          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
179          CALL DIAGNOSTICS_ADD2LIST( diagNum,
180         I          diagName, diagCode, diagUnits, diagTitle, myThid )
181    
182          diagName  = 'VVELMASS'
183          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
184          diagUnits = 'm/sec           '
185          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
186          CALL DIAGNOSTICS_ADD2LIST( diagNum,
187         I          diagName, diagCode, diagUnits, diagTitle, myThid )
188    
189          diagName  = 'WVELMASS'
190          diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
191          diagUnits = 'm/sec           '
192          diagCode  = 'WM      MR      '
193          CALL DIAGNOSTICS_ADD2LIST( diagNum,
194         I          diagName, diagCode, diagUnits, diagTitle, myThid )
195    
196          diagName  = 'UTHMASS '
197          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'
198          diagUnits = 'm-K/sec         '
199          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
200          CALL DIAGNOSTICS_ADD2LIST( diagNum,
201         I          diagName, diagCode, diagUnits, diagTitle, myThid )
202    
203          diagName  = 'VTHMASS '
204          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'
205          diagUnits = 'm-K/sec         '
206          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
207          CALL DIAGNOSTICS_ADD2LIST( diagNum,
208         I          diagName, diagCode, diagUnits, diagTitle, myThid )
209    
210          diagName  = 'WTHMASS '
211          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'
212          diagUnits = 'm-K/sec         '
213          diagCode  = 'WM      MR      '
214          CALL DIAGNOSTICS_ADD2LIST( diagNum,
215         I          diagName, diagCode, diagUnits, diagTitle, myThid )
216    
217          diagName  = 'USLTMASS'
218          diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'
219          diagUnits = 'm-kg/s-kg       '
220          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
221          CALL DIAGNOSTICS_ADD2LIST( diagNum,
222         I          diagName, diagCode, diagUnits, diagTitle, myThid )
223    
224          diagName  = 'VSLTMASS'
225          diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'
226          diagUnits = 'm-kg/s-kg       '
227          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
228          CALL DIAGNOSTICS_ADD2LIST( diagNum,
229         I          diagName, diagCode, diagUnits, diagTitle, myThid )
230    
231          diagName  = 'WSLTMASS'
232          diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'
233          diagUnits = 'm-kg/s-kg       '
234          diagCode  = 'WM      MR      '
235          CALL DIAGNOSTICS_ADD2LIST( diagNum,
236         I          diagName, diagCode, diagUnits, diagTitle, myThid )
237    
238          diagName  = 'UVELTH  '
239          diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'
240          diagUnits = 'm-K/sec         '
241          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
242          CALL DIAGNOSTICS_ADD2LIST( diagNum,
243         I          diagName, diagCode, diagUnits, diagTitle, myThid )
244    
245          diagName  = 'VVELTH  '
246          diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'
247          diagUnits = 'm-K/sec         '
248          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
249          CALL DIAGNOSTICS_ADD2LIST( diagNum,
250         I          diagName, diagCode, diagUnits, diagTitle, myThid )
251    
252          diagName  = 'WVELTH  '
253          diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'
254          diagUnits = 'm-K/sec         '
255          diagCode  = 'WM      MR      '
256          CALL DIAGNOSTICS_ADD2LIST( diagNum,
257         I          diagName, diagCode, diagUnits, diagTitle, myThid )
258    
259          diagName  = 'UVELSLT '
260          diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'
261          diagUnits = 'm-kg/s-kg       '
262          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
263          CALL DIAGNOSTICS_ADD2LIST( diagNum,
264         I          diagName, diagCode, diagUnits, diagTitle, myThid )
265    
266          diagName  = 'VVELSLT '
267          diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'
268          diagUnits = 'm-kg/s-kg       '
269          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
270          CALL DIAGNOSTICS_ADD2LIST( diagNum,
271         I          diagName, diagCode, diagUnits, diagTitle, myThid )
272    
273          diagName  = 'WVELSLT '
274          diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'
275          diagUnits = 'm-kg/s-kg       '
276          diagCode  = 'WM      MR      '
277          CALL DIAGNOSTICS_ADD2LIST( diagNum,
278         I          diagName, diagCode, diagUnits, diagTitle, myThid )
279    
280    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
281    
282          diagName  = 'RHOAnoma'
283          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
284          diagUnits = 'kg/m^3          '
285          diagCode  = 'SM      MR      '
286          CALL DIAGNOSTICS_ADD2LIST( diagNum,
287         I          diagName, diagCode, diagUnits, diagTitle, myThid )
288    
289          diagName  = 'PHIHYD  '
290          diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
291          diagUnits = 'm2/s2           '
292          diagCode  = 'SM      MR      '
293          CALL DIAGNOSTICS_ADD2LIST( diagNum,
294         I          diagName, diagCode, diagUnits, diagTitle, myThid )
295    
296          diagName  = 'PHIBOT  '
297          diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
298          diagUnits = 'm2/s2           '
299          diagCode  = 'SM      M1      '
300          CALL DIAGNOSTICS_ADD2LIST( diagNum,
301         I          diagName, diagCode, diagUnits, diagTitle, myThid )
302    
303          diagName  = 'PHIBOTSQ'
304          diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
305          diagUnits = 'm4/s4           '
306          diagCode  = 'SM      M1      '
307          CALL DIAGNOSTICS_ADD2LIST( diagNum,
308         I          diagName, diagCode, diagUnits, diagTitle, myThid )
309    
310          diagName  = 'DRHODR  '
311          diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '
312          diagUnits = 'kg/m4           '
313          diagCode  = 'SM      LR      '
314          CALL DIAGNOSTICS_ADD2LIST( diagNum,
315         I          diagName, diagCode, diagUnits, diagTitle, myThid )
316    
317    C--   Momentum terms:
318          diagName  = 'VISCA4  '
319          diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
320          diagUnits = 'm4/sec          '
321          diagCode  = 'SM      MR      '
322          CALL DIAGNOSTICS_ADD2LIST( diagNum,
323         I          diagName, diagCode, diagUnits, diagTitle, myThid )
324    
325          diagName  = 'VISCAH  '
326          diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
327          diagUnits = 'm2/sec          '
328          diagCode  = 'SM      MR      '
329          CALL DIAGNOSTICS_ADD2LIST( diagNum,
330         I          diagName, diagCode, diagUnits, diagTitle, myThid )
331    
332    C--   surface fluxes:
333          diagName  = 'TAUX    '
334          diagTitle = 'zonal surface wind stress, >0 increases uVel    '
335          diagUnits = 'N/m2            '
336          diagCode  = 'SU      U1      '
337          CALL DIAGNOSTICS_ADD2LIST( diagNum,
338         I          diagName, diagCode, diagUnits, diagTitle, myThid )
339    
340          diagName  = 'TAUY    '
341          diagTitle = 'meridional surf. wind stress, >0 increases vVel '
342          diagUnits = 'N/m2            '
343          diagCode  = 'SV      U1      '
344          CALL DIAGNOSTICS_ADD2LIST( diagNum,
345         I          diagName, diagCode, diagUnits, diagTitle, myThid )
346    
347          diagName  = 'TFLUX   '
348          diagTitle = 'net surface heat flux, >0 increases theta       '
349          diagUnits = 'W/m2            '
350          diagCode  = 'SM      U1      '
351          CALL DIAGNOSTICS_ADD2LIST( diagNum,
352         I          diagName, diagCode, diagUnits, diagTitle, myThid )
353    
354          diagName  = 'TRELAX  '
355          diagTitle = 'surface temperature relaxation, >0 increases theta'
356          diagUnits = 'W/m2            '
357          diagCode  = 'SM      U1      '
358          CALL DIAGNOSTICS_ADD2LIST( diagNum,
359         I          diagName, diagCode, diagUnits, diagTitle, myThid )
360    
361          diagName  = 'TICE    '
362          diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
363          diagUnits = 'W/m2            '
364          diagCode  = 'SM      U1      '
365          CALL DIAGNOSTICS_ADD2LIST( diagNum,
366         I          diagName, diagCode, diagUnits, diagTitle, myThid )
367    
368          diagName  = 'SFLUX   '
369          diagTitle = 'net surface salt flux, >0 increases salt        '
370          diagUnits = 'g/m2/s          '
371          diagCode  = 'SM      U1      '
372          CALL DIAGNOSTICS_ADD2LIST( diagNum,
373         I          diagName, diagCode, diagUnits, diagTitle, myThid )
374    
375          diagName  = 'SRELAX  '
376          diagTitle = 'surface salinity relaxation, >0 increases salt  '
377          diagUnits = 'g/m2/s          '
378          diagCode  = 'SM      U1      '
379          CALL DIAGNOSTICS_ADD2LIST( diagNum,
380       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
381    
382  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22