COMMON/surfdat/albedodat(ngridmx),                                &
     &   phisfi(ngridmx),                                               &
     &   zmea(ngridmx),zstd(ngridmx),phitop(ngridmx),                   &
     &   zsig(ngridmx),zgam(ngridmx),zthe(ngridmx),ttop,kp(ngridmx),    &
     &   p00,alb_n2b,alb_n2a,alb_ch4,alb_co,alb_tho,emis_n2b,emis_n2a,  &
     &   emis_ch4,emis_co,emis_tho,emis_tho_eq,alb_tho_eq,alb_ch4_eq,   &
     &   ITN2,ITCH4,ITH2O,ITN2d,ITCH4d,ITH2Od,alb_ch4_s,albspe,emispe,  &            
     &   alb_tho_spe,emis_tho_spe
  

      real albedodat ! albedo of bare ground
      real phisfi ! geopotential at ground level
      real phitop,ttop
      real zmea,zstd,zsig,zgam,zthe
      real kp ! TB ref pressure 
      real p00  ! globalaverage(kp)
! surface properties ! TB16
      real alb_n2b
      real alb_n2a
      real alb_ch4
      real alb_co
      real alb_tho
      real alb_tho_eq
      real alb_tho_spe
      real emis_n2b
      real emis_n2a
      real emis_ch4
      real emis_co
      real emis_tho
      real emis_tho_eq
      real emis_tho_spe
      real alb_ch4_eq
      real alb_ch4_s
      real albspe
      real emispe
      real ITN2,ITN2d
      real ITCH4,ITCH4d
      real ITH2O,ITH2Od