COMMON/surfdat/albedodat(ngridmx),                                &
     &   phisfi(ngridmx),albedice(2),emisice(2),emissiv,                &
     &   TESice_Ncoef,TESice_Scoef,                                     &
     &   iceradius(2) , dtemisice(2),                                   &
     &   zmea(ngridmx),zstd(ngridmx),                                   &
     &   zsig(ngridmx),zgam(ngridmx),zthe(ngridmx)

      COMMON/surfdatl/TESicealbedo

      real albedodat ! albedo of bare ground
      real phisfi ! geopotential at ground level
      real albedice ! default albedo for ice (1: North H. 2: South H.)
      real emisice ! ice emissivity; 1:Northern hemisphere 2:Southern hemisphere
      real emissiv ! emissivity of bare ground
      logical TESicealbedo ! use TES ice cap albedoes (if set to .true.)
      real TESice_Ncoef ! coefficient for TES ice albedo in Northern hemisphere
      real TESice_Scoef ! coefficient for TES ice albedo in Southern hemisphere
      real iceradius , dtemisice
      real zmea,zstd,zsig,zgam,zthe