CHARACTER*(80) diagTitle

      CHARACTER*(80) diagTitle

      diagTitle = 'Pressure Thickness of Layers on Fizhi Grid'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Sea-Ice enthalphy, 2nd  layer (mass weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Sea-Ice Albedo [0-1] (area weighted average)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'snow age over Sea-Ice'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'snow precip. (+=dw) over Sea-Ice (area weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'net heat flux from the Atmosphere (+=dw)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'fresh-water flux to the Atmosphere (+=up)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'heat flux out of the ocean (+=up)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'fresh-water flux out of the ocean (+=up)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'salt flux out of the ocean (+=up)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'ocean mixed layer temperature'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'ocean mixed layer salinity'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      CHARACTER*(80) diagTitle

        diagTitle = 'Sea-Ice fraction  [0-1]'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Sea-Ice thickness (area weighted average)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Snow thickness over Sea-Ice (area weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Surface Temperature over Sea-Ice (area weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Sea-Ice Temperature, 1srt layer (area weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Sea-Ice Temperature, 2nd  layer (area weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Sea-Ice enthalphy, 1srt layer (mass weighted)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Humidity Change due to Convection (g/kg/sec)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Relative Humidity (fraction)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Surface Evaporation  (mm/day)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Convective Precipitation  (mm/day)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Large Scale Precipitation (g/m2/sec)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Cloud Top Pressure Counter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Surface Air Temperature  (deg K)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Energy of Precipitation (snow, rain Temp) in (J/g)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Surf. Temp Change after 1 implicit time step'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      CHARACTER*(80) diagTitle

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Temp Change due to Convection (K/sec)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I              diagName, diagCode, diagUnits, diagTitle, myThid )

       WRITE(diagTitle,'(2A)') 'Merid Mass-Weighted Transp of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid )

       WRITE(diagTitle,'(2A)') 'Vert  Mass-Weighted Transp of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid )

       WRITE(diagTitle,'(2A)') 'Vertical   Advective Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       WRITE(diagTitle,'(2A)') 'Zonal      Advective Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       WRITE(diagTitle,'(2A)') 'Meridional Advective Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       WRITE(diagTitle,'(2A)') 'Vertical Diffusive Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       WRITE(diagTitle,'(2A)') 'Zonal      Diffusive Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       WRITE(diagTitle,'(2A)') 'Meridional Diffusive Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       WRITE(diagTitle,'(2A)') 'Vertical Diffusive Flux of ',

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

      CHARACTER*(80) diagTitle

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

     I              diagName, diagCode, diagUnits, diagTitle, myThid )

       WRITE(diagTitle,'(2A)') 'Zonal Mass-Weighted Transp of ',

      CHARACTER*(80) diagTitle

        diagTitle = 'KPP vertical eddy viscosity coefficient'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Vertical diffusion coefficient for salt & tracers'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Vertical diffusion coefficient for heat'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Nonlocal transport coefficient'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'KPP boundary layer depth, bulk Ri criterion'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Mixed layer depth, dT=.8degC density criterion'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Short-wave flux fraction penetrating mixing layer'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I     diagName, diagCode, diagUnits, diagTitle,

      CHARACTER*(*) diagTitle

          tdiag(diagNum) = diagTitle

        tdiag(diagNum) = diagTitle

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Perturbation of Sfc (Pa^2,m^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Potential Temperature (degC,K)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Sea Surface Temperature (degC,K)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Sea Surface Salinity (g/kg)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Salt anomaly (=SALT-35; g/kg)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Component of Velocity (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Component of Velocity (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Component of Velocity (r_units/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Potential Temperature (K^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle ='Product of horizontal Comp of velocity (cell center)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Density Anomaly (=Rho-rhoConst)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Zonal Transport of Density'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Meridional Transport of Density'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Vertical Transport of Potential Density'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'meridional surf. wind stress, >0 increases vVel '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'net surface heat flux, >0 increases theta       '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'surface temperature relaxation, >0 increases theta'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'net surface salt flux, >0 increases salt        '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'surface salinity relaxation, >0 increases salt  '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Atmospheric Pressure (Pa)     '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Surface (Pressure, Height) (Pa,m)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      CHARACTER*(80) diagTitle

      diagTitle = 'Perturbation of Surface (pressure, height) (Pa,m)'

      CHARACTER*(80) diagTitle

        diagTitle = 'Temperature Tendency due to Shapiro Filter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Specific Humidity Tendency due to Shapiro Filter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Zonal Wind Tendency due to Shapiro Filter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

        diagTitle = 'Meridional Wind Tendency due to Shapiro Filter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Vertical   Advective Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Zonal      Advective Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Meridional Advective Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Vertical Diffusive Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Zonal      Diffusive Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Meridional Diffusive Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Vertical Diffusive Flux of'//locName

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

      CHARACTER*(80) diagTitle

       diagTitle = 'Vertical   Advective Flux of Pot.Temperature'

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Zonal      Advective Flux of Pot.Temperature'

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Meridional Advective Flux of Pot.Temperature'

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Vertical Diffusive Flux of Pot.Temperature'

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Zonal      Diffusive Flux of Pot.Temperature'

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Meridional Diffusive Flux of Pot.Temperature'

     I              diagName, diagCode, diagUnits, diagTitle, myThid ) 

       diagTitle = 'Vertical Diffusive Flux of Pot.Temperature'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Canopy Interception Reservoir (mm)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Snow Depth (mm water equivalent)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Leaf Area Index  (%)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Green-ness Index  (%)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Lifting Condensation Level (mb)                 '

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Surface Albedo Counter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Convective Available Potential Energy'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Moist Static Energy in RAS subcloud layer'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Saturation Moist Static Energy'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Planetary Albedo (Fraction)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Pressures for Fizhi (mb)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle = 'Relative Humidity (fraction)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  = 'Total Precipitation  (mm/day)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  = 'Convective Precipitation  (mm/day)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  = 'Ground Temperature  (deg K)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  = 'Surface Air Temperature  (deg K)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  = 'Ground Temperature Adjustment (deg K)'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

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

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  =

     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagTitle  = 'Cloud Optical Depth Counter'

     I          diagName, diagCode, diagUnits, diagTitle, myThid )