& chlon(1,bi,bj),igrd(1,bi,bj),nchptot(bi,bj),nchpland(bi,bj),

      INTEGER nchptot(nSxIn,nSyIn),nchpland(nSxIn,nSyIn)

     & turbStart, nchp,nchptot,nchpland,

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

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

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

     .   tgz, landtype, nchptot, nchpland, ityp

      integer nchptot(Nsx,Nsy),nchpland(Nsx,Nsy)

     &         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,

         nchptot(bi,bj) = nchpland(bi,bj)

               nchptot(bi,bj)  = nchptot(bi,bj) + 1

               igrd (nchptot(bi,bj),bi,bj) = i + (j-1)*im

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

               chfr (nchptot(bi,bj),bi,bj) = tilefrac(i,j,k,bi,bj)

               chlon(nchptot(bi,bj),bi,bj) = lons(i,j,bi,bj)

               chlt (nchptot(bi,bj),bi,bj) = lats(i,j,bi,bj)

     &                 ', Total # of Tiles=', nchptot(bi,bj)

      INTEGER nchptot(nSx,nSy), nchpland(nSx,nSy)

     & maxtyp,nchp,nchptot,nchpland,lons,lats,surftype,tilefrac,

     .                 temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

     .                    temptile,nchp,nchptot(bi,bj))

      do i =1,nchptot

      do i =1,nchptot

      call MSC2GRD(igrd,chfr,dumsc(1,L),nchp,nchptot,fracland,

      call MSC2GRD(igrd,chfr,dvmsc(1,L),nchp,nchptot,fracland,

      call MSC2GRD(igrd,chfr,dtmsc(1,L),nchp,nchptot,fracland,

      call MSC2GRD(igrd,chfr,dqmsc(1,L,nt),nchp,nchptot,fracland,

      call MSC2GRD(igrd,chfr,  tke(1,L),nchp,nchptot,fracland,

      call MSC2GRD(igrd,chfr, fccmsc(1,L),nchp,nchptot,fracland,

      call MSC2GRD(igrd,chfr,qliqmsc(1,L),nchp,nchptot,fracland,

      integer nchp,nchptot,nchplnd

      call MSC2GRD(igrd,chfr,tground ,nchp,nchptot,fracland,tgz ,im,jm)

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

      numstrips   = ((nchptot-1)/istrip) + 1

      call GRD2MSC(pz(1,1),im,jm,igrd,pmsc,nchp,nchptot)

      call GRD2MSC(tgz,im,jm,igrd,tground,nchp,nchptot)

      call GRD2MSC(tmpdiag,im,jm,igrd,alwcoeff,nchp,nchptot)

      call GRD2MSC(tmpdiag,im,jm,igrd,blwcoeff,nchp,nchptot)

      call GRD2MSC(tmpdiag,im,jm,igrd,pardiff,nchp,nchptot)

      call GRD2MSC(tmpdiag,im,jm,igrd,pardirct,nchp,nchptot)

      call GRD2MSC(tmpdiag,im,jm,igrd,netsw,nchp,nchptot)

      call GRD2MSC(tmpdiag,im,jm,igrd,netlw,nchp,nchptot)

      call GRD2MSC(thkz,im,jm,igrd,icethk,nchp,nchptot)

      call GRD2MSC(rainlsp,im,jm,igrd,lsprec,nchp,nchptot)

      call GRD2MSC(rainconv,im,jm,igrd,cnvprec,nchp,nchptot)

      call GRD2MSC(snowfall,im,jm,igrd,snowprec,nchp,nchptot)

      do i = 1,nchptot

      do i =1,nchptot

      do I =1,nchptot

      do i = nchplnd+1,nchptot

      do i = nchplnd+1,nchptot

      do i =1,nchptot

       call MSC2GRD (igrd,chfr,qice,nchp,nchptot,fracland,tmpdiag,im,jm)

       call STRIPIT  (z0,stz0,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (tground,th(1,nlay+1),nchptot,nchp,istrip,1,nn)

       call STRIPIT  (pmsc,pe(1,nlay+1),nchptot,nchp,istrip,1,nn)

       call STRIPIT  (tke,qq,nchptot,nchp,istrip,nlay-1,nn)

       call STRIPIT  (ctmt,ctsave,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (xxmt,xxsave,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (yymt,yysave,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (zetamt,zetasave,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (xlmt,xlsave,nchptot,nchp,istrip,nlay,nn)

       call STRIPIT  (khmt,khsave,nchptot,nchp,istrip,nlay,nn)

       call STRIPITINT (water,stwatr,nchptot,nchp,istrip,1,nn)

       call STRIPITINT (igrd,igrdstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (chfr,chfrstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (icethk,icest,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (pardiff,pardf,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (pardirct,pardr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (chlt,lats,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (chlon,lons,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (lsprec,rainls,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (cnvprec,raincon,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (snowprec,newsnow,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (netsw,swnet,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (netlw,lwstrip,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (alwcoeff,alwrad,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (blwcoeff,blwrad,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (alai,laistrip,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (agrn,grnstrip,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (z2ch,z2str,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (sqscat,scatstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (rsoil1,rs1str,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (rsoil2,rs2str,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (rdc,rdcstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (u2fac,u2fstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (shg,shgstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (dqsdt,dqsdtstr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  ( qice, qicestr,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (dqice,dqicestr,nchptot,nchp,istrip,1,nn)

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

       call STRIPIT  (tground,tc,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (tdeep,td,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (qground,qa,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (swetshal,swet1,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (swetroot,swet2,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (swetdeep,swet3,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (snodep,snowdepth,nchptot,nchp,istrip,1,nn)

       call STRIPIT  (capac,capacity,nchptot,nchp,istrip,1,nn)

      call PASTIT (qq,tke,istrip,nchp,nchptot,nlay,nn)

      call PASTIT (ctsave,ctmt,istrip,nchp,nchptot,1,nn)

      call PASTIT (xxsave,xxmt,istrip,nchp,nchptot,1,nn)

      call PASTIT (yysave,yymt,istrip,nchp,nchptot,1,nn)

      call PASTIT (zetasave,zetamt,istrip,nchp,nchptot,1,nn)

      call PASTIT (xlsave,xlmt,istrip,nchp,nchptot,nlay,nn)

      call PASTIT (khsave,khmt,istrip,nchp,nchptot,nlay,nn)

      call PASTIT (qliq  ,qliqmsc,istrip,nchp,nchptot,nlay,nn)

      call PASTIT (turbfcc,fccmsc,istrip,nchp,nchptot,nlay,nn)

      call PASTIT (tc,tground,istrip,nchp,nchptot,1,nn)

      call PASTIT (td,tdeep,istrip,nchp,nchptot,1,nn)

      call PASTIT (qa,qground,istrip,nchp,nchptot,1,nn)

      call PASTIT (swet1,swetshal,istrip,nchp,nchptot,1,nn)

      call PASTIT (swet2,swetroot,istrip,nchp,nchptot,1,nn)

      call PASTIT (swet3,swetdeep,istrip,nchp,nchptot,1,nn)

      call PASTIT (capacity,capac,istrip,nchp,nchptot,1,nn)

      call PASTIT (snowdepth,snodep,istrip,nchp,nchptot,1,nn)

      call PASTIT (tends,dumsc(1,l),istrip,nchp,nchptot,1,nn)

      call PASTIT (tends,dvmsc(1,l),istrip,nchp,nchptot,1,nn)

      call PASTIT (tends,dtmsc(1,l),istrip,nchp,nchptot,1,nn)

      call PASTIT (tends,dqmsc(1,l,1),istrip,nchp,nchptot,1,nn)

     &  turbStart(bi,bj), nchp,nchptot,nchpland,

      do n = 1,nchptot

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

     &  nchptot,fracg,tmpij,im,jm)

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

     &               nchp,nchptot(bi,bj))

        call GETEMISS(fracl,im2,jm2,nchp,nchptot(bi,bj),nSx,nSy,bi,bj,

       do i = 1,nchptot(bi,bj)

     &                    nchp,nchptot(bi,bj),fracl,tmpij,im2,jm2)

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

      INTEGER sec,month,day,im,jm,nchp,nchptot,nchpland,nSx,nSy,bi,bj

      SUBROUTINE GETEMISS(fracg,im,jm,nchp,nchptot,nSx,nSy,bi,bj,

      INTEGER im,jm,nchp,nchptot,nSx,nSy,bi,bj

     &         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,

         nchptot(bi,bj) = nchpland(bi,bj)

               nchptot(bi,bj)  = nchptot(bi,bj) + 1

               igrd (nchptot(bi,bj),bi,bj) = i + (j-1)*im

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

               chfr (nchptot(bi,bj),bi,bj) = tilefrac(i,j,k,bi,bj)

               chlon(nchptot(bi,bj),bi,bj) = lons(i,j,bi,bj)

               chlt (nchptot(bi,bj),bi,bj) = lats(i,j,bi,bj)

     &                 ', Total # of Tiles=', nchptot(bi,bj)

      INTEGER nchptot(nSx,nSy), nchpland(nSx,nSy)

     & maxtyp,nchp,nchptot,nchpland,lons,lats,surftype,tilefrac,

      INTEGER nchptot(nSxin,nSyin),nchpland(nSxin,nSyin)

     & turbStart, nchp,nchptot,nchpland,

     &         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,

         nchptot(bi,bj) = nchpland(bi,bj)

               nchptot(bi,bj)  = nchptot(bi,bj) + 1

               igrd (nchptot(bi,bj),bi,bj) = i + (j-1)*im

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

               chfr (nchptot(bi,bj),bi,bj) = tilefrac(i,j,k,bi,bj)

               chlon(nchptot(bi,bj),bi,bj) = lons(i,j,bi,bj)

               chlt (nchptot(bi,bj),bi,bj) = lats(i,j,bi,bj)

     &                 ', Total # of Tiles=', nchptot(bi,bj)

      INTEGER nchptot(nSx,nSy), nchpland(nSx,nSy)

     & maxtyp,nchp,nchptot,nchpland,lons,lats,surftype,tilefrac,