/[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.11 - (hide annotations) (download)
Fri Jun 9 16:05:28 2006 UTC (18 years ago) by molod
Branch: MAIN
CVS Tags: checkpoint58j_post, checkpoint58i_post
Changes since 1.10: +2 -2 lines
Fix definition and set of fizhi meridional wind diagnostic

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

  ViewVC Help
Powered by ViewVC 1.1.22