.       (CSOIL(ChNo)*67.73*deepfac(ityp(ChNo)))

        CALL MNC_CW_I_W('I',bnam,0,0,'ityp', ityp, myThid)

      subroutine CHPPRM(nymd,nhms,mxchps,nchp,chlt,ityp,alai,

       integer nymd,nhms,nchp,mxchps,ityp(mxchps)

      zoch(i) = vgz0(k2,ityp(i))*facp + vgz0(k1,ityp(i))*facm

      rdcs(i) = vgrd(k2,ityp(i))*facp + vgrd(k1,ityp(i))*facm

      rootl   = vgrt(k2,ityp(i))*facp + vgrt(k1,ityp(i))*facm

      vroot  = rootl * vgroca(ityp (i))

      rsl1(i) = vgrdrs (ityp (i)) / (rootl * vgrotd (ityp (i)))

      rsl2(i) = alphaf / vgrotd (ityp (i))

      scat =       agrn(i) *(vgtr11(ityp(i)) + vgrf11(ityp(i)))

     &     + (1. - agrn(i))*(vgtr12(ityp(i)) + vgrf12(ityp(i)))

      d       = vgdd(k2,ityp(i))*facp + vgdd(k1,ityp(i))*facm

      ufac(i) = log( (vgz2(ityp(i)) - d) / zoch(i) ) 

      z2ch(i) = vgz2(ityp (i))

      subroutine EMISSIVITY (snowdepth,fraci,nchp,numpts,ityp,newemis)

      integer   ityp(nchp)

     .   igrd,ityp,chfr,snodep,ficetile,emiss)

       if(ityp(i).le.10)then

        newemis(i, 1) = (emis( 1,ityp(i))+emis(2,ityp(i)))/2.

        newemis(i, 2) = (emis( 2,ityp(i))+emis(3,ityp(i)))/2.

        newemis(i, 3) = (emis( 4,ityp(i))+emis(5,ityp(i)))/2.

        newemis(i, 4) =  emis( 6,ityp(i))

        newemis(i, 5) =  emis( 7,ityp(i))

        newemis(i, 6) =  emis( 8,ityp(i))

        newemis(i, 7) =  emis( 9,ityp(i))

        newemis(i, 8) = (emis(10,ityp(i))+emis(11,ityp(i)))/2.

        newemis(i, 9) =  emis(12,ityp(i))

        newemis(i,10) =  emis( 4,ityp(i))

         fac = snowdepth(i) / (snowdepth(i) + snwmid(ityp(i)))

     .        nchp,nchptot,nchpland,nSx,nSy,bi,bj,igrd,ityp,chfr,chlt,

      integer igrd(nchp,nSx,nSy),ityp(nchp,nSx,nSy)

       call GETLGR (sec,month,day,chlt,ityp,nchpland(bi,bj),

     .  agrn(1,bi,bj),zenith,snodep(1,bi,bj),ityp(1,bi,bj),nchpland)

     .   igrd,ityp,chfr,snowdep,fraci,emiss)

      integer igrd(nchp,nSx,nSy), ityp(nchp,nSx,nSy)

      call EMISSIVITY(snowdep(1,bi,bj),fraci,nchp,nchptot,ityp(1,bi,bj),

     .    nchptot(bi,bj),nchpland(bi,bj),nSx,nSy,bi,bj,igrd,ityp,

     .  albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlt,chlon,

     . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,

     . chfr(1,bi,bj),ityp(1,bi,bj),tcanopy(1,bi,bj),tdeep(1,bi,bj),

      integer ityp(nchp,Nsxin,Nsyin),igrd(nchp,Nsxin,Nsyin)

                  ityp (nchpland(bi,bj),bi,bj) = surftype(i,j,k,bi,bj)

                  ityp (nchptot(bi,bj),bi,bj) = surftype(i,j,k,bi,bj)

      integer igrd(nchp,Nsx,Nsy),ityp(nchp,Nsx,Nsy)

     . igrd,ityp,chfr,chlt,chlon)

     .   tgz, landtype, nchptot, nchpland, ityp

      integer ityp(nchp,Nsx,Nsy)

     .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)

     .  nchp,nchptot,nchpland,chfr,chlat,chlon,igrd,ityp,

      integer nchp,igrd(nchp),nchptot,nchpland,ityp(nchp)

     . nchp,chlat,chlon,igrd,nchptot,nchpland,chfr,ityp,

     5 ,nchp,nchptot,nchplnd,chfr,chlt,chlon,igrd,ityp,alai,agrn,thkz

      integer igrd(nchp),ityp(nchp)

      call CHPPRM(nymd,nhms,nchp,nchplnd,chlt,ityp,alai,

       if((ityp(i).eq.100).and.(icethk(i).eq.0. ))water(i) = 1

       else if (ityp(i).eq.100) then

       call STRIPITINT (ityp,types,nchptot,nchp,istrip,1,nn)

     . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,

      integer ityp(nchp,Nsxin,Nsyin),igrd(nchp,Nsxin,Nsyin)

                  ityp (nchpland(bi,bj),bi,bj) = surftype(i,j,k,bi,bj)

                  ityp (nchptot(bi,bj),bi,bj) = surftype(i,j,k,bi,bj)

      integer igrd(nchp,Nsx,Nsy),ityp(nchp,Nsx,Nsy)

     . ityp,chfr,chlt,chlon)

     .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)

      subroutine EMISSIVITY (snowdepth,fraci,nchp,numpts,ityp,newemis)

      integer   ityp(nchp)

     .   igrd,ityp,chfr,snodep,ficetile,emiss)

       if(ityp(i).le.10)then

        newemis(i, 1) = (emis( 1,ityp(i))+emis(2,ityp(i)))/2.

        newemis(i, 2) = (emis( 2,ityp(i))+emis(3,ityp(i)))/2.

        newemis(i, 3) = (emis( 4,ityp(i))+emis(5,ityp(i)))/2.

        newemis(i, 4) =  emis( 6,ityp(i))

        newemis(i, 5) =  emis( 7,ityp(i))

        newemis(i, 6) =  emis( 8,ityp(i))

        newemis(i, 7) =  emis( 9,ityp(i))

        newemis(i, 8) = (emis(10,ityp(i))+emis(11,ityp(i)))/2.

        newemis(i, 9) =  emis(12,ityp(i))

        newemis(i,10) =  emis( 4,ityp(i))

         fac = snowdepth(i) / (snowdepth(i) + snwmid(ityp(i)))

     .        nchp,nchptot,nchpland,nSx,nSy,bi,bj,igrd,ityp,chfr,chlt,

      integer igrd(nchp,nSx,nSy),ityp(nchp,nSx,nSy)

       call GETLGR (sec,month,day,chlt,ityp,nchpland(bi,bj),

     .  agrn(1,bi,bj),zenith,snodep(1,bi,bj),ityp(1,bi,bj),nchpland)

     .   igrd,ityp,chfr,snowdep,fraci,emiss)

      integer igrd(nchp,nSx,nSy), ityp(nchp,nSx,nSy)

      call EMISSIVITY(snowdep(1,bi,bj),fraci,nchp,nchptot,ityp(1,bi,bj),

     .    nchptot(bi,bj),nchpland(bi,bj),nSx,nSy,bi,bj,igrd,ityp,

     . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,

     . chfr(1,bi,bj),ityp(1,bi,bj),tcanopy(1,bi,bj),tdeep(1,bi,bj),

      integer ityp(nchp,Nsxin,Nsyin),igrd(nchp,Nsxin,Nsyin)

                  ityp (nchptot(bi,bj),bi,bj) = surftype(i,j,k,bi,bj)

      integer igrd(nchp,Nsx,Nsy),ityp(nchp,Nsx,Nsy)

     . ityp,chfr,chlt,chlon)

     .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)

     5 ,nchp,nchptot,nchplnd,chfr,chlt,chlon,igrd,ityp,alai,agrn,thkz

      integer igrd(nchp),ityp(nchp)

      call CHPPRM(nymd,nhms,nchp,nchplnd,chlt,ityp,alai,

       if((ityp(i).eq.100).and.(icethk(i).eq.0. ))water(i) = 1

       else if (ityp(i).eq.100) then

       call STRIPITINT (ityp,types,nchptot,nchp,istrip,1,nn)