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

      integer nchptot(Nsxin,Nsyin),nchpland(Nsxin,Nsyin)

     . 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)

          print *,'No of Total Tiles for bi=',bi,': ',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))

      integer nchp,nchptot,nchplnd

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

      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,

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

      do n=1,nchptot

      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 PSTBITINT(types,qdiaglsm(1,1),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(chfrstr,qdiaglsm(1,2),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(lats,qdiaglsm(1,3),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(lons,qdiaglsm(1,4),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(tc,qdiaglsm(1,6),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(td,qdiaglsm(1,7),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(qa,qdiaglsm(1,8),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(swet1,qdiaglsm(1,9),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(swet2,qdiaglsm(1,10),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(swet3,qdiaglsm(1,11),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(capacity,qdiaglsm(1,12),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(snowdepth,qdiaglsm(1,13),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(eturb,qdiaglsm(1,14),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(hsturb,qdiaglsm(1,15),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(cd,qdiaglsm(1,16),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(laistrip,qdiaglsm(1,17),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(grnstrip,qdiaglsm(1,18),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(eint,qdiaglsm(1,19),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(esoi,qdiaglsm(1,20),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(eveg,qdiaglsm(1,21),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(esno,qdiaglsm(1,22),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg1,qdiaglsm(1,23),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg2,qdiaglsm(1,24),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg3,qdiaglsm(1,25),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg4,qdiaglsm(1,26),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg5,qdiaglsm(1,27),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg6,qdiaglsm(1,28),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg7,qdiaglsm(1,29),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg8,qdiaglsm(1,30),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(strdg9,qdiaglsm(1,31),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(smelt,qdiaglsm(1,32),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(gdrain,qdiaglsm(1,33),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(runsrf,qdiaglsm(1,34),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(fwsoil,qdiaglsm(1,35),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(evpot,qdiaglsm(1,36),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(stt2m,qdiaglsm(1,37),istrip,nchp,nchptot,1,nn)

      call PSTBMPIT(stq2m,qdiaglsm(1,38),istrip,nchp,nchptot,1,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)

     .  nchp,nchptot,nchpland,

      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,

        call GRD2MSC(fraci,im2,jm2,igrd,ficetile,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

      do n = 1,nchptot

     .         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)

         print *,'Number of Total Tiles for bi=',bi,': ',nchptot(bi,bj)

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

     . maxtyp,nchp,nchptot,nchpland,lons,lats,surftype,tilefrac,igrd,

      integer nchptot(Nsxin,Nsyin),nchpland(Nsxin,Nsyin)

     . 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)

         print *,'Number of Total Tiles for bi=',bi,': ',nchptot(bi,bj)

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

     . maxtyp,nchp,nchptot,nchpland,lons,lats,surftype,tilefrac,igrd,