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

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

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


Revision 1.6 - (hide annotations) (download)
Sat Dec 18 23:53:04 2004 UTC (19 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57b_post
Changes since 1.5: +222 -5 lines
move out off EQUIVALENCE all diagnostics that fizhi does not use

1 jmc 1.6 C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_main_init.F,v 1.5 2004/12/17 19:17:57 dimitri Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "DIAG_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP 0
8     C !ROUTINE: DIAGNOSTICS_MAIN_INIT
9    
10     C !INTERFACE:
11     SUBROUTINE DIAGNOSTICS_MAIN_INIT( myThid )
12    
13     C !DESCRIPTION:
14     C Initialize available diagnostics list for variables of the main code
15     C (not part of a package): set the following attributes:
16     C name (=cdiag), parsing code (=gdiag), units (=udiag), and title (=tdiag)
17     C Notes: 1) diagnostics defined here do not require any EQUIVALENCE
18     C since they get filled-in with S/R FILL_DIAGNOSTICS
19     C 2) GDIAG is defined as character*16 and can be to character*1
20     C parse(16) with the following codes currently defined:
21    
22     C \begin{center}
23     C \begin{tabular}[h]{|c|c|}\hline
24     C \textbf{Positions} & \textbf{Characters}
25     C & \textbf{Meanings} \\\hline
26     C parse(1) & S & scalar \\
27     C & U & U-vector component \\
28     C & V & V-vector component \\
29     C parse(2) & U & C-grid U-Point \\
30     C & V & C-grid V-Point \\
31     C & M & C-grid Mass Point \\
32     C & Z & C-grid Corner Point \\
33     C parse(3) & -- & --Unused-- \\
34     C parse(4) & P & positive definite \\
35     C parse(5 ) & C & with counter array \\
36     C & D & disable an array for output \\
37     C parse(6--8) & '123' & 3-digit mate number \\
38     C parse(9) & U & model-level plus 1/2 \\
39     C & M & model-level middle \\
40     C & L & model-level minus 1/2 \\
41     C parse(10) & 0 & levels = 0 \\
42     C & 1 & levels = 1 \\
43     C & R & levels = Nr \\
44     C & L & levels = MAX(Nr,NrPhys) \\
45     C & M & levels = MAX(Nr,NrPhys) - 1 \\
46     C \end{tabular}
47     C \end{center}
48    
49     C !USES:
50     IMPLICIT NONE
51     #include "SIZE.h"
52     #include "EEPARAMS.h"
53     c #include "EESUPPORT.h"
54     #include "PARAMS.h"
55    
56     C !INPUT PARAMETERS:
57     INTEGER myThid
58     CEOP
59    
60     C !LOCAL VARIABLES:
61 jmc 1.3 INTEGER diagNum
62 jmc 1.1 CHARACTER*8 diagName
63     CHARACTER*16 diagCode
64     CHARACTER*16 diagUnits
65     CHARACTER*(80) diagTitle
66    
67     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
68     C For each output variable,
69     C specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
70     C and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
71     C----------------------------------------------------------------------
72    
73 jmc 1.6 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 jmc 1.2 diagName = 'PHIHYD '
282     diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
283     diagUnits = 'm2/s2 '
284     diagCode = 'SM MR '
285 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
286 jmc 1.2 I diagName, diagCode, diagUnits, diagTitle, myThid )
287    
288     diagName = 'PHIBOT '
289     diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
290     diagUnits = 'm2/s2 '
291 dimitri 1.4 diagCode = 'SM M1 '
292 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
293 jmc 1.2 I diagName, diagCode, diagUnits, diagTitle, myThid )
294    
295     diagName = 'PHIBOTSQ'
296     diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
297     diagUnits = 'm4/s4 '
298     diagCode = 'SM M1 '
299 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
300 jmc 1.2 I diagName, diagCode, diagUnits, diagTitle, myThid )
301    
302 jmc 1.1 diagName = 'DRHODR '
303     diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit) '
304     diagUnits = 'kg/m4 '
305     diagCode = 'SM LR '
306 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
307 jmc 1.1 I diagName, diagCode, diagUnits, diagTitle, myThid )
308    
309 jmc 1.6 C-- Momentum terms:
310     diagName = 'VISCA4 '
311     diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s) '
312     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 jmc 1.3 CALL DIAGNOSTICS_ADD2LIST( diagNum,
322 jmc 1.1 I diagName, diagCode, diagUnits, diagTitle, myThid )
323    
324 jmc 1.6 C-- surface fluxes:
325 dimitri 1.5 diagName = 'TAUX '
326     diagTitle = 'zonal surface wind stress, >0 increases uVel '
327     diagUnits = 'N/m2 '
328     diagCode = 'SU U1 '
329     CALL DIAGNOSTICS_ADD2LIST( diagNum,
330     I diagName, diagCode, diagUnits, diagTitle, myThid )
331    
332     diagName = 'TAUY '
333     diagTitle = 'meridional surf. wind stress, >0 increases vVel '
334     diagUnits = 'N/m2 '
335     diagCode = 'SV U1 '
336     CALL DIAGNOSTICS_ADD2LIST( diagNum,
337     I diagName, diagCode, diagUnits, diagTitle, myThid )
338    
339     diagName = 'TFLUX '
340     diagTitle = 'net surface heat flux, >0 increases theta '
341     diagUnits = 'W/m2 '
342     diagCode = 'SM U1 '
343     CALL DIAGNOSTICS_ADD2LIST( diagNum,
344     I diagName, diagCode, diagUnits, diagTitle, myThid )
345    
346     diagName = 'TRELAX '
347     diagTitle = 'surface temperature relaxation, >0 increases theta'
348     diagUnits = 'W/m2 '
349     diagCode = 'SM U1 '
350     CALL DIAGNOSTICS_ADD2LIST( diagNum,
351     I diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353     diagName = 'TICE '
354     diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
355     diagUnits = 'W/m2 '
356     diagCode = 'SM U1 '
357     CALL DIAGNOSTICS_ADD2LIST( diagNum,
358     I diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360     diagName = 'SFLUX '
361     diagTitle = 'net surface salt flux, >0 increases salt '
362     diagUnits = 'g/m2/s '
363     diagCode = 'SM U1 '
364     CALL DIAGNOSTICS_ADD2LIST( diagNum,
365     I diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367     diagName = 'SRELAX '
368     diagTitle = 'surface salinity relaxation, >0 increases salt '
369     diagUnits = 'g/m2/s '
370     diagCode = 'SM U1 '
371     CALL DIAGNOSTICS_ADD2LIST( diagNum,
372     I diagName, diagCode, diagUnits, diagTitle, myThid )
373    
374 jmc 1.1 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
375    
376     RETURN
377     END

  ViewVC Help
Powered by ViewVC 1.1.22