/[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.9 - (hide annotations) (download)
Thu Mar 11 23:56:45 2004 UTC (20 years, 3 months ago) by molod
Branch: MAIN
Changes since 1.8: +30 -0 lines
Finish diagnostics that are in timeave

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

  ViewVC Help
Powered by ViewVC 1.1.22