call STRIPITINT (landtype,lwi,im*jm,im*jm,istrip,1,nn)

      integer iw,ip,nn

        nn=mod(i,nt)+1

        pa = coef1(ip,iw-1,nn)*(1.-fp)+coef1(ip+1,iw-1,nn)*fp

        pb = coef1(ip,iw,  nn)*(1.-fp)+coef1(ip+1,iw,  nn)*fp

        pc = coef1(ip,iw+1,nn)*(1.-fp)+coef1(ip+1,iw+1,nn)*fp

        ba = coef2(ip,iw,  nn)*(1.-fp)+coef2(ip+1,iw,  nn)*fp

        bb = coef2(ip,iw+1,nn)*(1.-fp)+coef2(ip+1,iw+1,nn)*fp

        ca = coef3(ip,iw,  nn)*(1.-fp)+coef3(ip+1,iw,  nn)*fp

        cb = coef3(ip,iw+1,nn)*(1.-fp)+coef3(ip+1,iw+1,nn)*fp

      integer i,j,l,n,nn

      integer nn,lng,lenc,k,nltop,nlayr

       CALL CLOUD(nn,lng, LENC, K, NLTOP, nlayr, IB, RASBLF,SETRAS,FRAC

      integer i,j,iloop,indx,indgath,l,nn,num,numdeps,nt

      SUBROUTINE CLOUD(nn,lng, LENC, K, NLTOP, nlayr, IC, RASALF

      integer nn,lng,LENC,K,NLTOP,nlayr,ic,ntracedim, ntracer

          nn = 3*3600.0/tmstp + 1

      cldmin = rasmax*(1.0-tmstp/3600.)**nn   

      integer nn,irun,nlay

        call STRIP ( ugather(1,1,nt), ul(1,1,nt),im*jm,istrip,lm,nn )

      call PASTE (  srcld,   cldsr,istrip,im*jm,lm,nn )

      call PASTE (  cldlz,cldwater,istrip,im*jm,lm,nn )

      call PASTE ( clsbth,   cldls,istrip,im*jm,lm,nn )

      call PASTE ( clboth,   cpen ,istrip,im*jm,lm,nn )

      integer i,j,n,nn,iyear,iyearbc,nd,ndby3,imbc,jmbc

         write(6,1000) (nn,nymdbc(nn),nhmsbc(nn),nn=n,n+2)

        write(6,1000) (nn,nymdbc(nn),nhmsbc(nn),nn=ndby3*3+1,ndatebc)

      integer i,j,n,nn,iyear,iyearbc,nd,ndby3,imbc,jmbc

         write(6,1000) (nn,nymdbc(nn),nhmsbc(nn),nn=n,n+2)

        write(6,1000) (nn,nymdbc(nn),nhmsbc(nn),nn=ndby3*3+1,ndatebc)

      do 1000 nn = 1,npcs

      call STRIP ( cldtot,clro,im*jm,istrip,lm,nn )

      integer   i,j,L,nn,nsecf

      call PASTE ( dtsw , dtradsw ,istrip,im*jm,lm,nn )

      call PASTE ( dtswc, dtswclr ,istrip,im*jm,lm,nn )

      call PASTE ( flux   (1,1),osr   ,istrip,im*jm,1,nn )

      call PASTE ( fluxclr(1,1),osrclr,istrip,im*jm,1,nn )

      call PASTE ( flux   (1,lm+1),radswg,istrip,im*jm,1,nn )

      call PASTE ( fluxclr(1,lm+1),swgclr,istrip,im*jm,1,nn )

      call PASTE ( dirpar,fdirpar,istrip,im*jm,1,nn )

      call PASTE ( difpar,fdifpar,istrip,im*jm,1,nn )

      call PASTE ( tstrip,albedo,istrip,im*jm,1,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 (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,

      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,

     . nn,.false., 'UFLUX   ', 1, 1, bi, bj, myid)

     . nn,.false., 'VFLUX   ', 1, 1, bi, bj, myid)

      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,

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

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

      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,

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

      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,

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

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

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

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

      integer nn,irun,nlev,nltop,ntrace,ntracedim,itrtrb,nhms,nymd

      integer istnp1,istnlay,itrtrb,i,j,L,nn,nt

      integer nn,n,irun

      integer nn,irun,itype

      do 2000 nn = 1, numstrips

       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 TRBFLX(nn,th,thv,sh,u,v,qq,p,pe,pk,pke,dpstr,stwatr,stz0,

      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)