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

      diagTitle = 'meridional surface wind stress, >0 increases vVel'

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

      diagTitle = 'zonal 10-m wind speed, >0 eastward'

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

      diagTitle = 'meridional 10-m wind speed, >0 northward'

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

      diagTitle = '10-m wind speed modulus ( >= 0 )'

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

      diagTitle = 'surface (2-m) air temperature'

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

      diagTitle = 'surface (2-m) specific humidity'

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

      diagTitle = 'evaporation, > 0 increases salinity'

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

      diagTitle = 'precipitation, > 0 decreases salinity'

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

      diagTitle = 'snow precipitation, > 0 decreases salinity'

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

      diagTitle = 'net upward freshwater flux, > 0 increases salinity'

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

      diagTitle = 'atmospheric pressure field'

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

      diagTitle = 'river runoff, > 0 decreases salinity'

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

      diagTitle = 'river runoff temperature'

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

      diagTitle = 'upward salt flux, > 0 decreases salinity'

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

      CHARACTER*(80) diagTitle

      diagTitle = 'Sensible heat flux into ocean, >0 increases theta'

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

      diagTitle = 'Latent heat flux into ocean, >0 increases theta'

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

      diagTitle = 'Net upward longwave radiation, >0 decreases theta'

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

      diagTitle = 'Net upward shortwave radiation, >0 decreases theta'

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

      diagTitle = 'Downward longwave radiation, >0 increases theta'

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

      diagTitle = 'Downward shortwave radiation, >0 increases theta'

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

      diagTitle = 'Net upward heat flux (turb+rad), >0 decreases theta'

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

      diagTitle = 'zonal surface wind stress, >0 increases uVel    '

      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

      diagTitle = 'Frazil temperature forcing, > 0 increases theta'

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

      CHARACTER*(80) diagTitle

       WRITE(diagTitle,'(2A)') locName(1:ilnb),

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

     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 )