_RL pl(sNx,sNy,Nrphys,nSx,nSy),pkl(sNx,sNy,Nrphys,nSx,nSy)

      COMMON /saver/ p,pl,ple,dpres,pkle,pkl

       pl(i,j,L,bi,bj) = (pephy(i,j,L,bi,bj)+pephy(i,j,L+1,bi,bj))/200.

     & p(1,1,bi,bj),u,v,t,q,pl(1,1,1,bi,bj),ple(1,1,1,bi,bj),

     &   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,

     &   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,

     &  ptop,p,u,v,t,q,ntracer,ptracer,pl,ple,dpres,pkht,pkl,

      call GWDRAG (myid,p,pl,ple,dpres,pkl,u,v,t,q,phis_var,

      call RAYLEIGH(myid,pl,pkl,pkht,surfz,u,v,t,q,im,jm,lm,

     &    p,pl,ple,dpres,pkht,pkl,u,v,t,q,bi,bj,ntracer,ptracer,qq,

      _RL pl(im,jm,lm),ple(im,jm,lm+1),dpres(im,jm,lm)

     & p,u,v,t,q,pl,ple,dpres,pkht,pkl,surfz,fracland,landtype,radswt,

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

     .                  std,pl,ple,dpres,

      _RL ple(irun,Lm+1), pl(irun,Lm), dpres(irun,Lm)

       ro(i,L) = pl(i,Lm+1-L)/(rgas*t(i,Lm+1-L))

       dz(i,L) = agrav*roiave*( pl(i,Lm+2-L)-pl(i,Lm+1-L) )

       nbar(i) = nbar(i) + bv(i,L)*(pl(i,Lm+2-L)-pl(i,Lm+1-L))

       nbar(i) = nbar(i)/(pl(i,Lm)-pl(i,Lm+1-nbase(i)))

      _RL pl(im,jm,Lm)

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

       CALL STRIP (PLZ ,pl     ,im*jm,ISTRIP,lm   ,NN)

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

      subroutine IRRAD (m,n,ndim,np,pl,ta,wa,oa,ts,co2,

      _RL pl(istrip,lm)

      _RL pl(m,ndim,np+1),ta(m,ndim,np),wa(m,ndim,np),oa(m,ndim,np),

          if (pl(i,j,1) .eq. 0.0) then

              pl(i,j,1)=1.0e-4

         pa(i,j,k)=0.5*(pl(i,j,k)+pl(i,j,k+1))

         dp           = pl(i,j,k+1)-pl(i,j,k)

             dp          =pl(i,j,k+1)-pl(i,j,k)

       plev = pl(indx,L)

         if( pl(i).ge.p(L)   .and.

     .       pl(i).le.p(L+1)) then

         temp = ( pl(i)-p(L) )/( p(L+1)-p(L) )

      if( pl(i).lt.p(1)   ) acr(i) = a(1)

      if( pl(i).gt.p(lma) ) acr(i) = a(lma)

      acr(i) = acr(i) * (plb(i)-pl(i))

       call QSAT ( temp1(i,3),pl(i,L),temp1(i,2),temp1(i,6),.true. )

           call QSAT( tl(i,L),pl(i,L),rh,dum,.false.)

           if( pl(i,L).ge.10.0 ) CLSBTH(I,L) = fact

      subroutine SRCLOUDS (th,q,plk,pl,plke,cloud,cldwat,irun,irise,

      _RL   pl(irun,irise)

      call QSAT ( temp,pl(i,L),qs(i,L),dum,.false. )

      call CTEI ( th,q,cloud,cldwat,pl,plk,plke,irun,irise )

      subroutine CTEI ( th,q,cldfrc,cldwat,pl,plk,plke,im,lm )

      _RL plk(im,lm),pl(im,lm),cldfrc(im,lm)

      p = pl(i,L)

        pcheck = pl(i,L)

         ppbl = pl(i,nsubcl)

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

      subroutine OPTHK ( tl,pl,ple,lz,cf,cfm,lwi,im,jm,lm,tau )

      _RL  pl(im,jm,lm)

      subroutine SORAD(m,n,ndim,np,pl,ta,wa,oa,co2,

      _RL pl(m,ndim,np+1),ta(m,ndim,np),wa(m,ndim,np),oa(m,ndim,np)

          dp(i,j,k)=pl(i,j,k+1)-pl(i,j,k)

          scal(i,j,k)=dp(i,j,k)*(.5*(pl(i,j,k)+pl(i,j,k+1))/300.)**.8

       call QSAT ( tl(i,nlev),pl(i,nlev),shsat(i,nlev),dum,.false. )

       call QSAT ( tl(i,LL),pl(i,LL),shsat(i,LL),dum,.false. )

         if(pl(i,L).le.150.) then

         if(pl(i,L).le.150.) then

      _RL pl(sNx,sNy,Nrphys,nSx,nSy),pkl(sNx,sNy,Nrphys,nSx,nSy)