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

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

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

Revision 1.20 - (hide annotations) (download)
Tue Jul 6 22:49:13 2004 UTC (19 years, 11 months ago) by edhill
Branch: MAIN
Changes since 1.19: +2 -64 lines 
 o diagnostics are now producing NetCDF output
   - tested with hs94.cs-32x32x5 and seems to work well
   - still need some variable attributes:
     - list the levels
     - describe the time averaging periods where applicable
   - need the N+1 points for U and V values (or do we care?)
1 edhill 1.20 C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_init_vals.F,v 1.19 2004/07/06 05:22:48 edhill Exp $
2 edhill 1.16 C $Name: $
3    
4 edhill 1.18 #include "DIAG_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP 0
8 edhill 1.19 C !ROUTINE: DIAGNOSTICS_INIT_VALS
9 edhill 1.18
10     C !INTERFACE:
11 edhill 1.19 SUBROUTINE DIAGNOSTICS_INIT_VALS( myThid )
12 edhill 1.18
13     C !DESCRIPTION:
14     C Initialize Diagnostic Names, Levels, and GPARMS. Note that GDIAG
15     C is defined as character*16 and can be to character*1 parse(16)
16     C with the following codes currently defined:
17    
18 molod 1.1 C ***********************************************************************
19     C ***** ****
20     C ***** PURPOSE ****
21     C ***** Initialize Diagnostic Names, Levels, and GPARMS ****
22     C ***** ****
23     C ***** ****
24 molod 1.7 C ***** Note: GDIAG is defined as character*16 and can be ****
25 molod 1.1 C ***** equivalenced to character*1 parse(16) ****
26     C ***** ****
27     C ***** parse(1) => 'S' for a scalar diagnostic ****
28     C ***** => 'U' for a U-vector component diagnostic ****
29     C ***** => 'V' for a V-vector component diagnostic ****
30     C ***** ****
31     C ***** parse(2) => 'U' for a C-grid U-Point ****
32     C ***** => 'V' for a C-grid V-Point ****
33     C ***** => 'M' for a C-grid Mass Point ****
34 molod 1.7 C ***** => 'Z' for a C-grid Corner Point ****
35 molod 1.1 C ***** ****
36 molod 1.7 C ***** parse(3) => not used at present ****
37 molod 1.1 C ***** ****
38     C ***** parse(4) => 'P' if diagnostic is positive definite ****
39     C ***** ****
40     C ***** parse(5) => 'C' for a diagnostic with a counter array ****
41     C ***** => 'D' to disable an array for output ****
42     C ***** ****
43     C ***** parse(6-8) => 3-digit integer corresponding to vector ****
44     C ***** component mate if parse(1) = 'U' or 'V' ****
45 molod 1.7 C ***** or to the counter if parse(5)='C' ****
46 molod 1.1 C ***** ****
47 molod 1.7 C ***** parse(9) => 'U' if defined at Model-level plus 1/2 ****
48 molod 1.1 C ***** => 'M' if defined at Model-level middle ****
49 molod 1.7 C ***** => 'L' if defined at Model-level minus 1/2 ****
50 molod 1.1 C ***** ****
51     C ***** ****
52     C ***********************************************************************
53    
54 edhill 1.18 C \begin{center}
55     C \begin{tabular}[h]{|c|c|}\hline
56     C \textbf{Positions} & \textbf{Characters}
57     C & \textbf{Meanings} \\\hline
58     C parse(1) & S & scalar \\
59     C & U & U-vector component \\
60     C & V & V-vector component \\
61     C parse(2) & U & C-grid U-Point \\
62     C & V & C-grid V-Point \\
63     C & M & C-grid Mass Point \\
64     C & Z & C-grid Corner Point \\
65     C parse(3) & -- & --Unused-- \\
66     C parse(4) & P & positive definite \\
67     C parse(5 ) & C & with counter array \\
68     C & D & disable an array for output \\
69     C parse(6--8) & '123' & 3-digit mate number \\
70     C parse(9) & U & model-level plus 1/2 \\
71     C & M & model-level middle \\
72     C & L & model-level minus 1/2 \\
73     C parse(10) & 0 & levels = 0 \\
74     C & 1 & levels = 1 \\
75     C & R & levels = Nr \\
76     C & L & levels = MAX(Nr,NrPhys) \\
77     C & M & levels = MAX(Nr,NrPhys) - 1 \\
78     C \end{tabular}
79     C \end{center}
80    
81     C !USES:
82 molod 1.1 implicit none
83 molod 1.2 #include "SIZE.h"
84 edhill 1.19 #include "EEPARAMS.h"
85     #include "EESUPPORT.h"
86     #include "PARAMS.h"
87 molod 1.2
88 molod 1.4 #ifdef ALLOW_FIZHI
89 molod 1.2 #include "fizhi_SIZE.h"
90     #else
91     integer Nrphys
92 molod 1.6 parameter (Nrphys=0)
93 molod 1.2 #endif
94    
95 molod 1.1 #include "diagnostics_SIZE.h"
96     #include "diagnostics.h"
97    
98 edhill 1.18 C !INPUT PARAMETERS:
99 molod 1.2 integer myThid
100 edhill 1.18 CEOP
101 molod 1.17
102 edhill 1.18 integer n, nlevs
103    
104     nlevs = max(Nr,Nrphys)
105 molod 1.1
106 molod 1.5 do n=1,ndiagt
107 edhill 1.18 kdiag(n) = 0
108     cdiag(n) = ' '
109     gdiag(n) = ' '
110     udiag(n) = ' '
111     tdiag(n) = ' '
112 molod 1.5 enddo
113    
114 molod 1.1 c Diagnostic Names
115     c ----------------
116     cdiag( 1) = 'UFLUX '
117     cdiag( 2) = 'VFLUX '
118     cdiag( 3) = 'HFLUX '
119     cdiag( 4) = 'EFLUX '
120     cdiag( 5) = 'QICE '
121     cdiag( 6) = 'RADLWG '
122     cdiag( 7) = 'RADSWG '
123     cdiag( 8) = 'RI '
124     cdiag( 9) = 'CT '
125     cdiag(10) = 'CU '
126     cdiag(11) = 'KM '
127     cdiag(12) = 'KH '
128     cdiag(13) = 'TURBU '
129     cdiag(14) = 'TURBV '
130     cdiag(15) = 'TURBT '
131     cdiag(16) = 'TURBQ '
132     cdiag(17) = 'MOISTT '
133     cdiag(18) = 'MOISTQ '
134     cdiag(19) = 'RADLW '
135     cdiag(20) = 'RADSW '
136     cdiag(21) = 'PREACC '
137     cdiag(22) = 'PRECON '
138     cdiag(23) = 'TUFLUX '
139     cdiag(24) = 'TVFLUX '
140     cdiag(25) = 'TTFLUX '
141     cdiag(26) = 'TQFLUX '
142     cdiag(27) = 'CN '
143     cdiag(28) = 'WINDS '
144     cdiag(29) = 'DTSRF '
145     cdiag(30) = 'TGROUND '
146     cdiag(31) = 'TS '
147     cdiag(32) = 'DTG '
148     cdiag(33) = 'QG '
149     cdiag(34) = 'QS '
150     cdiag(35) = 'TGRLW '
151     cdiag(36) = 'LWGUP '
152     cdiag(37) = 'OLR '
153     cdiag(38) = 'OLRCLR '
154     cdiag(39) = 'LWGCLR '
155     cdiag(40) = 'LWCLR '
156     cdiag(41) = 'TLW '
157     cdiag(42) = 'SHRAD '
158     cdiag(43) = 'OZLW '
159     cdiag(44) = 'CLDRAS '
160     cdiag(45) = 'CLDTOT '
161     cdiag(46) = 'LWGDOWN '
162     cdiag(47) = 'GWDT '
163     cdiag(48) = 'RADSWT '
164     cdiag(49) = 'TAUCLD '
165     cdiag(50) = 'TAUCLDC '
166     cdiag(51) = 'CLDLOW '
167     cdiag(52) = 'EVAP '
168     cdiag(53) = 'DPDT '
169 molod 1.5 cdiag(54) = 'UAVE '
170     cdiag(55) = 'VAVE '
171     cdiag(56) = 'TAVE '
172     cdiag(57) = 'QAVE '
173 molod 1.1 cdiag(58) = 'OMEGA '
174     cdiag(59) = 'DUDT '
175     cdiag(60) = 'DVDT '
176     cdiag(61) = 'DTDT '
177     cdiag(62) = 'DQDT '
178     cdiag(63) = 'VORT '
179     cdiag(64) = 'PSIGDOT '
180     cdiag(65) = 'DTLS '
181     cdiag(66) = 'DQLS '
182     cdiag(67) = 'USTAR '
183     cdiag(68) = 'Z0 '
184     cdiag(69) = 'FRQTRB '
185     cdiag(70) = 'PBL '
186     cdiag(71) = 'SWCLR '
187     cdiag(72) = 'OSR '
188     cdiag(73) = 'OSRCLR '
189     cdiag(74) = 'CLDMAS '
190     cdiag(75) = 'UWND '
191     cdiag(76) = 'VWND '
192     cdiag(77) = 'TMPU '
193     cdiag(78) = 'SPHU '
194     cdiag(79) = 'RFT '
195     cdiag(80) = 'PS '
196     cdiag(81) = 'TKE '
197     cdiag(82) = 'SWGCLR '
198 molod 1.5 cdiag(83) = 'PAVE '
199 molod 1.1 cdiag(84) = 'SDIAG1 '
200     cdiag(85) = 'SDIAG2 '
201     cdiag(86) = 'UDIAG1 '
202     cdiag(87) = 'UDIAG2 '
203     cdiag(88) = 'DIABU '
204     cdiag(89) = 'DIABV '
205     cdiag(90) = 'DIABT '
206     cdiag(91) = 'DIABQ '
207     cdiag(92) = 'RFU '
208     cdiag(93) = 'RFV '
209     cdiag(94) = 'GWDU '
210     cdiag(95) = 'GWDV '
211     cdiag(96) = 'GWDUS '
212     cdiag(97) = 'GWDVS '
213     cdiag(98) = 'GWDUT '
214     cdiag(99) = 'GWDVT '
215     cdiag(100) = 'LZRAD '
216     cdiag(101) = 'SLP '
217     cdiag(102) = 'VAVEUQ '
218     cdiag(103) = 'VAVEVQ '
219     cdiag(104) = 'VAVEUT '
220     cdiag(105) = 'VAVEVT '
221     cdiag(106) = 'CLDFRC '
222     cdiag(107) = 'TPW '
223     cdiag(108) = 'U2M '
224     cdiag(109) = 'V2M '
225     cdiag(110) = 'T2M '
226     cdiag(111) = 'Q2M '
227     cdiag(112) = 'U10M '
228     cdiag(113) = 'V10M '
229     cdiag(114) = 'T10M '
230     cdiag(115) = 'Q10M '
231     cdiag(116) = 'DTRAIN '
232     cdiag(117) = 'QFILL '
233 molod 1.5 cdiag(118) = 'VAVEQAVE'
234 molod 1.1 cdiag(119) = 'VAVEQFIL'
235     cdiag(120) = 'SHAPU '
236     cdiag(121) = 'SHAPV '
237     cdiag(122) = 'SHAPT '
238     cdiag(123) = 'SHAPQ '
239     cdiag(124) = 'SDIAG3 '
240     cdiag(125) = 'SDIAG4 '
241     cdiag(126) = 'SDIAG5 '
242     cdiag(127) = 'SDIAG6 '
243     cdiag(128) = 'SDIAG7 '
244     cdiag(129) = 'SDIAG8 '
245     cdiag(130) = 'SDIAG9 '
246     cdiag(131) = 'SDIAG10 '
247     cdiag(132) = 'UDIAG3 '
248     cdiag(133) = 'UDIAG4 '
249     cdiag(134) = 'UDIAG5 '
250     cdiag(135) = 'UDIAG6 '
251     cdiag(136) = 'UDIAG7 '
252     cdiag(137) = 'UDIAG8 '
253     cdiag(138) = 'UDIAG9 '
254     cdiag(139) = 'UDIAG10 '
255     cdiag(140) = 'CLDMID '
256     cdiag(141) = 'CLDHI '
257     cdiag(142) = 'TAULOW '
258     cdiag(143) = 'TAUMID '
259     cdiag(144) = 'TAUHI '
260     cdiag(145) = 'CLDNP '
261     cdiag(146) = 'CLDPRS '
262     cdiag(147) = 'CTPCNT '
263     cdiag(148) = 'CLDTMP '
264     cdiag(149) = 'CTTCNT '
265     cdiag(150) = 'TAULOWC '
266     cdiag(151) = 'TAUMIDC '
267     cdiag(152) = 'TAUHIC '
268     cdiag(153) = 'TCANOPY '
269     cdiag(154) = 'TDEEP '
270     cdiag(155) = 'QCANOPY '
271     cdiag(156) = 'SMSHAL '
272     cdiag(157) = 'SMROOT '
273     cdiag(158) = 'SMDEEP '
274     cdiag(159) = 'CAPACITY'
275     cdiag(160) = 'SNOW '
276     cdiag(161) = 'RAINCON '
277     cdiag(162) = 'RAINLSP '
278     cdiag(163) = 'SNOWFALL'
279     cdiag(164) = 'RUNOFF '
280     cdiag(165) = 'FWSOIL '
281     cdiag(166) = 'GDRAIN '
282     cdiag(167) = 'SNOWMELT'
283     cdiag(168) = 'ERESV '
284     cdiag(169) = 'ESOIL '
285     cdiag(170) = 'EVEG '
286     cdiag(171) = 'ESNOW '
287     cdiag(172) = 'PARDF '
288     cdiag(173) = 'PARDR '
289     cdiag(174) = 'LAI '
290     cdiag(175) = 'GREEN '
291     cdiag(176) = 'DLWDTC '
292     cdiag(177) = 'DHDTC '
293     cdiag(178) = 'DEDTC '
294     cdiag(179) = 'VDTMOIST'
295     cdiag(180) = 'VDTTURB '
296     cdiag(181) = 'VDTLWRAD'
297     cdiag(182) = 'VDTSWRAD'
298 molod 1.5 cdiag(183) = 'VAVETAVE'
299 molod 1.1 cdiag(184) = 'PSUBCLD '
300     cdiag(185) = 'PSUBCLDC'
301     cdiag(186) = 'LCL '
302     cdiag(187) = 'SDIAGC '
303     cdiag(188) = 'SDIAGCC '
304     cdiag(189) = ' '
305     cdiag(190) = 'EVPOT '
306     cdiag(191) = 'RHCHECK '
307     cdiag(192) = 'DHDQA '
308     cdiag(193) = 'DEDQA '
309     cdiag(194) = 'DTC '
310     cdiag(195) = 'DQC '
311     cdiag(196) = 'TCDTC '
312     cdiag(197) = 'RADDTC '
313     cdiag(198) = 'SENDDTC '
314     cdiag(199) = 'LATDTC '
315     cdiag(200) = 'TDDTC '
316     cdiag(201) = 'QCDTC '
317     cdiag(202) = 'TRBQLIQ '
318     cdiag(203) = 'TRBFCC '
319     cdiag(204) = 'ALBEDO '
320     cdiag(205) = 'ALBVISDR'
321     cdiag(206) = 'ALBVISDF'
322     cdiag(207) = 'ALBNIRDR'
323     cdiag(208) = 'ALBNIRDF'
324     cdiag(209) = 'TAUAVE '
325     cdiag(210) = 'ALBEDOC '
326     cdiag(211) = 'RHCHECKT'
327     cdiag(212) = 'RHCHECKQ'
328     cdiag(213) = 'CLDLSP '
329     cdiag(214) = 'LZLOW '
330     cdiag(215) = 'LZMID '
331     cdiag(216) = 'LZHI '
332     cdiag(217) = 'LZTOT '
333     cdiag(218) = 'CAPE '
334     cdiag(219) = 'HSUBCL '
335     cdiag(220) = 'HSTAR '
336     cdiag(221) = 'CONT '
337     cdiag(222) = 'CONQ '
338     cdiag(223) = 'LST '
339     cdiag(224) = 'LSQ '
340     cdiag(225) = 'CONEVPT '
341     cdiag(226) = 'CONEVPQ '
342     cdiag(227) = 'LSEVPT '
343     cdiag(228) = 'LSEVPQ '
344     cdiag(229) = 'GWCUMID '
345     cdiag(230) = 'GWCVMID '
346     cdiag(231) = 'CLDSTD '
347     cdiag(232) = 'GWCUBAR '
348     cdiag(233) = 'GWCVBAR '
349     cdiag(234) = 'GWCUS '
350     cdiag(235) = 'GWCVS '
351     cdiag(236) = 'GWCU '
352     cdiag(237) = 'GWCV '
353 molod 1.7 cdiag(238) = 'ETAN '
354     cdiag(239) = 'ETANSQ '
355     cdiag(240) = 'THETA '
356     cdiag(241) = 'SALT '
357     cdiag(242) = 'UVEL '
358     cdiag(243) = 'VVEL '
359     cdiag(244) = 'WVEL '
360     cdiag(245) = 'THETASQ '
361     cdiag(246) = 'SALTSQ '
362     cdiag(247) = 'UVELSQ '
363     cdiag(248) = 'VVELSQ '
364     cdiag(249) = 'WVELSQ '
365     cdiag(250) = 'UVELVVEL'
366 molod 1.8 cdiag(251) = 'UVELMASS'
367     cdiag(252) = 'VVELMASS'
368     cdiag(253) = 'WVELMASS'
369     cdiag(254) = 'UTHMASS '
370     cdiag(255) = 'VTHMASS '
371     cdiag(256) = 'WTHMASS '
372     cdiag(257) = 'USLTMASS'
373     cdiag(258) = 'VSLTMASS'
374     cdiag(259) = 'WSLTMASS'
375 molod 1.9 cdiag(260) = 'UVELTH '
376     cdiag(261) = 'VVELTH '
377     cdiag(262) = 'WVELTH '
378     cdiag(263) = 'UVELSLT '
379     cdiag(264) = 'VVELSLT '
380     cdiag(265) = 'WVELSLT '
381 molod 1.10 cdiag(266) = 'UTRAC1 '
382     cdiag(267) = 'VTRAC1 '
383     cdiag(268) = 'WTRAC1 '
384     cdiag(269) = 'UTRAC2 '
385     cdiag(270) = 'VTRAC2 '
386     cdiag(271) = 'WTRAC2 '
387     cdiag(272) = 'UTRAC3 '
388     cdiag(273) = 'VTRAC3 '
389     cdiag(274) = 'WTRAC3 '
390     cdiag(275) = 'UTRAC4 '
391     cdiag(276) = 'VTRAC4 '
392     cdiag(277) = 'WTRAC4 '
393     cdiag(278) = 'UTRAC5 '
394     cdiag(279) = 'VTRAC5 '
395     cdiag(280) = 'WTRAC5 '
396 molod 1.12 cdiag(281) = 'TRAC1 '
397     cdiag(282) = 'TRAC2 '
398     cdiag(283) = 'TRAC3 '
399     cdiag(284) = 'TRAC4 '
400     cdiag(285) = 'TRAC5 '
401 molod 1.13 cdiag(286) = 'DICBIOA '
402     cdiag(287) = 'DICCARB '
403     cdiag(288) = 'DICTFLX '
404     cdiag(289) = 'DICOFLX '
405     cdiag(290) = 'DICCFLX '
406     cdiag(291) = 'DICPCO2 '
407     cdiag(292) = 'DICPHAV '
408 molod 1.14 cdiag(293) = 'DTCONV '
409     cdiag(294) = 'DQCONV '
410     cdiag(295) = 'RELHUM '
411     cdiag(296) = 'PRECLS '
412     cdiag(297) = 'ENPREC '
413 molod 1.1
414     c Diagnostic Descriptions
415 molod 1.2 c ----------------------------------------------------------------------
416 molod 1.1 tdiag( 1) = 'Zonal Wind Surface Stress (N/m^2)'
417     tdiag( 2) = 'Meridional Wind Surface Stress (N/m^2)'
418     tdiag( 3) = 'Sensible Heat Flux (pos.upwrd) (W/m^2)'
419     tdiag( 4) = 'Latent Heat Flux (pos.upwrd) (W/m^2)'
420     tdiag( 5) = 'Heat Conduction through Sea-Ice (W/m^2)'
421     tdiag( 6) = 'Net Upward Longwave Flux at the Ground (W/m^2)'
422     tdiag( 7) = 'Net Downward Shortwave Flux at the Ground (W/m^2)'
423     tdiag( 8) = 'Richardson Number (non-dimensional)'
424     tdiag( 9) = 'Surface Drag Coef for Scalars (non-dimensional)'
425     tdiag(10) = 'Surface Drag Coef for Winds (non-dimensional)'
426     tdiag(11) = 'Eddy Diffusivity Coefficient for Momentum (m^2/s)'
427     tdiag(12) = 'Eddy Diffusivity Coefficient for Scalars (m^2/s)'
428 molod 1.2 tdiag(13) = 'Zonal Wind Tendency from Turbulence (m/s/day)'
429     tdiag(14) = 'Meridional Wind Tendency from Turbulence (m/s/day)'
430     tdiag(15) =
431     . 'Temperature Tendency from Turbulence (K/day)'
432     tdiag(16) =
433     . 'Specific Humidity Tendency from Turbulence (g/kg/day)'
434     tdiag(17) =
435     . 'Temperature Tendency due to Moist Processes (K/day)'
436     tdiag(18) =
437     . 'Specific Humidity Tendency due to Moist Processes (g/kg/day)'
438     tdiag(19) =
439     . 'Temperature Tendency due to Longwave Radiation (K/day)'
440     tdiag(20) =
441     . 'Temperature Tendency due to Shortwave Radiation (K/day)'
442 molod 1.1 tdiag(21) = 'Total Precipitation (mm/day)'
443     tdiag(22) = 'Convective Precipitation (mm/day)'
444     tdiag(23) = 'Turbulent Flux of Zonal Momentum (N/m^2)'
445     tdiag(24) = 'Turbulent Flux of Meridional Momentum (N/m^2)'
446     tdiag(25) = 'Turbulent Flux of Sensible Heat (W/m^2)'
447     tdiag(26) = 'Turbulent Flux of Latent Heat (W/m^2)'
448     tdiag(27) = 'Neutral Drag Coefficient (non-dimensional)'
449     tdiag(28) = 'Surface Wind Speed (m/s)'
450     tdiag(29) = 'Air/Surface Virtual Temperature Difference (deg K)'
451     tdiag(30) = 'Ground Temperature (deg K)'
452     tdiag(31) = 'Surface Air Temperature (deg K)'
453     tdiag(32) = 'Ground Temperature Adjustment (deg K)'
454     tdiag(33) = 'Ground Specific Humidity (g/kg)'
455     tdiag(34) = 'Saturation Surface Specific Humidity (g/kg)'
456     tdiag(35) = 'Ground Temperature used by LW Radiation (deg K)'
457     tdiag(36) = 'Upward Longwave Radiation at the Ground (W/m^2)'
458     tdiag(37) = 'Outgoing Longwave Radiation (W/m^2)'
459     tdiag(38) = 'Clear Sky Outgoing Longwave Radiation (W/m^2)'
460     tdiag(39) = 'Clear Sky Net Longwave Flux at the Ground (W/m^2)'
461     tdiag(40) = 'Clear Sky Longwave Radiation Heating Rates (K/day)'
462     tdiag(41) = 'Temperature used by Radiation (deg K)'
463     tdiag(42) = 'Specific Humidity used by Radiation (g/kg)'
464 molod 1.2 tdiag(43) =
465     . 'Ozone mixing ratio used by Longwave Radiation (non-dimensional)'
466 molod 1.1 tdiag(44) = 'Convective Cloud Fraction (0-1)'
467     tdiag(45) = '3-D Total Cloud Fraction (0-1)'
468     tdiag(46) = 'Downward Longwave Radiation at the Ground (W/m^2)'
469 molod 1.2 tdiag(47) =
470     . 'Temperature Tendency due to Gravity Wave Drag (K/day)'
471 molod 1.1 tdiag(48) = 'Incident Shortwave Radiation at TOA (W/m^2)'
472 molod 1.2 tdiag(49) =
473     . 'Counted Cloud Optical Depth (non-dimensional) per 100 mb'
474 molod 1.1 tdiag(50) = 'Cloud Optical Depth Counter'
475     tdiag(51) = 'Low-Level ( 1000-700 hPa) Cloud Fraction (0-1)'
476     tdiag(52) = 'Surface Evaporation (mm/day)'
477     tdiag(53) = 'Surface Pressure Tendency (hPa/day)'
478 molod 1.5 tdiag(54) = 'Average Zonal Wind (m/s)'
479     tdiag(55) = 'Average Meridional Wind (m/s)'
480     tdiag(56) = 'Average Potential Temperature (K)'
481 molod 1.2 tdiag(57) =
482 molod 1.5 . 'Average Specific Humidity (g/kg)'
483 molod 1.1 tdiag(58) = 'Vertical Velocity (hPa/day)'
484     tdiag(59) = 'Total Zonal Wind Tendency (m/s/day)'
485     tdiag(60) = 'Total Meridional Wind Tendency (m/s/day)'
486     tdiag(61) = 'Total Temperature Tendency (K/day)'
487     tdiag(62) = 'Total Specific Humidity Tendency (g/kg/day)'
488     tdiag(63) = 'Relative Vorticity (x10^-4 sec^-1)'
489     tdiag(64) = 'Pi-Weighted Vertical Sigma Velocity (hPa/day)'
490 molod 1.2 tdiag(65) =
491     . 'Temperature Tendency due to Stratiform Processes (K/day)'
492     tdiag(66) =
493     . 'Specific Humidity Tendency from Stratiform Processes (g/kg/day)'
494 molod 1.1 tdiag(67) = 'Friction Velocity (m/s)'
495     tdiag(68) = 'Surface Roughness Z0 (m)'
496     tdiag(69) = 'Frequency of Turbulence (0-1)'
497     tdiag(70) = 'Estimated PBL Depth (hPa)'
498     tdiag(71) = 'Clear Sky Shortwave Heating Rates (K/day)'
499     tdiag(72) = 'Outgoing Shortwave Radiation (W/m^2)'
500     tdiag(73) = 'Clear Sky Outgoing Shortwave Radiation (W/m^2)'
501     tdiag(74) = 'Cloud Mass Flux (g/m^2/sec)'
502     tdiag(75) = 'Time-Averaged Zonal Wind (m/s)'
503     tdiag(76) = 'Time-Averaged Meridional Wind (m/s)'
504     tdiag(77) = 'Time-Averaged Temperature (deg K)'
505     tdiag(78) = 'Time-Averaged Specific Humidity (g/kg)'
506 molod 1.2 tdiag(79) =
507     . 'Temperature Tendency due to Rayleigh Friction (K/day)'
508 molod 1.1 tdiag(80) = 'Time-Averaged Surface Pressure (hPa)'
509     tdiag(81) = 'Time-Averaged Turbulent Kinetic Energy (m/s)^2'
510 molod 1.2 tdiag(82) =
511     . 'Clear Sky Net Downward SW Radiation at the Ground (W/m^2)'
512 molod 1.5 tdiag(83) = 'Average Surface Pressure (hPa)'
513 molod 1.1 tdiag(84) = 'User-Defined Surface Diagnostic #1'
514     tdiag(85) = 'User-Defined Surface Diagnostic #2'
515     tdiag(86) = 'User-Defined Upper-Air Diagnostic #1'
516     tdiag(87) = 'User-Defined Upper-Air Diagnostic #2'
517 molod 1.2 tdiag(88) =
518     . 'Zonal Wind Tendency due to Diabatic Forcing (m/s/day)'
519     tdiag(89) =
520     . 'Meridional Wind Tendency due to Diabatic Forcing (m/s/day)'
521     tdiag(90) =
522     . 'Temperature Tendency due to Diabatic Forcing (K/day)'
523     tdiag(91) =
524     . 'Specific Humidity Tendency due to Diabatic Forcing (g/kg/day)'
525     tdiag(92) =
526     . 'Zonal Wind Tendency due to Rayleigh Friction (m/s/day)'
527     tdiag(93) =
528     . 'Meridional Wind Tendency due to Rayleigh Friction (m/s/day)'
529     tdiag(94) =
530     . 'Zonal Wind Tendency due to Gravity Wave Drag (m/s/day)'
531     tdiag(95) =
532     . 'Meridional Wind Tendency due to Gravity Wave Drag (m/s/day)'
533 molod 1.1 tdiag(96) = 'Zonal Wind Gravity Wave Surface Stress (N/m^2)'
534     tdiag(97) = 'Meridional Wind Gravity Wave Surface Stress (N/m^2)'
535     tdiag(98) = 'Zonal Wind Gravity Wave Stress at Model Top (N/m^2)'
536 molod 1.2 tdiag(99) =
537     . 'Meridional Wind Gravity Wave Stress at Model Top (N/m^2)'
538     tdiag(100) =
539     . 'Estimated Cloud Liquid Water used in Radiation (mg/kg)'
540 molod 1.1 tdiag(101) = 'Time-Averaged Sea-Level Pressure (mb)'
541     tdiag(102) = 'Vertically Averaged UWND*SPHU (m/s g/kg)'
542     tdiag(103) = 'Vertically Averaged VWND*SPHU (m/s g/kg)'
543     tdiag(104) = 'Vertically Averaged UWND*TMPU (m/s deg K)'
544     tdiag(105) = 'Vertically Averaged VWND*TMPU (m/s deg K)'
545     tdiag(106) = '2-D Total Cloud Fraction (0-1)'
546     tdiag(107) = 'Total Precipitable Water (g/cm^2)'
547     tdiag(108) = 'Zonal Wind Interpolated to 2 Meters (m/s)'
548     tdiag(109) = 'Meridional Wind Interpolated to 2 Meters (m/s)'
549     tdiag(110) = 'Temperature Interpolated to 2 Meters (deg K)'
550     tdiag(111) = 'Specific Humidity Interpolated to 2 Meters (g/kg)'
551     tdiag(112) = 'Zonal Wind Interpolated to 10 Meters (m/s)'
552     tdiag(113) = 'Meridional Wind Interpolated to 10 Meters (m/s)'
553     tdiag(114) = 'Temperature Interpolated to 10 Meters (deg K)'
554     tdiag(115) = 'Specific Humidity Interpolated to 10 Meters (g/kg)'
555     tdiag(116) = 'Detrainment Cloud Mass Flux (g/m^2/sec)'
556     tdiag(117) = 'Filling of Negative Specific Humidity (g/kg/day)'
557 molod 1.5 tdiag(118) = 'Vertically Averaged QAVE (mm/day)'
558 molod 1.1 tdiag(119) = 'Vertically Averaged QFILL (mm/day)'
559 molod 1.2 tdiag(120) =
560     . 'Zonal Wind Tendency due to Shapiro Filter (m/s/day)'
561     tdiag(121) =
562     . 'Meridional Wind Tendency due to Shapiro Filter (m/s/day)'
563 molod 1.1 tdiag(122) = 'Temperature Tendency due to Shapiro Filter (K/day)'
564 molod 1.2 tdiag(123) =
565     . 'Specific Humidity Tendency due to Shapiro Filter (g/kg/day)'
566 molod 1.1 tdiag(124) = 'User-Defined Surface Diagnostic #3 '
567     tdiag(125) = 'User-Defined Surface Diagnostic #4 '
568     tdiag(126) = 'User-Defined Surface Diagnostic #5 '
569     tdiag(127) = 'User-Defined Surface Diagnostic #6 '
570     tdiag(128) = 'User-Defined Surface Diagnostic #7 '
571     tdiag(129) = 'User-Defined Surface Diagnostic #8 '
572     tdiag(130) = 'User-Defined Surface Diagnostic #9 '
573     tdiag(131) = 'User-Defined Surface Diagnostic #10 '
574     tdiag(132) = 'User-Defined Upper-Air Diagnostic #3 '
575     tdiag(133) = 'User-Defined Upper-Air Diagnostic #4 '
576     tdiag(134) = 'User-Defined Upper-Air Diagnostic #5 '
577     tdiag(135) = 'User-Defined Upper-Air Diagnostic #6 '
578     tdiag(136) = 'User-Defined Upper-Air Diagnostic #7 '
579     tdiag(137) = 'User-Defined Upper-Air Diagnostic #8 '
580     tdiag(138) = 'User-Defined Upper-Air Diagnostic #9 '
581     tdiag(139) = 'User-Defined Upper-Air Diagnostic #10 '
582     tdiag(140) = 'Mid-Level (700-400 hPa) Cloud Fraction (0-1)'
583     tdiag(141) = 'High-Level (above 400 hPa) Cloud Fraction (0-1)'
584     tdiag(142) = 'Low-Level (1000-700 hPa) Optical Depth'
585     tdiag(143) = 'Mid-Level (700-400 hPa) Optical Depth'
586     tdiag(144) = 'High-Level (above 400 hPa) Optical Depth'
587     tdiag(145) = 'Non-Precipitating Cloud Fraction (0-1)'
588     tdiag(146) = 'Cloud Top Pressure (when cloudy) (hPa)'
589     tdiag(147) = 'Cloud Top Pressure Counter'
590     tdiag(148) = 'Cloud Top Temperature (when cloudy) (deg K)'
591     tdiag(149) = 'Cloud Top Temperature Counter'
592     tdiag(150) = 'Low-Level Optical Depth Counter'
593     tdiag(151) = 'Mid-Level Optical Depth Counter'
594     tdiag(152) = 'High-Level Optical Depth Counter'
595 molod 1.2 tdiag(153) =
596     . 'Skin Temperature of the Canopy/Soil (SST over water) (deg K)'
597 molod 1.1 tdiag(154) = 'Deep Soil Temperature (deg K)'
598     tdiag(155) = 'Specific Humidity of the Canopy (g/kg)'
599 molod 1.2 tdiag(156) =
600     . 'Soil Moisture Field Capacity Fraction in Shallow (5 mm) Layer'
601     tdiag(157) = 'Soil Moisture Field Capacity Fraction in Root Zone'
602     tdiag(158) =
603     . 'Soil Moisture Field Capacity Fraction in Recharge Layer '
604 molod 1.1 tdiag(159) = 'Canopy Interception Reservoir (mm)'
605     tdiag(160) = 'Snow Depth (mm water equivalent)'
606     tdiag(161) = 'Convective Rainfall (liquid precipitate) (mm/day)'
607     tdiag(162) = 'Large-scale Rainfall (liquid precipitate) (mm/day)'
608     tdiag(163) = 'Total Snowfall (solid precipitate) (mm/day)'
609 molod 1.2 tdiag(164) =
610     . 'Water from precipitation not infiltrated into soil (mm/s)'
611 molod 1.1 tdiag(165) = 'Rainwater Infiltration into top soil layer (mm/s)'
612 molod 1.2 tdiag(166) =
613     . 'Diffusion of Moisture across bottom of root zone (mm/s)'
614 molod 1.1 tdiag(167) = 'Rate of Snow Melt (mm/s)'
615 molod 1.2 tdiag(168) =
616     . 'Latent Heat Flux from Canopy Interception Reservoir (W/m^2)'
617 molod 1.1 tdiag(169) = 'Latent Heat Flux from Bare Soil (W/m^2)'
618 molod 1.2 tdiag(170) =
619     . 'Latent Heat Flux (Transpiration) from Vegetated Surface (W/m^2)'
620 molod 1.1 tdiag(171) = 'Latent Heat Flux from Snow Pack (W/m^2)'
621 molod 1.2 tdiag(172) =
622     . 'Diffuse-beam Photosynthetically Active Radiation (W/m^2)'
623     tdiag(173) =
624     . 'Direct-beam Photosynthetically Active Radiation (W/m^2)'
625 molod 1.1 tdiag(174) = 'Leaf Area Index (%)'
626     tdiag(175) = 'Green-ness Index (%)'
627 molod 1.2 tdiag(176) =
628     . 'Derivative of Upward LW Radiation wrt TCANOPY (W/m^2 /deg K)'
629     tdiag(177) =
630     . 'Derivative of Sensible Heat Flux wrt TCANOPY (W/m^2 /deg K)'
631     tdiag(178) =
632     . 'Derivative of Latent Heat Flux wrt TCANOPY (W/m^2 /deg K)'
633 molod 1.1 tdiag(179) = 'Vertically Averaged MOISTT (K/day)'
634     tdiag(180) = 'Vertically Averaged TURBT (K/day)'
635     tdiag(181) = 'Vertically Averaged RADLW (K/day)'
636     tdiag(182) = 'Vertically Averaged RADSW (K/day)'
637 molod 1.5 tdiag(183) = 'Vertically Averaged TAVE (K/day)'
638 molod 1.1 tdiag(184) = 'Pressure Thickness of the RAS Subcloud Layer (mb)'
639 molod 1.2 tdiag(185) =
640     . 'Pressure Thickness of the RAS Subcloud Layer Counter'
641     tdiag(186) = 'Lifting Condensation Level (mb) '
642     tdiag(187) = 'User-Defined Counted Surface Diagnostic '
643     tdiag(188) = 'User-Defined Counted Surface Diagnostic Counter'
644     tdiag(189) = ' '
645 molod 1.1 tdiag(190) = 'Potential Evapotranspiration (W/m^2)'
646 molod 1.2 tdiag(191) =
647     . 'Net precip due to supersaturation above critical RH (mm/day)'
648     tdiag(192) =
649     . 'Derivative of Sensible Heat wrt Canopy Humidity (W/m^2)'
650     tdiag(193) =
651     . 'Derivative of Latent Heat wrt Canopy Humidity (W/m^2)'
652 molod 1.1 tdiag(194) = 'Total Change in Canopy Temperature (deg/sec)'
653     tdiag(195) = 'Change in Canopy Humidity (kg/kg/sec)'
654 molod 1.2 tdiag(196) =
655     . 'Change in Canopy Temperature - Energy Term (deg/sec)'
656     tdiag(197) =
657     . 'Change in Canopy Temperature - Net Radiation Term (deg/sec)'
658     tdiag(198) =
659     . 'Change in Canopy Temperature - Sensible Heat Term (deg/sec)'
660     tdiag(199) =
661     . 'Change in Canopy Temperature - Latent Heat Term (deg/sec)'
662     tdiag(200) =
663     . 'Change in Canopy Temperature - Ground Temp Term (deg/sec)'
664     tdiag(201) =
665     . 'Change in Canopy Temperature - Humidity Term (deg/sec)'
666 molod 1.1 tdiag(202) = 'Cloud Liquid Water Generated in Turbulence (mg/kg)'
667     tdiag(203) = 'Cloud Fraction Generated in Turbulence (fraction)'
668     tdiag(204) = 'Mean Surface Albedo (SW_Up/SW_Down) 0-1 '
669     tdiag(205) = 'Direct Beam VIS Surface Albedo (0-1) '
670     tdiag(206) = 'Diffuse Beam VIS Surface Albedo (0-1) '
671     tdiag(207) = 'Direct Beam NIR Surface Albedo (0-1) '
672     tdiag(208) = 'Diffuse Beam NIR Surface Albedo (0-1) '
673 molod 1.2 tdiag(209) =
674     . 'Averaged Cloud Optical Depth (non-dimensional) per 100 mb'
675 molod 1.1 tdiag(210) = 'Surface Albedo Counter'
676 molod 1.2 tdiag(211) =
677     . 'Temperature Change due to humidity above critical RH (K/day)'
678     tdiag(212) =
679     . 'Spec.Humid. Change due to humidity above critical RH (g/kg/day)'
680 molod 1.1 tdiag(213) = 'Large-Scale Cloud Fraction (0-1)'
681 molod 1.2 tdiag(214) =
682     . 'Vert integrated low-level (below 700mb) cloud water (mg/cm**2)'
683     tdiag(215) =
684     . 'Vert integrated mid-level (700-400mb) cloud water (mg/cm**2)'
685     tdiag(216) =
686     . 'Vert integrated high-level (above 400mb) cloud water (mg/cm**2)'
687     tdiag(217) =
688     . 'Vertically integrated total suspended cloud water (mg/cm**2)'
689 molod 1.1 tdiag(218) = 'Convective Available Potential Energy'
690     tdiag(219) = 'Moist Static Energy in RAS subcloud layer'
691     tdiag(220) = 'Saturation Moist Static Energy'
692 molod 1.2 tdiag(221) =
693     . 'Temperature Tendency due to Convective Precipitation (K/day)'
694     tdiag(222) =
695     . 'Moisture Tendency due to Convective Precipitation (g/kg/day)'
696     tdiag(223) =
697     . 'Temperature Tendency due to Large Scale Precipitation (K/day)'
698     tdiag(224) =
699     . 'Moisture Tendency due to Large Scale Precipitation (g/kg/day)'
700     tdiag(225) =
701     . 'Temperature Tendency due to Conv Rain Re-evaporation (K/day)'
702     tdiag(226) =
703     . 'Moisture Tendency due to Conv Rain Re-evaporation (g/kg/day)'
704     tdiag(227) =
705     . 'Temperature Tendency due to Large Scale Rain Re-evap (K/day)'
706     tdiag(228) =
707     . 'Moisture Tendency due to Large Scale Rain Re-evap (g/kg/day)'
708     tdiag(229) =
709     . 'Mid-Lev Mean U-Wind for Gravity Wave Drag w/Clouds (m/sec)'
710     tdiag(230) =
711     . 'Mid-Lev Mean V-Wind for Gravity Wave Drag w/Clouds (m/sec)'
712     tdiag(231) =
713     . 'Mid-Lev Cloud Perturb Length for Gravity Wave Drag w/Clouds (m)'
714     tdiag(232) =
715     . 'Mid-Lev Perturb U-Wind for Gravity Wave Drag w/Clouds (m/sec)'
716     tdiag(233) =
717     . 'Mid-Lev Perturb V-Wind for Gravity Wave Drag w/Clouds (m/sec)'
718     tdiag(234) =
719     . 'Mid-Lev U-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)'
720     tdiag(235) =
721     . 'Mid-Lev V-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)'
722     tdiag(236) =
723     . 'U-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
724     tdiag(237) =
725     . 'V-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
726 molod 1.7 tdiag(238) = 'Perturbation of Surface (pressure, height) (hPa,m)'
727     tdiag(239) = 'Square of Perturbation of Sfc (hPa^2,m^2)'
728     tdiag(240) = 'Potential Temperature (K)'
729     tdiag(241) = 'Salt or Water Vapor mixing ratio (g/kg)'
730     tdiag(242) = 'Zonal Component of Velocity (m/s)'
731     tdiag(243) = 'Meridional Component of Velocity (m/s)'
732     tdiag(244) = 'Vertical Component of Velocity (m/s)'
733     tdiag(245) = 'Square of Potential Temperature (K^2)'
734     tdiag(246) = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
735     tdiag(247) = 'Square of Zonal Comp of Velocity (m^2/s^2)'
736     tdiag(248) = 'Square of Meridional Comp of Velocity (m^2/s^2)'
737     tdiag(249) = 'Square of Vertical Comp of Velocity (m^2/s^2)'
738     tdiag(250) = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
739 molod 1.8 tdiag(251) = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
740     tdiag(252) = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
741     tdiag(253) = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
742     tdiag(254) = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'
743     tdiag(255) = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'
744     tdiag(256) = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'
745     tdiag(257) = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'
746     tdiag(258) = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'
747     tdiag(259) = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'
748 molod 1.9 tdiag(260) = 'Zonal Transp of Pot Temp (m-K/s)'
749     tdiag(261) = 'Meridional Transp of Pot Temp (m-K/s)'
750     tdiag(262) = 'Vertical Transp of Pot Temp (m-K/s)'
751     tdiag(263) = 'Zonal Transp of Salt (m-kg/s-kg)'
752     tdiag(264) = 'Meridional Transp of Salt (m-kg/s-kg)'
753     tdiag(265) = 'Vertical Transp of Salt (m-kg/s-kg)'
754 molod 1.10 tdiag(266) = 'Zonal Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
755     tdiag(267) = 'Merid Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
756     tdiag(268) = 'Vert Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
757     tdiag(269) = 'Zonal Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
758     tdiag(270) = 'Merid Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
759     tdiag(271) = 'Vert Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
760     tdiag(272) = 'Zonal Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
761     tdiag(273) = 'Merid Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
762     tdiag(274) = 'Vert Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
763     tdiag(275) = 'Zonal Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
764     tdiag(276) = 'Merid Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
765     tdiag(277) = 'Vert Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
766     tdiag(278) = 'Zonal Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
767     tdiag(279) = 'Merid Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
768     tdiag(280) = 'Vert Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
769 molod 1.13 tdiag(281) = 'Mass-Weighted Tracer 1 ()'
770     tdiag(282) = 'Mass-Weighted Tracer 2 ()'
771     tdiag(283) = 'Mass-Weighted Tracer 3 ()'
772     tdiag(284) = 'Mass-Weighted Tracer 4 ()'
773     tdiag(285) = 'Mass-Weighted Tracer 5 ()'
774     tdiag(286) = 'Biological Productivity (mol/m3/s)'
775     tdiag(287) = 'Carbonate chg-biol prod and remin (mol eq/m3/s)'
776     tdiag(288) = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
777     tdiag(289) = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
778     tdiag(290) = 'Flux of CO2 - air-sea exch (mol/m2/s)'
779     tdiag(291) = 'Partial Pressure of CO2 (atm)'
780     tdiag(292) = 'pH (dimensionless)'
781 molod 1.14 tdiag(293) = 'Temp Change due to Convection (K/sec)'
782     tdiag(294) = 'Humidity Change due to Convection (g/kg/sec)'
783     tdiag(295) = 'Relative Humidity (fraction)'
784     tdiag(296) = 'Large Scale Precipitation (g/m2/sec)'
785     tdiag(297) = 'Energy of Precipitation (snow, rain Temp) in (J/g)'
786 molod 1.1
787     c Diagnostic Units
788     c ----------------
789     udiag( 1) = 'N/m2 '
790     udiag( 2) = 'N/m2 '
791     udiag( 3) = 'watt/m2 '
792     udiag( 4) = 'watt/m2 '
793     udiag( 5) = 'watt/m2 '
794     udiag( 6) = 'watt/m2 '
795     udiag( 7) = 'watt/m2 '
796     udiag( 8) = 'non-dimensional '
797     udiag( 9) = 'non-dimensional '
798     udiag(10) = 'non-dimensional '
799     udiag(11) = 'm2/sec '
800     udiag(12) = 'm2/sec '
801     udiag(13) = 'm/s/day '
802     udiag(14) = 'm/s/day '
803     udiag(15) = 'K/day '
804     udiag(16) = 'gram/kg/day '
805     udiag(17) = 'K/day '
806     udiag(18) = 'gram/kg/day '
807     udiag(19) = 'K/day '
808     udiag(20) = 'K/day '
809     udiag(21) = 'mm/day '
810     udiag(22) = 'mm/day '
811     udiag(23) = 'N/m2 '
812     udiag(24) = 'N/m2 '
813     udiag(25) = 'watt/m2 '
814     udiag(26) = 'watt/m2 '
815     udiag(27) = 'non-dimensional '
816     udiag(28) = 'm/sec '
817     udiag(29) = 'K '
818     udiag(30) = 'K '
819     udiag(31) = 'K '
820     udiag(32) = 'K '
821     udiag(33) = 'gram/kg '
822     udiag(34) = 'gram/kg '
823     udiag(35) = 'K '
824     udiag(36) = 'watt/m2 '
825     udiag(37) = 'watt/m2 '
826     udiag(38) = 'watt/m2 '
827     udiag(39) = 'watt/m2 '
828     udiag(40) = 'K/day '
829     udiag(41) = 'K '
830     udiag(42) = 'gram/kg '
831     udiag(43) = 'non-dimensional '
832     udiag(44) = '0-1 '
833     udiag(45) = '0-1 '
834     udiag(46) = 'watt/m2 '
835     udiag(47) = 'K/day '
836     udiag(48) = 'watt/m2 '
837     udiag(49) = 'non-dimensional '
838     udiag(50) = 'count '
839     udiag(51) = '0-1 '
840     udiag(52) = 'mm/day '
841     udiag(53) = 'hPa/day '
842 molod 1.5 udiag(54) = 'm/s '
843     udiag(55) = 'm/s '
844     udiag(56) = 'K '
845     udiag(57) = 'gram/kg '
846 molod 1.1 udiag(58) = 'hPa/day '
847     udiag(59) = 'm/s/day '
848     udiag(60) = 'm/s/day '
849     udiag(61) = 'K/day '
850     udiag(62) = 'gram/kg/day '
851     udiag(63) = '(e-4 sec-1) '
852     udiag(64) = 'hPa/day '
853     udiag(65) = 'K/day '
854     udiag(66) = 'gram/kg/day '
855     udiag(67) = 'm/sec '
856     udiag(68) = 'm '
857     udiag(69) = '0-1 '
858     udiag(70) = 'hPa '
859     udiag(71) = 'K/day '
860     udiag(72) = 'watt/m2 '
861     udiag(73) = 'watt/m2 '
862     udiag(74) = 'g/m2/sec '
863     udiag(75) = 'm/sec '
864     udiag(76) = 'm/sec '
865     udiag(77) = 'K '
866     udiag(78) = 'gram/kg '
867     udiag(79) = 'K/day '
868     udiag(80) = 'hPa '
869     udiag(81) = '(m/sec)2 '
870     udiag(82) = 'watt/m2 '
871     udiag(83) = 'hPa/day '
872     udiag(84) = 'user-defined '
873     udiag(85) = 'user-defined '
874     udiag(86) = 'user-defined '
875     udiag(87) = 'user-defined '
876     udiag(88) = 'm/s/day '
877     udiag(89) = 'm/s/day '
878     udiag(90) = 'K/day '
879     udiag(91) = 'gram/kg/day '
880     udiag(92) = 'm/s/day '
881     udiag(93) = 'm/s/day '
882     udiag(94) = 'm/s/day '
883     udiag(95) = 'm/s/day '
884     udiag(96) = 'N/m2 '
885     udiag(97) = 'N/m2 '
886     udiag(98) = 'N/m2 '
887     udiag(99) = 'N/m2 '
888     udiag(100) = 'mg/kg '
889     udiag(101) = 'hPa '
890     udiag(102) = 'm/s gram/kg '
891     udiag(103) = 'm/s gram/kg '
892     udiag(104) = 'm/s K '
893     udiag(105) = 'm/s K '
894     udiag(106) = '0-1 '
895     udiag(107) = 'gram/cm2 '
896     udiag(108) = 'm/sec '
897     udiag(109) = 'm/sec '
898     udiag(110) = 'K '
899     udiag(111) = 'gram/kg '
900     udiag(112) = 'm/sec '
901     udiag(113) = 'm/sec '
902     udiag(114) = 'K '
903     udiag(115) = 'gram/kg '
904     udiag(116) = 'g/m2/sec '
905     udiag(117) = 'gram/kg/day '
906     udiag(118) = 'mm/day '
907     udiag(119) = 'mm/day '
908     udiag(120) = 'm/s/day '
909     udiag(121) = 'm/s/day '
910     udiag(122) = 'K/day '
911     udiag(123) = 'gram/kg/day '
912     udiag(124) = 'user-defined '
913     udiag(125) = 'user-defined '
914     udiag(126) = 'user-defined '
915     udiag(127) = 'user-defined '
916     udiag(128) = 'user-defined '
917     udiag(129) = 'user-defined '
918     udiag(130) = 'user-defined '
919     udiag(131) = 'user-defined '
920     udiag(132) = 'user-defined '
921     udiag(133) = 'user-defined '
922     udiag(134) = 'user-defined '
923     udiag(135) = 'user-defined '
924     udiag(136) = 'user-defined '
925     udiag(137) = 'user-defined '
926     udiag(138) = 'user-defined '
927     udiag(139) = 'user-defined '
928     udiag(140) = '0-1 '
929     udiag(141) = '0-1 '
930     udiag(142) = 'non-dimensional '
931     udiag(143) = 'non-dimensional '
932     udiag(144) = 'non-dimensional '
933     udiag(145) = '0-1 '
934     udiag(146) = 'hPa '
935     udiag(147) = 'count '
936     udiag(148) = 'K '
937     udiag(149) = 'count '
938     udiag(150) = 'count '
939     udiag(151) = 'count '
940     udiag(152) = 'count '
941     udiag(153) = 'K '
942     udiag(154) = 'K '
943     udiag(155) = 'gram/kg '
944     udiag(156) = '0-1 '
945     udiag(157) = '0-1 '
946     udiag(158) = '0-1 '
947     udiag(159) = 'mm '
948     udiag(160) = 'mm water equiv '
949     udiag(161) = 'mm/day '
950     udiag(162) = 'mm/day '
951     udiag(163) = 'mm/day '
952     udiag(164) = 'mm/sec '
953     udiag(165) = 'mm/sec '
954     udiag(166) = 'mm/sec '
955     udiag(167) = 'mm/sec '
956     udiag(168) = 'watt/m2 '
957     udiag(169) = 'watt/m2 '
958     udiag(170) = 'watt/m2 '
959     udiag(171) = 'watt/m2 '
960     udiag(172) = 'watt/m2 '
961     udiag(173) = 'watt/m2 '
962     udiag(174) = 'percent '
963     udiag(175) = 'percent '
964     udiag(176) = 'watt/m2/K '
965     udiag(177) = 'watt/m2/K '
966     udiag(178) = 'watt/m2/K '
967     udiag(179) = 'K/day '
968     udiag(180) = 'K/day '
969     udiag(181) = 'K/day '
970     udiag(182) = 'K/day '
971     udiag(183) = 'K/day '
972     udiag(184) = 'hPa '
973     udiag(185) = 'count '
974     udiag(186) = 'hPa '
975     udiag(187) = 'user-defined '
976     udiag(188) = 'count '
977     udiag(189) = ' '
978     udiag(190) = 'watt/m2 '
979     udiag(191) = 'mm/day '
980     udiag(192) = 'watt/m2 '
981     udiag(193) = 'watt/m2 '
982     udiag(194) = 'K/sec '
983     udiag(195) = 'sec-1 '
984     udiag(196) = 'K/sec '
985     udiag(197) = 'K/sec '
986     udiag(198) = 'K/sec '
987     udiag(199) = 'K/sec '
988     udiag(200) = 'K/sec '
989     udiag(201) = 'K/sec '
990     udiag(202) = 'mg/kg '
991     udiag(203) = '0-1 '
992     udiag(204) = '0-1 '
993     udiag(205) = '0-1 '
994     udiag(206) = '0-1 '
995     udiag(207) = '0-1 '
996     udiag(208) = '0-1 '
997     udiag(209) = 'non-dimensional '
998     udiag(210) = 'count '
999     udiag(211) = 'K/day '
1000     udiag(212) = 'g/kg/day '
1001     udiag(213) = '0-1 '
1002     udiag(214) = 'mg/cm2 '
1003     udiag(215) = 'mg/cm2 '
1004     udiag(216) = 'mg/cm2 '
1005     udiag(217) = 'mg/cm2 '
1006     udiag(218) = 'mb m2/sec2 '
1007     udiag(219) = 'm2/sec2 '
1008     udiag(220) = 'm2/sec2 '
1009     udiag(221) = 'K/day '
1010     udiag(222) = 'g/kg/day '
1011     udiag(223) = 'K/day '
1012     udiag(224) = 'g/kg/day '
1013     udiag(225) = 'K/day '
1014     udiag(226) = 'g/kg/day '
1015     udiag(227) = 'K/day '
1016     udiag(228) = 'g/kg/day '
1017     udiag(229) = 'm/sec '
1018     udiag(230) = 'm/sec '
1019     udiag(231) = 'm '
1020     udiag(232) = 'm/sec '
1021     udiag(233) = 'm/sec '
1022     udiag(234) = 'N/m2 '
1023     udiag(235) = 'N/m2 '
1024     udiag(236) = 'm/sec/day '
1025     udiag(237) = 'm/sec/day '
1026 molod 1.7 udiag(238) = 'm or hPa '
1027     udiag(239) = 'm2 or hPa2 '
1028     udiag(240) = 'K '
1029     udiag(241) = 'g/kg '
1030     udiag(242) = 'm/sec '
1031     udiag(243) = 'm/sec '
1032     udiag(244) = 'm/sec '
1033     udiag(245) = 'K2 '
1034     udiag(246) = 'g2/kg2 '
1035     udiag(247) = 'm2/sec2 '
1036     udiag(248) = 'm2/sec2 '
1037     udiag(249) = 'm2/sec2 '
1038     udiag(250) = 'm2/sec2 '
1039 molod 1.8 udiag(251) = 'm/sec '
1040     udiag(252) = 'm/sec '
1041     udiag(253) = 'm/sec '
1042     udiag(254) = 'm-K/sec '
1043     udiag(255) = 'm-K/sec '
1044     udiag(256) = 'm-K/sec '
1045     udiag(257) = 'm-kg/s-kg '
1046     udiag(258) = 'm-kg/s-kg '
1047     udiag(259) = 'm-kg/s-kg '
1048 molod 1.9 udiag(260) = 'm-K/sec '
1049     udiag(261) = 'm-K/sec '
1050     udiag(262) = 'm-K/sec '
1051     udiag(263) = 'm-kg/s-kg '
1052     udiag(264) = 'm-kg/s-kg '
1053     udiag(265) = 'm-kg/s-kg '
1054 molod 1.10 udiag(266) = 'm-kg/s-kg '
1055     udiag(267) = 'm-kg/s-kg '
1056     udiag(268) = 'm-kg/s-kg '
1057     udiag(269) = 'm-kg/s-kg '
1058     udiag(270) = 'm-kg/s-kg '
1059     udiag(271) = 'm-kg/s-kg '
1060     udiag(272) = 'm-kg/s-kg '
1061     udiag(273) = 'm-kg/s-kg '
1062     udiag(274) = 'm-kg/s-kg '
1063     udiag(275) = 'm-kg/s-kg '
1064     udiag(276) = 'm-kg/s-kg '
1065     udiag(277) = 'm-kg/s-kg '
1066     udiag(278) = 'm-kg/s-kg '
1067     udiag(279) = 'm-kg/s-kg '
1068     udiag(280) = 'm-kg/s-kg '
1069 molod 1.12 udiag(281) = 'kg/kg '
1070     udiag(282) = 'kg/kg '
1071     udiag(283) = 'kg/kg '
1072     udiag(284) = 'kg/kg '
1073     udiag(285) = 'kg/kg '
1074 molod 1.13 udiag(286) = 'mol/m3/sec '
1075     udiag(287) = 'mol eq/m3/sec '
1076     udiag(288) = 'mol/m3/sec '
1077     udiag(289) = 'mol/m3/sec '
1078     udiag(290) = 'mol/m2/sec '
1079     udiag(291) = 'atm '
1080     udiag(292) = 'dimensionless '
1081 molod 1.14 udiag(293) = 'k/sec '
1082     udiag(294) = 'g/kg/sec '
1083     udiag(295) = 'Fraction '
1084     udiag(296) = 'g/m2/sec '
1085     udiag(297) = 'J/g '
1086 molod 1.1
1087     c Diagnostic Parms
1088 edhill 1.18 c ---------------- 123456789
1089     gdiag( 1) = 'UM 002L1 '
1090     gdiag( 2) = 'VM 001L1 '
1091     gdiag( 3) = 'SM L1 '
1092     gdiag( 4) = 'SM L1 '
1093     gdiag( 5) = 'SM L1 '
1094     gdiag( 6) = 'SM L1 '
1095     gdiag( 7) = 'SM P L1 '
1096     gdiag( 8) = 'SM LL '
1097     gdiag( 9) = 'SM L1 '
1098     gdiag(10) = 'SM L1 '
1099     gdiag(11) = 'SM LL '
1100     gdiag(12) = 'SM LL '
1101     gdiag(13) = 'UM 014ML '
1102     gdiag(14) = 'VM 013ML '
1103     gdiag(15) = 'SM ML '
1104     gdiag(16) = 'SM ML '
1105     gdiag(17) = 'SM ML '
1106     gdiag(18) = 'SM ML '
1107     gdiag(19) = 'SM ML '
1108     gdiag(20) = 'SM P ML '
1109     gdiag(21) = 'SM P L1 '
1110     gdiag(22) = 'SM P L1 '
1111     gdiag(23) = 'UM 024LL '
1112     gdiag(24) = 'VM 023LL '
1113     gdiag(25) = 'SM LL '
1114     gdiag(26) = 'SM LL '
1115     gdiag(27) = 'SM L1 '
1116     gdiag(28) = 'SM L1 '
1117     gdiag(29) = 'SM L1 '
1118     gdiag(30) = 'SM L1 '
1119     gdiag(31) = 'SM L1 '
1120     gdiag(32) = 'SM L1 '
1121     gdiag(33) = 'SM P L1 '
1122     gdiag(34) = 'SM P L1 '
1123     gdiag(35) = 'SM L1 '
1124     gdiag(36) = 'SM L1 '
1125     gdiag(37) = 'SM U1 '
1126     gdiag(38) = 'SM U1 '
1127     gdiag(39) = 'SM L1 '
1128     gdiag(40) = 'SM ML '
1129     gdiag(41) = 'SM ML '
1130     gdiag(42) = 'SM P ML '
1131     gdiag(43) = 'SM ML '
1132     gdiag(44) = 'SM P ML '
1133     gdiag(45) = 'SM P ML '
1134     gdiag(46) = 'SM L1 '
1135     gdiag(47) = 'SM ML '
1136     gdiag(48) = 'SM P U1 '
1137     gdiag(49) = 'SM C050ML '
1138     gdiag(50) = 'SM D ML '
1139     gdiag(51) = 'SM P M1 '
1140     gdiag(52) = 'SM L1 '
1141     gdiag(53) = 'SM L1 '
1142     gdiag(54) = 'UU 055MR '
1143     gdiag(55) = 'VV 054MR '
1144     gdiag(56) = 'SM MR '
1145     gdiag(57) = 'SM MR '
1146     gdiag(58) = 'SM ML '
1147     gdiag(59) = 'UU 060ML '
1148     gdiag(60) = 'VV 059ML '
1149     gdiag(61) = 'SM ML '
1150     gdiag(62) = 'SM ML '
1151     gdiag(63) = 'SM ML '
1152     gdiag(64) = 'SM LL '
1153     gdiag(65) = 'SM ML '
1154     gdiag(66) = 'SM ML '
1155     gdiag(67) = 'SM L1 '
1156     gdiag(68) = 'SM L1 '
1157     gdiag(69) = 'SM P MM '
1158     gdiag(70) = 'SM P M1 '
1159     gdiag(71) = 'SM P ML '
1160     gdiag(72) = 'SM P U1 '
1161     gdiag(73) = 'SM P U1 '
1162     gdiag(74) = 'SM P LL '
1163     gdiag(75) = 'UU 076ML '
1164     gdiag(76) = 'VV 075ML '
1165     gdiag(77) = 'SM ML '
1166     gdiag(78) = 'SM P ML '
1167     gdiag(79) = 'SM ML '
1168     gdiag(80) = 'SM L1 '
1169     gdiag(81) = 'SM P LL '
1170     gdiag(82) = 'SM P L1 '
1171     gdiag(83) = 'SM L1 '
1172     gdiag(84) = 'SM L1 '
1173     gdiag(85) = 'SM L1 '
1174     gdiag(86) = 'SM ML '
1175     gdiag(87) = 'SM ML '
1176     gdiag(88) = 'UU 089ML '
1177     gdiag(89) = 'VV 088ML '
1178     gdiag(90) = 'SM ML '
1179     gdiag(91) = 'SM ML '
1180     gdiag(92) = 'UU 093ML '
1181     gdiag(93) = 'VV 092ML '
1182     gdiag(94) = 'UM 095ML '
1183     gdiag(95) = 'VM 094ML '
1184     gdiag(96) = 'UM 097L1 '
1185     gdiag(97) = 'VM 096L1 '
1186     gdiag(98) = 'UM 099U1 '
1187     gdiag(99) = 'VM 098U1 '
1188     gdiag(100) = 'SM ML '
1189     gdiag(101) = 'SM L1 '
1190     gdiag(102) = 'UM 103M1 '
1191     gdiag(103) = 'VM 102M1 '
1192     gdiag(104) = 'UM 105M1 '
1193     gdiag(105) = 'VM 104M1 '
1194     gdiag(106) = 'SM P M1 '
1195     gdiag(107) = 'SM P M1 '
1196     gdiag(108) = 'UM 109M1 '
1197     gdiag(109) = 'VM 108M1 '
1198     gdiag(110) = 'SM M1 '
1199     gdiag(111) = 'SM P M1 '
1200     gdiag(112) = 'UM 113M1 '
1201     gdiag(113) = 'VM 112M1 '
1202     gdiag(114) = 'SM M1 '
1203     gdiag(115) = 'SM P M1 '
1204     gdiag(116) = 'SM P LL '
1205     gdiag(117) = 'SM ML '
1206     gdiag(118) = 'SM M1 '
1207     gdiag(119) = 'SM P L1 '
1208     gdiag(120) = 'UU 121ML '
1209     gdiag(121) = 'VV 120ML '
1210     gdiag(122) = 'SM ML '
1211     gdiag(123) = 'SM ML '
1212     gdiag(124) = 'SM L1 '
1213     gdiag(125) = 'SM L1 '
1214     gdiag(126) = 'SM L1 '
1215     gdiag(127) = 'SM L1 '
1216     gdiag(128) = 'SM L1 '
1217     gdiag(129) = 'SM L1 '
1218     gdiag(130) = 'SM L1 '
1219     gdiag(131) = 'SM L1 '
1220     gdiag(132) = 'SM ML '
1221     gdiag(133) = 'SM ML '
1222     gdiag(134) = 'SM ML '
1223     gdiag(135) = 'SM ML '
1224     gdiag(136) = 'SM ML '
1225     gdiag(137) = 'SM ML '
1226     gdiag(138) = 'SM ML '
1227     gdiag(139) = 'SM ML '
1228     gdiag(140) = 'SM P M1 '
1229     gdiag(141) = 'SM P M1 '
1230     gdiag(142) = 'SM PC150M1 '
1231     gdiag(143) = 'SM PC151M1 '
1232     gdiag(144) = 'SM PC152M1 '
1233     gdiag(145) = 'SM ML '
1234     gdiag(146) = 'SM PC147M1 '
1235     gdiag(147) = 'SM D M1 '
1236     gdiag(148) = 'SM PC149M1 '
1237     gdiag(149) = 'SM D M1 '
1238     gdiag(150) = 'SM D M1 '
1239     gdiag(151) = 'SM D M1 '
1240     gdiag(152) = 'SM D M1 '
1241     gdiag(153) = 'SM L1 '
1242     gdiag(154) = 'SM L1 '
1243     gdiag(155) = 'SM L1 '
1244     gdiag(156) = 'SM L1 '
1245     gdiag(157) = 'SM L1 '
1246     gdiag(158) = 'SM L1 '
1247     gdiag(159) = 'SM L1 '
1248     gdiag(160) = 'SM L1 '
1249     gdiag(161) = 'SM L1 '
1250     gdiag(162) = 'SM L1 '
1251     gdiag(163) = 'SM P L1 '
1252     gdiag(164) = 'SM L1 '
1253     gdiag(165) = 'SM L1 '
1254     gdiag(166) = 'SM L1 '
1255     gdiag(167) = 'SM L1 '
1256     gdiag(168) = 'SM L1 '
1257     gdiag(169) = 'SM L1 '
1258     gdiag(170) = 'SM L1 '
1259     gdiag(171) = 'SM L1 '
1260     gdiag(172) = 'SM P L1 '
1261     gdiag(173) = 'SM P L1 '
1262     gdiag(174) = 'SM L1 '
1263     gdiag(175) = 'SM L1 '
1264     gdiag(176) = 'SM L1 '
1265     gdiag(177) = 'SM L1 '
1266     gdiag(178) = 'SM L1 '
1267     gdiag(179) = 'SM M1 '
1268     gdiag(180) = 'SM M1 '
1269     gdiag(181) = 'SM M1 '
1270     gdiag(182) = 'SM P M1 '
1271     gdiag(183) = 'SM M1 '
1272     gdiag(184) = 'SM PC185M1 '
1273     gdiag(185) = 'SM PD M1 '
1274     gdiag(186) = 'SM L1 '
1275     gdiag(187) = 'SM C188L1 '
1276     gdiag(188) = 'SM D L1 '
1277     gdiag(189) = ' 0 '
1278     gdiag(190) = 'SM L1 '
1279     gdiag(191) = 'SM P L1 '
1280     gdiag(192) = 'SM L1 '
1281     gdiag(193) = 'SM L1 '
1282     gdiag(194) = 'SM L1 '
1283     gdiag(195) = 'SM L1 '
1284     gdiag(196) = 'SM L1 '
1285     gdiag(197) = 'SM L1 '
1286     gdiag(198) = 'SM L1 '
1287     gdiag(199) = 'SM L1 '
1288     gdiag(200) = 'SM L1 '
1289     gdiag(201) = 'SM L1 '
1290     gdiag(202) = 'SM ML '
1291     gdiag(203) = 'SM ML '
1292     gdiag(204) = 'SM C210L1 '
1293     gdiag(205) = 'SM L1 '
1294     gdiag(206) = 'SM L1 '
1295     gdiag(207) = 'SM L1 '
1296     gdiag(208) = 'SM L1 '
1297     gdiag(209) = 'SM ML '
1298     gdiag(210) = 'SM D L1 '
1299     gdiag(211) = 'SM ML '
1300     gdiag(212) = 'SM ML '
1301     gdiag(213) = 'SM P ML '
1302     gdiag(214) = 'SM P L1 '
1303     gdiag(215) = 'SM P L1 '
1304     gdiag(216) = 'SM P L1 '
1305     gdiag(217) = 'SM P L1 '
1306     gdiag(218) = 'SM P L1 '
1307     gdiag(219) = 'SM P L1 '
1308     gdiag(220) = 'SM P ML '
1309     gdiag(221) = 'SM ML '
1310     gdiag(222) = 'SM ML '
1311     gdiag(223) = 'SM ML '
1312     gdiag(224) = 'SM ML '
1313     gdiag(225) = 'SM ML '
1314     gdiag(226) = 'SM ML '
1315     gdiag(227) = 'SM ML '
1316     gdiag(228) = 'SM ML '
1317     gdiag(229) = 'UM 230L1 '
1318     gdiag(230) = 'VM 229L1 '
1319     gdiag(231) = 'SM L1 '
1320     gdiag(232) = 'UM 233L1 '
1321     gdiag(233) = 'VM 232L1 '
1322     gdiag(234) = 'UM 235L1 '
1323     gdiag(235) = 'VM 234L1 '
1324     gdiag(236) = 'UM 237ML '
1325     gdiag(237) = 'VM 236ML '
1326     gdiag(238) = 'SM M1 '
1327     gdiag(239) = 'SM M1 '
1328     gdiag(240) = 'SM MR '
1329     gdiag(241) = 'SM MR '
1330     gdiag(242) = 'UU 243MR '
1331     gdiag(243) = 'VV 242MR '
1332     gdiag(244) = 'WM LR '
1333     gdiag(245) = 'SM MR '
1334     gdiag(246) = 'SM MR '
1335     gdiag(247) = 'UU 248MR '
1336     gdiag(248) = 'VV 247MR '
1337     gdiag(249) = 'WM LR '
1338     gdiag(250) = 'UZ 250MR '
1339     gdiag(251) = 'UU 252MR '
1340     gdiag(252) = 'VV 251MR '
1341     gdiag(253) = 'WM MR '
1342     gdiag(254) = 'UU 255MR '
1343     gdiag(255) = 'VV 254MR '
1344     gdiag(256) = 'WM MR '
1345     gdiag(257) = 'UU 258MR '
1346     gdiag(258) = 'VV 257MR '
1347     gdiag(259) = 'WM MR '
1348     gdiag(260) = 'UU 261MR '
1349     gdiag(261) = 'VV 260MR '
1350     gdiag(262) = 'WM MR '
1351     gdiag(263) = 'UU 264MR '
1352     gdiag(264) = 'VV 263MR '
1353     gdiag(265) = 'WM MR '
1354     gdiag(266) = 'UU 267MR '
1355     gdiag(267) = 'VV 266MR '
1356     gdiag(268) = 'WM MR '
1357     gdiag(269) = 'UU 270MR '
1358     gdiag(270) = 'VV 269MR '
1359     gdiag(271) = 'WM MR '
1360     gdiag(272) = 'UU 273MR '
1361     gdiag(273) = 'VV 272MR '
1362     gdiag(274) = 'WM MR '
1363     gdiag(275) = 'UU 276MR '
1364     gdiag(276) = 'VV 275MR '
1365     gdiag(277) = 'WM MR '
1366     gdiag(278) = 'UU 279MR '
1367     gdiag(279) = 'VV 278MR '
1368     gdiag(280) = 'WM MR '
1369     gdiag(281) = 'SM P MR '
1370     gdiag(282) = 'SM P MR '
1371     gdiag(283) = 'SM P MR '
1372     gdiag(284) = 'SM P MR '
1373     gdiag(285) = 'SM P MR '
1374     gdiag(286) = 'SM P MR '
1375     gdiag(287) = 'SM P MR '
1376     gdiag(288) = 'SM P L1 '
1377     gdiag(289) = 'SM P L1 '
1378     gdiag(290) = 'SM P L1 '
1379     gdiag(291) = 'SM P M1 '
1380     gdiag(292) = 'SM P M1 '
1381     gdiag(293) = 'SM MR '
1382     gdiag(294) = 'SM MR '
1383     gdiag(295) = 'SM MR '
1384     gdiag(296) = 'SM M1 '
1385     gdiag(297) = 'SM M1 '
1386    
1387     c Diagnostic Levels
1388     c -----------------
1389 edhill 1.20 do n = 1,ndiagt
1390 edhill 1.18 IF (gdiag(n)(10:10) .eq. '0') KDIAG(n) = 0
1391     IF (gdiag(n)(10:10) .eq. '1') KDIAG(n) = 1
1392     IF (gdiag(n)(10:10) .eq. 'R') KDIAG(n) = Nr
1393     IF (gdiag(n)(10:10) .eq. 'L') KDIAG(n) = nlevs
1394     IF (gdiag(n)(10:10) .eq. 'M') KDIAG(n) = nlevs - 1
1395     enddo
1396 molod 1.1
1397     do n=1,ndiagt
1398 edhill 1.18 idiag(n) = 0
1399     ndiag(n) = 0
1400 molod 1.1 enddo
1401    
1402     RETURN
1403     END
1404 edhill 1.19
1405    
1406     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
  ViewVC Help
Powered by ViewVC 1.1.22