subroutine diagnostics_init_vals (myThid) C *********************************************************************** C ***** **** C ***** PURPOSE **** C ***** Initialize Diagnostic Names, Levels, and GPARMS **** C ***** **** C ***** **** C ***** Note: GDIAG is defined as character*16 and can be **** C ***** equivalenced to character*1 parse(16) **** C ***** **** C ***** parse(1) => 'S' for a scalar diagnostic **** C ***** => 'U' for a U-vector component diagnostic **** C ***** => 'V' for a V-vector component diagnostic **** C ***** **** C ***** parse(2) => 'U' for a C-grid U-Point **** C ***** => 'V' for a C-grid V-Point **** C ***** => 'M' for a C-grid Mass Point **** C ***** => 'Z' for a C-grid Corner Point **** C ***** **** C ***** parse(3) => not used at present **** C ***** **** C ***** parse(4) => 'P' if diagnostic is positive definite **** C ***** **** C ***** parse(5) => 'C' for a diagnostic with a counter array **** C ***** => 'D' to disable an array for output **** C ***** **** C ***** parse(6-8) => 3-digit integer corresponding to vector **** C ***** component mate if parse(1) = 'U' or 'V' **** C ***** or to the counter if parse(5)='C' **** C ***** **** C ***** parse(9) => 'U' if defined at Model-level plus 1/2 **** C ***** => 'M' if defined at Model-level middle **** C ***** => 'L' if defined at Model-level minus 1/2 **** C ***** **** C ***** **** C *********************************************************************** implicit none #include "CPP_OPTIONS.h" #include "SIZE.h" #ifdef ALLOW_FIZHI #include "fizhi_SIZE.h" #else integer Nrphys parameter (Nrphys=0) #endif #include "diagnostics_SIZE.h" #include "diagnostics.h" integer myThid integer n do n=1,ndiagt kdiag (n) = 0 cdiag(n) = ' ' gdiag(n) = ' ' udiag(n) = ' ' tdiag(n) = ' ' enddo c Diagnostic Names c ---------------- cdiag( 1) = 'UFLUX ' cdiag( 2) = 'VFLUX ' cdiag( 3) = 'HFLUX ' cdiag( 4) = 'EFLUX ' cdiag( 5) = 'QICE ' cdiag( 6) = 'RADLWG ' cdiag( 7) = 'RADSWG ' cdiag( 8) = 'RI ' cdiag( 9) = 'CT ' cdiag(10) = 'CU ' cdiag(11) = 'KM ' cdiag(12) = 'KH ' cdiag(13) = 'TURBU ' cdiag(14) = 'TURBV ' cdiag(15) = 'TURBT ' cdiag(16) = 'TURBQ ' cdiag(17) = 'MOISTT ' cdiag(18) = 'MOISTQ ' cdiag(19) = 'RADLW ' cdiag(20) = 'RADSW ' cdiag(21) = 'PREACC ' cdiag(22) = 'PRECON ' cdiag(23) = 'TUFLUX ' cdiag(24) = 'TVFLUX ' cdiag(25) = 'TTFLUX ' cdiag(26) = 'TQFLUX ' cdiag(27) = 'CN ' cdiag(28) = 'WINDS ' cdiag(29) = 'DTSRF ' cdiag(30) = 'TGROUND ' cdiag(31) = 'TS ' cdiag(32) = 'DTG ' cdiag(33) = 'QG ' cdiag(34) = 'QS ' cdiag(35) = 'TGRLW ' cdiag(36) = 'LWGUP ' cdiag(37) = 'OLR ' cdiag(38) = 'OLRCLR ' cdiag(39) = 'LWGCLR ' cdiag(40) = 'LWCLR ' cdiag(41) = 'TLW ' cdiag(42) = 'SHRAD ' cdiag(43) = 'OZLW ' cdiag(44) = 'CLDRAS ' cdiag(45) = 'CLDTOT ' cdiag(46) = 'LWGDOWN ' cdiag(47) = 'GWDT ' cdiag(48) = 'RADSWT ' cdiag(49) = 'TAUCLD ' cdiag(50) = 'TAUCLDC ' cdiag(51) = 'CLDLOW ' cdiag(52) = 'EVAP ' cdiag(53) = 'DPDT ' cdiag(54) = 'UAVE ' cdiag(55) = 'VAVE ' cdiag(56) = 'TAVE ' cdiag(57) = 'QAVE ' cdiag(58) = 'OMEGA ' cdiag(59) = 'DUDT ' cdiag(60) = 'DVDT ' cdiag(61) = 'DTDT ' cdiag(62) = 'DQDT ' cdiag(63) = 'VORT ' cdiag(64) = 'PSIGDOT ' cdiag(65) = 'DTLS ' cdiag(66) = 'DQLS ' cdiag(67) = 'USTAR ' cdiag(68) = 'Z0 ' cdiag(69) = 'FRQTRB ' cdiag(70) = 'PBL ' cdiag(71) = 'SWCLR ' cdiag(72) = 'OSR ' cdiag(73) = 'OSRCLR ' cdiag(74) = 'CLDMAS ' cdiag(75) = 'UWND ' cdiag(76) = 'VWND ' cdiag(77) = 'TMPU ' cdiag(78) = 'SPHU ' cdiag(79) = 'RFT ' cdiag(80) = 'PS ' cdiag(81) = 'TKE ' cdiag(82) = 'SWGCLR ' cdiag(83) = 'PAVE ' cdiag(84) = 'SDIAG1 ' cdiag(85) = 'SDIAG2 ' cdiag(86) = 'UDIAG1 ' cdiag(87) = 'UDIAG2 ' cdiag(88) = 'DIABU ' cdiag(89) = 'DIABV ' cdiag(90) = 'DIABT ' cdiag(91) = 'DIABQ ' cdiag(92) = 'RFU ' cdiag(93) = 'RFV ' cdiag(94) = 'GWDU ' cdiag(95) = 'GWDV ' cdiag(96) = 'GWDUS ' cdiag(97) = 'GWDVS ' cdiag(98) = 'GWDUT ' cdiag(99) = 'GWDVT ' cdiag(100) = 'LZRAD ' cdiag(101) = 'SLP ' cdiag(102) = 'VAVEUQ ' cdiag(103) = 'VAVEVQ ' cdiag(104) = 'VAVEUT ' cdiag(105) = 'VAVEVT ' cdiag(106) = 'CLDFRC ' cdiag(107) = 'TPW ' cdiag(108) = 'U2M ' cdiag(109) = 'V2M ' cdiag(110) = 'T2M ' cdiag(111) = 'Q2M ' cdiag(112) = 'U10M ' cdiag(113) = 'V10M ' cdiag(114) = 'T10M ' cdiag(115) = 'Q10M ' cdiag(116) = 'DTRAIN ' cdiag(117) = 'QFILL ' cdiag(118) = 'VAVEQAVE' cdiag(119) = 'VAVEQFIL' cdiag(120) = 'SHAPU ' cdiag(121) = 'SHAPV ' cdiag(122) = 'SHAPT ' cdiag(123) = 'SHAPQ ' cdiag(124) = 'SDIAG3 ' cdiag(125) = 'SDIAG4 ' cdiag(126) = 'SDIAG5 ' cdiag(127) = 'SDIAG6 ' cdiag(128) = 'SDIAG7 ' cdiag(129) = 'SDIAG8 ' cdiag(130) = 'SDIAG9 ' cdiag(131) = 'SDIAG10 ' cdiag(132) = 'UDIAG3 ' cdiag(133) = 'UDIAG4 ' cdiag(134) = 'UDIAG5 ' cdiag(135) = 'UDIAG6 ' cdiag(136) = 'UDIAG7 ' cdiag(137) = 'UDIAG8 ' cdiag(138) = 'UDIAG9 ' cdiag(139) = 'UDIAG10 ' cdiag(140) = 'CLDMID ' cdiag(141) = 'CLDHI ' cdiag(142) = 'TAULOW ' cdiag(143) = 'TAUMID ' cdiag(144) = 'TAUHI ' cdiag(145) = 'CLDNP ' cdiag(146) = 'CLDPRS ' cdiag(147) = 'CTPCNT ' cdiag(148) = 'CLDTMP ' cdiag(149) = 'CTTCNT ' cdiag(150) = 'TAULOWC ' cdiag(151) = 'TAUMIDC ' cdiag(152) = 'TAUHIC ' cdiag(153) = 'TCANOPY ' cdiag(154) = 'TDEEP ' cdiag(155) = 'QCANOPY ' cdiag(156) = 'SMSHAL ' cdiag(157) = 'SMROOT ' cdiag(158) = 'SMDEEP ' cdiag(159) = 'CAPACITY' cdiag(160) = 'SNOW ' cdiag(161) = 'RAINCON ' cdiag(162) = 'RAINLSP ' cdiag(163) = 'SNOWFALL' cdiag(164) = 'RUNOFF ' cdiag(165) = 'FWSOIL ' cdiag(166) = 'GDRAIN ' cdiag(167) = 'SNOWMELT' cdiag(168) = 'ERESV ' cdiag(169) = 'ESOIL ' cdiag(170) = 'EVEG ' cdiag(171) = 'ESNOW ' cdiag(172) = 'PARDF ' cdiag(173) = 'PARDR ' cdiag(174) = 'LAI ' cdiag(175) = 'GREEN ' cdiag(176) = 'DLWDTC ' cdiag(177) = 'DHDTC ' cdiag(178) = 'DEDTC ' cdiag(179) = 'VDTMOIST' cdiag(180) = 'VDTTURB ' cdiag(181) = 'VDTLWRAD' cdiag(182) = 'VDTSWRAD' cdiag(183) = 'VAVETAVE' cdiag(184) = 'PSUBCLD ' cdiag(185) = 'PSUBCLDC' cdiag(186) = 'LCL ' cdiag(187) = 'SDIAGC ' cdiag(188) = 'SDIAGCC ' cdiag(189) = ' ' cdiag(190) = 'EVPOT ' cdiag(191) = 'RHCHECK ' cdiag(192) = 'DHDQA ' cdiag(193) = 'DEDQA ' cdiag(194) = 'DTC ' cdiag(195) = 'DQC ' cdiag(196) = 'TCDTC ' cdiag(197) = 'RADDTC ' cdiag(198) = 'SENDDTC ' cdiag(199) = 'LATDTC ' cdiag(200) = 'TDDTC ' cdiag(201) = 'QCDTC ' cdiag(202) = 'TRBQLIQ ' cdiag(203) = 'TRBFCC ' cdiag(204) = 'ALBEDO ' cdiag(205) = 'ALBVISDR' cdiag(206) = 'ALBVISDF' cdiag(207) = 'ALBNIRDR' cdiag(208) = 'ALBNIRDF' cdiag(209) = 'TAUAVE ' cdiag(210) = 'ALBEDOC ' cdiag(211) = 'RHCHECKT' cdiag(212) = 'RHCHECKQ' cdiag(213) = 'CLDLSP ' cdiag(214) = 'LZLOW ' cdiag(215) = 'LZMID ' cdiag(216) = 'LZHI ' cdiag(217) = 'LZTOT ' cdiag(218) = 'CAPE ' cdiag(219) = 'HSUBCL ' cdiag(220) = 'HSTAR ' cdiag(221) = 'CONT ' cdiag(222) = 'CONQ ' cdiag(223) = 'LST ' cdiag(224) = 'LSQ ' cdiag(225) = 'CONEVPT ' cdiag(226) = 'CONEVPQ ' cdiag(227) = 'LSEVPT ' cdiag(228) = 'LSEVPQ ' cdiag(229) = 'GWCUMID ' cdiag(230) = 'GWCVMID ' cdiag(231) = 'CLDSTD ' cdiag(232) = 'GWCUBAR ' cdiag(233) = 'GWCVBAR ' cdiag(234) = 'GWCUS ' cdiag(235) = 'GWCVS ' cdiag(236) = 'GWCU ' cdiag(237) = 'GWCV ' cdiag(238) = 'ETAN ' cdiag(239) = 'ETANSQ ' cdiag(240) = 'THETA ' cdiag(241) = 'SALT ' cdiag(242) = 'UVEL ' cdiag(243) = 'VVEL ' cdiag(244) = 'WVEL ' cdiag(245) = 'THETASQ ' cdiag(246) = 'SALTSQ ' cdiag(247) = 'UVELSQ ' cdiag(248) = 'VVELSQ ' cdiag(249) = 'WVELSQ ' cdiag(250) = 'UVELVVEL' cdiag(251) = 'UVELMASS' cdiag(252) = 'VVELMASS' cdiag(253) = 'WVELMASS' cdiag(254) = 'UTHMASS ' cdiag(255) = 'VTHMASS ' cdiag(256) = 'WTHMASS ' cdiag(257) = 'USLTMASS' cdiag(258) = 'VSLTMASS' cdiag(259) = 'WSLTMASS' cdiag(260) = 'UVELTH ' cdiag(261) = 'VVELTH ' cdiag(262) = 'WVELTH ' cdiag(263) = 'UVELSLT ' cdiag(264) = 'VVELSLT ' cdiag(265) = 'WVELSLT ' cdiag(266) = 'UTRAC1 ' cdiag(267) = 'VTRAC1 ' cdiag(268) = 'WTRAC1 ' cdiag(269) = 'UTRAC2 ' cdiag(270) = 'VTRAC2 ' cdiag(271) = 'WTRAC2 ' cdiag(272) = 'UTRAC3 ' cdiag(273) = 'VTRAC3 ' cdiag(274) = 'WTRAC3 ' cdiag(275) = 'UTRAC4 ' cdiag(276) = 'VTRAC4 ' cdiag(277) = 'WTRAC4 ' cdiag(278) = 'UTRAC5 ' cdiag(279) = 'VTRAC5 ' cdiag(280) = 'WTRAC5 ' c Diagnostic Descriptions c ---------------------------------------------------------------------- tdiag( 1) = 'Zonal Wind Surface Stress (N/m^2)' tdiag( 2) = 'Meridional Wind Surface Stress (N/m^2)' tdiag( 3) = 'Sensible Heat Flux (pos.upwrd) (W/m^2)' tdiag( 4) = 'Latent Heat Flux (pos.upwrd) (W/m^2)' tdiag( 5) = 'Heat Conduction through Sea-Ice (W/m^2)' tdiag( 6) = 'Net Upward Longwave Flux at the Ground (W/m^2)' tdiag( 7) = 'Net Downward Shortwave Flux at the Ground (W/m^2)' tdiag( 8) = 'Richardson Number (non-dimensional)' tdiag( 9) = 'Surface Drag Coef for Scalars (non-dimensional)' tdiag(10) = 'Surface Drag Coef for Winds (non-dimensional)' tdiag(11) = 'Eddy Diffusivity Coefficient for Momentum (m^2/s)' tdiag(12) = 'Eddy Diffusivity Coefficient for Scalars (m^2/s)' tdiag(13) = 'Zonal Wind Tendency from Turbulence (m/s/day)' tdiag(14) = 'Meridional Wind Tendency from Turbulence (m/s/day)' tdiag(15) = . 'Temperature Tendency from Turbulence (K/day)' tdiag(16) = . 'Specific Humidity Tendency from Turbulence (g/kg/day)' tdiag(17) = . 'Temperature Tendency due to Moist Processes (K/day)' tdiag(18) = . 'Specific Humidity Tendency due to Moist Processes (g/kg/day)' tdiag(19) = . 'Temperature Tendency due to Longwave Radiation (K/day)' tdiag(20) = . 'Temperature Tendency due to Shortwave Radiation (K/day)' tdiag(21) = 'Total Precipitation (mm/day)' tdiag(22) = 'Convective Precipitation (mm/day)' tdiag(23) = 'Turbulent Flux of Zonal Momentum (N/m^2)' tdiag(24) = 'Turbulent Flux of Meridional Momentum (N/m^2)' tdiag(25) = 'Turbulent Flux of Sensible Heat (W/m^2)' tdiag(26) = 'Turbulent Flux of Latent Heat (W/m^2)' tdiag(27) = 'Neutral Drag Coefficient (non-dimensional)' tdiag(28) = 'Surface Wind Speed (m/s)' tdiag(29) = 'Air/Surface Virtual Temperature Difference (deg K)' tdiag(30) = 'Ground Temperature (deg K)' tdiag(31) = 'Surface Air Temperature (deg K)' tdiag(32) = 'Ground Temperature Adjustment (deg K)' tdiag(33) = 'Ground Specific Humidity (g/kg)' tdiag(34) = 'Saturation Surface Specific Humidity (g/kg)' tdiag(35) = 'Ground Temperature used by LW Radiation (deg K)' tdiag(36) = 'Upward Longwave Radiation at the Ground (W/m^2)' tdiag(37) = 'Outgoing Longwave Radiation (W/m^2)' tdiag(38) = 'Clear Sky Outgoing Longwave Radiation (W/m^2)' tdiag(39) = 'Clear Sky Net Longwave Flux at the Ground (W/m^2)' tdiag(40) = 'Clear Sky Longwave Radiation Heating Rates (K/day)' tdiag(41) = 'Temperature used by Radiation (deg K)' tdiag(42) = 'Specific Humidity used by Radiation (g/kg)' tdiag(43) = . 'Ozone mixing ratio used by Longwave Radiation (non-dimensional)' tdiag(44) = 'Convective Cloud Fraction (0-1)' tdiag(45) = '3-D Total Cloud Fraction (0-1)' tdiag(46) = 'Downward Longwave Radiation at the Ground (W/m^2)' tdiag(47) = . 'Temperature Tendency due to Gravity Wave Drag (K/day)' tdiag(48) = 'Incident Shortwave Radiation at TOA (W/m^2)' tdiag(49) = . 'Counted Cloud Optical Depth (non-dimensional) per 100 mb' tdiag(50) = 'Cloud Optical Depth Counter' tdiag(51) = 'Low-Level ( 1000-700 hPa) Cloud Fraction (0-1)' tdiag(52) = 'Surface Evaporation (mm/day)' tdiag(53) = 'Surface Pressure Tendency (hPa/day)' tdiag(54) = 'Average Zonal Wind (m/s)' tdiag(55) = 'Average Meridional Wind (m/s)' tdiag(56) = 'Average Potential Temperature (K)' tdiag(57) = . 'Average Specific Humidity (g/kg)' tdiag(58) = 'Vertical Velocity (hPa/day)' tdiag(59) = 'Total Zonal Wind Tendency (m/s/day)' tdiag(60) = 'Total Meridional Wind Tendency (m/s/day)' tdiag(61) = 'Total Temperature Tendency (K/day)' tdiag(62) = 'Total Specific Humidity Tendency (g/kg/day)' tdiag(63) = 'Relative Vorticity (x10^-4 sec^-1)' tdiag(64) = 'Pi-Weighted Vertical Sigma Velocity (hPa/day)' tdiag(65) = . 'Temperature Tendency due to Stratiform Processes (K/day)' tdiag(66) = . 'Specific Humidity Tendency from Stratiform Processes (g/kg/day)' tdiag(67) = 'Friction Velocity (m/s)' tdiag(68) = 'Surface Roughness Z0 (m)' tdiag(69) = 'Frequency of Turbulence (0-1)' tdiag(70) = 'Estimated PBL Depth (hPa)' tdiag(71) = 'Clear Sky Shortwave Heating Rates (K/day)' tdiag(72) = 'Outgoing Shortwave Radiation (W/m^2)' tdiag(73) = 'Clear Sky Outgoing Shortwave Radiation (W/m^2)' tdiag(74) = 'Cloud Mass Flux (g/m^2/sec)' tdiag(75) = 'Time-Averaged Zonal Wind (m/s)' tdiag(76) = 'Time-Averaged Meridional Wind (m/s)' tdiag(77) = 'Time-Averaged Temperature (deg K)' tdiag(78) = 'Time-Averaged Specific Humidity (g/kg)' tdiag(79) = . 'Temperature Tendency due to Rayleigh Friction (K/day)' tdiag(80) = 'Time-Averaged Surface Pressure (hPa)' tdiag(81) = 'Time-Averaged Turbulent Kinetic Energy (m/s)^2' tdiag(82) = . 'Clear Sky Net Downward SW Radiation at the Ground (W/m^2)' tdiag(83) = 'Average Surface Pressure (hPa)' tdiag(84) = 'User-Defined Surface Diagnostic #1' tdiag(85) = 'User-Defined Surface Diagnostic #2' tdiag(86) = 'User-Defined Upper-Air Diagnostic #1' tdiag(87) = 'User-Defined Upper-Air Diagnostic #2' tdiag(88) = . 'Zonal Wind Tendency due to Diabatic Forcing (m/s/day)' tdiag(89) = . 'Meridional Wind Tendency due to Diabatic Forcing (m/s/day)' tdiag(90) = . 'Temperature Tendency due to Diabatic Forcing (K/day)' tdiag(91) = . 'Specific Humidity Tendency due to Diabatic Forcing (g/kg/day)' tdiag(92) = . 'Zonal Wind Tendency due to Rayleigh Friction (m/s/day)' tdiag(93) = . 'Meridional Wind Tendency due to Rayleigh Friction (m/s/day)' tdiag(94) = . 'Zonal Wind Tendency due to Gravity Wave Drag (m/s/day)' tdiag(95) = . 'Meridional Wind Tendency due to Gravity Wave Drag (m/s/day)' tdiag(96) = 'Zonal Wind Gravity Wave Surface Stress (N/m^2)' tdiag(97) = 'Meridional Wind Gravity Wave Surface Stress (N/m^2)' tdiag(98) = 'Zonal Wind Gravity Wave Stress at Model Top (N/m^2)' tdiag(99) = . 'Meridional Wind Gravity Wave Stress at Model Top (N/m^2)' tdiag(100) = . 'Estimated Cloud Liquid Water used in Radiation (mg/kg)' tdiag(101) = 'Time-Averaged Sea-Level Pressure (mb)' tdiag(102) = 'Vertically Averaged UWND*SPHU (m/s g/kg)' tdiag(103) = 'Vertically Averaged VWND*SPHU (m/s g/kg)' tdiag(104) = 'Vertically Averaged UWND*TMPU (m/s deg K)' tdiag(105) = 'Vertically Averaged VWND*TMPU (m/s deg K)' tdiag(106) = '2-D Total Cloud Fraction (0-1)' tdiag(107) = 'Total Precipitable Water (g/cm^2)' tdiag(108) = 'Zonal Wind Interpolated to 2 Meters (m/s)' tdiag(109) = 'Meridional Wind Interpolated to 2 Meters (m/s)' tdiag(110) = 'Temperature Interpolated to 2 Meters (deg K)' tdiag(111) = 'Specific Humidity Interpolated to 2 Meters (g/kg)' tdiag(112) = 'Zonal Wind Interpolated to 10 Meters (m/s)' tdiag(113) = 'Meridional Wind Interpolated to 10 Meters (m/s)' tdiag(114) = 'Temperature Interpolated to 10 Meters (deg K)' tdiag(115) = 'Specific Humidity Interpolated to 10 Meters (g/kg)' tdiag(116) = 'Detrainment Cloud Mass Flux (g/m^2/sec)' tdiag(117) = 'Filling of Negative Specific Humidity (g/kg/day)' tdiag(118) = 'Vertically Averaged QAVE (mm/day)' tdiag(119) = 'Vertically Averaged QFILL (mm/day)' tdiag(120) = . 'Zonal Wind Tendency due to Shapiro Filter (m/s/day)' tdiag(121) = . 'Meridional Wind Tendency due to Shapiro Filter (m/s/day)' tdiag(122) = 'Temperature Tendency due to Shapiro Filter (K/day)' tdiag(123) = . 'Specific Humidity Tendency due to Shapiro Filter (g/kg/day)' tdiag(124) = 'User-Defined Surface Diagnostic #3 ' tdiag(125) = 'User-Defined Surface Diagnostic #4 ' tdiag(126) = 'User-Defined Surface Diagnostic #5 ' tdiag(127) = 'User-Defined Surface Diagnostic #6 ' tdiag(128) = 'User-Defined Surface Diagnostic #7 ' tdiag(129) = 'User-Defined Surface Diagnostic #8 ' tdiag(130) = 'User-Defined Surface Diagnostic #9 ' tdiag(131) = 'User-Defined Surface Diagnostic #10 ' tdiag(132) = 'User-Defined Upper-Air Diagnostic #3 ' tdiag(133) = 'User-Defined Upper-Air Diagnostic #4 ' tdiag(134) = 'User-Defined Upper-Air Diagnostic #5 ' tdiag(135) = 'User-Defined Upper-Air Diagnostic #6 ' tdiag(136) = 'User-Defined Upper-Air Diagnostic #7 ' tdiag(137) = 'User-Defined Upper-Air Diagnostic #8 ' tdiag(138) = 'User-Defined Upper-Air Diagnostic #9 ' tdiag(139) = 'User-Defined Upper-Air Diagnostic #10 ' tdiag(140) = 'Mid-Level (700-400 hPa) Cloud Fraction (0-1)' tdiag(141) = 'High-Level (above 400 hPa) Cloud Fraction (0-1)' tdiag(142) = 'Low-Level (1000-700 hPa) Optical Depth' tdiag(143) = 'Mid-Level (700-400 hPa) Optical Depth' tdiag(144) = 'High-Level (above 400 hPa) Optical Depth' tdiag(145) = 'Non-Precipitating Cloud Fraction (0-1)' tdiag(146) = 'Cloud Top Pressure (when cloudy) (hPa)' tdiag(147) = 'Cloud Top Pressure Counter' tdiag(148) = 'Cloud Top Temperature (when cloudy) (deg K)' tdiag(149) = 'Cloud Top Temperature Counter' tdiag(150) = 'Low-Level Optical Depth Counter' tdiag(151) = 'Mid-Level Optical Depth Counter' tdiag(152) = 'High-Level Optical Depth Counter' tdiag(153) = . 'Skin Temperature of the Canopy/Soil (SST over water) (deg K)' tdiag(154) = 'Deep Soil Temperature (deg K)' tdiag(155) = 'Specific Humidity of the Canopy (g/kg)' tdiag(156) = . 'Soil Moisture Field Capacity Fraction in Shallow (5 mm) Layer' tdiag(157) = 'Soil Moisture Field Capacity Fraction in Root Zone' tdiag(158) = . 'Soil Moisture Field Capacity Fraction in Recharge Layer ' tdiag(159) = 'Canopy Interception Reservoir (mm)' tdiag(160) = 'Snow Depth (mm water equivalent)' tdiag(161) = 'Convective Rainfall (liquid precipitate) (mm/day)' tdiag(162) = 'Large-scale Rainfall (liquid precipitate) (mm/day)' tdiag(163) = 'Total Snowfall (solid precipitate) (mm/day)' tdiag(164) = . 'Water from precipitation not infiltrated into soil (mm/s)' tdiag(165) = 'Rainwater Infiltration into top soil layer (mm/s)' tdiag(166) = . 'Diffusion of Moisture across bottom of root zone (mm/s)' tdiag(167) = 'Rate of Snow Melt (mm/s)' tdiag(168) = . 'Latent Heat Flux from Canopy Interception Reservoir (W/m^2)' tdiag(169) = 'Latent Heat Flux from Bare Soil (W/m^2)' tdiag(170) = . 'Latent Heat Flux (Transpiration) from Vegetated Surface (W/m^2)' tdiag(171) = 'Latent Heat Flux from Snow Pack (W/m^2)' tdiag(172) = . 'Diffuse-beam Photosynthetically Active Radiation (W/m^2)' tdiag(173) = . 'Direct-beam Photosynthetically Active Radiation (W/m^2)' tdiag(174) = 'Leaf Area Index (%)' tdiag(175) = 'Green-ness Index (%)' tdiag(176) = . 'Derivative of Upward LW Radiation wrt TCANOPY (W/m^2 /deg K)' tdiag(177) = . 'Derivative of Sensible Heat Flux wrt TCANOPY (W/m^2 /deg K)' tdiag(178) = . 'Derivative of Latent Heat Flux wrt TCANOPY (W/m^2 /deg K)' tdiag(179) = 'Vertically Averaged MOISTT (K/day)' tdiag(180) = 'Vertically Averaged TURBT (K/day)' tdiag(181) = 'Vertically Averaged RADLW (K/day)' tdiag(182) = 'Vertically Averaged RADSW (K/day)' tdiag(183) = 'Vertically Averaged TAVE (K/day)' tdiag(184) = 'Pressure Thickness of the RAS Subcloud Layer (mb)' tdiag(185) = . 'Pressure Thickness of the RAS Subcloud Layer Counter' tdiag(186) = 'Lifting Condensation Level (mb) ' tdiag(187) = 'User-Defined Counted Surface Diagnostic ' tdiag(188) = 'User-Defined Counted Surface Diagnostic Counter' tdiag(189) = ' ' tdiag(190) = 'Potential Evapotranspiration (W/m^2)' tdiag(191) = . 'Net precip due to supersaturation above critical RH (mm/day)' tdiag(192) = . 'Derivative of Sensible Heat wrt Canopy Humidity (W/m^2)' tdiag(193) = . 'Derivative of Latent Heat wrt Canopy Humidity (W/m^2)' tdiag(194) = 'Total Change in Canopy Temperature (deg/sec)' tdiag(195) = 'Change in Canopy Humidity (kg/kg/sec)' tdiag(196) = . 'Change in Canopy Temperature - Energy Term (deg/sec)' tdiag(197) = . 'Change in Canopy Temperature - Net Radiation Term (deg/sec)' tdiag(198) = . 'Change in Canopy Temperature - Sensible Heat Term (deg/sec)' tdiag(199) = . 'Change in Canopy Temperature - Latent Heat Term (deg/sec)' tdiag(200) = . 'Change in Canopy Temperature - Ground Temp Term (deg/sec)' tdiag(201) = . 'Change in Canopy Temperature - Humidity Term (deg/sec)' tdiag(202) = 'Cloud Liquid Water Generated in Turbulence (mg/kg)' tdiag(203) = 'Cloud Fraction Generated in Turbulence (fraction)' tdiag(204) = 'Mean Surface Albedo (SW_Up/SW_Down) 0-1 ' tdiag(205) = 'Direct Beam VIS Surface Albedo (0-1) ' tdiag(206) = 'Diffuse Beam VIS Surface Albedo (0-1) ' tdiag(207) = 'Direct Beam NIR Surface Albedo (0-1) ' tdiag(208) = 'Diffuse Beam NIR Surface Albedo (0-1) ' tdiag(209) = . 'Averaged Cloud Optical Depth (non-dimensional) per 100 mb' tdiag(210) = 'Surface Albedo Counter' tdiag(211) = . 'Temperature Change due to humidity above critical RH (K/day)' tdiag(212) = . 'Spec.Humid. Change due to humidity above critical RH (g/kg/day)' tdiag(213) = 'Large-Scale Cloud Fraction (0-1)' tdiag(214) = . 'Vert integrated low-level (below 700mb) cloud water (mg/cm**2)' tdiag(215) = . 'Vert integrated mid-level (700-400mb) cloud water (mg/cm**2)' tdiag(216) = . 'Vert integrated high-level (above 400mb) cloud water (mg/cm**2)' tdiag(217) = . 'Vertically integrated total suspended cloud water (mg/cm**2)' tdiag(218) = 'Convective Available Potential Energy' tdiag(219) = 'Moist Static Energy in RAS subcloud layer' tdiag(220) = 'Saturation Moist Static Energy' tdiag(221) = . 'Temperature Tendency due to Convective Precipitation (K/day)' tdiag(222) = . 'Moisture Tendency due to Convective Precipitation (g/kg/day)' tdiag(223) = . 'Temperature Tendency due to Large Scale Precipitation (K/day)' tdiag(224) = . 'Moisture Tendency due to Large Scale Precipitation (g/kg/day)' tdiag(225) = . 'Temperature Tendency due to Conv Rain Re-evaporation (K/day)' tdiag(226) = . 'Moisture Tendency due to Conv Rain Re-evaporation (g/kg/day)' tdiag(227) = . 'Temperature Tendency due to Large Scale Rain Re-evap (K/day)' tdiag(228) = . 'Moisture Tendency due to Large Scale Rain Re-evap (g/kg/day)' tdiag(229) = . 'Mid-Lev Mean U-Wind for Gravity Wave Drag w/Clouds (m/sec)' tdiag(230) = . 'Mid-Lev Mean V-Wind for Gravity Wave Drag w/Clouds (m/sec)' tdiag(231) = . 'Mid-Lev Cloud Perturb Length for Gravity Wave Drag w/Clouds (m)' tdiag(232) = . 'Mid-Lev Perturb U-Wind for Gravity Wave Drag w/Clouds (m/sec)' tdiag(233) = . 'Mid-Lev Perturb V-Wind for Gravity Wave Drag w/Clouds (m/sec)' tdiag(234) = . 'Mid-Lev U-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)' tdiag(235) = . 'Mid-Lev V-Wind Stress for Gravity Wave Drag w/Clouds (N/m^2)' tdiag(236) = . 'U-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)' tdiag(237) = . 'V-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)' tdiag(238) = 'Perturbation of Surface (pressure, height) (hPa,m)' tdiag(239) = 'Square of Perturbation of Sfc (hPa^2,m^2)' tdiag(240) = 'Potential Temperature (K)' tdiag(241) = 'Salt or Water Vapor mixing ratio (g/kg)' tdiag(242) = 'Zonal Component of Velocity (m/s)' tdiag(243) = 'Meridional Component of Velocity (m/s)' tdiag(244) = 'Vertical Component of Velocity (m/s)' tdiag(245) = 'Square of Potential Temperature (K^2)' tdiag(246) = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)' tdiag(247) = 'Square of Zonal Comp of Velocity (m^2/s^2)' tdiag(248) = 'Square of Meridional Comp of Velocity (m^2/s^2)' tdiag(249) = 'Square of Vertical Comp of Velocity (m^2/s^2)' tdiag(250) = 'Meridional Transport of Zonal Momentum (m^2/s^2)' tdiag(251) = 'Zonal Mass-Weighted Comp of Velocity (m/s)' tdiag(252) = 'Meridional Mass-Weighted Comp of Velocity (m/s)' tdiag(253) = 'Vertical Mass-Weighted Comp of Velocity (m/s)' tdiag(254) = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)' tdiag(255) = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)' tdiag(256) = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)' tdiag(257) = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)' tdiag(258) = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)' tdiag(259) = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)' tdiag(260) = 'Zonal Transp of Pot Temp (m-K/s)' tdiag(261) = 'Meridional Transp of Pot Temp (m-K/s)' tdiag(262) = 'Vertical Transp of Pot Temp (m-K/s)' tdiag(263) = 'Zonal Transp of Salt (m-kg/s-kg)' tdiag(264) = 'Meridional Transp of Salt (m-kg/s-kg)' tdiag(265) = 'Vertical Transp of Salt (m-kg/s-kg)' tdiag(266) = 'Zonal Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)' tdiag(267) = 'Merid Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)' tdiag(268) = 'Vert Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)' tdiag(269) = 'Zonal Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)' tdiag(270) = 'Merid Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)' tdiag(271) = 'Vert Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)' tdiag(272) = 'Zonal Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)' tdiag(273) = 'Merid Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)' tdiag(274) = 'Vert Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)' tdiag(275) = 'Zonal Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)' tdiag(276) = 'Merid Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)' tdiag(277) = 'Vert Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)' tdiag(278) = 'Zonal Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' tdiag(279) = 'Merid Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' tdiag(280) = 'Vert Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' c Diagnostic Units c ---------------- udiag( 1) = 'N/m2 ' udiag( 2) = 'N/m2 ' udiag( 3) = 'watt/m2 ' udiag( 4) = 'watt/m2 ' udiag( 5) = 'watt/m2 ' udiag( 6) = 'watt/m2 ' udiag( 7) = 'watt/m2 ' udiag( 8) = 'non-dimensional ' udiag( 9) = 'non-dimensional ' udiag(10) = 'non-dimensional ' udiag(11) = 'm2/sec ' udiag(12) = 'm2/sec ' udiag(13) = 'm/s/day ' udiag(14) = 'm/s/day ' udiag(15) = 'K/day ' udiag(16) = 'gram/kg/day ' udiag(17) = 'K/day ' udiag(18) = 'gram/kg/day ' udiag(19) = 'K/day ' udiag(20) = 'K/day ' udiag(21) = 'mm/day ' udiag(22) = 'mm/day ' udiag(23) = 'N/m2 ' udiag(24) = 'N/m2 ' udiag(25) = 'watt/m2 ' udiag(26) = 'watt/m2 ' udiag(27) = 'non-dimensional ' udiag(28) = 'm/sec ' udiag(29) = 'K ' udiag(30) = 'K ' udiag(31) = 'K ' udiag(32) = 'K ' udiag(33) = 'gram/kg ' udiag(34) = 'gram/kg ' udiag(35) = 'K ' udiag(36) = 'watt/m2 ' udiag(37) = 'watt/m2 ' udiag(38) = 'watt/m2 ' udiag(39) = 'watt/m2 ' udiag(40) = 'K/day ' udiag(41) = 'K ' udiag(42) = 'gram/kg ' udiag(43) = 'non-dimensional ' udiag(44) = '0-1 ' udiag(45) = '0-1 ' udiag(46) = 'watt/m2 ' udiag(47) = 'K/day ' udiag(48) = 'watt/m2 ' udiag(49) = 'non-dimensional ' udiag(50) = 'count ' udiag(51) = '0-1 ' udiag(52) = 'mm/day ' udiag(53) = 'hPa/day ' udiag(54) = 'm/s ' udiag(55) = 'm/s ' udiag(56) = 'K ' udiag(57) = 'gram/kg ' udiag(58) = 'hPa/day ' udiag(59) = 'm/s/day ' udiag(60) = 'm/s/day ' udiag(61) = 'K/day ' udiag(62) = 'gram/kg/day ' udiag(63) = '(e-4 sec-1) ' udiag(64) = 'hPa/day ' udiag(65) = 'K/day ' udiag(66) = 'gram/kg/day ' udiag(67) = 'm/sec ' udiag(68) = 'm ' udiag(69) = '0-1 ' udiag(70) = 'hPa ' udiag(71) = 'K/day ' udiag(72) = 'watt/m2 ' udiag(73) = 'watt/m2 ' udiag(74) = 'g/m2/sec ' udiag(75) = 'm/sec ' udiag(76) = 'm/sec ' udiag(77) = 'K ' udiag(78) = 'gram/kg ' udiag(79) = 'K/day ' udiag(80) = 'hPa ' udiag(81) = '(m/sec)2 ' udiag(82) = 'watt/m2 ' udiag(83) = 'hPa/day ' udiag(84) = 'user-defined ' udiag(85) = 'user-defined ' udiag(86) = 'user-defined ' udiag(87) = 'user-defined ' udiag(88) = 'm/s/day ' udiag(89) = 'm/s/day ' udiag(90) = 'K/day ' udiag(91) = 'gram/kg/day ' udiag(92) = 'm/s/day ' udiag(93) = 'm/s/day ' udiag(94) = 'm/s/day ' udiag(95) = 'm/s/day ' udiag(96) = 'N/m2 ' udiag(97) = 'N/m2 ' udiag(98) = 'N/m2 ' udiag(99) = 'N/m2 ' udiag(100) = 'mg/kg ' udiag(101) = 'hPa ' udiag(102) = 'm/s gram/kg ' udiag(103) = 'm/s gram/kg ' udiag(104) = 'm/s K ' udiag(105) = 'm/s K ' udiag(106) = '0-1 ' udiag(107) = 'gram/cm2 ' udiag(108) = 'm/sec ' udiag(109) = 'm/sec ' udiag(110) = 'K ' udiag(111) = 'gram/kg ' udiag(112) = 'm/sec ' udiag(113) = 'm/sec ' udiag(114) = 'K ' udiag(115) = 'gram/kg ' udiag(116) = 'g/m2/sec ' udiag(117) = 'gram/kg/day ' udiag(118) = 'mm/day ' udiag(119) = 'mm/day ' udiag(120) = 'm/s/day ' udiag(121) = 'm/s/day ' udiag(122) = 'K/day ' udiag(123) = 'gram/kg/day ' udiag(124) = 'user-defined ' udiag(125) = 'user-defined ' udiag(126) = 'user-defined ' udiag(127) = 'user-defined ' udiag(128) = 'user-defined ' udiag(129) = 'user-defined ' udiag(130) = 'user-defined ' udiag(131) = 'user-defined ' udiag(132) = 'user-defined ' udiag(133) = 'user-defined ' udiag(134) = 'user-defined ' udiag(135) = 'user-defined ' udiag(136) = 'user-defined ' udiag(137) = 'user-defined ' udiag(138) = 'user-defined ' udiag(139) = 'user-defined ' udiag(140) = '0-1 ' udiag(141) = '0-1 ' udiag(142) = 'non-dimensional ' udiag(143) = 'non-dimensional ' udiag(144) = 'non-dimensional ' udiag(145) = '0-1 ' udiag(146) = 'hPa ' udiag(147) = 'count ' udiag(148) = 'K ' udiag(149) = 'count ' udiag(150) = 'count ' udiag(151) = 'count ' udiag(152) = 'count ' udiag(153) = 'K ' udiag(154) = 'K ' udiag(155) = 'gram/kg ' udiag(156) = '0-1 ' udiag(157) = '0-1 ' udiag(158) = '0-1 ' udiag(159) = 'mm ' udiag(160) = 'mm water equiv ' udiag(161) = 'mm/day ' udiag(162) = 'mm/day ' udiag(163) = 'mm/day ' udiag(164) = 'mm/sec ' udiag(165) = 'mm/sec ' udiag(166) = 'mm/sec ' udiag(167) = 'mm/sec ' udiag(168) = 'watt/m2 ' udiag(169) = 'watt/m2 ' udiag(170) = 'watt/m2 ' udiag(171) = 'watt/m2 ' udiag(172) = 'watt/m2 ' udiag(173) = 'watt/m2 ' udiag(174) = 'percent ' udiag(175) = 'percent ' udiag(176) = 'watt/m2/K ' udiag(177) = 'watt/m2/K ' udiag(178) = 'watt/m2/K ' udiag(179) = 'K/day ' udiag(180) = 'K/day ' udiag(181) = 'K/day ' udiag(182) = 'K/day ' udiag(183) = 'K/day ' udiag(184) = 'hPa ' udiag(185) = 'count ' udiag(186) = 'hPa ' udiag(187) = 'user-defined ' udiag(188) = 'count ' udiag(189) = ' ' udiag(190) = 'watt/m2 ' udiag(191) = 'mm/day ' udiag(192) = 'watt/m2 ' udiag(193) = 'watt/m2 ' udiag(194) = 'K/sec ' udiag(195) = 'sec-1 ' udiag(196) = 'K/sec ' udiag(197) = 'K/sec ' udiag(198) = 'K/sec ' udiag(199) = 'K/sec ' udiag(200) = 'K/sec ' udiag(201) = 'K/sec ' udiag(202) = 'mg/kg ' udiag(203) = '0-1 ' udiag(204) = '0-1 ' udiag(205) = '0-1 ' udiag(206) = '0-1 ' udiag(207) = '0-1 ' udiag(208) = '0-1 ' udiag(209) = 'non-dimensional ' udiag(210) = 'count ' udiag(211) = 'K/day ' udiag(212) = 'g/kg/day ' udiag(213) = '0-1 ' udiag(214) = 'mg/cm2 ' udiag(215) = 'mg/cm2 ' udiag(216) = 'mg/cm2 ' udiag(217) = 'mg/cm2 ' udiag(218) = 'mb m2/sec2 ' udiag(219) = 'm2/sec2 ' udiag(220) = 'm2/sec2 ' udiag(221) = 'K/day ' udiag(222) = 'g/kg/day ' udiag(223) = 'K/day ' udiag(224) = 'g/kg/day ' udiag(225) = 'K/day ' udiag(226) = 'g/kg/day ' udiag(227) = 'K/day ' udiag(228) = 'g/kg/day ' udiag(229) = 'm/sec ' udiag(230) = 'm/sec ' udiag(231) = 'm ' udiag(232) = 'm/sec ' udiag(233) = 'm/sec ' udiag(234) = 'N/m2 ' udiag(235) = 'N/m2 ' udiag(236) = 'm/sec/day ' udiag(237) = 'm/sec/day ' udiag(238) = 'm or hPa ' udiag(239) = 'm2 or hPa2 ' udiag(240) = 'K ' udiag(241) = 'g/kg ' udiag(242) = 'm/sec ' udiag(243) = 'm/sec ' udiag(244) = 'm/sec ' udiag(245) = 'K2 ' udiag(246) = 'g2/kg2 ' udiag(247) = 'm2/sec2 ' udiag(248) = 'm2/sec2 ' udiag(249) = 'm2/sec2 ' udiag(250) = 'm2/sec2 ' udiag(251) = 'm/sec ' udiag(252) = 'm/sec ' udiag(253) = 'm/sec ' udiag(254) = 'm-K/sec ' udiag(255) = 'm-K/sec ' udiag(256) = 'm-K/sec ' udiag(257) = 'm-kg/s-kg ' udiag(258) = 'm-kg/s-kg ' udiag(259) = 'm-kg/s-kg ' udiag(260) = 'm-K/sec ' udiag(261) = 'm-K/sec ' udiag(262) = 'm-K/sec ' udiag(263) = 'm-kg/s-kg ' udiag(264) = 'm-kg/s-kg ' udiag(265) = 'm-kg/s-kg ' udiag(266) = 'm-kg/s-kg ' udiag(267) = 'm-kg/s-kg ' udiag(268) = 'm-kg/s-kg ' udiag(269) = 'm-kg/s-kg ' udiag(270) = 'm-kg/s-kg ' udiag(271) = 'm-kg/s-kg ' udiag(272) = 'm-kg/s-kg ' udiag(273) = 'm-kg/s-kg ' udiag(274) = 'm-kg/s-kg ' udiag(275) = 'm-kg/s-kg ' udiag(276) = 'm-kg/s-kg ' udiag(277) = 'm-kg/s-kg ' udiag(278) = 'm-kg/s-kg ' udiag(279) = 'm-kg/s-kg ' udiag(280) = 'm-kg/s-kg ' c Diagnostic Parms c ---------------- gdiag( 1) = 'UM 002L ' gdiag( 2) = 'VM 001L ' gdiag( 3) = 'SM L ' gdiag( 4) = 'SM L ' gdiag( 5) = 'SM L ' gdiag( 6) = 'SM L ' gdiag( 7) = 'SM P L ' gdiag( 8) = 'SM L ' gdiag( 9) = 'SM L ' gdiag(10) = 'SM L ' gdiag(11) = 'SM L ' gdiag(12) = 'SM L ' gdiag(13) = 'UM 014M ' gdiag(14) = 'VM 013M ' gdiag(15) = 'SM M ' gdiag(16) = 'SM M ' gdiag(17) = 'SM M ' gdiag(18) = 'SM M ' gdiag(19) = 'SM M ' gdiag(20) = 'SM P M ' gdiag(21) = 'SM P L ' gdiag(22) = 'SM P L ' gdiag(23) = 'UM 024L ' gdiag(24) = 'VM 023L ' gdiag(25) = 'SM L ' gdiag(26) = 'SM L ' gdiag(27) = 'SM L ' gdiag(28) = 'SM L ' gdiag(29) = 'SM L ' gdiag(30) = 'SM L ' gdiag(31) = 'SM L ' gdiag(32) = 'SM L ' gdiag(33) = 'SM P L ' gdiag(34) = 'SM P L ' gdiag(35) = 'SM L ' gdiag(36) = 'SM L ' gdiag(37) = 'SM U ' gdiag(38) = 'SM U ' gdiag(39) = 'SM L ' gdiag(40) = 'SM M ' gdiag(41) = 'SM M ' gdiag(42) = 'SM P M ' gdiag(43) = 'SM M ' gdiag(44) = 'SM P M ' gdiag(45) = 'SM P M ' gdiag(46) = 'SM L ' gdiag(47) = 'SM M ' gdiag(48) = 'SM P U ' gdiag(49) = 'SM C050M ' gdiag(50) = 'SM D M ' gdiag(51) = 'SM P M ' gdiag(52) = 'SM L ' gdiag(53) = 'SM L ' gdiag(54) = 'UU 055M ' gdiag(55) = 'VV 054M ' gdiag(56) = 'SM M ' gdiag(57) = 'SM M ' gdiag(58) = 'SM M ' gdiag(59) = 'UU 060M ' gdiag(60) = 'VV 059M ' gdiag(61) = 'SM M ' gdiag(62) = 'SM M ' gdiag(63) = 'SM M ' gdiag(64) = 'SM L ' gdiag(65) = 'SM M ' gdiag(66) = 'SM M ' gdiag(67) = 'SM L ' gdiag(68) = 'SM L ' gdiag(69) = 'SM P M ' gdiag(70) = 'SM P M ' gdiag(71) = 'SM P M ' gdiag(72) = 'SM P U ' gdiag(73) = 'SM P U ' gdiag(74) = 'SM P L ' gdiag(75) = 'UU 076M ' gdiag(76) = 'VV 075M ' gdiag(77) = 'SM M ' gdiag(78) = 'SM P M ' gdiag(79) = 'SM M ' gdiag(80) = 'SM L ' gdiag(81) = 'SM P L ' gdiag(82) = 'SM P L ' gdiag(83) = 'SM L ' gdiag(84) = 'SM L ' gdiag(85) = 'SM L ' gdiag(86) = 'SM M ' gdiag(87) = 'SM M ' gdiag(88) = 'UU 089M ' gdiag(89) = 'VV 088M ' gdiag(90) = 'SM M ' gdiag(91) = 'SM M ' gdiag(92) = 'UU 093M ' gdiag(93) = 'VV 092M ' gdiag(94) = 'UM 095M ' gdiag(95) = 'VM 094M ' gdiag(96) = 'UM 097L ' gdiag(97) = 'VM 096L ' gdiag(98) = 'UM 099U ' gdiag(99) = 'VM 098U ' gdiag(100) = 'SM M ' gdiag(101) = 'SM L ' gdiag(102) = 'UM 103M ' gdiag(103) = 'VM 102M ' gdiag(104) = 'UM 105M ' gdiag(105) = 'VM 104M ' gdiag(106) = 'SM P M ' gdiag(107) = 'SM P M ' gdiag(108) = 'UM 109M ' gdiag(109) = 'VM 108M ' gdiag(110) = 'SM M ' gdiag(111) = 'SM P M ' gdiag(112) = 'UM 113M ' gdiag(113) = 'VM 112M ' gdiag(114) = 'SM M ' gdiag(115) = 'SM P M ' gdiag(116) = 'SM P L ' gdiag(117) = 'SM M ' gdiag(118) = 'SM M ' gdiag(119) = 'SM P L ' gdiag(120) = 'UU 121M ' gdiag(121) = 'VV 120M ' gdiag(122) = 'SM M ' gdiag(123) = 'SM M ' gdiag(124) = 'SM L ' gdiag(125) = 'SM L ' gdiag(126) = 'SM L ' gdiag(127) = 'SM L ' gdiag(128) = 'SM L ' gdiag(129) = 'SM L ' gdiag(130) = 'SM L ' gdiag(131) = 'SM L ' gdiag(132) = 'SM M ' gdiag(133) = 'SM M ' gdiag(134) = 'SM M ' gdiag(135) = 'SM M ' gdiag(136) = 'SM M ' gdiag(137) = 'SM M ' gdiag(138) = 'SM M ' gdiag(139) = 'SM M ' gdiag(140) = 'SM P M ' gdiag(141) = 'SM P M ' gdiag(142) = 'SM PC150M ' gdiag(143) = 'SM PC151M ' gdiag(144) = 'SM PC152M ' gdiag(145) = 'SM M ' gdiag(146) = 'SM PC147M ' gdiag(147) = 'SM D M ' gdiag(148) = 'SM PC149M ' gdiag(149) = 'SM D M ' gdiag(150) = 'SM D M ' gdiag(151) = 'SM D M ' gdiag(152) = 'SM D M ' gdiag(153) = 'SM L ' gdiag(154) = 'SM L ' gdiag(155) = 'SM L ' gdiag(156) = 'SM L ' gdiag(157) = 'SM L ' gdiag(158) = 'SM L ' gdiag(159) = 'SM L ' gdiag(160) = 'SM L ' gdiag(161) = 'SM L ' gdiag(162) = 'SM L ' gdiag(163) = 'SM P L ' gdiag(164) = 'SM L ' gdiag(165) = 'SM L ' gdiag(166) = 'SM L ' gdiag(167) = 'SM L ' gdiag(168) = 'SM L ' gdiag(169) = 'SM L ' gdiag(170) = 'SM L ' gdiag(171) = 'SM L ' gdiag(172) = 'SM P L ' gdiag(173) = 'SM P L ' gdiag(174) = 'SM L ' gdiag(175) = 'SM L ' gdiag(176) = 'SM L ' gdiag(177) = 'SM L ' gdiag(178) = 'SM L ' gdiag(179) = 'SM M ' gdiag(180) = 'SM M ' gdiag(181) = 'SM M ' gdiag(182) = 'SM P M ' gdiag(183) = 'SM M ' gdiag(184) = 'SM PC185M ' gdiag(185) = 'SM PD M ' gdiag(186) = 'SM L ' gdiag(187) = 'SM C188L ' gdiag(188) = 'SM D L ' gdiag(189) = ' ' gdiag(190) = 'SM L ' gdiag(191) = 'SM P L ' gdiag(192) = 'SM L ' gdiag(193) = 'SM L ' gdiag(194) = 'SM L ' gdiag(195) = 'SM L ' gdiag(196) = 'SM L ' gdiag(197) = 'SM L ' gdiag(198) = 'SM L ' gdiag(199) = 'SM L ' gdiag(200) = 'SM L ' gdiag(201) = 'SM L ' gdiag(202) = 'SM M ' gdiag(203) = 'SM M ' gdiag(204) = 'SM C210L ' gdiag(205) = 'SM L ' gdiag(206) = 'SM L ' gdiag(207) = 'SM L ' gdiag(208) = 'SM L ' gdiag(209) = 'SM M ' gdiag(210) = 'SM D L ' gdiag(211) = 'SM M ' gdiag(212) = 'SM M ' gdiag(213) = 'SM P M ' gdiag(214) = 'SM P L ' gdiag(215) = 'SM P L ' gdiag(216) = 'SM P L ' gdiag(217) = 'SM P L ' gdiag(218) = 'SM P L ' gdiag(219) = 'SM P L ' gdiag(220) = 'SM P M ' gdiag(221) = 'SM M ' gdiag(222) = 'SM M ' gdiag(223) = 'SM M ' gdiag(224) = 'SM M ' gdiag(225) = 'SM M ' gdiag(226) = 'SM M ' gdiag(227) = 'SM M ' gdiag(228) = 'SM M ' gdiag(229) = 'UM 230L ' gdiag(230) = 'VM 229L ' gdiag(231) = 'SM L ' gdiag(232) = 'UM 233L ' gdiag(233) = 'VM 232L ' gdiag(234) = 'UM 235L ' gdiag(235) = 'VM 234L ' gdiag(236) = 'UM 237M ' gdiag(237) = 'VM 236M ' gdiag(238) = 'SM M ' gdiag(239) = 'SM M ' gdiag(240) = 'SM M ' gdiag(241) = 'SM M ' gdiag(242) = 'UU 243M ' gdiag(243) = 'VV 242M ' gdiag(244) = 'WM L ' gdiag(245) = 'SM M ' gdiag(246) = 'SM M ' gdiag(247) = 'UU 243M ' gdiag(248) = 'VV 242M ' gdiag(249) = 'WM L ' gdiag(250) = 'UZ M ' gdiag(251) = 'UU 252M ' gdiag(252) = 'VV 251M ' gdiag(253) = 'WM M ' gdiag(254) = 'UU 255M ' gdiag(255) = 'VV 254M ' gdiag(256) = 'WM M ' gdiag(257) = 'UU 258M ' gdiag(258) = 'VV 257M ' gdiag(259) = 'WM M ' gdiag(260) = 'UU 261M ' gdiag(261) = 'VV 260M ' gdiag(262) = 'WM M ' gdiag(263) = 'UU 264M ' gdiag(264) = 'VV 263M ' gdiag(265) = 'WM M ' gdiag(266) = 'UU 267M ' gdiag(267) = 'VV 266M ' gdiag(268) = 'WM M ' gdiag(269) = 'UU 270M ' gdiag(270) = 'VV 269M ' gdiag(271) = 'WM M ' gdiag(272) = 'UU 273M ' gdiag(273) = 'VV 272M ' gdiag(274) = 'WM M ' gdiag(275) = 'UU 276M ' gdiag(276) = 'VV 275M ' gdiag(277) = 'WM M ' gdiag(278) = 'UU 279M ' gdiag(279) = 'VV 278M ' gdiag(280) = 'WM M ' c Diagnostic Levels c ----------------- KDIAG( 1) = 1 KDIAG( 2) = 1 KDIAG( 3) = 1 KDIAG( 4) = 1 KDIAG( 5) = 1 KDIAG( 6) = 1 KDIAG( 7) = 1 KDIAG( 8) = nrphys KDIAG( 9) = 1 KDIAG(10) = 1 KDIAG(11) = nrphys KDIAG(12) = nrphys KDIAG(13) = nrphys KDIAG(14) = nrphys KDIAG(15) = nrphys KDIAG(16) = nrphys KDIAG(17) = nrphys KDIAG(18) = nrphys KDIAG(19) = nrphys KDIAG(20) = nrphys KDIAG(21) = 1 KDIAG(22) = 1 KDIAG(23) = nrphys KDIAG(24) = nrphys KDIAG(25) = nrphys KDIAG(26) = nrphys KDIAG(27) = 1 KDIAG(28) = 1 KDIAG(29) = 1 KDIAG(30) = 1 KDIAG(31) = 1 KDIAG(32) = 1 KDIAG(33) = 1 KDIAG(34) = 1 KDIAG(35) = 1 KDIAG(36) = 1 KDIAG(37) = 1 KDIAG(38) = 1 KDIAG(39) = 1 KDIAG(40) = nrphys KDIAG(41) = nrphys KDIAG(42) = nrphys KDIAG(43) = nrphys KDIAG(44) = nrphys KDIAG(45) = nrphys KDIAG(46) = 1 KDIAG(47) = nrphys KDIAG(48) = 1 KDIAG(49) = nrphys KDIAG(50) = nrphys KDIAG(51) = 1 KDIAG(52) = 1 KDIAG(53) = 1 KDIAG(54) = Nr KDIAG(55) = Nr KDIAG(56) = Nr KDIAG(57) = Nr KDIAG(58) = nrphys KDIAG(59) = nrphys KDIAG(60) = nrphys KDIAG(61) = nrphys KDIAG(62) = nrphys KDIAG(63) = nrphys KDIAG(64) = nrphys KDIAG(65) = nrphys KDIAG(66) = nrphys KDIAG(67) = 1 KDIAG(68) = 1 KDIAG(69) = nrphys-1 KDIAG(70) = 1 KDIAG(71) = nrphys KDIAG(72) = 1 KDIAG(73) = 1 KDIAG(74) = nrphys KDIAG(75) = nrphys KDIAG(76) = nrphys KDIAG(77) = nrphys KDIAG(78) = nrphys KDIAG(79) = nrphys KDIAG(80) = 1 KDIAG(81) = nrphys KDIAG(82) = 1 KDIAG(83) = 1 KDIAG(84) = 1 KDIAG(85) = 1 KDIAG(86) = nrphys KDIAG(87) = nrphys KDIAG(88) = nrphys KDIAG(89) = nrphys KDIAG(90) = nrphys KDIAG(91) = nrphys KDIAG(92) = nrphys KDIAG(93) = nrphys KDIAG(94) = nrphys KDIAG(95) = nrphys KDIAG(96) = 1 KDIAG(97) = 1 KDIAG(98) = 1 KDIAG(99) = 1 KDIAG(100)= nrphys KDIAG(101)= 1 KDIAG(102)= 1 KDIAG(103)= 1 KDIAG(104)= 1 KDIAG(105)= 1 KDIAG(106)= 1 KDIAG(107)= 1 KDIAG(108)= 1 KDIAG(109)= 1 KDIAG(110)= 1 KDIAG(111)= 1 KDIAG(112)= 1 KDIAG(113)= 1 KDIAG(114)= 1 KDIAG(115)= 1 KDIAG(116)= nrphys KDIAG(117)= nrphys KDIAG(118)= 1 KDIAG(119)= 1 KDIAG(120)= nrphys KDIAG(121)= nrphys KDIAG(122)= nrphys KDIAG(123)= nrphys KDIAG(124)= 1 KDIAG(125)= 1 KDIAG(126)= 1 KDIAG(127)= 1 KDIAG(128)= 1 KDIAG(129)= 1 KDIAG(130)= 1 KDIAG(131)= 1 KDIAG(132)= nrphys KDIAG(133)= nrphys KDIAG(134)= nrphys KDIAG(135)= nrphys KDIAG(136)= nrphys KDIAG(137)= nrphys KDIAG(138)= nrphys KDIAG(139)= nrphys KDIAG(140)= 1 KDIAG(141)= 1 KDIAG(142)= 1 KDIAG(143)= 1 KDIAG(144)= 1 KDIAG(145)= nrphys KDIAG(146)= 1 KDIAG(147)= 1 KDIAG(148)= 1 KDIAG(149)= 1 KDIAG(150)= 1 KDIAG(151)= 1 KDIAG(152)= 1 KDIAG(153)= 1 KDIAG(154)= 1 KDIAG(155)= 1 KDIAG(156)= 1 KDIAG(157)= 1 KDIAG(158)= 1 KDIAG(159)= 1 KDIAG(160)= 1 KDIAG(161)= 1 KDIAG(162)= 1 KDIAG(163)= 1 KDIAG(164)= 1 KDIAG(165)= 1 KDIAG(166)= 1 KDIAG(167)= 1 KDIAG(168)= 1 KDIAG(169)= 1 KDIAG(170)= 1 KDIAG(171)= 1 KDIAG(172)= 1 KDIAG(173)= 1 KDIAG(174)= 1 KDIAG(175)= 1 KDIAG(176)= 1 KDIAG(177)= 1 KDIAG(178)= 1 KDIAG(179)= 1 KDIAG(180)= 1 KDIAG(181)= 1 KDIAG(182)= 1 KDIAG(183)= 1 KDIAG(184)= 1 KDIAG(185)= 1 KDIAG(186)= 1 KDIAG(187)= 1 KDIAG(188)= 1 KDIAG(189)= 0 KDIAG(190)= 1 KDIAG(191)= 1 KDIAG(192)= 1 KDIAG(193)= 1 KDIAG(194)= 1 KDIAG(195)= 1 KDIAG(196)= 1 KDIAG(197)= 1 KDIAG(198)= 1 KDIAG(199)= 1 KDIAG(200)= 1 KDIAG(201)= 1 KDIAG(202)= nrphys KDIAG(203)= nrphys KDIAG(204)= 1 KDIAG(205)= 1 KDIAG(206)= 1 KDIAG(207)= 1 KDIAG(208)= 1 KDIAG(209)= nrphys KDIAG(210)= 1 KDIAG(211)= nrphys KDIAG(212)= nrphys KDIAG(213)= nrphys KDIAG(214)= 1 KDIAG(215)= 1 KDIAG(216)= 1 KDIAG(217)= 1 KDIAG(218)= 1 KDIAG(219)= 1 KDIAG(220)= nrphys KDIAG(221)= nrphys KDIAG(222)= nrphys KDIAG(223)= nrphys KDIAG(224)= nrphys KDIAG(225)= nrphys KDIAG(226)= nrphys KDIAG(227)= nrphys KDIAG(228)= nrphys KDIAG(229)= 1 KDIAG(230)= 1 KDIAG(231)= 1 KDIAG(232)= 1 KDIAG(233)= 1 KDIAG(234)= 1 KDIAG(235)= 1 KDIAG(236)= nrphys KDIAG(237)= nrphys KDIAG(238)= 1 KDIAG(239)= 1 KDIAG(240)= Nr KDIAG(241)= Nr KDIAG(242)= Nr KDIAG(243)= Nr KDIAG(244)= Nr KDIAG(245)= Nr KDIAG(246)= Nr KDIAG(247)= Nr KDIAG(248)= Nr KDIAG(249)= Nr KDIAG(250)= Nr KDIAG(251)= Nr KDIAG(252)= Nr KDIAG(253)= Nr KDIAG(254)= Nr KDIAG(255)= Nr KDIAG(256)= Nr KDIAG(257)= Nr KDIAG(258)= Nr KDIAG(259)= Nr KDIAG(260)= Nr KDIAG(261)= Nr KDIAG(262)= Nr KDIAG(263)= Nr KDIAG(264)= Nr KDIAG(265)= Nr KDIAG(266)= Nr KDIAG(267)= Nr KDIAG(268)= Nr KDIAG(269)= Nr KDIAG(270)= Nr KDIAG(271)= Nr KDIAG(272)= Nr KDIAG(273)= Nr KDIAG(274)= Nr KDIAG(275)= Nr KDIAG(276)= Nr KDIAG(277)= Nr KDIAG(278)= Nr KDIAG(279)= Nr KDIAG(280)= Nr do n=1,ndiagt idiag (n) = 0 ndiag(n) = 0 enddo RETURN END