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

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

      do i=1,nchp

         do i = 1,nchp

         do i = 1,nchp

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stuflux,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stvflux,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (sttflux,igrd,chfrstr,istrip,nchp,nn,

      integer nchp,nchptot,nchplnd

      call DIAG_VEGTILE_FILL (stqflux,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (sri,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (skh,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (skm,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (sct,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (scu,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (swinds,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stuflux(1,nlay),igrd,chfrstr,istrip,nchp,

      call DIAG_VEGTILE_FILL (stvflux(1,nlay),igrd,chfrstr,istrip,nchp,

      call DIAG_VEGTILE_FILL (sustar,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (sz0,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (frqtrb,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (pbldpth,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stu2m,igrd,chfrstr,istrip,nchp,nn,

      _RL ctmt(nchp),xxmt(nchp),yymt(nchp),zetamt(nchp)

      call DIAG_VEGTILE_FILL (stv2m,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stt2m,igrd,chfrstr,istrip,nchp,nn,

      _RL xlmt(nchp,nlay),khmt(nchp,nlay),tke(nchp,nlay)

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

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stu10m,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stv10m,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (stt10m,igrd,chfrstr,istrip,nchp,nn,

      _RL tcanopy(nchp),tdeep(nchp),ecanopy(nchp),swetshal(nchp)

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      _RL swetroot(nchp),swetdeep(nchp),snodep(nchp),capac(nchp)

      call DIAG_VEGTILE_FILL (td,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (qa,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (swet1,igrd,chfrstr,istrip,nchp,nn,

      _RL chfr(nchp),chlt(nchp),chlon(nchp)

      call DIAG_VEGTILE_FILL (swet2,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (swet3,igrd,chfrstr,istrip,nchp,nn,

      integer igrd(nchp),ityp(nchp)

      call DIAG_VEGTILE_FILL (capacity,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (snowdepth,igrd,chfrstr,istrip,nchp,nn,

      _RL alai(nchp),agrn(nchp),thkz(im,jm)

      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)

      _RL alwcoeff(nchp),blwcoeff(nchp)

      _RL netsw(nchp)

      _RL cnvprec(nchp),lsprec(nchp)

      _RL snowprec(nchp)

      _RL pardiff(nchp),pardirct(nchp)

      _RL pmsc(nchp)

      _RL netlw(nchp)

      _RL sqscat(nchp), rsoil1(nchp)

      _RL rsoil2(nchp)

      _RL rdc(nchp),u2fac(nchp)

      _RL z2ch(nchp)

      _RL zoch(nchp),cdrc(nchp)

      _RL cdsc(nchp)

      _RL dqsdt(nchp)

      _RL tground(nchp),qground(nchp)

      _RL utility(nchp)

      _RL    qice(nchp)

      _RL   dqice(nchp)

      _RL dumsc(nchp,nlay),dvmsc(nchp,nlay)

      _RL dtmsc(nchp,nlay),dqmsc(nchp,nlay,ntracers)

      _RL shg(nchp),z0(nchp),icethk(nchp)

      integer water(nchp)

      _RL qliqmsc(nchp,nlay),fccmsc(nchp,nlay)

      _RL qdiaglsm(nchp,ndlsm)

      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)

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

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

       call STRIP2TILE(uz,igrd,u,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(vz,igrd,v,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(tz,igrd,th,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(qz(1,1,1,1),igrd,sh,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(dpres,igrd,dpstr,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(plz,igrd,p,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(plze,igrd,pe,nchp,ijall,istrip,nlay+1,nn)

       call STRIP2TILE(pkz,igrd,pk,nchp,ijall,istrip,nlay,nn)

       call STRIP2TILE(pkht,igrd,pke,nchp,ijall,istrip,nlay+1,nn)

       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 DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

       call DIAG_VEGTILE_FILL (stdiag,igrd,chfrstr,istrip,nchp,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 DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (grnstrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (laistrip,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (pardr,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (pardf,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (blwrad,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (dhsdtc,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (dedtc,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (dhsdqa,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (dedqa,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (hlwdwn,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (runoff,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (fwsoil,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (gdrain,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (smelt,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (eveg,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (esno,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (esoi,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (eint,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (evpot,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg1,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg2,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg3,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg4,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg5,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg6,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg7,igrd,chfrstr,istrip,nchp,nn,

      call DIAG_VEGTILE_FILL (strdg8,igrd,chfrstr,istrip,nchp,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 STRIP2TILE(uz(1,1,l),igrd,tmpstrip,nchp,ijall,

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

      call STRIP2TILE(vz(1,1,l),igrd,tmpstrip,nchp,ijall,

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

      call STRIP2TILE(tz(1,1,l),igrd,tmpstrip,nchp,ijall,

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

      call STRIP2TILE(qz(1,1,l,1),igrd,tmpstrip,nchp,ijall,

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

      call DIAG_VEGTILE_FILL (tmpstrip,igrd,chfrstr,istrip,nchp,nn,

      _RL ecanopy(nchp,Nsx,Nsy)

      _RL swetshal(nchp,Nsx,Nsy)

      _RL swetroot(nchp,Nsx,Nsy)

      _RL swetdeep(nchp,Nsx,Nsy)

      _RL snodep(nchp,Nsx,Nsy)

      _RL capac(nchp,Nsx,Nsy)

      _RL chlt(nchp,Nsx,Nsy)

      _RL chlon(nchp,Nsx,Nsy)

      integer igrd(nchp,Nsx,Nsy)

      _RL tcanopy(nchp,Nsx,Nsy)

      _RL tdeep(nchp,Nsx,Nsy)

      dim(1)      = nchp

      ie(1)       = nchp

      dim(1)      = nchp

      ie(1)       = nchp

     .  nchp,nchptot,nchpland,

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

      _RL temptile(nchp)

       call GRD2MSC(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))

       do n = 1,nchp

     . phis_var(1,1,bi,bj),tgz(1,1,bi,bj),sea_ice,nchp,chlat(1,bi,bj),

      integer Nrphin,Nsxin,Nsyin,bi,bj,nchp

      _RL ctmt(nchp,Nsxin,Nsyin),xxmt(nchp,Nsxin,Nsyin)

      _RL yymt(nchp,Nsxin,Nsyin)

      _RL zetamt(nchp,Nsxin,Nsyin)

      _RL xlmt(nchp,Nrphin,Nsxin,Nsyin),khmt(nchp,Nrphin,Nsxin,Nsyin)

      _RL tke(nchp,Nrphin,Nsxin,Nsyin)

      _RL chfr(nchp,Nsxin,Nsyin),alai(nchp,Nsxin,Nsyin)

      _RL agrn(nchp,Nsxin,Nsyin)

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

      _RL chlat(nchp,Nsxin,Nsyin),chlon(nchp,Nsxin,Nsyin)

      _RL tcanopy(nchp,Nsxin,Nsyin),tdeep(nchp,Nsxin,Nsyin)

      _RL ecanopy(nchp,Nsxin,Nsyin),swetshal(nchp,Nsxin,Nsyin)

      _RL swetroot(nchp,Nsxin,Nsyin),swetdeep(nchp,Nsxin,Nsyin)

      _RL snodep(nchp,Nsxin,Nsyin),capac(nchp,Nsxin,Nsyin)

     . nchp,nchptot,nchpland,