INTEGER istrip,npcs

      istrip = im*jm

      call SWRIO ( nymd,nhms,bi,bj,ndsw,myid,istrip,npcs,

      call LWRIO ( nymd,nhms,bi,bj,myid,istrip,npcs,

      call TURBIO (im,jm,lm,istrip,nymd,nhms,bi,bj,

     &   turbu,turbv,turbt,im,jm,lm,bi,bj,istrip,npcs,imglobal)

      call MOISTIO (ndmoist,istrip,npcs,

      call STRIPIT ( phis_std,std,im*jm,im*jm,istrip,1,n )

      call STRIPIT ( pz,ps,im*jm,im*jm,istrip,1 ,n )

      call STRIPIT ( uz,us,im*jm,im*jm,istrip,Lm,n )

      call STRIPIT ( vz,vs,im*jm,im*jm,istrip,Lm,n )

      call STRIPIT ( tv,ts,im*jm,im*jm,istrip,Lm,n )

      call STRIPIT ( pl,plstr,im*jm,im*jm,istrip,Lm,n )

      call STRIPIT ( ple,plestr,im*jm,im*jm,istrip,Lm+1,n )

      call STRIPIT ( dpres,dpresstr,im*jm,im*jm,istrip,Lm,n )

      call STRIPITINT ( nthin,nthinstr,im*jm,im*jm,istrip,1,n )

      call STRIPITINT ( nbase,nbasestr,im*jm,im*jm,istrip,1,n )

     .            istrip,Lm,nthinstr,nbasestr,lstar )

      call PASTIT( dragus,dragu,istrip,im*jm,im*jm,Lm,n )

      call PASTIT( dragvs,dragv,istrip,im*jm,im*jm,Lm,n )

      call PASTIT( dragxs,dragx,istrip,im*jm,im*jm,1 ,n )

      call PASTIT( dragys,dragy,istrip,im*jm,im*jm,1 ,n )

      integer myid,im,jm,Lm,bi,bj,istrip,npcs,imglobal

     .         dudt,dvdt,dtdt,im,jm,Lm,bi,bj,istrip,npcs,imglobal)

      _RL std(istrip), ps(istrip)

      _RL us(istrip,Lm), vs(istrip,Lm), ts(istrip,Lm)

      _RL dragus(istrip,Lm), dragvs(istrip,Lm) 

      _RL dragxs(istrip), dragys(istrip)

      _RL plstr(istrip,Lm),plestr(istrip,Lm+1),dpresstr(istrip,Lm)

      integer nthinstr(istrip),nbasestr(istrip)

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

      call OPTHK ( tzl,pl,ple,cldlz,clro,clmo,lwi,istrip,1,lm,tauc )

      do i = 1,istrip

      do i = 1,istrip

      integer nymd,nhms,istrip,npcs,bi,bj,myid

      do i = 1,istrip

      call IRRAD ( istrip,1,1,lm,ple,tzl,shzl,ozl,tsurf,co2,

       do i = 1,istrip

      _RL pl(istrip,lm)

      _RL pk(istrip,lm)

      _RL pke(istrip,lm)

      _RL ple(istrip,lm+1)

      _RL dtrad(istrip,lm),dtradc(istrip,lm)

      _RL cldlz(istrip,lm)

      _RL dfdts(istrip,lm+1),dtdtg(istrip,lm)

      _RL emiss(istrip,10)

      _RL taual(istrip,lm,10)

      subroutine LWRIO (nymd,nhms,bi,bj,myid,istrip,npcs,

      _RL ssaal(istrip,lm,10)

      _RL asyal(istrip,lm,10)

      _RL cwc(istrip,lm,3)

      _RL reff(istrip,lm,3)

      _RL tauc(istrip,lm,3)

      integer lwi(istrip)

      _RL tmpstrip(istrip,lm)

      _RL saveu  (istrip,lm,ntracerin+2-ptracer)

      _RL usubcl (istrip,   ntracerin+2-ptracer)

      _RL     ple(istrip,lm+1)

      _RL      dp(istrip,lm)

      _RL     cvu(istrip,lm,ntracerin+2-ptracer)

      integer pbl(istrip),depths(lm)

      _RL   cldlz(istrip,lm), cldwater(im,jm,lm)

      _RL   rhfrac(istrip), rhmin, pup, ppbl, rhcrit(istrip,lm)

      _RL    clfrac (istrip,lm)

      _RL    cldmas (istrip,lm)

      _RL    detrain(istrip,lm)

      _RL    psubcld    (istrip), psubcldg (im,jm)

      _RL    psubcld_cnt(istrip), psubcldgc(im,jm)

      integer ndmoist,istrip,npcs

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

       do i = 1,istrip

      CALL RAS ( NN,istrip,nindeces(nsubcl),NLRAS,NLTOP,lm,TMSTP

      _RL   srcld(istrip,lm)

      subroutine MOISTIO (ndmoist,istrip,npcs,

      call PASTE( CVTH,deltgather,istrip,im*jm,lm,NN )

      call PASTE(  CVQ,delqgather,istrip,im*jm,lm,NN )

      call PASTE( CVU(1,1,nt),delugather(1,1,nt),istrip,im*jm,lm,NN )

      call PASTE( cldmas,tmpgather,istrip,im*jm,lm,NN)

      call PASTE(detrain,pkegather,istrip,im*jm,lm,NN)

      call PASTE(psubcld    ,psubcldg ,istrip,im*jm,1,NN)

      call PASTE(psubcld_cnt,psubcldgc,istrip,im*jm,1,NN)

       call PASTE(tmp1(1,L),deltrnev(1,L),istrip,im*jm,1,NN)

       call PASTE(tmp1(1,L),delqrnev(1,L),istrip,im*jm,1,NN)

      do i=1,istrip

      call SRCLOUDS (saveth,saveq,plk,pl,plke,clsbth,cldlz,istrip,lm,

      do i=1,istrip

      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 )

      do i = 1,istrip

      do i = 1,istrip

        do i = 1,istrip

      _RL saveth (istrip,lm)

      _RL saveq  (istrip,lm)

      integer nymd,nhms,bi,bj,ndswr,myid,istrip,npcs

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

      do i= 1,istrip

      call SORAD ( istrip,1,1,lm,ple,tzl,qzl,ozl,co2,

      do i=1,istrip

      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 )

      do i=1,istrip

      call PASTE ( tstrip,albedo,istrip,im*jm,1,nn )

      subroutine SWRIO (nymd,nhms,bi,bj,ndswr,myid,istrip,npcs,

      _RL fdirir(istrip),fdifir(istrip)

      _RL fdiruv(istrip),fdifuv(istrip)

      _RL flux(istrip,lm+1)

      _RL fluxclr(istrip,lm+1)

      _RL dtsw(istrip,lm)

      _RL dtswc(istrip,lm)

      _RL taul   (istrip,lm)

      _RL reff   (istrip,lm,2)

      _RL tauc   (istrip,lm,2)

      _RL taua   (istrip,lm)

      _RL tstrip (istrip)

      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,

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

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

      do i=1,istrip

      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,

      do i=1,istrip

      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 im,jm,nlay,istrip,nymd,nhms,bi,bj,ndturb,nltop

      _RL lats(istrip),lons(istrip),cosz(istrip),icest(istrip)

      _RL rainls(istrip),raincon(istrip),newsnow(istrip)

      _RL pardf(istrip),pardr(istrip),swnet(istrip)

      _RL hlwdwn(istrip),alwrad(istrip),blwrad(istrip)

      _RL tmpnlay(istrip)

      _RL laistrip(istrip),grnstrip(istrip),z2str(istrip),cd(istrip)

      _RL scatstr(istrip), rs1str(istrip), rs2str(istrip)

      _RL rdcstr(istrip),u2fstr(istrip),dqsdtstr(istrip)

      _RL eturb(istrip),dedqa(istrip),dedtc(istrip)

      _RL hsturb(istrip),dhsdqa(istrip),dhsdtc(istrip)

      _RL savetc(istrip),saveqa(istrip),lwstrip(istrip)

      _RL chfrstr(istrip),psurf(istrip),shgstr(istrip)

      integer types(istrip),igrdstr(istrip)

      _RL evap(istrip),shflux(istrip),runoff(istrip),bomb(istrip)

      _RL eint(istrip),esoi(istrip),eveg(istrip),esno(istrip)

      _RL smelt(istrip),hlatn(istrip),hlwup(istrip),gdrain(istrip)

      _RL runsrf(istrip),fwsoil(istrip),evpot(istrip)

      _RL strdg1(istrip),strdg2(istrip),strdg3(istrip),strdg4(istrip)

      _RL strdg5(istrip),strdg6(istrip),strdg7(istrip),strdg8(istrip)

      _RL strdg9(istrip),tmpstrip(istrip),qicestr(istrip)

      _RL dqicestr(istrip)

      _RL u(istrip,nlay+1), v(istrip,nlay+1), th(istrip,nlay+1)

      _RL sh(istrip,nlay+1), thv(istrip,nlay+1), pe(istrip,nlay+1)

      _RL tracers(istrip,nlay+1,ntracers)

      _RL dpstr(istrip,nlay),pke(istrip,nlay+1)

      _RL pk(istrip,nlay), qq(istrip,nlay),   p(istrip,nlay)

      _RL sri(istrip,nlay), skh(istrip,nlay), skm(istrip,nlay)

      _RL stuflux(istrip,nlay), stvflux(istrip,nlay)

      _RL sttflux(istrip,nlay), stqflux(istrip,nlay)

      _RL frqtrb(istrip,nlay-1)

      _RL dshdthg(istrip,nlay),dthdthg(istrip,nlay)

      _RL dshdshg(istrip,nlay),dthdshg(istrip,nlay)

      _RL transth(istrip,nlay), transsh(istrip,nlay)

      _RL tc(istrip),td(istrip),qa(istrip)

      _RL swet1(istrip),swet2(istrip),swet3(istrip)

      _RL capacity(istrip),snowdepth(istrip)

      _RL stz0(istrip)

      _RL stdiag(istrip)

      _RL tends(istrip),sustar(istrip), sz0(istrip),pbldpth(istrip)

      _RL sct(istrip), scu(istrip), swinds(istrip)

      _RL stu2m(istrip),stv2m(istrip),stt2m(istrip),stq2m(istrip)

      _RL stu10m(istrip),stv10m(istrip),stt10m(istrip),stq10m(istrip)

      integer  stwatr(istrip)

      _RL  wspeed(istrip)

      _RL ctsave(istrip),xxsave(istrip),yysave(istrip)

      _RL zetasave(istrip)

      _RL xlsave(istrip,nlay),khsave(istrip,nlay)

      _RL qliq(istrip,nlay),turbfcc(istrip,nlay)

      SUBROUTINE TURBIO (im,jm,nlay,istrip,nymd,nhms,bi,bj, qbeg

       istnp1   =   istrip * (nlay+1)

       istnlay  =   istrip * nlay

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

       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 ASTRO ( 20040321,nhms,lats,lons,istrip,cosz,ra )

       call ASTRO ( nymd,nhms,lats,lons,istrip,cosz,ra )

      do i = 1,istrip

      do i =1,istrip

      do i =1,istrip

      do i=1,istrip

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

      do i =1,istrip

      do i =1,istrip

      do i =1,istrip

      do i = 1,istrip

       do i=1,istrip

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

     4 stq10m,istrip,nlay,nltop,nymd,nhms,grav,cp,rgas,faceps,virtcon,

      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)

      do i = 1,istrip

      do i = 1,istrip

      do i = 1,istrip

      do i = 1,istrip

      do i = 1,istrip

      do i = 1,istrip

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

      do i = 1,istrip

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

      do i = 1,istrip

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

      do i =1,istrip

      do i =1,istrip

      do i =1,istrip

      do i =1,istrip

      do i =1,istrip

     1                                                 istrip,1,nn)

      do i =1,istrip

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

     1                                                 istrip,1,nn)

      do i =1,istrip

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

     1                                                 istrip,1,nn)

      do i =1,istrip

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

     1                                                 istrip,1,nn)

      do i =1,istrip

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

      do i=1,istrip

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

      do i=1,istrip

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

      do i=1,istrip

      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,

      do i=1,istrip

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

      do i=1,istrip

      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,