/[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.12 - (hide annotations) (download)
Wed Mar 24 23:04:03 2004 UTC (20 years, 3 months ago) by molod
Branch: MAIN
Changes since 1.11: +25 -0 lines
More code to finish tracer diagnostics

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.10 cdiag(266) = 'UTRAC1 '
331     cdiag(267) = 'VTRAC1 '
332     cdiag(268) = 'WTRAC1 '
333     cdiag(269) = 'UTRAC2 '
334     cdiag(270) = 'VTRAC2 '
335     cdiag(271) = 'WTRAC2 '
336     cdiag(272) = 'UTRAC3 '
337     cdiag(273) = 'VTRAC3 '
338     cdiag(274) = 'WTRAC3 '
339     cdiag(275) = 'UTRAC4 '
340     cdiag(276) = 'VTRAC4 '
341     cdiag(277) = 'WTRAC4 '
342     cdiag(278) = 'UTRAC5 '
343     cdiag(279) = 'VTRAC5 '
344     cdiag(280) = 'WTRAC5 '
345 molod 1.12 cdiag(281) = 'TRAC1 '
346     cdiag(282) = 'TRAC2 '
347     cdiag(283) = 'TRAC3 '
348     cdiag(284) = 'TRAC4 '
349     cdiag(285) = 'TRAC5 '
350 molod 1.1
351     c Diagnostic Descriptions
352 molod 1.2 c ----------------------------------------------------------------------
353 molod 1.1 tdiag( 1) = 'Zonal Wind Surface Stress (N/m^2)'
354     tdiag( 2) = 'Meridional Wind Surface Stress (N/m^2)'
355     tdiag( 3) = 'Sensible Heat Flux (pos.upwrd) (W/m^2)'
356     tdiag( 4) = 'Latent Heat Flux (pos.upwrd) (W/m^2)'
357     tdiag( 5) = 'Heat Conduction through Sea-Ice (W/m^2)'
358     tdiag( 6) = 'Net Upward Longwave Flux at the Ground (W/m^2)'
359     tdiag( 7) = 'Net Downward Shortwave Flux at the Ground (W/m^2)'
360     tdiag( 8) = 'Richardson Number (non-dimensional)'
361     tdiag( 9) = 'Surface Drag Coef for Scalars (non-dimensional)'
362     tdiag(10) = 'Surface Drag Coef for Winds (non-dimensional)'
363     tdiag(11) = 'Eddy Diffusivity Coefficient for Momentum (m^2/s)'
364     tdiag(12) = 'Eddy Diffusivity Coefficient for Scalars (m^2/s)'
365 molod 1.2 tdiag(13) = 'Zonal Wind Tendency from Turbulence (m/s/day)'
366     tdiag(14) = 'Meridional Wind Tendency from Turbulence (m/s/day)'
367     tdiag(15) =
368     . 'Temperature Tendency from Turbulence (K/day)'
369     tdiag(16) =
370     . 'Specific Humidity Tendency from Turbulence (g/kg/day)'
371     tdiag(17) =
372     . 'Temperature Tendency due to Moist Processes (K/day)'
373     tdiag(18) =
374     . 'Specific Humidity Tendency due to Moist Processes (g/kg/day)'
375     tdiag(19) =
376     . 'Temperature Tendency due to Longwave Radiation (K/day)'
377     tdiag(20) =
378     . 'Temperature Tendency due to Shortwave Radiation (K/day)'
379 molod 1.1 tdiag(21) = 'Total Precipitation (mm/day)'
380     tdiag(22) = 'Convective Precipitation (mm/day)'
381     tdiag(23) = 'Turbulent Flux of Zonal Momentum (N/m^2)'
382     tdiag(24) = 'Turbulent Flux of Meridional Momentum (N/m^2)'
383     tdiag(25) = 'Turbulent Flux of Sensible Heat (W/m^2)'
384     tdiag(26) = 'Turbulent Flux of Latent Heat (W/m^2)'
385     tdiag(27) = 'Neutral Drag Coefficient (non-dimensional)'
386     tdiag(28) = 'Surface Wind Speed (m/s)'
387     tdiag(29) = 'Air/Surface Virtual Temperature Difference (deg K)'
388     tdiag(30) = 'Ground Temperature (deg K)'
389     tdiag(31) = 'Surface Air Temperature (deg K)'
390     tdiag(32) = 'Ground Temperature Adjustment (deg K)'
391     tdiag(33) = 'Ground Specific Humidity (g/kg)'
392     tdiag(34) = 'Saturation Surface Specific Humidity (g/kg)'
393     tdiag(35) = 'Ground Temperature used by LW Radiation (deg K)'
394     tdiag(36) = 'Upward Longwave Radiation at the Ground (W/m^2)'
395     tdiag(37) = 'Outgoing Longwave Radiation (W/m^2)'
396     tdiag(38) = 'Clear Sky Outgoing Longwave Radiation (W/m^2)'
397     tdiag(39) = 'Clear Sky Net Longwave Flux at the Ground (W/m^2)'
398     tdiag(40) = 'Clear Sky Longwave Radiation Heating Rates (K/day)'
399     tdiag(41) = 'Temperature used by Radiation (deg K)'
400     tdiag(42) = 'Specific Humidity used by Radiation (g/kg)'
401 molod 1.2 tdiag(43) =
402     . 'Ozone mixing ratio used by Longwave Radiation (non-dimensional)'
403 molod 1.1 tdiag(44) = 'Convective Cloud Fraction (0-1)'
404     tdiag(45) = '3-D Total Cloud Fraction (0-1)'
405     tdiag(46) = 'Downward Longwave Radiation at the Ground (W/m^2)'
406 molod 1.2 tdiag(47) =
407     . 'Temperature Tendency due to Gravity Wave Drag (K/day)'
408 molod 1.1 tdiag(48) = 'Incident Shortwave Radiation at TOA (W/m^2)'
409 molod 1.2 tdiag(49) =
410     . 'Counted Cloud Optical Depth (non-dimensional) per 100 mb'
411 molod 1.1 tdiag(50) = 'Cloud Optical Depth Counter'
412     tdiag(51) = 'Low-Level ( 1000-700 hPa) Cloud Fraction (0-1)'
413     tdiag(52) = 'Surface Evaporation (mm/day)'
414     tdiag(53) = 'Surface Pressure Tendency (hPa/day)'
415 molod 1.5 tdiag(54) = 'Average Zonal Wind (m/s)'
416     tdiag(55) = 'Average Meridional Wind (m/s)'
417     tdiag(56) = 'Average Potential Temperature (K)'
418 molod 1.2 tdiag(57) =
419 molod 1.5 . 'Average Specific Humidity (g/kg)'
420 molod 1.1 tdiag(58) = 'Vertical Velocity (hPa/day)'
421     tdiag(59) = 'Total Zonal Wind Tendency (m/s/day)'
422     tdiag(60) = 'Total Meridional Wind Tendency (m/s/day)'
423     tdiag(61) = 'Total Temperature Tendency (K/day)'
424     tdiag(62) = 'Total Specific Humidity Tendency (g/kg/day)'
425     tdiag(63) = 'Relative Vorticity (x10^-4 sec^-1)'
426     tdiag(64) = 'Pi-Weighted Vertical Sigma Velocity (hPa/day)'
427 molod 1.2 tdiag(65) =
428     . 'Temperature Tendency due to Stratiform Processes (K/day)'
429     tdiag(66) =
430     . 'Specific Humidity Tendency from Stratiform Processes (g/kg/day)'
431 molod 1.1 tdiag(67) = 'Friction Velocity (m/s)'
432     tdiag(68) = 'Surface Roughness Z0 (m)'
433     tdiag(69) = 'Frequency of Turbulence (0-1)'
434     tdiag(70) = 'Estimated PBL Depth (hPa)'
435     tdiag(71) = 'Clear Sky Shortwave Heating Rates (K/day)'
436     tdiag(72) = 'Outgoing Shortwave Radiation (W/m^2)'
437     tdiag(73) = 'Clear Sky Outgoing Shortwave Radiation (W/m^2)'
438     tdiag(74) = 'Cloud Mass Flux (g/m^2/sec)'
439     tdiag(75) = 'Time-Averaged Zonal Wind (m/s)'
440     tdiag(76) = 'Time-Averaged Meridional Wind (m/s)'
441     tdiag(77) = 'Time-Averaged Temperature (deg K)'
442     tdiag(78) = 'Time-Averaged Specific Humidity (g/kg)'
443 molod 1.2 tdiag(79) =
444     . 'Temperature Tendency due to Rayleigh Friction (K/day)'
445 molod 1.1 tdiag(80) = 'Time-Averaged Surface Pressure (hPa)'
446     tdiag(81) = 'Time-Averaged Turbulent Kinetic Energy (m/s)^2'
447 molod 1.2 tdiag(82) =
448     . 'Clear Sky Net Downward SW Radiation at the Ground (W/m^2)'
449 molod 1.5 tdiag(83) = 'Average Surface Pressure (hPa)'
450 molod 1.1 tdiag(84) = 'User-Defined Surface Diagnostic #1'
451     tdiag(85) = 'User-Defined Surface Diagnostic #2'
452     tdiag(86) = 'User-Defined Upper-Air Diagnostic #1'
453     tdiag(87) = 'User-Defined Upper-Air Diagnostic #2'
454 molod 1.2 tdiag(88) =
455     . 'Zonal Wind Tendency due to Diabatic Forcing (m/s/day)'
456     tdiag(89) =
457     . 'Meridional Wind Tendency due to Diabatic Forcing (m/s/day)'
458     tdiag(90) =
459     . 'Temperature Tendency due to Diabatic Forcing (K/day)'
460     tdiag(91) =
461     . 'Specific Humidity Tendency due to Diabatic Forcing (g/kg/day)'
462     tdiag(92) =
463     . 'Zonal Wind Tendency due to Rayleigh Friction (m/s/day)'
464     tdiag(93) =
465     . 'Meridional Wind Tendency due to Rayleigh Friction (m/s/day)'
466     tdiag(94) =
467     . 'Zonal Wind Tendency due to Gravity Wave Drag (m/s/day)'
468     tdiag(95) =
469     . 'Meridional Wind Tendency due to Gravity Wave Drag (m/s/day)'
470 molod 1.1 tdiag(96) = 'Zonal Wind Gravity Wave Surface Stress (N/m^2)'
471     tdiag(97) = 'Meridional Wind Gravity Wave Surface Stress (N/m^2)'
472     tdiag(98) = 'Zonal Wind Gravity Wave Stress at Model Top (N/m^2)'
473 molod 1.2 tdiag(99) =
474     . 'Meridional Wind Gravity Wave Stress at Model Top (N/m^2)'
475     tdiag(100) =
476     . 'Estimated Cloud Liquid Water used in Radiation (mg/kg)'
477 molod 1.1 tdiag(101) = 'Time-Averaged Sea-Level Pressure (mb)'
478     tdiag(102) = 'Vertically Averaged UWND*SPHU (m/s g/kg)'
479     tdiag(103) = 'Vertically Averaged VWND*SPHU (m/s g/kg)'
480     tdiag(104) = 'Vertically Averaged UWND*TMPU (m/s deg K)'
481     tdiag(105) = 'Vertically Averaged VWND*TMPU (m/s deg K)'
482     tdiag(106) = '2-D Total Cloud Fraction (0-1)'
483     tdiag(107) = 'Total Precipitable Water (g/cm^2)'
484     tdiag(108) = 'Zonal Wind Interpolated to 2 Meters (m/s)'
485     tdiag(109) = 'Meridional Wind Interpolated to 2 Meters (m/s)'
486     tdiag(110) = 'Temperature Interpolated to 2 Meters (deg K)'
487     tdiag(111) = 'Specific Humidity Interpolated to 2 Meters (g/kg)'
488     tdiag(112) = 'Zonal Wind Interpolated to 10 Meters (m/s)'
489     tdiag(113) = 'Meridional Wind Interpolated to 10 Meters (m/s)'
490     tdiag(114) = 'Temperature Interpolated to 10 Meters (deg K)'
491     tdiag(115) = 'Specific Humidity Interpolated to 10 Meters (g/kg)'
492     tdiag(116) = 'Detrainment Cloud Mass Flux (g/m^2/sec)'
493     tdiag(117) = 'Filling of Negative Specific Humidity (g/kg/day)'
494 molod 1.5 tdiag(118) = 'Vertically Averaged QAVE (mm/day)'
495 molod 1.1 tdiag(119) = 'Vertically Averaged QFILL (mm/day)'
496 molod 1.2 tdiag(120) =
497     . 'Zonal Wind Tendency due to Shapiro Filter (m/s/day)'
498     tdiag(121) =
499     . 'Meridional Wind Tendency due to Shapiro Filter (m/s/day)'
500 molod 1.1 tdiag(122) = 'Temperature Tendency due to Shapiro Filter (K/day)'
501 molod 1.2 tdiag(123) =
502     . 'Specific Humidity Tendency due to Shapiro Filter (g/kg/day)'
503 molod 1.1 tdiag(124) = 'User-Defined Surface Diagnostic #3 '
504     tdiag(125) = 'User-Defined Surface Diagnostic #4 '
505     tdiag(126) = 'User-Defined Surface Diagnostic #5 '
506     tdiag(127) = 'User-Defined Surface Diagnostic #6 '
507     tdiag(128) = 'User-Defined Surface Diagnostic #7 '
508     tdiag(129) = 'User-Defined Surface Diagnostic #8 '
509     tdiag(130) = 'User-Defined Surface Diagnostic #9 '
510     tdiag(131) = 'User-Defined Surface Diagnostic #10 '
511     tdiag(132) = 'User-Defined Upper-Air Diagnostic #3 '
512     tdiag(133) = 'User-Defined Upper-Air Diagnostic #4 '
513     tdiag(134) = 'User-Defined Upper-Air Diagnostic #5 '
514     tdiag(135) = 'User-Defined Upper-Air Diagnostic #6 '
515     tdiag(136) = 'User-Defined Upper-Air Diagnostic #7 '
516     tdiag(137) = 'User-Defined Upper-Air Diagnostic #8 '
517     tdiag(138) = 'User-Defined Upper-Air Diagnostic #9 '
518     tdiag(139) = 'User-Defined Upper-Air Diagnostic #10 '
519     tdiag(140) = 'Mid-Level (700-400 hPa) Cloud Fraction (0-1)'
520     tdiag(141) = 'High-Level (above 400 hPa) Cloud Fraction (0-1)'
521     tdiag(142) = 'Low-Level (1000-700 hPa) Optical Depth'
522     tdiag(143) = 'Mid-Level (700-400 hPa) Optical Depth'
523     tdiag(144) = 'High-Level (above 400 hPa) Optical Depth'
524     tdiag(145) = 'Non-Precipitating Cloud Fraction (0-1)'
525     tdiag(146) = 'Cloud Top Pressure (when cloudy) (hPa)'
526     tdiag(147) = 'Cloud Top Pressure Counter'
527     tdiag(148) = 'Cloud Top Temperature (when cloudy) (deg K)'
528     tdiag(149) = 'Cloud Top Temperature Counter'
529     tdiag(150) = 'Low-Level Optical Depth Counter'
530     tdiag(151) = 'Mid-Level Optical Depth Counter'
531     tdiag(152) = 'High-Level Optical Depth Counter'
532 molod 1.2 tdiag(153) =
533     . 'Skin Temperature of the Canopy/Soil (SST over water) (deg K)'
534 molod 1.1 tdiag(154) = 'Deep Soil Temperature (deg K)'
535     tdiag(155) = 'Specific Humidity of the Canopy (g/kg)'
536 molod 1.2 tdiag(156) =
537     . 'Soil Moisture Field Capacity Fraction in Shallow (5 mm) Layer'
538     tdiag(157) = 'Soil Moisture Field Capacity Fraction in Root Zone'
539     tdiag(158) =
540     . 'Soil Moisture Field Capacity Fraction in Recharge Layer '
541 molod 1.1 tdiag(159) = 'Canopy Interception Reservoir (mm)'
542     tdiag(160) = 'Snow Depth (mm water equivalent)'
543     tdiag(161) = 'Convective Rainfall (liquid precipitate) (mm/day)'
544     tdiag(162) = 'Large-scale Rainfall (liquid precipitate) (mm/day)'
545     tdiag(163) = 'Total Snowfall (solid precipitate) (mm/day)'
546 molod 1.2 tdiag(164) =
547     . 'Water from precipitation not infiltrated into soil (mm/s)'
548 molod 1.1 tdiag(165) = 'Rainwater Infiltration into top soil layer (mm/s)'
549 molod 1.2 tdiag(166) =
550     . 'Diffusion of Moisture across bottom of root zone (mm/s)'
551 molod 1.1 tdiag(167) = 'Rate of Snow Melt (mm/s)'
552 molod 1.2 tdiag(168) =
553     . 'Latent Heat Flux from Canopy Interception Reservoir (W/m^2)'
554 molod 1.1 tdiag(169) = 'Latent Heat Flux from Bare Soil (W/m^2)'
555 molod 1.2 tdiag(170) =
556     . 'Latent Heat Flux (Transpiration) from Vegetated Surface (W/m^2)'
557 molod 1.1 tdiag(171) = 'Latent Heat Flux from Snow Pack (W/m^2)'
558 molod 1.2 tdiag(172) =
559     . 'Diffuse-beam Photosynthetically Active Radiation (W/m^2)'
560     tdiag(173) =
561     . 'Direct-beam Photosynthetically Active Radiation (W/m^2)'
562 molod 1.1 tdiag(174) = 'Leaf Area Index (%)'
563     tdiag(175) = 'Green-ness Index (%)'
564 molod 1.2 tdiag(176) =
565     . 'Derivative of Upward LW Radiation wrt TCANOPY (W/m^2 /deg K)'
566     tdiag(177) =
567     . 'Derivative of Sensible Heat Flux wrt TCANOPY (W/m^2 /deg K)'
568     tdiag(178) =
569     . 'Derivative of Latent Heat Flux wrt TCANOPY (W/m^2 /deg K)'
570 molod 1.1 tdiag(179) = 'Vertically Averaged MOISTT (K/day)'
571     tdiag(180) = 'Vertically Averaged TURBT (K/day)'
572     tdiag(181) = 'Vertically Averaged RADLW (K/day)'
573     tdiag(182) = 'Vertically Averaged RADSW (K/day)'
574 molod 1.5 tdiag(183) = 'Vertically Averaged TAVE (K/day)'
575 molod 1.1 tdiag(184) = 'Pressure Thickness of the RAS Subcloud Layer (mb)'
576 molod 1.2 tdiag(185) =
577     . 'Pressure Thickness of the RAS Subcloud Layer Counter'
578     tdiag(186) = 'Lifting Condensation Level (mb) '
579     tdiag(187) = 'User-Defined Counted Surface Diagnostic '
580     tdiag(188) = 'User-Defined Counted Surface Diagnostic Counter'
581     tdiag(189) = ' '
582 molod 1.1 tdiag(190) = 'Potential Evapotranspiration (W/m^2)'
583 molod 1.2 tdiag(191) =
584     . 'Net precip due to supersaturation above critical RH (mm/day)'
585     tdiag(192) =
586     . 'Derivative of Sensible Heat wrt Canopy Humidity (W/m^2)'
587     tdiag(193) =
588     . 'Derivative of Latent Heat wrt Canopy Humidity (W/m^2)'
589 molod 1.1 tdiag(194) = 'Total Change in Canopy Temperature (deg/sec)'
590     tdiag(195) = 'Change in Canopy Humidity (kg/kg/sec)'
591 molod 1.2 tdiag(196) =
592     . 'Change in Canopy Temperature - Energy Term (deg/sec)'
593     tdiag(197) =
594     . 'Change in Canopy Temperature - Net Radiation Term (deg/sec)'
595     tdiag(198) =
596     . 'Change in Canopy Temperature - Sensible Heat Term (deg/sec)'
597     tdiag(199) =
598     . 'Change in Canopy Temperature - Latent Heat Term (deg/sec)'
599     tdiag(200) =
600     . 'Change in Canopy Temperature - Ground Temp Term (deg/sec)'
601     tdiag(201) =
602     . 'Change in Canopy Temperature - Humidity Term (deg/sec)'
603 molod 1.1 tdiag(202) = 'Cloud Liquid Water Generated in Turbulence (mg/kg)'
604     tdiag(203) = 'Cloud Fraction Generated in Turbulence (fraction)'
605     tdiag(204) = 'Mean Surface Albedo (SW_Up/SW_Down) 0-1 '
606     tdiag(205) = 'Direct Beam VIS Surface Albedo (0-1) '
607     tdiag(206) = 'Diffuse Beam VIS Surface Albedo (0-1) '
608     tdiag(207) = 'Direct Beam NIR Surface Albedo (0-1) '
609     tdiag(208) = 'Diffuse Beam NIR Surface Albedo (0-1) '
610 molod 1.2 tdiag(209) =
611     . 'Averaged Cloud Optical Depth (non-dimensional) per 100 mb'
612 molod 1.1 tdiag(210) = 'Surface Albedo Counter'
613 molod 1.2 tdiag(211) =
614     . 'Temperature Change due to humidity above critical RH (K/day)'
615     tdiag(212) =
616     . 'Spec.Humid. Change due to humidity above critical RH (g/kg/day)'
617 molod 1.1 tdiag(213) = 'Large-Scale Cloud Fraction (0-1)'
618 molod 1.2 tdiag(214) =
619     . 'Vert integrated low-level (below 700mb) cloud water (mg/cm**2)'
620     tdiag(215) =
621     . 'Vert integrated mid-level (700-400mb) cloud water (mg/cm**2)'
622     tdiag(216) =
623     . 'Vert integrated high-level (above 400mb) cloud water (mg/cm**2)'
624     tdiag(217) =
625     . 'Vertically integrated total suspended cloud water (mg/cm**2)'
626 molod 1.1 tdiag(218) = 'Convective Available Potential Energy'
627     tdiag(219) = 'Moist Static Energy in RAS subcloud layer'
628     tdiag(220) = 'Saturation Moist Static Energy'
629 molod 1.2 tdiag(221) =
630     . 'Temperature Tendency due to Convective Precipitation (K/day)'
631     tdiag(222) =
632     . 'Moisture Tendency due to Convective Precipitation (g/kg/day)'
633     tdiag(223) =
634     . 'Temperature Tendency due to Large Scale Precipitation (K/day)'
635     tdiag(224) =
636     . 'Moisture Tendency due to Large Scale Precipitation (g/kg/day)'
637     tdiag(225) =
638     . 'Temperature Tendency due to Conv Rain Re-evaporation (K/day)'
639     tdiag(226) =
640     . 'Moisture Tendency due to Conv Rain Re-evaporation (g/kg/day)'
641     tdiag(227) =
642     . 'Temperature Tendency due to Large Scale Rain Re-evap (K/day)'
643     tdiag(228) =
644     . 'Moisture Tendency due to Large Scale Rain Re-evap (g/kg/day)'
645     tdiag(229) =
646     . 'Mid-Lev Mean U-Wind for Gravity Wave Drag w/Clouds (m/sec)'
647     tdiag(230) =
648     . 'Mid-Lev Mean V-Wind for Gravity Wave Drag w/Clouds (m/sec)'
649     tdiag(231) =
650     . 'Mid-Lev Cloud Perturb Length for Gravity Wave Drag w/Clouds (m)'
651     tdiag(232) =
652     . 'Mid-Lev Perturb U-Wind for Gravity Wave Drag w/Clouds (m/sec)'
653     tdiag(233) =
654     . 'Mid-Lev Perturb V-Wind for Gravity Wave Drag w/Clouds (m/sec)'
655     tdiag(234) =
656     . 'Mid-Lev U-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)'
657     tdiag(235) =
658     . 'Mid-Lev V-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)'
659     tdiag(236) =
660     . 'U-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
661     tdiag(237) =
662     . 'V-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
663 molod 1.7 tdiag(238) = 'Perturbation of Surface (pressure, height) (hPa,m)'
664     tdiag(239) = 'Square of Perturbation of Sfc (hPa^2,m^2)'
665     tdiag(240) = 'Potential Temperature (K)'
666     tdiag(241) = 'Salt or Water Vapor mixing ratio (g/kg)'
667     tdiag(242) = 'Zonal Component of Velocity (m/s)'
668     tdiag(243) = 'Meridional Component of Velocity (m/s)'
669     tdiag(244) = 'Vertical Component of Velocity (m/s)'
670     tdiag(245) = 'Square of Potential Temperature (K^2)'
671     tdiag(246) = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
672     tdiag(247) = 'Square of Zonal Comp of Velocity (m^2/s^2)'
673     tdiag(248) = 'Square of Meridional Comp of Velocity (m^2/s^2)'
674     tdiag(249) = 'Square of Vertical Comp of Velocity (m^2/s^2)'
675     tdiag(250) = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
676 molod 1.8 tdiag(251) = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
677     tdiag(252) = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
678     tdiag(253) = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
679     tdiag(254) = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'
680     tdiag(255) = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'
681     tdiag(256) = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'
682     tdiag(257) = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'
683     tdiag(258) = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'
684     tdiag(259) = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'
685 molod 1.9 tdiag(260) = 'Zonal Transp of Pot Temp (m-K/s)'
686     tdiag(261) = 'Meridional Transp of Pot Temp (m-K/s)'
687     tdiag(262) = 'Vertical Transp of Pot Temp (m-K/s)'
688     tdiag(263) = 'Zonal Transp of Salt (m-kg/s-kg)'
689     tdiag(264) = 'Meridional Transp of Salt (m-kg/s-kg)'
690     tdiag(265) = 'Vertical Transp of Salt (m-kg/s-kg)'
691 molod 1.10 tdiag(266) = 'Zonal Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
692     tdiag(267) = 'Merid Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
693     tdiag(268) = 'Vert Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
694     tdiag(269) = 'Zonal Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
695     tdiag(270) = 'Merid Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
696     tdiag(271) = 'Vert Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
697     tdiag(272) = 'Zonal Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
698     tdiag(273) = 'Merid Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
699     tdiag(274) = 'Vert Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
700     tdiag(275) = 'Zonal Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
701     tdiag(276) = 'Merid Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
702     tdiag(277) = 'Vert Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
703     tdiag(278) = 'Zonal Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
704     tdiag(279) = 'Merid Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
705     tdiag(280) = 'Vert Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
706 molod 1.12 tdiag(281) = 'Mass-Weighted Tracer 1 (kgkg)'
707     tdiag(282) = 'Mass-Weighted Tracer 2 (kgkg)'
708     tdiag(283) = 'Mass-Weighted Tracer 3 (kgkg)'
709     tdiag(284) = 'Mass-Weighted Tracer 4 (kgkg)'
710     tdiag(285) = 'Mass-Weighted Tracer 5 (kgkg)'
711 molod 1.1
712     c Diagnostic Units
713     c ----------------
714     udiag( 1) = 'N/m2 '
715     udiag( 2) = 'N/m2 '
716     udiag( 3) = 'watt/m2 '
717     udiag( 4) = 'watt/m2 '
718     udiag( 5) = 'watt/m2 '
719     udiag( 6) = 'watt/m2 '
720     udiag( 7) = 'watt/m2 '
721     udiag( 8) = 'non-dimensional '
722     udiag( 9) = 'non-dimensional '
723     udiag(10) = 'non-dimensional '
724     udiag(11) = 'm2/sec '
725     udiag(12) = 'm2/sec '
726     udiag(13) = 'm/s/day '
727     udiag(14) = 'm/s/day '
728     udiag(15) = 'K/day '
729     udiag(16) = 'gram/kg/day '
730     udiag(17) = 'K/day '
731     udiag(18) = 'gram/kg/day '
732     udiag(19) = 'K/day '
733     udiag(20) = 'K/day '
734     udiag(21) = 'mm/day '
735     udiag(22) = 'mm/day '
736     udiag(23) = 'N/m2 '
737     udiag(24) = 'N/m2 '
738     udiag(25) = 'watt/m2 '
739     udiag(26) = 'watt/m2 '
740     udiag(27) = 'non-dimensional '
741     udiag(28) = 'm/sec '
742     udiag(29) = 'K '
743     udiag(30) = 'K '
744     udiag(31) = 'K '
745     udiag(32) = 'K '
746     udiag(33) = 'gram/kg '
747     udiag(34) = 'gram/kg '
748     udiag(35) = 'K '
749     udiag(36) = 'watt/m2 '
750     udiag(37) = 'watt/m2 '
751     udiag(38) = 'watt/m2 '
752     udiag(39) = 'watt/m2 '
753     udiag(40) = 'K/day '
754     udiag(41) = 'K '
755     udiag(42) = 'gram/kg '
756     udiag(43) = 'non-dimensional '
757     udiag(44) = '0-1 '
758     udiag(45) = '0-1 '
759     udiag(46) = 'watt/m2 '
760     udiag(47) = 'K/day '
761     udiag(48) = 'watt/m2 '
762     udiag(49) = 'non-dimensional '
763     udiag(50) = 'count '
764     udiag(51) = '0-1 '
765     udiag(52) = 'mm/day '
766     udiag(53) = 'hPa/day '
767 molod 1.5 udiag(54) = 'm/s '
768     udiag(55) = 'm/s '
769     udiag(56) = 'K '
770     udiag(57) = 'gram/kg '
771 molod 1.1 udiag(58) = 'hPa/day '
772     udiag(59) = 'm/s/day '
773     udiag(60) = 'm/s/day '
774     udiag(61) = 'K/day '
775     udiag(62) = 'gram/kg/day '
776     udiag(63) = '(e-4 sec-1) '
777     udiag(64) = 'hPa/day '
778     udiag(65) = 'K/day '
779     udiag(66) = 'gram/kg/day '
780     udiag(67) = 'm/sec '
781     udiag(68) = 'm '
782     udiag(69) = '0-1 '
783     udiag(70) = 'hPa '
784     udiag(71) = 'K/day '
785     udiag(72) = 'watt/m2 '
786     udiag(73) = 'watt/m2 '
787     udiag(74) = 'g/m2/sec '
788     udiag(75) = 'm/sec '
789     udiag(76) = 'm/sec '
790     udiag(77) = 'K '
791     udiag(78) = 'gram/kg '
792     udiag(79) = 'K/day '
793     udiag(80) = 'hPa '
794     udiag(81) = '(m/sec)2 '
795     udiag(82) = 'watt/m2 '
796     udiag(83) = 'hPa/day '
797     udiag(84) = 'user-defined '
798     udiag(85) = 'user-defined '
799     udiag(86) = 'user-defined '
800     udiag(87) = 'user-defined '
801     udiag(88) = 'm/s/day '
802     udiag(89) = 'm/s/day '
803     udiag(90) = 'K/day '
804     udiag(91) = 'gram/kg/day '
805     udiag(92) = 'm/s/day '
806     udiag(93) = 'm/s/day '
807     udiag(94) = 'm/s/day '
808     udiag(95) = 'm/s/day '
809     udiag(96) = 'N/m2 '
810     udiag(97) = 'N/m2 '
811     udiag(98) = 'N/m2 '
812     udiag(99) = 'N/m2 '
813     udiag(100) = 'mg/kg '
814     udiag(101) = 'hPa '
815     udiag(102) = 'm/s gram/kg '
816     udiag(103) = 'm/s gram/kg '
817     udiag(104) = 'm/s K '
818     udiag(105) = 'm/s K '
819     udiag(106) = '0-1 '
820     udiag(107) = 'gram/cm2 '
821     udiag(108) = 'm/sec '
822     udiag(109) = 'm/sec '
823     udiag(110) = 'K '
824     udiag(111) = 'gram/kg '
825     udiag(112) = 'm/sec '
826     udiag(113) = 'm/sec '
827     udiag(114) = 'K '
828     udiag(115) = 'gram/kg '
829     udiag(116) = 'g/m2/sec '
830     udiag(117) = 'gram/kg/day '
831     udiag(118) = 'mm/day '
832     udiag(119) = 'mm/day '
833     udiag(120) = 'm/s/day '
834     udiag(121) = 'm/s/day '
835     udiag(122) = 'K/day '
836     udiag(123) = 'gram/kg/day '
837     udiag(124) = 'user-defined '
838     udiag(125) = 'user-defined '
839     udiag(126) = 'user-defined '
840     udiag(127) = 'user-defined '
841     udiag(128) = 'user-defined '
842     udiag(129) = 'user-defined '
843     udiag(130) = 'user-defined '
844     udiag(131) = 'user-defined '
845     udiag(132) = 'user-defined '
846     udiag(133) = 'user-defined '
847     udiag(134) = 'user-defined '
848     udiag(135) = 'user-defined '
849     udiag(136) = 'user-defined '
850     udiag(137) = 'user-defined '
851     udiag(138) = 'user-defined '
852     udiag(139) = 'user-defined '
853     udiag(140) = '0-1 '
854     udiag(141) = '0-1 '
855     udiag(142) = 'non-dimensional '
856     udiag(143) = 'non-dimensional '
857     udiag(144) = 'non-dimensional '
858     udiag(145) = '0-1 '
859     udiag(146) = 'hPa '
860     udiag(147) = 'count '
861     udiag(148) = 'K '
862     udiag(149) = 'count '
863     udiag(150) = 'count '
864     udiag(151) = 'count '
865     udiag(152) = 'count '
866     udiag(153) = 'K '
867     udiag(154) = 'K '
868     udiag(155) = 'gram/kg '
869     udiag(156) = '0-1 '
870     udiag(157) = '0-1 '
871     udiag(158) = '0-1 '
872     udiag(159) = 'mm '
873     udiag(160) = 'mm water equiv '
874     udiag(161) = 'mm/day '
875     udiag(162) = 'mm/day '
876     udiag(163) = 'mm/day '
877     udiag(164) = 'mm/sec '
878     udiag(165) = 'mm/sec '
879     udiag(166) = 'mm/sec '
880     udiag(167) = 'mm/sec '
881     udiag(168) = 'watt/m2 '
882     udiag(169) = 'watt/m2 '
883     udiag(170) = 'watt/m2 '
884     udiag(171) = 'watt/m2 '
885     udiag(172) = 'watt/m2 '
886     udiag(173) = 'watt/m2 '
887     udiag(174) = 'percent '
888     udiag(175) = 'percent '
889     udiag(176) = 'watt/m2/K '
890     udiag(177) = 'watt/m2/K '
891     udiag(178) = 'watt/m2/K '
892     udiag(179) = 'K/day '
893     udiag(180) = 'K/day '
894     udiag(181) = 'K/day '
895     udiag(182) = 'K/day '
896     udiag(183) = 'K/day '
897     udiag(184) = 'hPa '
898     udiag(185) = 'count '
899     udiag(186) = 'hPa '
900     udiag(187) = 'user-defined '
901     udiag(188) = 'count '
902     udiag(189) = ' '
903     udiag(190) = 'watt/m2 '
904     udiag(191) = 'mm/day '
905     udiag(192) = 'watt/m2 '
906     udiag(193) = 'watt/m2 '
907     udiag(194) = 'K/sec '
908     udiag(195) = 'sec-1 '
909     udiag(196) = 'K/sec '
910     udiag(197) = 'K/sec '
911     udiag(198) = 'K/sec '
912     udiag(199) = 'K/sec '
913     udiag(200) = 'K/sec '
914     udiag(201) = 'K/sec '
915     udiag(202) = 'mg/kg '
916     udiag(203) = '0-1 '
917     udiag(204) = '0-1 '
918     udiag(205) = '0-1 '
919     udiag(206) = '0-1 '
920     udiag(207) = '0-1 '
921     udiag(208) = '0-1 '
922     udiag(209) = 'non-dimensional '
923     udiag(210) = 'count '
924     udiag(211) = 'K/day '
925     udiag(212) = 'g/kg/day '
926     udiag(213) = '0-1 '
927     udiag(214) = 'mg/cm2 '
928     udiag(215) = 'mg/cm2 '
929     udiag(216) = 'mg/cm2 '
930     udiag(217) = 'mg/cm2 '
931     udiag(218) = 'mb m2/sec2 '
932     udiag(219) = 'm2/sec2 '
933     udiag(220) = 'm2/sec2 '
934     udiag(221) = 'K/day '
935     udiag(222) = 'g/kg/day '
936     udiag(223) = 'K/day '
937     udiag(224) = 'g/kg/day '
938     udiag(225) = 'K/day '
939     udiag(226) = 'g/kg/day '
940     udiag(227) = 'K/day '
941     udiag(228) = 'g/kg/day '
942     udiag(229) = 'm/sec '
943     udiag(230) = 'm/sec '
944     udiag(231) = 'm '
945     udiag(232) = 'm/sec '
946     udiag(233) = 'm/sec '
947     udiag(234) = 'N/m2 '
948     udiag(235) = 'N/m2 '
949     udiag(236) = 'm/sec/day '
950     udiag(237) = 'm/sec/day '
951 molod 1.7 udiag(238) = 'm or hPa '
952     udiag(239) = 'm2 or hPa2 '
953     udiag(240) = 'K '
954     udiag(241) = 'g/kg '
955     udiag(242) = 'm/sec '
956     udiag(243) = 'm/sec '
957     udiag(244) = 'm/sec '
958     udiag(245) = 'K2 '
959     udiag(246) = 'g2/kg2 '
960     udiag(247) = 'm2/sec2 '
961     udiag(248) = 'm2/sec2 '
962     udiag(249) = 'm2/sec2 '
963     udiag(250) = 'm2/sec2 '
964 molod 1.8 udiag(251) = 'm/sec '
965     udiag(252) = 'm/sec '
966     udiag(253) = 'm/sec '
967     udiag(254) = 'm-K/sec '
968     udiag(255) = 'm-K/sec '
969     udiag(256) = 'm-K/sec '
970     udiag(257) = 'm-kg/s-kg '
971     udiag(258) = 'm-kg/s-kg '
972     udiag(259) = 'm-kg/s-kg '
973 molod 1.9 udiag(260) = 'm-K/sec '
974     udiag(261) = 'm-K/sec '
975     udiag(262) = 'm-K/sec '
976     udiag(263) = 'm-kg/s-kg '
977     udiag(264) = 'm-kg/s-kg '
978     udiag(265) = 'm-kg/s-kg '
979 molod 1.10 udiag(266) = 'm-kg/s-kg '
980     udiag(267) = 'm-kg/s-kg '
981     udiag(268) = 'm-kg/s-kg '
982     udiag(269) = 'm-kg/s-kg '
983     udiag(270) = 'm-kg/s-kg '
984     udiag(271) = 'm-kg/s-kg '
985     udiag(272) = 'm-kg/s-kg '
986     udiag(273) = 'm-kg/s-kg '
987     udiag(274) = 'm-kg/s-kg '
988     udiag(275) = 'm-kg/s-kg '
989     udiag(276) = 'm-kg/s-kg '
990     udiag(277) = 'm-kg/s-kg '
991     udiag(278) = 'm-kg/s-kg '
992     udiag(279) = 'm-kg/s-kg '
993     udiag(280) = 'm-kg/s-kg '
994 molod 1.12 udiag(281) = 'kg/kg '
995     udiag(282) = 'kg/kg '
996     udiag(283) = 'kg/kg '
997     udiag(284) = 'kg/kg '
998     udiag(285) = 'kg/kg '
999 molod 1.1
1000     c Diagnostic Parms
1001     c ----------------
1002 molod 1.7 gdiag( 1) = 'UM 002L '
1003     gdiag( 2) = 'VM 001L '
1004     gdiag( 3) = 'SM L '
1005     gdiag( 4) = 'SM L '
1006     gdiag( 5) = 'SM L '
1007     gdiag( 6) = 'SM L '
1008     gdiag( 7) = 'SM P L '
1009     gdiag( 8) = 'SM L '
1010     gdiag( 9) = 'SM L '
1011     gdiag(10) = 'SM L '
1012     gdiag(11) = 'SM L '
1013     gdiag(12) = 'SM L '
1014     gdiag(13) = 'UM 014M '
1015     gdiag(14) = 'VM 013M '
1016     gdiag(15) = 'SM M '
1017     gdiag(16) = 'SM M '
1018     gdiag(17) = 'SM M '
1019     gdiag(18) = 'SM M '
1020 molod 1.3 gdiag(19) = 'SM M '
1021     gdiag(20) = 'SM P M '
1022 molod 1.7 gdiag(21) = 'SM P L '
1023     gdiag(22) = 'SM P L '
1024     gdiag(23) = 'UM 024L '
1025     gdiag(24) = 'VM 023L '
1026     gdiag(25) = 'SM L '
1027     gdiag(26) = 'SM L '
1028     gdiag(27) = 'SM L '
1029     gdiag(28) = 'SM L '
1030     gdiag(29) = 'SM L '
1031     gdiag(30) = 'SM L '
1032     gdiag(31) = 'SM L '
1033     gdiag(32) = 'SM L '
1034     gdiag(33) = 'SM P L '
1035     gdiag(34) = 'SM P L '
1036     gdiag(35) = 'SM L '
1037     gdiag(36) = 'SM L '
1038     gdiag(37) = 'SM U '
1039     gdiag(38) = 'SM U '
1040     gdiag(39) = 'SM L '
1041 molod 1.3 gdiag(40) = 'SM M '
1042 molod 1.7 gdiag(41) = 'SM M '
1043     gdiag(42) = 'SM P M '
1044     gdiag(43) = 'SM M '
1045     gdiag(44) = 'SM P M '
1046     gdiag(45) = 'SM P M '
1047     gdiag(46) = 'SM L '
1048     gdiag(47) = 'SM M '
1049     gdiag(48) = 'SM P U '
1050     gdiag(49) = 'SM C050M '
1051     gdiag(50) = 'SM D M '
1052     gdiag(51) = 'SM P M '
1053     gdiag(52) = 'SM L '
1054 molod 1.3 gdiag(53) = 'SM L '
1055     gdiag(54) = 'UU 055M '
1056     gdiag(55) = 'VV 054M '
1057     gdiag(56) = 'SM M '
1058     gdiag(57) = 'SM M '
1059     gdiag(58) = 'SM M '
1060     gdiag(59) = 'UU 060M '
1061     gdiag(60) = 'VV 059M '
1062     gdiag(61) = 'SM M '
1063     gdiag(62) = 'SM M '
1064     gdiag(63) = 'SM M '
1065     gdiag(64) = 'SM L '
1066 molod 1.7 gdiag(65) = 'SM M '
1067     gdiag(66) = 'SM M '
1068     gdiag(67) = 'SM L '
1069     gdiag(68) = 'SM L '
1070     gdiag(69) = 'SM P M '
1071     gdiag(70) = 'SM P M '
1072 molod 1.3 gdiag(71) = 'SM P M '
1073 molod 1.7 gdiag(72) = 'SM P U '
1074     gdiag(73) = 'SM P U '
1075     gdiag(74) = 'SM P L '
1076 molod 1.3 gdiag(75) = 'UU 076M '
1077     gdiag(76) = 'VV 075M '
1078     gdiag(77) = 'SM M '
1079     gdiag(78) = 'SM P M '
1080     gdiag(79) = 'SM M '
1081     gdiag(80) = 'SM L '
1082     gdiag(81) = 'SM P L '
1083 molod 1.7 gdiag(82) = 'SM P L '
1084 molod 1.3 gdiag(83) = 'SM L '
1085 molod 1.7 gdiag(84) = 'SM L '
1086     gdiag(85) = 'SM L '
1087     gdiag(86) = 'SM M '
1088     gdiag(87) = 'SM M '
1089 molod 1.3 gdiag(88) = 'UU 089M '
1090     gdiag(89) = 'VV 088M '
1091     gdiag(90) = 'SM M '
1092     gdiag(91) = 'SM M '
1093     gdiag(92) = 'UU 093M '
1094     gdiag(93) = 'VV 092M '
1095 molod 1.7 gdiag(94) = 'UM 095M '
1096     gdiag(95) = 'VM 094M '
1097     gdiag(96) = 'UM 097L '
1098     gdiag(97) = 'VM 096L '
1099     gdiag(98) = 'UM 099U '
1100     gdiag(99) = 'VM 098U '
1101     gdiag(100) = 'SM M '
1102 molod 1.3 gdiag(101) = 'SM L '
1103     gdiag(102) = 'UM 103M '
1104     gdiag(103) = 'VM 102M '
1105     gdiag(104) = 'UM 105M '
1106     gdiag(105) = 'VM 104M '
1107 molod 1.7 gdiag(106) = 'SM P M '
1108 molod 1.3 gdiag(107) = 'SM P M '
1109 molod 1.7 gdiag(108) = 'UM 109M '
1110     gdiag(109) = 'VM 108M '
1111     gdiag(110) = 'SM M '
1112     gdiag(111) = 'SM P M '
1113     gdiag(112) = 'UM 113M '
1114     gdiag(113) = 'VM 112M '
1115     gdiag(114) = 'SM M '
1116     gdiag(115) = 'SM P M '
1117     gdiag(116) = 'SM P L '
1118 molod 1.3 gdiag(117) = 'SM M '
1119     gdiag(118) = 'SM M '
1120     gdiag(119) = 'SM P L '
1121     gdiag(120) = 'UU 121M '
1122     gdiag(121) = 'VV 120M '
1123     gdiag(122) = 'SM M '
1124     gdiag(123) = 'SM M '
1125 molod 1.7 gdiag(124) = 'SM L '
1126     gdiag(125) = 'SM L '
1127     gdiag(126) = 'SM L '
1128     gdiag(127) = 'SM L '
1129     gdiag(128) = 'SM L '
1130     gdiag(129) = 'SM L '
1131     gdiag(130) = 'SM L '
1132     gdiag(131) = 'SM L '
1133     gdiag(132) = 'SM M '
1134     gdiag(133) = 'SM M '
1135     gdiag(134) = 'SM M '
1136     gdiag(135) = 'SM M '
1137     gdiag(136) = 'SM M '
1138     gdiag(137) = 'SM M '
1139     gdiag(138) = 'SM M '
1140     gdiag(139) = 'SM M '
1141     gdiag(140) = 'SM P M '
1142     gdiag(141) = 'SM P M '
1143     gdiag(142) = 'SM PC150M '
1144     gdiag(143) = 'SM PC151M '
1145     gdiag(144) = 'SM PC152M '
1146     gdiag(145) = 'SM M '
1147     gdiag(146) = 'SM PC147M '
1148     gdiag(147) = 'SM D M '
1149     gdiag(148) = 'SM PC149M '
1150     gdiag(149) = 'SM D M '
1151     gdiag(150) = 'SM D M '
1152     gdiag(151) = 'SM D M '
1153     gdiag(152) = 'SM D M '
1154     gdiag(153) = 'SM L '
1155     gdiag(154) = 'SM L '
1156     gdiag(155) = 'SM L '
1157     gdiag(156) = 'SM L '
1158     gdiag(157) = 'SM L '
1159     gdiag(158) = 'SM L '
1160     gdiag(159) = 'SM L '
1161     gdiag(160) = 'SM L '
1162     gdiag(161) = 'SM L '
1163     gdiag(162) = 'SM L '
1164     gdiag(163) = 'SM P L '
1165     gdiag(164) = 'SM L '
1166     gdiag(165) = 'SM L '
1167     gdiag(166) = 'SM L '
1168     gdiag(167) = 'SM L '
1169     gdiag(168) = 'SM L '
1170     gdiag(169) = 'SM L '
1171     gdiag(170) = 'SM L '
1172     gdiag(171) = 'SM L '
1173     gdiag(172) = 'SM P L '
1174     gdiag(173) = 'SM P L '
1175     gdiag(174) = 'SM L '
1176     gdiag(175) = 'SM L '
1177     gdiag(176) = 'SM L '
1178     gdiag(177) = 'SM L '
1179     gdiag(178) = 'SM L '
1180 molod 1.3 gdiag(179) = 'SM M '
1181     gdiag(180) = 'SM M '
1182     gdiag(181) = 'SM M '
1183     gdiag(182) = 'SM P M '
1184     gdiag(183) = 'SM M '
1185 molod 1.7 gdiag(184) = 'SM PC185M '
1186     gdiag(185) = 'SM PD M '
1187     gdiag(186) = 'SM L '
1188     gdiag(187) = 'SM C188L '
1189     gdiag(188) = 'SM D L '
1190 molod 1.1 gdiag(189) = ' '
1191 molod 1.7 gdiag(190) = 'SM L '
1192 molod 1.3 gdiag(191) = 'SM P L '
1193 molod 1.7 gdiag(192) = 'SM L '
1194     gdiag(193) = 'SM L '
1195     gdiag(194) = 'SM L '
1196     gdiag(195) = 'SM L '
1197     gdiag(196) = 'SM L '
1198     gdiag(197) = 'SM L '
1199     gdiag(198) = 'SM L '
1200     gdiag(199) = 'SM L '
1201     gdiag(200) = 'SM L '
1202     gdiag(201) = 'SM L '
1203     gdiag(202) = 'SM M '
1204     gdiag(203) = 'SM M '
1205     gdiag(204) = 'SM C210L '
1206     gdiag(205) = 'SM L '
1207     gdiag(206) = 'SM L '
1208     gdiag(207) = 'SM L '
1209     gdiag(208) = 'SM L '
1210     gdiag(209) = 'SM M '
1211     gdiag(210) = 'SM D L '
1212 molod 1.3 gdiag(211) = 'SM M '
1213     gdiag(212) = 'SM M '
1214 molod 1.7 gdiag(213) = 'SM P M '
1215     gdiag(214) = 'SM P L '
1216     gdiag(215) = 'SM P L '
1217     gdiag(216) = 'SM P L '
1218     gdiag(217) = 'SM P L '
1219     gdiag(218) = 'SM P L '
1220     gdiag(219) = 'SM P L '
1221     gdiag(220) = 'SM P M '
1222     gdiag(221) = 'SM M '
1223     gdiag(222) = 'SM M '
1224     gdiag(223) = 'SM M '
1225     gdiag(224) = 'SM M '
1226     gdiag(225) = 'SM M '
1227     gdiag(226) = 'SM M '
1228     gdiag(227) = 'SM M '
1229     gdiag(228) = 'SM M '
1230     gdiag(229) = 'UM 230L '
1231     gdiag(230) = 'VM 229L '
1232     gdiag(231) = 'SM L '
1233     gdiag(232) = 'UM 233L '
1234     gdiag(233) = 'VM 232L '
1235     gdiag(234) = 'UM 235L '
1236     gdiag(235) = 'VM 234L '
1237     gdiag(236) = 'UM 237M '
1238     gdiag(237) = 'VM 236M '
1239     gdiag(238) = 'SM M '
1240     gdiag(239) = 'SM M '
1241     gdiag(240) = 'SM M '
1242     gdiag(241) = 'SM M '
1243     gdiag(242) = 'UU 243M '
1244     gdiag(243) = 'VV 242M '
1245     gdiag(244) = 'WM L '
1246     gdiag(245) = 'SM M '
1247     gdiag(246) = 'SM M '
1248 molod 1.11 gdiag(247) = 'UU 248M '
1249     gdiag(248) = 'VV 247M '
1250 molod 1.7 gdiag(249) = 'WM L '
1251     gdiag(250) = 'UZ M '
1252 molod 1.8 gdiag(251) = 'UU 252M '
1253     gdiag(252) = 'VV 251M '
1254     gdiag(253) = 'WM M '
1255     gdiag(254) = 'UU 255M '
1256     gdiag(255) = 'VV 254M '
1257     gdiag(256) = 'WM M '
1258     gdiag(257) = 'UU 258M '
1259     gdiag(258) = 'VV 257M '
1260     gdiag(259) = 'WM M '
1261 molod 1.9 gdiag(260) = 'UU 261M '
1262     gdiag(261) = 'VV 260M '
1263     gdiag(262) = 'WM M '
1264     gdiag(263) = 'UU 264M '
1265     gdiag(264) = 'VV 263M '
1266     gdiag(265) = 'WM M '
1267 molod 1.10 gdiag(266) = 'UU 267M '
1268     gdiag(267) = 'VV 266M '
1269     gdiag(268) = 'WM M '
1270     gdiag(269) = 'UU 270M '
1271     gdiag(270) = 'VV 269M '
1272     gdiag(271) = 'WM M '
1273     gdiag(272) = 'UU 273M '
1274     gdiag(273) = 'VV 272M '
1275     gdiag(274) = 'WM M '
1276     gdiag(275) = 'UU 276M '
1277     gdiag(276) = 'VV 275M '
1278     gdiag(277) = 'WM M '
1279     gdiag(278) = 'UU 279M '
1280     gdiag(279) = 'VV 278M '
1281     gdiag(280) = 'WM M '
1282 molod 1.12 gdiag(281) = 'SM P M '
1283     gdiag(282) = 'SM P M '
1284     gdiag(283) = 'SM P M '
1285     gdiag(284) = 'SM P M '
1286     gdiag(285) = 'SM P M '
1287 molod 1.1
1288     c Diagnostic Levels
1289     c -----------------
1290     KDIAG( 1) = 1
1291     KDIAG( 2) = 1
1292     KDIAG( 3) = 1
1293     KDIAG( 4) = 1
1294     KDIAG( 5) = 1
1295     KDIAG( 6) = 1
1296     KDIAG( 7) = 1
1297     KDIAG( 8) = nrphys
1298     KDIAG( 9) = 1
1299     KDIAG(10) = 1
1300     KDIAG(11) = nrphys
1301     KDIAG(12) = nrphys
1302     KDIAG(13) = nrphys
1303     KDIAG(14) = nrphys
1304     KDIAG(15) = nrphys
1305     KDIAG(16) = nrphys
1306     KDIAG(17) = nrphys
1307     KDIAG(18) = nrphys
1308     KDIAG(19) = nrphys
1309     KDIAG(20) = nrphys
1310     KDIAG(21) = 1
1311     KDIAG(22) = 1
1312     KDIAG(23) = nrphys
1313     KDIAG(24) = nrphys
1314     KDIAG(25) = nrphys
1315     KDIAG(26) = nrphys
1316     KDIAG(27) = 1
1317     KDIAG(28) = 1
1318     KDIAG(29) = 1
1319     KDIAG(30) = 1
1320     KDIAG(31) = 1
1321     KDIAG(32) = 1
1322     KDIAG(33) = 1
1323     KDIAG(34) = 1
1324     KDIAG(35) = 1
1325     KDIAG(36) = 1
1326     KDIAG(37) = 1
1327     KDIAG(38) = 1
1328     KDIAG(39) = 1
1329     KDIAG(40) = nrphys
1330     KDIAG(41) = nrphys
1331     KDIAG(42) = nrphys
1332     KDIAG(43) = nrphys
1333     KDIAG(44) = nrphys
1334     KDIAG(45) = nrphys
1335     KDIAG(46) = 1
1336     KDIAG(47) = nrphys
1337     KDIAG(48) = 1
1338     KDIAG(49) = nrphys
1339     KDIAG(50) = nrphys
1340     KDIAG(51) = 1
1341     KDIAG(52) = 1
1342     KDIAG(53) = 1
1343 molod 1.5 KDIAG(54) = Nr
1344     KDIAG(55) = Nr
1345     KDIAG(56) = Nr
1346     KDIAG(57) = Nr
1347 molod 1.1 KDIAG(58) = nrphys
1348     KDIAG(59) = nrphys
1349     KDIAG(60) = nrphys
1350     KDIAG(61) = nrphys
1351     KDIAG(62) = nrphys
1352     KDIAG(63) = nrphys
1353     KDIAG(64) = nrphys
1354     KDIAG(65) = nrphys
1355     KDIAG(66) = nrphys
1356     KDIAG(67) = 1
1357     KDIAG(68) = 1
1358     KDIAG(69) = nrphys-1
1359     KDIAG(70) = 1
1360     KDIAG(71) = nrphys
1361     KDIAG(72) = 1
1362     KDIAG(73) = 1
1363     KDIAG(74) = nrphys
1364     KDIAG(75) = nrphys
1365     KDIAG(76) = nrphys
1366     KDIAG(77) = nrphys
1367     KDIAG(78) = nrphys
1368     KDIAG(79) = nrphys
1369     KDIAG(80) = 1
1370     KDIAG(81) = nrphys
1371     KDIAG(82) = 1
1372     KDIAG(83) = 1
1373     KDIAG(84) = 1
1374     KDIAG(85) = 1
1375     KDIAG(86) = nrphys
1376     KDIAG(87) = nrphys
1377     KDIAG(88) = nrphys
1378     KDIAG(89) = nrphys
1379     KDIAG(90) = nrphys
1380     KDIAG(91) = nrphys
1381     KDIAG(92) = nrphys
1382     KDIAG(93) = nrphys
1383     KDIAG(94) = nrphys
1384     KDIAG(95) = nrphys
1385     KDIAG(96) = 1
1386     KDIAG(97) = 1
1387     KDIAG(98) = 1
1388     KDIAG(99) = 1
1389     KDIAG(100)= nrphys
1390     KDIAG(101)= 1
1391     KDIAG(102)= 1
1392     KDIAG(103)= 1
1393     KDIAG(104)= 1
1394     KDIAG(105)= 1
1395     KDIAG(106)= 1
1396     KDIAG(107)= 1
1397     KDIAG(108)= 1
1398     KDIAG(109)= 1
1399     KDIAG(110)= 1
1400     KDIAG(111)= 1
1401     KDIAG(112)= 1
1402     KDIAG(113)= 1
1403     KDIAG(114)= 1
1404     KDIAG(115)= 1
1405     KDIAG(116)= nrphys
1406     KDIAG(117)= nrphys
1407     KDIAG(118)= 1
1408     KDIAG(119)= 1
1409     KDIAG(120)= nrphys
1410     KDIAG(121)= nrphys
1411     KDIAG(122)= nrphys
1412     KDIAG(123)= nrphys
1413     KDIAG(124)= 1
1414     KDIAG(125)= 1
1415     KDIAG(126)= 1
1416     KDIAG(127)= 1
1417     KDIAG(128)= 1
1418     KDIAG(129)= 1
1419     KDIAG(130)= 1
1420     KDIAG(131)= 1
1421     KDIAG(132)= nrphys
1422     KDIAG(133)= nrphys
1423     KDIAG(134)= nrphys
1424     KDIAG(135)= nrphys
1425     KDIAG(136)= nrphys
1426     KDIAG(137)= nrphys
1427     KDIAG(138)= nrphys
1428     KDIAG(139)= nrphys
1429     KDIAG(140)= 1
1430     KDIAG(141)= 1
1431     KDIAG(142)= 1
1432     KDIAG(143)= 1
1433     KDIAG(144)= 1
1434     KDIAG(145)= nrphys
1435     KDIAG(146)= 1
1436     KDIAG(147)= 1
1437     KDIAG(148)= 1
1438     KDIAG(149)= 1
1439     KDIAG(150)= 1
1440     KDIAG(151)= 1
1441     KDIAG(152)= 1
1442     KDIAG(153)= 1
1443     KDIAG(154)= 1
1444     KDIAG(155)= 1
1445     KDIAG(156)= 1
1446     KDIAG(157)= 1
1447     KDIAG(158)= 1
1448     KDIAG(159)= 1
1449     KDIAG(160)= 1
1450     KDIAG(161)= 1
1451     KDIAG(162)= 1
1452     KDIAG(163)= 1
1453     KDIAG(164)= 1
1454     KDIAG(165)= 1
1455     KDIAG(166)= 1
1456     KDIAG(167)= 1
1457     KDIAG(168)= 1
1458     KDIAG(169)= 1
1459     KDIAG(170)= 1
1460     KDIAG(171)= 1
1461     KDIAG(172)= 1
1462     KDIAG(173)= 1
1463     KDIAG(174)= 1
1464     KDIAG(175)= 1
1465     KDIAG(176)= 1
1466     KDIAG(177)= 1
1467     KDIAG(178)= 1
1468     KDIAG(179)= 1
1469     KDIAG(180)= 1
1470     KDIAG(181)= 1
1471     KDIAG(182)= 1
1472     KDIAG(183)= 1
1473     KDIAG(184)= 1
1474     KDIAG(185)= 1
1475     KDIAG(186)= 1
1476     KDIAG(187)= 1
1477     KDIAG(188)= 1
1478     KDIAG(189)= 0
1479     KDIAG(190)= 1
1480     KDIAG(191)= 1
1481     KDIAG(192)= 1
1482     KDIAG(193)= 1
1483     KDIAG(194)= 1
1484     KDIAG(195)= 1
1485     KDIAG(196)= 1
1486     KDIAG(197)= 1
1487     KDIAG(198)= 1
1488     KDIAG(199)= 1
1489     KDIAG(200)= 1
1490     KDIAG(201)= 1
1491     KDIAG(202)= nrphys
1492     KDIAG(203)= nrphys
1493     KDIAG(204)= 1
1494     KDIAG(205)= 1
1495     KDIAG(206)= 1
1496     KDIAG(207)= 1
1497     KDIAG(208)= 1
1498     KDIAG(209)= nrphys
1499     KDIAG(210)= 1
1500     KDIAG(211)= nrphys
1501     KDIAG(212)= nrphys
1502     KDIAG(213)= nrphys
1503     KDIAG(214)= 1
1504     KDIAG(215)= 1
1505     KDIAG(216)= 1
1506     KDIAG(217)= 1
1507     KDIAG(218)= 1
1508     KDIAG(219)= 1
1509     KDIAG(220)= nrphys
1510     KDIAG(221)= nrphys
1511     KDIAG(222)= nrphys
1512     KDIAG(223)= nrphys
1513     KDIAG(224)= nrphys
1514     KDIAG(225)= nrphys
1515     KDIAG(226)= nrphys
1516     KDIAG(227)= nrphys
1517     KDIAG(228)= nrphys
1518     KDIAG(229)= 1
1519     KDIAG(230)= 1
1520     KDIAG(231)= 1
1521     KDIAG(232)= 1
1522     KDIAG(233)= 1
1523     KDIAG(234)= 1
1524     KDIAG(235)= 1
1525     KDIAG(236)= nrphys
1526     KDIAG(237)= nrphys
1527 molod 1.7 KDIAG(238)= 1
1528     KDIAG(239)= 1
1529     KDIAG(240)= Nr
1530     KDIAG(241)= Nr
1531     KDIAG(242)= Nr
1532     KDIAG(243)= Nr
1533     KDIAG(244)= Nr
1534     KDIAG(245)= Nr
1535     KDIAG(246)= Nr
1536     KDIAG(247)= Nr
1537     KDIAG(248)= Nr
1538     KDIAG(249)= Nr
1539     KDIAG(250)= Nr
1540 molod 1.8 KDIAG(251)= Nr
1541     KDIAG(252)= Nr
1542     KDIAG(253)= Nr
1543     KDIAG(254)= Nr
1544     KDIAG(255)= Nr
1545     KDIAG(256)= Nr
1546     KDIAG(257)= Nr
1547     KDIAG(258)= Nr
1548     KDIAG(259)= Nr
1549 molod 1.9 KDIAG(260)= Nr
1550     KDIAG(261)= Nr
1551     KDIAG(262)= Nr
1552     KDIAG(263)= Nr
1553     KDIAG(264)= Nr
1554     KDIAG(265)= Nr
1555 molod 1.10 KDIAG(266)= Nr
1556     KDIAG(267)= Nr
1557     KDIAG(268)= Nr
1558     KDIAG(269)= Nr
1559     KDIAG(270)= Nr
1560     KDIAG(271)= Nr
1561     KDIAG(272)= Nr
1562     KDIAG(273)= Nr
1563     KDIAG(274)= Nr
1564     KDIAG(275)= Nr
1565     KDIAG(276)= Nr
1566     KDIAG(277)= Nr
1567     KDIAG(278)= Nr
1568     KDIAG(279)= Nr
1569     KDIAG(280)= Nr
1570 molod 1.12 KDIAG(281)= Nr
1571     KDIAG(282)= Nr
1572     KDIAG(283)= Nr
1573     KDIAG(284)= Nr
1574     KDIAG(285)= Nr
1575 molod 1.1
1576     do n=1,ndiagt
1577     idiag (n) = 0
1578 molod 1.5 ndiag(n) = 0
1579 molod 1.1 enddo
1580    
1581     RETURN
1582     END

  ViewVC Help
Powered by ViewVC 1.1.22