/[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.3 - (hide annotations) (download)
Thu Feb 26 19:52:05 2004 UTC (20 years, 4 months ago) by molod
Branch: MAIN
Changes since 1.2: +52 -53 lines
Still fixing bugs

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

  ViewVC Help
Powered by ViewVC 1.1.22