if( pz(i,j)-ple(i,j,Lm+2-nthini).gt.25. ) then

       if( ple(i,j,Lm+2-nbasei).lt.(0.667*pz(i,j)) ) then

      if( (0.667*pz(i,j))-ple(i,j,Lm+2-nbase(i,j)) .gt. 

     .           ple(i,j,Lm+3-nbase(i,j))-(0.667*pz(i,j)) ) then

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

       dtdt(i,j,L) = dtdt(i,j,L) + dragt(i,j,L)*pz(i,j)/pkz(i,j,L)

      sumu(i,j) = sumu(i,j) + dragu(i,j,L)*dpres(i,j,L)/pz(i,j)

        tmpdiag(i,j) = dragx(i,j) + sumu(i,j)*pz(i,j)/grav*100

      sumu(i,j) = sumu(i,j) + dragv(i,j,L)*dpres(i,j,L)/pz(i,j)

        tmpdiag(i,j) = dragy(i,j) + sumu(i,j)*pz(i,j)/grav*100

      _RL pz(im,jm)

      subroutine GWDRAG (myid,pz,pl,ple,dpres,pkz,uz,vz,tz,qz,phis_var,

      _RL pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1)

     .                  pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,

     .             *  ( tz(i,j,L) + dtmoist(i,j,L)*tmstp/pz(i,j) )

      _RL pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1),dpres(im,jm,lm)

        pigather(indx) =     pz(i,j)

     .   pz,plz,plze,dpres,pkht,pkl,uz,vz,tz,qz,bi,bj,ntracerin,ptracer,

       tmpimjm(i,j) = 1./pz(i,j)

        tmpdiag(i,j)=(dtmoist(i,j,L)*sday*pkzgather(indgath,L)/pz(i,j))

        tmpdiag(i,j)=(dqmoist(i,j,L,1)*sday*1000./pz(i,j))

      _RL pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1),dpres(im,jm,lm)

     .        pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,co2,

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

      thtgz(i) = pz(i,1) * atimstp

     1 ,ndturb,nltop,ptop, pz, uz, vz, tz, qz, ntracers,ptracers

      call GRD2MSC(pz(1,1),im,jm,igrd,pmsc,nchp,nchptot)