/[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.6 by jmc, Sat Dec 18 23:53:04 2004 UTC
# Line 70  C     specify Name (cdiag, 8c), Descript Line 70  C     specify Name (cdiag, 8c), Descript
70  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
71  C----------------------------------------------------------------------  C----------------------------------------------------------------------
72    
73    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
74    C-    state variables of the main code (and related quadratic var):
75    
76          diagName  = 'ETAN    '
77          diagTitle = 'Perturbation of Surface (pressure, height) (hPa,m)'
78          diagUnits = 'm or hPa        '
79          diagCode  = 'SM      M1      '
80          CALL DIAGNOSTICS_ADD2LIST( diagNum,
81         I          diagName, diagCode, diagUnits, diagTitle, myThid )
82    
83          diagName  = 'ETANSQ  '
84          diagTitle = 'Square of Perturbation of Sfc (hPa^2,m^2)'
85          diagUnits = 'm2 or hPa2      '
86          diagCode  = 'SM      M1      '
87          CALL DIAGNOSTICS_ADD2LIST( diagNum,
88         I          diagName, diagCode, diagUnits, diagTitle, myThid )
89    
90          diagName  = 'DETADT2 '
91          diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
92          diagUnits = 'm2 or Pa2 /sec2 '
93          diagCode  = 'SM      M1      '
94          CALL DIAGNOSTICS_ADD2LIST( diagNum,
95         I          diagName, diagCode, diagUnits, diagTitle, myThid )
96    
97          diagName  = 'THETA   '
98          diagTitle = 'Potential Temperature (K)'
99          diagUnits = 'K               '
100          diagCode  = 'SM      MR      '
101          CALL DIAGNOSTICS_ADD2LIST( diagNum,
102         I          diagName, diagCode, diagUnits, diagTitle, myThid )
103    
104          diagName  = 'SALT    '
105          diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
106          diagUnits = 'g/kg            '
107          diagCode  = 'SM      MR      '
108          CALL DIAGNOSTICS_ADD2LIST( diagNum,
109         I          diagName, diagCode, diagUnits, diagTitle, myThid )
110    
111          diagName  = 'UVEL    '
112          diagTitle = 'Zonal Component of Velocity (m/s)'
113          diagUnits = 'm/sec           '
114          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
115          CALL DIAGNOSTICS_ADD2LIST( diagNum,
116         I          diagName, diagCode, diagUnits, diagTitle, myThid )
117    
118          diagName  = 'VVEL    '
119          diagTitle = 'Meridional Component of Velocity (m/s)'
120          diagUnits = 'm/sec           '
121          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
122          CALL DIAGNOSTICS_ADD2LIST( diagNum,
123         I          diagName, diagCode, diagUnits, diagTitle, myThid )
124    
125          diagName  = 'WVEL    '
126          diagTitle = 'Vertical Component of Velocity (m/s)'
127          diagUnits = 'm/sec           '
128          diagCode  = 'WM      LR      '
129          CALL DIAGNOSTICS_ADD2LIST( diagNum,
130         I          diagName, diagCode, diagUnits, diagTitle, myThid )
131    
132          diagName  = 'THETASQ '
133          diagTitle = 'Square of Potential Temperature (K^2)'
134          diagUnits = 'K2              '
135          diagCode  = 'SM      MR      '
136          CALL DIAGNOSTICS_ADD2LIST( diagNum,
137         I          diagName, diagCode, diagUnits, diagTitle, myThid )
138    
139          diagName  = 'SALTSQ  '
140          diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
141          diagUnits = 'g2/kg2          '
142          diagCode  = 'SM      MR      '
143          CALL DIAGNOSTICS_ADD2LIST( diagNum,
144         I          diagName, diagCode, diagUnits, diagTitle, myThid )
145    
146          diagName  = 'UVELSQ  '
147          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
148          diagUnits = 'm2/sec2         '
149          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
150          CALL DIAGNOSTICS_ADD2LIST( diagNum,
151         I          diagName, diagCode, diagUnits, diagTitle, myThid )
152    
153          diagName  = 'VVELSQ  '
154          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
155          diagUnits = 'm2/sec2         '
156          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
157          CALL DIAGNOSTICS_ADD2LIST( diagNum,
158         I          diagName, diagCode, diagUnits, diagTitle, myThid )
159    
160          diagName  = 'WVELSQ  '
161          diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
162          diagUnits = 'm2/sec2         '
163          diagCode  = 'WM      LR      '
164          CALL DIAGNOSTICS_ADD2LIST( diagNum,
165         I          diagName, diagCode, diagUnits, diagTitle, myThid )
166    
167          diagName  = 'UVELVVEL'
168          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
169          diagUnits = 'm2/sec2         '
170          WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
171          CALL DIAGNOSTICS_ADD2LIST( diagNum,
172         I          diagName, diagCode, diagUnits, diagTitle, myThid )
173    
174          diagName  = 'UVELMASS'
175          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
176          diagUnits = 'm/sec           '
177          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
178          CALL DIAGNOSTICS_ADD2LIST( diagNum,
179         I          diagName, diagCode, diagUnits, diagTitle, myThid )
180    
181          diagName  = 'VVELMASS'
182          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
183          diagUnits = 'm/sec           '
184          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
185          CALL DIAGNOSTICS_ADD2LIST( diagNum,
186         I          diagName, diagCode, diagUnits, diagTitle, myThid )
187    
188          diagName  = 'WVELMASS'
189          diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
190          diagUnits = 'm/sec           '
191          diagCode  = 'WM      MR      '
192          CALL DIAGNOSTICS_ADD2LIST( diagNum,
193         I          diagName, diagCode, diagUnits, diagTitle, myThid )
194    
195          diagName  = 'UTHMASS '
196          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'
197          diagUnits = 'm-K/sec         '
198          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
199          CALL DIAGNOSTICS_ADD2LIST( diagNum,
200         I          diagName, diagCode, diagUnits, diagTitle, myThid )
201    
202          diagName  = 'VTHMASS '
203          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'
204          diagUnits = 'm-K/sec         '
205          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
206          CALL DIAGNOSTICS_ADD2LIST( diagNum,
207         I          diagName, diagCode, diagUnits, diagTitle, myThid )
208    
209          diagName  = 'WTHMASS '
210          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'
211          diagUnits = 'm-K/sec         '
212          diagCode  = 'WM      MR      '
213          CALL DIAGNOSTICS_ADD2LIST( diagNum,
214         I          diagName, diagCode, diagUnits, diagTitle, myThid )
215    
216          diagName  = 'USLTMASS'
217          diagTitle = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'
218          diagUnits = 'm-kg/s-kg       '
219          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
220          CALL DIAGNOSTICS_ADD2LIST( diagNum,
221         I          diagName, diagCode, diagUnits, diagTitle, myThid )
222    
223          diagName  = 'VSLTMASS'
224          diagTitle = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'
225          diagUnits = 'm-kg/s-kg       '
226          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
227          CALL DIAGNOSTICS_ADD2LIST( diagNum,
228         I          diagName, diagCode, diagUnits, diagTitle, myThid )
229    
230          diagName  = 'WSLTMASS'
231          diagTitle = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'
232          diagUnits = 'm-kg/s-kg       '
233          diagCode  = 'WM      MR      '
234          CALL DIAGNOSTICS_ADD2LIST( diagNum,
235         I          diagName, diagCode, diagUnits, diagTitle, myThid )
236    
237          diagName  = 'UVELTH  '
238          diagTitle = 'Zonal Transp of Pot Temp (m-K/s)'
239          diagUnits = 'm-K/sec         '
240          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
241          CALL DIAGNOSTICS_ADD2LIST( diagNum,
242         I          diagName, diagCode, diagUnits, diagTitle, myThid )
243    
244          diagName  = 'VVELTH  '
245          diagTitle = 'Meridional Transp of Pot Temp (m-K/s)'
246          diagUnits = 'm-K/sec         '
247          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
248          CALL DIAGNOSTICS_ADD2LIST( diagNum,
249         I          diagName, diagCode, diagUnits, diagTitle, myThid )
250    
251          diagName  = 'WVELTH  '
252          diagTitle = 'Vertical Transp of Pot Temp (m-K/s)'
253          diagUnits = 'm-K/sec         '
254          diagCode  = 'WM      MR      '
255          CALL DIAGNOSTICS_ADD2LIST( diagNum,
256         I          diagName, diagCode, diagUnits, diagTitle, myThid )
257    
258          diagName  = 'UVELSLT '
259          diagTitle = 'Zonal Transp of Salt (m-kg/s-kg)'
260          diagUnits = 'm-kg/s-kg       '
261          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
262          CALL DIAGNOSTICS_ADD2LIST( diagNum,
263         I          diagName, diagCode, diagUnits, diagTitle, myThid )
264    
265          diagName  = 'VVELSLT '
266          diagTitle = 'Meridional Transp of Salt (m-kg/s-kg)'
267          diagUnits = 'm-kg/s-kg       '
268          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
269          CALL DIAGNOSTICS_ADD2LIST( diagNum,
270         I          diagName, diagCode, diagUnits, diagTitle, myThid )
271    
272          diagName  = 'WVELSLT '
273          diagTitle = 'Vertical Transp of Salt (m-kg/s-kg)'
274          diagUnits = 'm-kg/s-kg       '
275          diagCode  = 'WM      MR      '
276          CALL DIAGNOSTICS_ADD2LIST( diagNum,
277         I          diagName, diagCode, diagUnits, diagTitle, myThid )
278    
279    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
280    
281        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
282        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
283        diagUnits = 'm2/s2           '        diagUnits = 'm2/s2           '
# Line 98  C--------------------------------------- Line 306  C---------------------------------------
306        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
307       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
308    
309        diagName  = 'DETADT2 '  C--   Momentum terms:
310        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagName  = 'VISCA4  '
311        diagUnits = 'm2 or Pa2 /sec2 '        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
312        diagCode  = 'SM      M1      '        diagUnits = 'm4/sec          '
313          diagCode  = 'SM      MR      '
314          CALL DIAGNOSTICS_ADD2LIST( diagNum,
315         I          diagName, diagCode, diagUnits, diagTitle, myThid )
316    
317          diagName  = 'VISCAH  '
318          diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
319          diagUnits = 'm2/sec          '
320          diagCode  = 'SM      MR      '
321        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
322       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
323    
324    C--   surface fluxes:
325        diagName  = 'TAUX    '        diagName  = 'TAUX    '
326        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
327        diagUnits = 'N/m2            '        diagUnits = 'N/m2            '
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22