/[MITgcm]/MITgcm/pkg/fizhi/fizhi_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/fizhi/fizhi_diagnostics_init.F

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


Revision 1.2 - (hide annotations) (download)
Thu Jul 28 00:09:21 2005 UTC (18 years, 11 months ago) by molod
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57n_post
Changes since 1.1: +89 -89 lines
Fix fizhi diagnostic initialization of counter, mate numbers and vector locations

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

  ViewVC Help
Powered by ViewVC 1.1.22