integer nn,irun,nlay

      _RL cldlz(irun,nlay)

      _RL rhcrit(irun,nlay)

       do L = 1,nlay

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

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

     & .false., 'FRQTRB  ', 0, nlay-1, bi, bj, myid)

      do L = 1,nlay

      do l = 1,nlay

      integer im,jm,nlay,istrip,nymd,nhms,bi,bj,ndturb,nltop

      call CTEI ( tz,qz,fcctmp,qliqtmp,plz,pkz,pkht,im*jm,nlay )

      do l = 1,nlay

      CALL DIAGNOSTICS_FILL(fcctmp,'TRBFCC  ',0,nlay,3,bi,bj,myid)

     &                             'TRBQLIQ ',0,nlay,3,bi,bj,myid)

        tmpdiag(i,1) = tz(i,1,nlay) * pkht(i,1,nlay)

      _RL pz(im*jm,1),uz(im*jm,1,nlay),vz(im*jm,1,nlay)

     &                            'TURBV   ',0,nlay,3,bi,bj,myid)

     &                            'TURBU   ',0,nlay,3,bi,bj,myid)

      _RL tz(im*jm,1,nlay),qz(im*jm,1,nlay,ntracers)

     &                            'TURBQ   ',0,nlay,3,bi,bj,myid)

       do L = 1,nlay

      _RL plz(im*jm,1,nlay),plze(im*jm,1,nlay+1),dpres(im*jm,1,nlay)

      do l =1,nlay

      _RL pkht(im*jm,1,nlay+1),pkz(im*jm,1,nlay)

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

      _RL duturb(im*jm,1,nlay),dvturb(im*jm,1,nlay)

      _RL dtturb(im*jm,1,nlay),dqturb(im*jm,1,nlay,ntracers)

      _RL qliqavesw(im*jm,1,nlay),qliqavelw(im*jm,1,nlay)

      _RL fccavelw (im*jm,1,nlay),fccavesw (im*jm,1,nlay)

      _RL qqgrid   (im*jm,1,nlay)

      _RL qliqtmp(im*jm,1,nlay)

      _RL  fcctmp(im*jm,1,nlay)

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

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

      _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 xlsave(istrip,nlay),khsave(istrip,nlay)

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

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

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

       istnp1   =   istrip * (nlay+1)

       istnlay  =   istrip * nlay

     1                        nlay,irun)

      integer nlev,nlay,irun

      _RL EE(irun*(nlay-1),1),RF(irun*(nlay-1),1)

     1                      nlay,irun)

      integer nlev,nlay,irun

      INTEGER INTSTB(irun*nlay,1), INTQ(irun*nlay,1)

      _RL F2(irun*(nlay-1),1),F3(irun*(nlay-1),1)

      _RL F4(irun*(nlay-1),1),XQ(irun*(nlay-1),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  (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  (xlmt,xlsave,nchptot,nchp,istrip,nlay,nn)

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

       u(i,nlay+1) = 0.

       v(i,nlay+1) = 0.

      th(i,nlay+1) = th(i,nlay+1) / pke(i,nlay+1)

      sh(i,nlay+1) = qa(i)

      tmpstrip(i) = sh(i,nlay+1)*1000

      do L = 1,nlay+1

      sh(i,nlay+1) = qa(i)

      do L =1,nlay

        stdiag(i) = ( thv(i,nlay+1)-thv(i,nlay) ) / pke(i,nlay+1)

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

      call PASTIT (qq,tke,istrip,nchp,nchptot,nlay,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)

       evpot(i) = transsh(i,nlay) * (shgstr(i) - sh(i,nlay))

       tmpnlay(i) = th(i,nlay)*pk(i,nlay)

       psurf(i) = pe(i,nlay+1)

       wspeed(i) = sqrt(u(i,nlay)*u(i,nlay) + v(i,nlay)*v(i,nlay))

       eturb(i) = eturb(i) * pke(i,nlay+1)

       dedqa(i) = dedqa(i) * pke(i,nlay+1)

       hsturb(i) = hsturb(i) * pke(i,nlay+1)

     I   tmpnlay, sh(1,nlay), cd, cosz, pardr, pardf,

     &          pke(nland+1,nlay+1),    icest(nland+1),

      th(i,nlay+1) = tc(i) / pke(i,nlay+1)

      do L = 1,nlay

      th(i,L) = th(i,L) + dthdthg(i,L)*(tc(i)-savetc(i))/pke(i,nlay+1)

      sh(i,nlay+1) = qa(i)

      do L = 1,nlay

      do L = 1,nlay

      do  l=1,nlay

      tmpstrip(i) = stqflux(i,nlay) * 86400

      tmpstrip(i) = stqflux(i,nlay) * alhl

      tmpstrip(i) = sttflux(i,nlay) * cp

     & .false., 'TUFLUX  ', 0, nlay, bi, bj, myid)

     & .false., 'TVFLUX  ', 0, nlay, bi, bj, myid)

      do l=1,nlay

     & .false., 'TTFLUX  ', 0, nlay, bi, bj, myid)

      do l=1,nlay

     & .false., 'TQFLUX  ', 0, nlay, bi, bj, myid)

     & .false., 'RI      ', 0, nlay, bi, bj, myid)

     & .false., 'KH      ', 0, nlay, bi, bj, myid)

     & .false., 'KM      ', 0, nlay, bi, bj, myid)