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

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

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


Revision 1.7 - (show annotations) (download)
Wed Feb 15 20:51:36 2006 UTC (18 years, 4 months ago) by molod
Branch: MAIN
CVS Tags: checkpoint58e_post, checkpoint58d_post, checkpoint58c_post, checkpoint58b_post
Changes since 1.6: +33 -1 lines
Add total tend on fizhi grid

1 C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_diagnostics_init.F,v 1.6 2006/01/13 23:18:48 molod Exp $
2 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 = '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 diagName = 'UFLUX '
62 diagTitle = 'Zonal Wind Surface Stress (N/m^2)'
63 diagUnits = 'N/m2 '
64 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'L1 '
65 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 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'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 diagCode = 'SM L1 '
93 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 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
135 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 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
142 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 diagCode = 'SM ML '
190 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 diagCode = 'SM L1 '
197 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 diagCode = 'SM L1 '
204 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 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'LL '
211 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 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'LL '
218 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 diagCode = 'SM U1 '
388 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 WRITE(diagCode,'(A,I3.3,A)') 'SM C',diagNum+2,'ML '
396 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 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'MR '
431 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 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'MR '
438 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 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
467 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 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
474 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 = 'DTLS '
499 diagTitle =
500 . 'Temperature Tendency due to Stratiform Processes (K/day)'
501 diagUnits = 'K/day '
502 diagCode = 'SM ML '
503 CALL DIAGNOSTICS_ADD2LIST( diagNum,
504 I diagName, diagCode, diagUnits, diagTitle, myThid )
505
506 diagName = 'DQLS '
507 diagTitle =
508 . 'Specific Humidity Tendency from Stratiform Processes (g/kg/day)'
509 diagUnits = 'gram/kg/day '
510 diagCode = 'SM ML '
511 CALL DIAGNOSTICS_ADD2LIST( diagNum,
512 I diagName, diagCode, diagUnits, diagTitle, myThid )
513
514 diagName = 'USTAR '
515 diagTitle = 'Friction Velocity (m/s)'
516 diagUnits = 'm/sec '
517 diagCode = 'SM L1 '
518 CALL DIAGNOSTICS_ADD2LIST( diagNum,
519 I diagName, diagCode, diagUnits, diagTitle, myThid )
520
521 diagName = 'Z0 '
522 diagTitle = 'Surface Roughness Z0 (m)'
523 diagUnits = 'm '
524 diagCode = 'SM L1 '
525 CALL DIAGNOSTICS_ADD2LIST( diagNum,
526 I diagName, diagCode, diagUnits, diagTitle, myThid )
527
528 diagName = 'FRQTRB '
529 diagTitle = 'Frequency of Turbulence (0-1)'
530 diagUnits = '0-1 '
531 diagCode = 'SM P MM '
532 CALL DIAGNOSTICS_ADD2LIST( diagNum,
533 I diagName, diagCode, diagUnits, diagTitle, myThid )
534
535 diagName = 'PBL '
536 diagTitle = 'Estimated PBL Depth (hPa)'
537 diagUnits = 'hPa '
538 diagCode = 'SM P M1 '
539 CALL DIAGNOSTICS_ADD2LIST( diagNum,
540 I diagName, diagCode, diagUnits, diagTitle, myThid )
541
542 diagName = 'SWCLR '
543 diagTitle = 'Clear Sky Shortwave Heating Rates (K/day)'
544 diagUnits = 'K/day '
545 diagCode = 'SM ML '
546 CALL DIAGNOSTICS_ADD2LIST( diagNum,
547 I diagName, diagCode, diagUnits, diagTitle, myThid )
548
549 diagName = 'OSR '
550 diagTitle = 'Outgoing Shortwave Radiation (W/m^2)'
551 diagUnits = 'watt/m2 '
552 diagCode = 'SM U1 '
553 CALL DIAGNOSTICS_ADD2LIST( diagNum,
554 I diagName, diagCode, diagUnits, diagTitle, myThid )
555
556 diagName = 'OSRCLR '
557 diagTitle = 'Clear Sky Outgoing Shortwave Radiation (W/m^2)'
558 diagUnits = 'watt/m2 '
559 diagCode = 'SM U1 '
560 CALL DIAGNOSTICS_ADD2LIST( diagNum,
561 I diagName, diagCode, diagUnits, diagTitle, myThid )
562
563 diagName = 'CLDMAS '
564 diagTitle = 'Cloud Mass Flux (g/m^2/sec)'
565 diagUnits = 'g/m2/sec '
566 diagCode = 'SM P LL '
567 CALL DIAGNOSTICS_ADD2LIST( diagNum,
568 I diagName, diagCode, diagUnits, diagTitle, myThid )
569
570 diagName = 'UWND '
571 diagTitle = 'Time-Averaged Zonal Wind (m/s)'
572 diagUnits = 'm/sec '
573 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
574 CALL DIAGNOSTICS_ADD2LIST( diagNum,
575 I diagName, diagCode, diagUnits, diagTitle, myThid )
576
577 diagName = 'VWND '
578 diagTitle = 'Time-Averaged Meridional Wind (m/s)'
579 diagUnits = 'm/sec '
580 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
581 CALL DIAGNOSTICS_ADD2LIST( diagNum,
582 I diagName, diagCode, diagUnits, diagTitle, myThid )
583
584 diagName = 'TMPU '
585 diagTitle = 'Time-Averaged Temperature (deg K)'
586 diagUnits = 'K '
587 diagCode = 'SM ML '
588 CALL DIAGNOSTICS_ADD2LIST( diagNum,
589 I diagName, diagCode, diagUnits, diagTitle, myThid )
590
591 diagName = 'SPHU '
592 diagTitle = 'Time-Averaged Specific Humidity (g/kg)'
593 diagUnits = 'gram/kg '
594 diagCode = 'SM P ML '
595 CALL DIAGNOSTICS_ADD2LIST( diagNum,
596 I diagName, diagCode, diagUnits, diagTitle, myThid )
597
598 diagName = 'RFT '
599 diagTitle =
600 . 'Temperature Tendency due to Rayleigh Friction (K/day)'
601 diagUnits = 'K/day '
602 diagCode = 'SM ML '
603 CALL DIAGNOSTICS_ADD2LIST( diagNum,
604 I diagName, diagCode, diagUnits, diagTitle, myThid )
605
606 diagName = 'PS '
607 diagTitle = 'Time-Averaged Surface Pressure (hPa)'
608 diagUnits = 'hPa '
609 diagCode = 'SM L1 '
610 CALL DIAGNOSTICS_ADD2LIST( diagNum,
611 I diagName, diagCode, diagUnits, diagTitle, myThid )
612
613 diagName = 'TKE '
614 diagTitle = 'Time-Averaged Turbulent Kinetic Energy (m/s)^2'
615 diagUnits = '(m/sec)2 '
616 diagCode = 'SM P LL '
617 CALL DIAGNOSTICS_ADD2LIST( diagNum,
618 I diagName, diagCode, diagUnits, diagTitle, myThid )
619
620 diagName = 'SWGCLR '
621 diagTitle =
622 . 'Clear Sky Net Downward SW Radiation at the Ground (W/m^2)'
623 diagUnits = 'watt/m2 '
624 diagCode = 'SM L1 '
625 CALL DIAGNOSTICS_ADD2LIST( diagNum,
626 I diagName, diagCode, diagUnits, diagTitle, myThid )
627
628 diagName = 'PAVE '
629 diagTitle = 'Average Surface Pressure (hPa)'
630 diagUnits = 'hPa/day '
631 diagCode = 'SM L1 '
632 CALL DIAGNOSTICS_ADD2LIST( diagNum,
633 I diagName, diagCode, diagUnits, diagTitle, myThid )
634
635 diagName = 'DIABU '
636 diagTitle =
637 . 'Zonal Wind Tendency due to Diabatic Forcing (m/s/day)'
638 diagUnits = 'm/s/day '
639 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
640 CALL DIAGNOSTICS_ADD2LIST( diagNum,
641 I diagName, diagCode, diagUnits, diagTitle, myThid )
642
643 diagName = 'DIABV '
644 diagTitle =
645 . 'Meridional Wind Tendency due to Diabatic Forcing (m/s/day)'
646 diagUnits = 'm/s/day '
647 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
648 CALL DIAGNOSTICS_ADD2LIST( diagNum,
649 I diagName, diagCode, diagUnits, diagTitle, myThid )
650
651 diagName = 'DIABT '
652 diagTitle =
653 . 'Temperature Tendency due to Diabatic Forcing (K/day)'
654 diagUnits = 'K/day '
655 diagCode = 'SM ML '
656 CALL DIAGNOSTICS_ADD2LIST( diagNum,
657 I diagName, diagCode, diagUnits, diagTitle, myThid )
658
659 diagName = 'DIABQ '
660 diagTitle =
661 . 'Specific Humidity Tendency due to Diabatic Forcing (g/kg/day)'
662 diagUnits = 'gram/kg/day '
663 diagCode = 'SM ML '
664 CALL DIAGNOSTICS_ADD2LIST( diagNum,
665 I diagName, diagCode, diagUnits, diagTitle, myThid )
666
667 diagName = 'RFU '
668 diagTitle =
669 . 'Zonal Wind Tendency due to Rayleigh Friction (m/s/day)'
670 diagUnits = 'm/s/day '
671 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
672 CALL DIAGNOSTICS_ADD2LIST( diagNum,
673 I diagName, diagCode, diagUnits, diagTitle, myThid )
674
675 diagName = 'RFV '
676 diagTitle =
677 . 'Meridional Wind Tendency due to Rayleigh Friction (m/s/day)'
678 diagUnits = 'm/s/day '
679 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
680 CALL DIAGNOSTICS_ADD2LIST( diagNum,
681 I diagName, diagCode, diagUnits, diagTitle, myThid )
682
683 diagName = 'GWDU '
684 diagTitle =
685 . 'Zonal Wind Tendency due to Gravity Wave Drag (m/s/day)'
686 diagUnits = 'm/s/day '
687 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
688 CALL DIAGNOSTICS_ADD2LIST( diagNum,
689 I diagName, diagCode, diagUnits, diagTitle, myThid )
690
691 diagName = 'GWDV '
692 diagTitle =
693 . 'Meridional Wind Tendency due to Gravity Wave Drag (m/s/day)'
694 diagUnits = 'm/s/day '
695 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
696 CALL DIAGNOSTICS_ADD2LIST( diagNum,
697 I diagName, diagCode, diagUnits, diagTitle, myThid )
698
699 diagName = 'GWDUS '
700 diagTitle = 'Zonal Wind Gravity Wave Surface Stress (N/m^2)'
701 diagUnits = 'N/m2 '
702 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'L1 '
703 CALL DIAGNOSTICS_ADD2LIST( diagNum,
704 I diagName, diagCode, diagUnits, diagTitle, myThid )
705
706 diagName = 'GWDVS '
707 diagTitle = 'Meridional Wind Gravity Wave Surface Stress (N/m^2)'
708 diagUnits = 'N/m2 '
709 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'L1 '
710 CALL DIAGNOSTICS_ADD2LIST( diagNum,
711 I diagName, diagCode, diagUnits, diagTitle, myThid )
712
713 diagName = 'GWDUT '
714 diagTitle = 'Zonal Wind Gravity Wave Stress at Model Top (N/m^2)'
715 diagUnits = 'N/m2 '
716 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'U1 '
717 CALL DIAGNOSTICS_ADD2LIST( diagNum,
718 I diagName, diagCode, diagUnits, diagTitle, myThid )
719
720 diagName = 'GWDVT '
721 diagTitle =
722 . 'Meridional Wind Gravity Wave Stress at Model Top (N/m^2)'
723 diagUnits = 'N/m2 '
724 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'U1 '
725 CALL DIAGNOSTICS_ADD2LIST( diagNum,
726 I diagName, diagCode, diagUnits, diagTitle, myThid )
727
728 diagName = 'LZRAD '
729 diagTitle =
730 . 'Estimated Cloud Liquid Water used in Radiation (mg/kg)'
731 diagUnits = 'mg/kg '
732 diagCode = 'SM ML '
733 CALL DIAGNOSTICS_ADD2LIST( diagNum,
734 I diagName, diagCode, diagUnits, diagTitle, myThid )
735
736 diagName = 'SLP '
737 diagTitle = 'Time-Averaged Sea-Level Pressure (mb)'
738 diagUnits = 'hPa '
739 diagCode = 'SM L1 '
740 CALL DIAGNOSTICS_ADD2LIST( diagNum,
741 I diagName, diagCode, diagUnits, diagTitle, myThid )
742
743 diagName = 'VAVEUQ '
744 diagTitle = 'Vertically Averaged UWND*SPHU (m/s g/kg)'
745 diagUnits = 'm/s gram/kg '
746 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'M1 '
747 CALL DIAGNOSTICS_ADD2LIST( diagNum,
748 I diagName, diagCode, diagUnits, diagTitle, myThid )
749
750 diagName = 'VAVEVQ '
751 diagTitle = 'Vertically Averaged VWND*SPHU (m/s g/kg)'
752 diagUnits = 'm/s gram/kg '
753 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'M1 '
754 CALL DIAGNOSTICS_ADD2LIST( diagNum,
755 I diagName, diagCode, diagUnits, diagTitle, myThid )
756
757 diagName = 'VAVEUT '
758 diagTitle = 'Vertically Averaged UWND*TMPU (m/s deg K)'
759 diagUnits = 'm/s K '
760 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'M1 '
761 CALL DIAGNOSTICS_ADD2LIST( diagNum,
762 I diagName, diagCode, diagUnits, diagTitle, myThid )
763
764 diagName = 'VAVEVT '
765 diagTitle = 'Vertically Averaged VWND*TMPU (m/s deg K)'
766 diagUnits = 'm/s K '
767 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'M1 '
768 CALL DIAGNOSTICS_ADD2LIST( diagNum,
769 I diagName, diagCode, diagUnits, diagTitle, myThid )
770
771 diagName = 'CLDFRC '
772 diagTitle = '2-D Total Cloud Fraction (0-1)'
773 diagUnits = '0-1 '
774 diagCode = 'SM P M1 '
775 CALL DIAGNOSTICS_ADD2LIST( diagNum,
776 I diagName, diagCode, diagUnits, diagTitle, myThid )
777
778 diagName = 'TPW '
779 diagTitle = 'Total Precipitable Water (g/cm^2)'
780 diagUnits = 'gram/cm2 '
781 diagCode = 'SM P M1 '
782 CALL DIAGNOSTICS_ADD2LIST( diagNum,
783 I diagName, diagCode, diagUnits, diagTitle, myThid )
784
785 diagName = 'U2M '
786 diagTitle = 'Zonal Wind Interpolated to 2 Meters (m/s)'
787 diagUnits = 'm/sec '
788 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'M1 '
789 CALL DIAGNOSTICS_ADD2LIST( diagNum,
790 I diagName, diagCode, diagUnits, diagTitle, myThid )
791
792 diagName = 'V2M '
793 diagTitle = 'Meridional Wind Interpolated to 2 Meters (m/s)'
794 diagUnits = 'm/sec '
795 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'M1 '
796 CALL DIAGNOSTICS_ADD2LIST( diagNum,
797 I diagName, diagCode, diagUnits, diagTitle, myThid )
798
799 diagName = 'T2M '
800 diagTitle = 'Temperature Interpolated to 2 Meters (deg K)'
801 diagUnits = 'K '
802 diagCode = 'SM M1 '
803 CALL DIAGNOSTICS_ADD2LIST( diagNum,
804 I diagName, diagCode, diagUnits, diagTitle, myThid )
805
806 diagName = 'Q2M '
807 diagTitle = 'Specific Humidity Interpolated to 2 Meters (g/kg)'
808 diagUnits = 'gram/kg '
809 diagCode = 'SM P M1 '
810 CALL DIAGNOSTICS_ADD2LIST( diagNum,
811 I diagName, diagCode, diagUnits, diagTitle, myThid )
812
813 diagName = 'U10M '
814 diagTitle = 'Zonal Wind Interpolated to 10 Meters (m/s)'
815 diagUnits = 'm/sec '
816 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'M1 '
817 CALL DIAGNOSTICS_ADD2LIST( diagNum,
818 I diagName, diagCode, diagUnits, diagTitle, myThid )
819
820 diagName = 'V10M '
821 diagTitle = 'Meridional Wind Interpolated to 10 Meters (m/s)'
822 diagUnits = 'm/sec '
823 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'M1 '
824 CALL DIAGNOSTICS_ADD2LIST( diagNum,
825 I diagName, diagCode, diagUnits, diagTitle, myThid )
826
827 diagName = 'T10M '
828 diagTitle = 'Temperature Interpolated to 10 Meters (deg K)'
829 diagUnits = 'K '
830 diagCode = 'SM M1 '
831 CALL DIAGNOSTICS_ADD2LIST( diagNum,
832 I diagName, diagCode, diagUnits, diagTitle, myThid )
833
834 diagName = 'Q10M '
835 diagTitle = 'Specific Humidity Interpolated to 10 Meters (g/kg)'
836 diagUnits = 'gram/kg '
837 diagCode = 'SM P M1 '
838 CALL DIAGNOSTICS_ADD2LIST( diagNum,
839 I diagName, diagCode, diagUnits, diagTitle, myThid )
840
841 diagName = 'DTRAIN '
842 diagTitle = 'Detrainment Cloud Mass Flux (g/m^2/sec)'
843 diagUnits = 'g/m2/sec '
844 diagCode = 'SM P LL '
845 CALL DIAGNOSTICS_ADD2LIST( diagNum,
846 I diagName, diagCode, diagUnits, diagTitle, myThid )
847
848 diagName = 'QFILL '
849 diagTitle = 'Filling of Negative Specific Humidity (g/kg/day)'
850 diagUnits = 'gram/kg/day '
851 diagCode = 'SM ML '
852 CALL DIAGNOSTICS_ADD2LIST( diagNum,
853 I diagName, diagCode, diagUnits, diagTitle, myThid )
854
855 diagName = 'VAVEQAVE'
856 diagTitle = 'Vertically Averaged QAVE (mm/day)'
857 diagUnits = 'mm/day '
858 diagCode = 'SM M1 '
859 CALL DIAGNOSTICS_ADD2LIST( diagNum,
860 I diagName, diagCode, diagUnits, diagTitle, myThid )
861
862 diagName = 'VAVEQFIL'
863 diagTitle = 'Vertically Averaged QFILL (mm/day)'
864 diagUnits = 'mm/day '
865 diagCode = 'SM P L1 '
866 CALL DIAGNOSTICS_ADD2LIST( diagNum,
867 I diagName, diagCode, diagUnits, diagTitle, myThid )
868
869 diagName = 'CLDMID '
870 diagTitle = 'Mid-Level (700-400 hPa) Cloud Fraction (0-1)'
871 diagUnits = '0-1 '
872 diagCode = 'SM P M1 '
873 CALL DIAGNOSTICS_ADD2LIST( diagNum,
874 I diagName, diagCode, diagUnits, diagTitle, myThid )
875
876 diagName = 'CLDHI '
877 diagTitle = 'High-Level (above 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 = 'TAULOW '
884 diagTitle = 'Low-Level (1000-700 hPa) Optical Depth'
885 diagUnits = 'non-dimensional '
886 WRITE(diagCode,'(A,I3.3,A)') 'SM PC',diagNum+2,'M1 '
887 CALL DIAGNOSTICS_ADD2LIST( diagNum,
888 I diagName, diagCode, diagUnits, diagTitle, myThid )
889
890 diagName = 'TAULOWC '
891 diagTitle = 'Low-Level Optical Depth Counter'
892 diagUnits = 'count '
893 diagCode = 'SM D M1 '
894 CALL DIAGNOSTICS_ADD2LIST( diagNum,
895 I diagName, diagCode, diagUnits, diagTitle, myThid )
896
897 diagName = 'TAUMID '
898 diagTitle = 'Mid-Level (700-400 hPa) Optical Depth'
899 diagUnits = 'non-dimensional '
900 WRITE(diagCode,'(A,I3.3,A)') 'SM PC',diagNum+2,'M1 '
901 CALL DIAGNOSTICS_ADD2LIST( diagNum,
902 I diagName, diagCode, diagUnits, diagTitle, myThid )
903
904 diagName = 'TAUMIDC '
905 diagTitle = 'Mid-Level Optical Depth Counter'
906 diagUnits = 'count '
907 diagCode = 'SM D M1 '
908 CALL DIAGNOSTICS_ADD2LIST( diagNum,
909 I diagName, diagCode, diagUnits, diagTitle, myThid )
910
911 diagName = 'TAUHI '
912 diagTitle = 'High-Level (above 400 hPa) Optical Depth'
913 diagUnits = 'non-dimensional '
914 WRITE(diagCode,'(A,I3.3,A)') 'SM PC',diagNum+2,'M1 '
915 CALL DIAGNOSTICS_ADD2LIST( diagNum,
916 I diagName, diagCode, diagUnits, diagTitle, myThid )
917
918 diagName = 'TAUHIC '
919 diagTitle = 'High-Level Optical Depth Counter'
920 diagUnits = 'count '
921 diagCode = 'SM D M1 '
922 CALL DIAGNOSTICS_ADD2LIST( diagNum,
923 I diagName, diagCode, diagUnits, diagTitle, myThid )
924
925 diagName = 'CLDNP '
926 diagTitle = 'Non-Precipitating Cloud Fraction (0-1)'
927 diagUnits = '0-1 '
928 diagCode = 'SM ML '
929 CALL DIAGNOSTICS_ADD2LIST( diagNum,
930 I diagName, diagCode, diagUnits, diagTitle, myThid )
931
932 diagName = 'CLDPRS '
933 diagTitle = 'Cloud Top Pressure (when cloudy) (hPa)'
934 diagUnits = 'hPa '
935 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 = 'CTPCNT '
940 diagTitle = 'Cloud Top Pressure Counter'
941 diagUnits = 'count '
942 diagCode = 'SM D M1 '
943 CALL DIAGNOSTICS_ADD2LIST( diagNum,
944 I diagName, diagCode, diagUnits, diagTitle, myThid )
945
946 diagName = 'CLDTMP '
947 diagTitle = 'Cloud Top Temperature (when cloudy) (deg K)'
948 diagUnits = 'K '
949 WRITE(diagCode,'(A,I3.3,A)') 'SM PC',diagNum+2,'M1 '
950 CALL DIAGNOSTICS_ADD2LIST( diagNum,
951 I diagName, diagCode, diagUnits, diagTitle, myThid )
952
953 diagName = 'CTTCNT '
954 diagTitle = 'Cloud Top Temperature Counter'
955 diagUnits = 'count '
956 diagCode = 'SM D M1 '
957 CALL DIAGNOSTICS_ADD2LIST( diagNum,
958 I diagName, diagCode, diagUnits, diagTitle, myThid )
959
960 diagName = 'TCANOPY '
961 diagTitle =
962 . 'Skin Temperature of the Canopy/Soil (SST over water) (deg K)'
963 diagUnits = 'K '
964 diagCode = 'SM L1 '
965 CALL DIAGNOSTICS_ADD2LIST( diagNum,
966 I diagName, diagCode, diagUnits, diagTitle, myThid )
967
968 diagName = 'TDEEP '
969 diagTitle = 'Deep Soil Temperature (deg K)'
970 diagUnits = 'K '
971 diagCode = 'SM L1 '
972 CALL DIAGNOSTICS_ADD2LIST( diagNum,
973 I diagName, diagCode, diagUnits, diagTitle, myThid )
974
975 diagName = 'QCANOPY '
976 diagTitle = 'Specific Humidity of the Canopy (g/kg)'
977 diagUnits = 'gram/kg '
978 diagCode = 'SM L1 '
979 CALL DIAGNOSTICS_ADD2LIST( diagNum,
980 I diagName, diagCode, diagUnits, diagTitle, myThid )
981
982 diagName = 'SMSHAL '
983 diagTitle =
984 . 'Soil Moisture Field Capacity Fraction in Shallow (5 mm) Layer'
985 diagUnits = '0-1 '
986 diagCode = 'SM L1 '
987 CALL DIAGNOSTICS_ADD2LIST( diagNum,
988 I diagName, diagCode, diagUnits, diagTitle, myThid )
989
990 diagName = 'SMROOT '
991 diagTitle = 'Soil Moisture Field Capacity Fraction in Root Zone'
992 diagUnits = '0-1 '
993 diagCode = 'SM L1 '
994 CALL DIAGNOSTICS_ADD2LIST( diagNum,
995 I diagName, diagCode, diagUnits, diagTitle, myThid )
996
997 diagName = 'SMDEEP '
998 diagTitle =
999 . 'Soil Moisture Field Capacity Fraction in Recharge Layer '
1000 diagUnits = '0-1 '
1001 diagCode = 'SM L1 '
1002 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1003 I diagName, diagCode, diagUnits, diagTitle, myThid )
1004
1005 diagName = 'CAPACITY'
1006 diagTitle = 'Canopy Interception Reservoir (mm)'
1007 diagUnits = 'mm '
1008 diagCode = 'SM L1 '
1009 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1010 I diagName, diagCode, diagUnits, diagTitle, myThid )
1011
1012 diagName = 'SNOW '
1013 diagTitle = 'Snow Depth (mm water equivalent)'
1014 diagUnits = 'mm water equiv '
1015 diagCode = 'SM L1 '
1016 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1017 I diagName, diagCode, diagUnits, diagTitle, myThid )
1018
1019 diagName = 'RAINCON '
1020 diagTitle = 'Convective Rainfall (liquid precipitate) (mm/day)'
1021 diagUnits = 'mm/day '
1022 diagCode = 'SM L1 '
1023 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1024 I diagName, diagCode, diagUnits, diagTitle, myThid )
1025
1026 diagName = 'RAINLSP '
1027 diagTitle = 'Large-scale 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 = 'SNOWFALL'
1034 diagTitle = 'Total Snowfall (solid 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 = 'RUNOFF '
1041 diagTitle =
1042 . 'Water from precipitation not infiltrated into soil (mm/s)'
1043 diagUnits = 'mm/sec '
1044 diagCode = 'SM L1 '
1045 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1046 I diagName, diagCode, diagUnits, diagTitle, myThid )
1047
1048 diagName = 'FWSOIL '
1049 diagTitle = 'Rainwater Infiltration into top soil layer (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 = 'GDRAIN '
1056 diagTitle =
1057 . 'Diffusion of Moisture across bottom of root zone (mm/s)'
1058 diagUnits = 'mm/sec '
1059 diagCode = 'SM L1 '
1060 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1061 I diagName, diagCode, diagUnits, diagTitle, myThid )
1062
1063 diagName = 'SNOWMELT'
1064 diagTitle = 'Rate of Snow Melt (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 = 'ERESV '
1071 diagTitle =
1072 . 'Latent Heat Flux from Canopy Interception Reservoir (W/m^2)'
1073 diagUnits = 'watt/m2 '
1074 diagCode = 'SM L1 '
1075 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1076 I diagName, diagCode, diagUnits, diagTitle, myThid )
1077
1078 diagName = 'ESOIL '
1079 diagTitle = 'Latent Heat Flux from Bare Soil (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 = 'EVEG '
1086 diagTitle =
1087 . 'Latent Heat Flux (Transpiration) from Vegetated Surface (W/m^2)'
1088 diagUnits = 'watt/m2 '
1089 diagCode = 'SM L1 '
1090 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1091 I diagName, diagCode, diagUnits, diagTitle, myThid )
1092
1093 diagName = 'ESNOW '
1094 diagTitle = 'Latent Heat Flux from Snow Pack (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 = 'PARDF '
1101 diagTitle =
1102 . 'Diffuse-beam Photosynthetically Active Radiation (W/m^2)'
1103 diagUnits = 'watt/m2 '
1104 diagCode = 'SM P L1 '
1105 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1106 I diagName, diagCode, diagUnits, diagTitle, myThid )
1107
1108 diagName = 'PARDR '
1109 diagTitle =
1110 . 'Direct-beam Photosynthetically Active Radiation (W/m^2)'
1111 diagUnits = 'watt/m2 '
1112 diagCode = 'SM P L1 '
1113 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1114 I diagName, diagCode, diagUnits, diagTitle, myThid )
1115
1116 diagName = 'LAI '
1117 diagTitle = 'Leaf Area Index (%)'
1118 diagUnits = 'percent '
1119 diagCode = 'SM L1 '
1120 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1121 I diagName, diagCode, diagUnits, diagTitle, myThid )
1122
1123 diagName = 'GREEN '
1124 diagTitle = 'Green-ness Index (%)'
1125 diagUnits = 'percent '
1126 diagCode = 'SM L1 '
1127 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1128 I diagName, diagCode, diagUnits, diagTitle, myThid )
1129
1130 diagName = 'DLWDTC '
1131 diagTitle =
1132 . 'Derivative of Upward LW Radiation wrt TCANOPY (W/m^2 /deg K)'
1133 diagUnits = 'watt/m2/K '
1134 diagCode = 'SM L1 '
1135 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1136 I diagName, diagCode, diagUnits, diagTitle, myThid )
1137
1138 diagName = 'DHDTC '
1139 diagTitle =
1140 . 'Derivative of Sensible Heat Flux wrt TCANOPY (W/m^2 /deg K)'
1141 diagUnits = 'watt/m2/K '
1142 diagCode = 'SM L1 '
1143 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1144 I diagName, diagCode, diagUnits, diagTitle, myThid )
1145
1146 diagName = 'DEDTC '
1147 diagTitle =
1148 . 'Derivative of Latent Heat Flux wrt TCANOPY (W/m^2 /deg K)'
1149 diagUnits = 'watt/m2/K '
1150 diagCode = 'SM L1 '
1151 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1152 I diagName, diagCode, diagUnits, diagTitle, myThid )
1153
1154 diagName = 'VDTMOIST'
1155 diagTitle = 'Vertically Averaged MOISTT (K/day)'
1156 diagUnits = 'K/day '
1157 diagCode = 'SM M1 '
1158 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1159 I diagName, diagCode, diagUnits, diagTitle, myThid )
1160
1161 diagName = 'VDTTURB '
1162 diagTitle = 'Vertically Averaged TURBT (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 = 'VDTLWRAD'
1169 diagTitle = 'Vertically Averaged RADLW (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 = 'VDTSWRAD'
1176 diagTitle = 'Vertically Averaged RADSW (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 = 'VAVETAVE'
1183 diagTitle = 'Vertically Averaged TAVE (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 = 'PSUBCLD '
1190 diagTitle = 'Pressure Thickness of the RAS Subcloud Layer (mb)'
1191 diagUnits = 'hPa '
1192 WRITE(diagCode,'(A,I3.3,A)') 'SM PC',diagNum+2,'M1 '
1193 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1194 I diagName, diagCode, diagUnits, diagTitle, myThid )
1195
1196 diagName = 'PSUBCLDC'
1197 diagTitle =
1198 . 'Pressure Thickness of the RAS Subcloud Layer Counter'
1199 diagUnits = 'count '
1200 diagCode = 'SM PD M1 '
1201 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1202 I diagName, diagCode, diagUnits, diagTitle, myThid )
1203
1204 diagName = 'LCL '
1205 diagTitle = 'Lifting Condensation Level (mb) '
1206 diagUnits = 'hPa '
1207 diagCode = 'SM L1 '
1208 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1209 I diagName, diagCode, diagUnits, diagTitle, myThid )
1210
1211 diagName = 'EVPOT '
1212 diagTitle = 'Potential Evapotranspiration (W/m^2)'
1213 diagUnits = 'watt/m2 '
1214 diagCode = 'SM L1 '
1215 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1216 I diagName, diagCode, diagUnits, diagTitle, myThid )
1217
1218 diagName = 'RHCHECK '
1219 diagTitle =
1220 . 'Net precip due to supersaturation above critical RH (mm/day)'
1221 diagUnits = 'mm/day '
1222 diagCode = 'SM L1 '
1223 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1224 I diagName, diagCode, diagUnits, diagTitle, myThid )
1225
1226 diagName = 'DHDQA '
1227 diagTitle =
1228 . 'Derivative of Sensible Heat wrt Canopy Humidity (W/m^2)'
1229 diagUnits = 'watt/m2 '
1230 diagCode = 'SM L1 '
1231 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1232 I diagName, diagCode, diagUnits, diagTitle, myThid )
1233
1234 diagName = 'DEDQA '
1235 diagTitle =
1236 . 'Derivative of Latent Heat wrt Canopy Humidity (W/m^2)'
1237 diagUnits = 'watt/m2 '
1238 diagCode = 'SM L1 '
1239 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1240 I diagName, diagCode, diagUnits, diagTitle, myThid )
1241
1242 diagName = 'DTC '
1243 diagTitle = 'Total Change in Canopy Temperature (deg/sec)'
1244 diagUnits = 'K/sec '
1245 diagCode = 'SM L1 '
1246 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1247 I diagName, diagCode, diagUnits, diagTitle, myThid )
1248
1249 diagName = 'DQC '
1250 diagTitle = 'Change in Canopy Humidity (kg/kg/sec)'
1251 diagUnits = 'sec-1 '
1252 diagCode = 'SM L1 '
1253 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1254 I diagName, diagCode, diagUnits, diagTitle, myThid )
1255
1256 diagName = 'TCDTC '
1257 diagTitle =
1258 . 'Change in Canopy Temperature - Energy Term (deg/sec)'
1259 diagUnits = 'K/sec '
1260 diagCode = 'SM L1 '
1261 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1262 I diagName, diagCode, diagUnits, diagTitle, myThid )
1263
1264 diagName = 'RADDTC '
1265 diagTitle =
1266 . 'Change in Canopy Temperature - Net Radiation Term (deg/sec)'
1267 diagUnits = 'K/sec '
1268 diagCode = 'SM L1 '
1269 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1270 I diagName, diagCode, diagUnits, diagTitle, myThid )
1271
1272 diagName = 'SENDDTC '
1273 diagTitle =
1274 . 'Change in Canopy Temperature - Sensible Heat Term (deg/sec)'
1275 diagUnits = 'K/sec '
1276 diagCode = 'SM L1 '
1277 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1278 I diagName, diagCode, diagUnits, diagTitle, myThid )
1279
1280 diagName = 'LATDTC '
1281 diagTitle =
1282 . 'Change in Canopy Temperature - Latent Heat Term (deg/sec)'
1283 diagUnits = 'K/sec '
1284 diagCode = 'SM L1 '
1285 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1286 I diagName, diagCode, diagUnits, diagTitle, myThid )
1287
1288 diagName = 'TDDTC '
1289 diagTitle =
1290 . 'Change in Canopy Temperature - Ground Temp Term (deg/sec)'
1291 diagUnits = 'K/sec '
1292 diagCode = 'SM L1 '
1293 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1294 I diagName, diagCode, diagUnits, diagTitle, myThid )
1295
1296 diagName = 'QCDTC '
1297 diagTitle =
1298 . 'Change in Canopy Temperature - Humidity Term (deg/sec)'
1299 diagUnits = 'K/sec '
1300 diagCode = 'SM L1 '
1301 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1302 I diagName, diagCode, diagUnits, diagTitle, myThid )
1303
1304 diagName = 'TRBQLIQ '
1305 diagTitle = 'Cloud Liquid Water Generated in Turbulence (mg/kg)'
1306 diagUnits = 'mg/kg '
1307 diagCode = 'SM ML '
1308 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1309 I diagName, diagCode, diagUnits, diagTitle, myThid )
1310
1311 diagName = 'TRBFCC '
1312 diagTitle = 'Cloud Fraction Generated in Turbulence (fraction)'
1313 diagUnits = '0-1 '
1314 diagCode = 'SM ML '
1315 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1316 I diagName, diagCode, diagUnits, diagTitle, myThid )
1317
1318 diagName = 'ALBEDO '
1319 diagTitle = 'Mean Surface Albedo (SW_Up/SW_Down) 0-1 '
1320 diagUnits = '0-1 '
1321 WRITE(diagCode,'(A,I3.3,A)') 'SM C',diagNum+2,'L1 '
1322 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1323 I diagName, diagCode, diagUnits, diagTitle, myThid )
1324
1325 diagName = 'ALBEDOC '
1326 diagTitle = 'Surface Albedo Counter'
1327 diagUnits = 'count '
1328 diagCode = 'SM D L1 '
1329 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1330 I diagName, diagCode, diagUnits, diagTitle, myThid )
1331
1332 diagName = 'ALBVISDR'
1333 diagTitle = 'Direct Beam VIS Surface Albedo (0-1) '
1334 diagUnits = '0-1 '
1335 diagCode = 'SM L1 '
1336 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1337 I diagName, diagCode, diagUnits, diagTitle, myThid )
1338
1339 diagName = 'ALBVISDF'
1340 diagTitle = 'Diffuse 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 = 'ALBNIRDR'
1347 diagTitle = 'Direct Beam NIR 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 = 'ALBNIRDF'
1354 diagTitle = 'Diffuse 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 = 'TAUAVE '
1361 diagTitle =
1362 . 'Averaged Cloud Optical Depth (non-dimensional) per 100 mb'
1363 diagUnits = 'non-dimensional '
1364 diagCode = 'SM ML '
1365 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1366 I diagName, diagCode, diagUnits, diagTitle, myThid )
1367
1368 diagName = 'RHCHECKT'
1369 diagTitle =
1370 . 'Temperature Change due to humidity above critical RH (K/day)'
1371 diagUnits = 'K/day '
1372 diagCode = 'SM ML '
1373 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1374 I diagName, diagCode, diagUnits, diagTitle, myThid )
1375
1376 diagName = 'RHCHECKQ'
1377 diagTitle =
1378 . 'Spec.Humid. Change due to humidity above critical RH (g/kg/day)'
1379 diagUnits = 'g/kg/day '
1380 diagCode = 'SM ML '
1381 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1382 I diagName, diagCode, diagUnits, diagTitle, myThid )
1383
1384 diagName = 'CLDLSP '
1385 diagTitle = 'Large-Scale Cloud Fraction (0-1)'
1386 diagUnits = '0-1 '
1387 diagCode = 'SM P ML '
1388 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1389 I diagName, diagCode, diagUnits, diagTitle, myThid )
1390
1391 diagName = 'LZLOW '
1392 diagTitle =
1393 . 'Vert integrated low-level (below 700mb) cloud water (mg/cm**2)'
1394 diagUnits = 'mg/cm2 '
1395 diagCode = 'SM P L1 '
1396 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1397 I diagName, diagCode, diagUnits, diagTitle, myThid )
1398
1399 diagName = 'LZMID '
1400 diagTitle =
1401 . 'Vert integrated mid-level (700-400mb) cloud water (mg/cm**2)'
1402 diagUnits = 'mg/cm2 '
1403 diagCode = 'SM P L1 '
1404 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1405 I diagName, diagCode, diagUnits, diagTitle, myThid )
1406
1407 diagName = 'LZHI '
1408 diagTitle =
1409 . 'Vert integrated high-level (above 400mb) cloud water (mg/cm**2)'
1410 diagUnits = 'mg/cm2 '
1411 diagCode = 'SM P L1 '
1412 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1413 I diagName, diagCode, diagUnits, diagTitle, myThid )
1414
1415 diagName = 'LZTOT '
1416 diagTitle =
1417 . 'Vertically integrated total suspended cloud water (mg/cm**2)'
1418 diagUnits = 'mg/cm2 '
1419 diagCode = 'SM P L1 '
1420 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1421 I diagName, diagCode, diagUnits, diagTitle, myThid )
1422
1423 diagName = 'CAPE '
1424 diagTitle = 'Convective Available Potential Energy'
1425 diagUnits = 'mb m2/sec2 '
1426 diagCode = 'SM P L1 '
1427 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1428 I diagName, diagCode, diagUnits, diagTitle, myThid )
1429
1430 diagName = 'HSUBCL '
1431 diagTitle = 'Moist Static Energy in RAS subcloud layer'
1432 diagUnits = 'm2/sec2 '
1433 diagCode = 'SM L1 '
1434 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1435 I diagName, diagCode, diagUnits, diagTitle, myThid )
1436
1437 diagName = 'HSTAR '
1438 diagTitle = 'Saturation Moist Static Energy'
1439 diagUnits = 'm2/sec2 '
1440 diagCode = 'SM ML '
1441 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1442 I diagName, diagCode, diagUnits, diagTitle, myThid )
1443
1444 diagName = 'CONT '
1445 diagTitle =
1446 . 'Temperature Tendency due to Convective Precipitation (K/day)'
1447 diagUnits = 'K/day '
1448 diagCode = 'SM ML '
1449 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1450 I diagName, diagCode, diagUnits, diagTitle, myThid )
1451
1452 diagName = 'CONQ '
1453 diagTitle =
1454 . 'Moisture Tendency due to Convective Precipitation (g/kg/day)'
1455 diagUnits = 'g/kg/day '
1456 diagCode = 'SM ML '
1457 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1458 I diagName, diagCode, diagUnits, diagTitle, myThid )
1459
1460 diagName = 'LST '
1461 diagTitle =
1462 . 'Temperature Tendency due to Large Scale Precipitation (K/day)'
1463 diagUnits = 'K/day '
1464 diagCode = 'SM ML '
1465 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1466 I diagName, diagCode, diagUnits, diagTitle, myThid )
1467
1468 diagName = 'LSQ '
1469 diagTitle =
1470 . 'Moisture Tendency due to Large Scale Precipitation (g/kg/day)'
1471 diagUnits = 'g/kg/day '
1472 diagCode = 'SM ML '
1473 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1474 I diagName, diagCode, diagUnits, diagTitle, myThid )
1475
1476 diagName = 'CONEVPT '
1477 diagTitle =
1478 . 'Temperature Tendency due to Conv Rain Re-evaporation (K/day)'
1479 diagUnits = 'K/day '
1480 diagCode = 'SM ML '
1481 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1482 I diagName, diagCode, diagUnits, diagTitle, myThid )
1483
1484 diagName = 'CONEVPQ '
1485 diagTitle =
1486 . 'Moisture Tendency due to Conv Rain Re-evaporation (g/kg/day)'
1487 diagUnits = 'g/kg/day '
1488 diagCode = 'SM ML '
1489 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1490 I diagName, diagCode, diagUnits, diagTitle, myThid )
1491
1492 diagName = 'LSEVPT '
1493 diagTitle =
1494 . 'Temperature Tendency due to Large Scale Rain Re-evap (K/day)'
1495 diagUnits = 'K/day '
1496 diagCode = 'SM ML '
1497 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1498 I diagName, diagCode, diagUnits, diagTitle, myThid )
1499
1500 diagName = 'LSEVPQ '
1501 diagTitle =
1502 . 'Moisture Tendency due to Large Scale Rain Re-evap (g/kg/day)'
1503 diagUnits = 'g/kg/day '
1504 diagCode = 'SM ML '
1505 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1506 I diagName, diagCode, diagUnits, diagTitle, myThid )
1507
1508 diagName = 'GWCUMID '
1509 diagTitle =
1510 . 'Mid-Lev Mean U-Wind for Gravity Wave Drag w/Clouds (m/sec)'
1511 diagUnits = 'm/sec '
1512 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'L1 '
1513 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1514 I diagName, diagCode, diagUnits, diagTitle, myThid )
1515
1516 diagName = 'GWCVMID '
1517 diagTitle =
1518 . 'Mid-Lev Mean V-Wind for Gravity Wave Drag w/Clouds (m/sec)'
1519 diagUnits = 'm/sec '
1520 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'L1 '
1521 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1522 I diagName, diagCode, diagUnits, diagTitle, myThid )
1523
1524 diagName = 'CLDSTD '
1525 diagTitle =
1526 . 'Mid-Lev Cloud Perturb Length for Gravity Wave Drag w/Clouds (m)'
1527 diagUnits = 'm '
1528 diagCode = 'SM L1 '
1529 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1530 I diagName, diagCode, diagUnits, diagTitle, myThid )
1531
1532 diagName = 'GWCUBAR '
1533 diagTitle =
1534 . 'Mid-Lev Perturb U-Wind for Gravity Wave Drag w/Clouds (m/sec)'
1535 diagUnits = 'm/sec '
1536 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'L1 '
1537 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1538 I diagName, diagCode, diagUnits, diagTitle, myThid )
1539
1540 diagName = 'GWCVBAR '
1541 diagTitle =
1542 . 'Mid-Lev Perturb V-Wind for Gravity Wave Drag w/Clouds (m/sec)'
1543 diagUnits = 'm/sec '
1544 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'L1 '
1545 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1546 I diagName, diagCode, diagUnits, diagTitle, myThid )
1547
1548 diagName = 'GWCUS '
1549 diagTitle =
1550 . 'Mid-Lev U-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)'
1551 diagUnits = 'N/m2 '
1552 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'L1 '
1553 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1554 I diagName, diagCode, diagUnits, diagTitle, myThid )
1555
1556 diagName = 'GWCVS '
1557 diagTitle =
1558 . 'Mid-Lev V-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)'
1559 diagUnits = 'N/m2 '
1560 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'L1 '
1561 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1562 I diagName, diagCode, diagUnits, diagTitle, myThid )
1563
1564 diagName = 'GWCU '
1565 diagTitle =
1566 . 'U-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
1567 diagUnits = 'm/sec/day '
1568 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
1569 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1570 I diagName, diagCode, diagUnits, diagTitle, myThid )
1571
1572 diagName = 'GWCV '
1573 diagTitle =
1574 . 'V-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
1575 diagUnits = 'm/sec/day '
1576 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
1577 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1578 I diagName, diagCode, diagUnits, diagTitle, myThid )
1579
1580 diagName = 'PLALBEDO'
1581 diagTitle = 'Planetary Albedo (Fraction)'
1582 diagUnits = 'Fraction '
1583 diagCode = 'SM L1 '
1584 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1585 I diagName, diagCode, diagUnits, diagTitle, myThid )
1586
1587 diagName = 'FIZPRES '
1588 diagTitle = 'Pressures for Fizhi (mb)'
1589 diagUnits = 'mb (hPa) '
1590 diagCode = 'SM ML '
1591 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1592 I diagName, diagCode, diagUnits, diagTitle, myThid )
1593
1594 diagName = 'CORRDU '
1595 diagTitle =
1596 . 'Zonal Wind Tendency Gridalt Correction Term (m/s/day)'
1597 diagUnits = 'm/s/day '
1598 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
1599 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1600 I diagName, diagCode, diagUnits, diagTitle, myThid )
1601
1602 diagName = 'CORRDV '
1603 diagTitle =
1604 . 'Meridional Wind Tendency Gridalt Correction Term (m/s/day)'
1605 diagUnits = 'm/s/day '
1606 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
1607 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1608 I diagName, diagCode, diagUnits, diagTitle, myThid )
1609
1610 diagName = 'CORRDT '
1611 diagTitle =
1612 . 'Temperature Tendency Gridalt Correction Term (K/day)'
1613 diagUnits = 'K/day '
1614 diagCode = 'SM ML '
1615 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1616 I diagName, diagCode, diagUnits, diagTitle, myThid )
1617
1618 diagName = 'CORRDQ '
1619 diagTitle =
1620 . 'Specific Humidity Tendency Gridalt Correction Term (g/kg/day)'
1621 diagUnits = 'gram/kg/day '
1622 diagCode = 'SM ML '
1623 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1624 I diagName, diagCode, diagUnits, diagTitle, myThid )
1625
1626 diagName = 'DIABUDYN'
1627 diagTitle =
1628 . 'Zonal Wind Tend due to Diab Forcing - Dyn Grid (m/s/day)'
1629 diagUnits = 'm/s/day '
1630 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'MR '
1631 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1632 I diagName, diagCode, diagUnits, diagTitle, myThid )
1633
1634 diagName = 'DIABVDYN'
1635 diagTitle =
1636 . 'Meridional Wind Tend due to Diab Forcing - Dyn Grid (m/s/day)'
1637 diagUnits = 'm/s/day '
1638 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'MR '
1639 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1640 I diagName, diagCode, diagUnits, diagTitle, myThid )
1641
1642 diagName = 'DIABTDYN'
1643 diagTitle =
1644 . 'Temperature Tend due to Diab Forcing - Dyn Grid (K/day)'
1645 diagUnits = 'K/day '
1646 diagCode = 'SM MR '
1647 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1648 I diagName, diagCode, diagUnits, diagTitle, myThid )
1649
1650 diagName = 'DIABQDYN'
1651 diagTitle =
1652 . 'Specific Hum Tend due to Diab Forcing - Dyn Grid (g/kg/day)'
1653 diagUnits = 'gram/kg/day '
1654 diagCode = 'SM MR '
1655 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1656 I diagName, diagCode, diagUnits, diagTitle, myThid )
1657
1658 diagName = 'TENDUFIZ'
1659 diagTitle =
1660 . 'Total Zonal Wind Tend - Fizhi Grid (m/s/day)'
1661 diagUnits = 'm/s/day '
1662 WRITE(diagCode,'(A,I3.3,A)') 'UM ',diagNum+2,'ML '
1663 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1664 I diagName, diagCode, diagUnits, diagTitle, myThid )
1665
1666 diagName = 'TENDVFIZ'
1667 diagTitle =
1668 . 'Total Meridional Wind Tend - Fizhi Grid (m/s/day)'
1669 diagUnits = 'm/s/day '
1670 WRITE(diagCode,'(A,I3.3,A)') 'VM ', diagNum ,'ML '
1671 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1672 I diagName, diagCode, diagUnits, diagTitle, myThid )
1673
1674 diagName = 'TENDTFIZ'
1675 diagTitle =
1676 . 'Total Temperature Tend - Fizhi Grid (K/day)'
1677 diagUnits = 'K/day '
1678 diagCode = 'SM ML '
1679 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1680 I diagName, diagCode, diagUnits, diagTitle, myThid )
1681
1682 diagName = 'TENDQFIZ'
1683 diagTitle =
1684 . 'Total Specific Hum Tend - Fizhi Grid (g/kg/day)'
1685 diagUnits = 'gram/kg/day '
1686 diagCode = 'SM ML '
1687 CALL DIAGNOSTICS_ADD2LIST( diagNum,
1688 I diagName, diagCode, diagUnits, diagTitle, myThid )
1689
1690 #endif /* ALLOW_DIAGNOSTICS */
1691
1692 RETURN
1693 END

  ViewVC Help
Powered by ViewVC 1.1.22