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

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

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


Revision 1.6 - (hide annotations) (download)
Mon May 16 15:07:45 2005 UTC (19 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57h_done
Changes since 1.5: +1 -2 lines
move initialization of idiag from _init_early to _set_pointers

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

  ViewVC Help
Powered by ViewVC 1.1.22