diagTitle = 'Cloud Top Temperature (when cloudy)  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Cloud Top Temperature Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

      diagTitle = 'Richardson Number  (non-dimensional)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Deep Soil Temperature  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Specific Humidity of the Canopy  (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Soil Moisture Field Capacity Fraction in Root Zone'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Canopy Interception Reservoir (mm)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Snow Depth (mm water equivalent)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Convective Rainfall (liquid precipitate)  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Large-scale Rainfall (liquid precipitate)  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Total Snowfall (solid precipitate)  (mm/day)'

      diagTitle = 'Surface Drag Coef for Scalars (non-dimensional)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Rainwater Infiltration into top soil layer  (mm/s)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Rate of Snow Melt  (mm/s)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Latent Heat Flux from Bare Soil  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Latent Heat Flux from Snow Pack  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

      diagTitle = 'Surface Drag Coef for Winds (non-dimensional)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Leaf Area Index  (%)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Green-ness Index  (%)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged MOISTT  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged TURBT  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged RADLW  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Eddy Diffusivity Coefficient for Momentum  (m^2/s)'

      diagTitle = 'Vertically Averaged RADSW  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged TAVE  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Pressure Thickness of the RAS Subcloud Layer (mb)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Lifting Condensation Level (mb)                 '

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Potential Evapotranspiration (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Total Change in Canopy Temperature (deg/sec)'

      diagTitle = 'Eddy Diffusivity Coefficient for Scalars  (m^2/s)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Change in Canopy Humidity (kg/kg/sec)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Cloud Liquid Water Generated in Turbulence (mg/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Cloud Fraction Generated in Turbulence (fraction)'

      diagTitle = 'Zonal Wind Tendency from Turbulence  (m/s/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Mean Surface Albedo (SW_Up/SW_Down)  0-1 '

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Surface Albedo Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Direct  Beam VIS Surface Albedo (0-1)   '

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Diffuse Beam VIS Surface Albedo (0-1)   '

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Direct  Beam NIR Surface Albedo (0-1)   '

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Diffuse Beam NIR Surface Albedo (0-1)   '

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Large-Scale Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

      diagTitle = 'Meridional Wind Tendency from Turbulence  (m/s/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Convective Available Potential Energy'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Moist Static Energy in RAS subcloud layer'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Saturation Moist Static Energy'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Planetary Albedo (Fraction)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Pressures for Fizhi (mb)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

      diagTitle = 'Zonal Wind Tendency from Moist Proc  (m/s/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Meridional Wind Tendency from Moist Proc  (m/s/day)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Total Precipitation  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Convective Precipitation  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Turbulent Flux of Zonal Momentum  (N/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Turbulent Flux of Meridional Momentum  (N/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Turbulent Flux of Sensible Heat  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Turbulent Flux of Latent Heat  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Neutral Drag Coefficient  (non-dimensional)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Surface Wind Speed  (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Air/Surface Virtual Temperature Difference  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Ground Temperature  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Surface Air Temperature  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Ground Temperature Adjustment (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Ground Specific Humidity  (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Saturation Surface Specific Humidity  (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Ground Temperature used by LW Radiation  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Upward Longwave Radiation at the Ground  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Outgoing Longwave Radiation   (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Clear Sky Outgoing Longwave Radiation  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Clear Sky Net Longwave Flux at the Ground  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Clear Sky Longwave Radiation Heating Rates  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Temperature used by Radiation  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Specific Humidity used by Radiation  (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Convective Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = '3-D Total Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Downward Longwave Radiation at the Ground  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Incident Shortwave Radiation at TOA (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Cloud Optical Depth Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Low-Level ( 1000-700 hPa) Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Surface Evaporation  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Surface Pressure Tendency  (hPa/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Average Zonal Wind (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      CHARACTER*(80) diagTitle

      diagTitle = 'Average Meridional Wind (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Average Potential Temperature (K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertical Velocity  (hPa/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Total Zonal Wind Tendency  (m/s/day)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Total Meridional Wind Tendency  (m/s/day)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Total Temperature Tendency  (K/day)'

      diagTitle = 'Sensible Heat Flux (pos.upwrd)  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Total Specific Humidity Tendency  (g/kg/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Relative Vorticity  (x10^-4  sec^-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Friction Velocity  (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Surface Roughness Z0  (m)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Frequency of Turbulence  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Estimated PBL Depth  (hPa)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Clear Sky Shortwave Heating Rates  (K/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Outgoing Shortwave Radiation  (W/m^2)'

      diagTitle = 'Latent   Heat Flux (pos.upwrd)  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Clear Sky Outgoing Shortwave Radiation  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Cloud Mass Flux  (g/m^2/sec)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Time-Averaged Zonal Wind  (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Time-Averaged Meridional Wind  (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Time-Averaged Temperature  (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Time-Averaged Specific Humidity  (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Time-Averaged Surface Pressure  (hPa)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Sea Level Pressure (mb)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Zonal Wind Surface Stress  (N/m^2)'

      diagTitle = 'Time-Averaged Turbulent Kinetic Energy  (m/s)^2'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Average Surface Pressure (hPa)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

      diagTitle = 'Meridional Wind Surface Stress  (N/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Zonal Wind Gravity Wave Surface Stress (N/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Meridional Wind Gravity Wave Surface Stress (N/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Zonal Wind Gravity Wave Stress at Model Top (N/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle =

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged UWND*SPHU  (m/s g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Vertically Averaged VWND*SPHU  (m/s g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Vertically Averaged UWND*TMPU  (m/s deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Vertically Averaged VWND*TMPU  (m/s deg K)'

      diagTitle = 'Heat Conduction through Sea-Ice  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = '2-D Total Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Total Precipitable Water  (g/cm^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Zonal Wind Interpolated to 2 Meters (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Meridional Wind Interpolated to 2 Meters (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Temperature Interpolated to 2 Meters (deg K)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Specific Humidity Interpolated to 2 Meters (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Zonal Wind Interpolated to 10 Meters (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Meridional Wind Interpolated to 10 Meters (m/s)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Temperature Interpolated to 10 Meters (deg K)'

      diagTitle = 'Net Upward   Longwave Flux at the Ground  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Specific Humidity Interpolated to 10 Meters (g/kg)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Detrainment Cloud Mass Flux  (g/m^2/sec)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Filling of Negative Specific Humidity (g/kg/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged QAVE  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Vertically Averaged QFILL  (mm/day)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Mid-Level (700-400 hPa) Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'High-Level (above 400 hPa) Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Low-Level (1000-700 hPa) Optical Depth'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Low-Level Optical Depth Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Mid-Level (700-400 hPa) Optical Depth'

      diagTitle = 'Net Downward Shortwave Flux at the Ground  (W/m^2)'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Mid-Level Optical Depth Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'High-Level (above 400 hPa) Optical Depth'

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'High-Level Optical Depth Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Non-Precipitating Cloud Fraction  (0-1)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagTitle = 'Cloud Top Pressure (when cloudy)  (hPa)'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

      diagTitle = 'Cloud Top Pressure Counter'

     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      CHARACTER*(80) diagTitle

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Meridional Diffusive Flux of Pot.Temperature'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Vertical Diffusive Flux of Pot.Temperature'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Pot.Temp.   1rst Order Moment Sx'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Pot.Temp.   1rst Order Moment Sy'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Pot.Temp.   1rst Order Moment Sz'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Pot.Temp.   2nd Order Moment Sxx'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Pot.Temp.   2nd Order Moment Syy'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Pot.Temp.   2nd Order Moment Szz'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Pot.Temp.   2nd Order Moment Sxy'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Pot.Temp.   2nd Order Moment Sxz'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Pot.Temp.   2nd Order Moment Syz'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Pot.Temp.   sub-grid variance'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Vertical   Advective Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Zonal      Advective Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Meridional Advective Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Vertical Diffusive Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Zonal      Diffusive Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Meridional Diffusive Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Vertical Diffusive Flux of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Filling of Negative Values of '//locName

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = locName//'1rst Order Moment Sx'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = locName//'1rst Order Moment Sy'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = locName//'1rst Order Moment Sz'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = locName//'2nd Order Moment Sxx'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = locName//'2nd Order Moment Syy'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = locName//'2nd Order Moment Szz'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = locName//'2nd Order Moment Sxy'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = locName//'2nd Order Moment Sxz'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = locName//'2nd Order Moment Syz'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = locName//'sub-grid variance'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      CHARACTER*(80) diagTitle

       diagTitle = 'Vertical   Advective Flux of Pot.Temperature'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Zonal      Advective Flux of Pot.Temperature'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Meridional Advective Flux of Pot.Temperature'

     I      diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

       diagTitle = 'Vertical Diffusive Flux of Pot.Temperature'

     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

       diagTitle = 'Zonal      Diffusive Flux of Pot.Temperature'

      CHARACTER*(80) diagTitle

        diagTitle = 'GGL90 sub-grid turbulent kinetic energy'

     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'Mixing length scale              '

     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'Prandtl number used in GGL90'

     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'GGL90 eddy viscosity at U-point'

     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'GGL90 eddy viscosity at V-point'

     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'GGL90 diffusion coefficient for temperature'

     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'K_13 element (U.point, Z.dir) of GM-Redi tensor'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'K_23 element (V.point, Z.dir) of GM-Redi tensor'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'K_31 element (W.point, X.dir) of GM-Redi tensor'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'K_32 element (W.point, Y.dir) of GM-Redi tensor'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'K_33 element (W.point, Z.dir) of GM-Redi tensor'

     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

        diagTitle = 'GM Bolus transport stream-function : X component'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'GM Bolus transport stream-function : Y component'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'Redi Off-diagonal Temperature flux: X component'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'Redi Off-diagonal Temperature flux: Y component'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'Zonal Mass-Weight Bolus Transp of Pot Temp'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'Meridional Mass-Weight Bolus Transp of Pot Temp'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'eddy-induced stress X-comp. estimated from Kwx'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )

        diagTitle = 'eddy-induced stress Y-comp. estimated from Kwy'

     I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )