do L=1,lm

      integer myid,im,jm,lm,bi,bj

      _RL zsurf(im,jm),pres(im,jm,lm),pkap(im,jm,lm)

      _RL pekap(im,jm,lm+1)

      _RL u(im,jm,lm),v(im,jm,lm),t(im,jm,lm),s(im,jm,lm)

      _RL rfu(im,jm,lm),rfv(im,jm,lm),rft(im,jm,lm)

      _RL rf(im,jm,lm)

      _RL z(im,jm,lm)

      _RL dz(im,jm,lm)

      subroutine RAYLEIGH(myid,pres,pkap,pekap,zsurf,u,v,t,s,im,jm,lm,

      do L=1,lm

       z(i,j,lm) = zsurf(i,j) +  0.5 * dz(i,j,lm)

      do L=lm-1,1,-1

      do L=1,lm

        do L =1,lm

       CALL STRIP (OZ,OZL      ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (PLZE,ple    ,im*jm,ISTRIP,lm+1 ,NN)

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

       CALL STRIP (PKZ ,pk     ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (PKHT,pke    ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (TZ,TZL      ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (qz,SHZL     ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (cldtot,CLRO ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (cldmxo,CLMO ,im*jm,ISTRIP,lm   ,NN)

       CALL STRIP (lwlz,cldlz  ,im*jm,ISTRIP,lm   ,NN)

       DO L = 1,lm

       DO L = 1,lm

      DO L = 1,lm

      TZL(I,lm) = ( 2*tzl(i,lm)+tsurf(i) )/3.0  

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

      do L = 1,lm

      do L = 1,lm

      do L = 1,lm

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

      integer im,jm,lm        

       do L = 1,lm

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

        flx   (i,lm+1) = -flx   (i,lm+1)

      _RL dpres(im,jm,lm),pkht(im,jm,lm+1),pkz(im,jm,lm)

        flxclr(i,lm+1) = -flxclr(i,lm+1)

             dsgmt4(i) = - dfdts(i,lm+1)

       CALL PASTE ( flx   (1,lm+1),RADLWG,ISTRIP,im*jm,1,NN )

       CALL PASTE ( flxclr(1,lm+1),LWGCLR,ISTRIP,im*jm,1,NN )

      _RL tz(im,jm,lm),qz(im,jm,lm),oz(im,jm,lm)

      _RL co2,cfc11,cfc12,cfc22,methane(lm),n2o(lm)    

      DO L = 1,lm

      DTRAD (I,L) = ple(i,lm+1) * DTRAD (I,L)/pk(I,L)

      DTRADC(I,L) = ple(i,lm+1) * DTRADC(I,L)/pk(I,L)

       dtdtg(I,L) = ple(i,lm+1) * dtdtg (I,L)/pk(I,L)

      CALL PASTE ( tmpstrip ,tmpimjm ,ISTRIP,im*jm,lm,NN )

      CALL PASTE ( DTRAD ,DTRADLW ,ISTRIP,im*jm,lm,NN )

      CALL PASTE ( DTRADC,DTRADLWC,ISTRIP,im*jm,lm,NN )

      CALL PASTE ( dtdtg ,dlwdtg  ,ISTRIP,im*jm,lm,NN )

      do L = 1,lm

      _RL dtradlw(im,jm,lm),dlwdtg (im,jm,lm) 

      _RL dtradlwc(im,jm,lm),lwgclr(im,jm)     

      _RL cldlw(im,jm,lm),clwmo(im,jm,lm),lwlz(im,jm,lm)

      do L = 1,lm

      _RL qliqave(im,jm,lm),fccave(im,jm,lm)

      _RL cldtot (im,jm,lm)

      _RL cldmxo (im,jm,lm)

      _RL pl(istrip,lm)

      _RL pk(istrip,lm)

      _RL pke(istrip,lm)

      _RL ple(istrip,lm+1)

      _RL ADELPL(ISTRIP,lm)

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

      _RL OZL(ISTRIP,lm),TZL(ISTRIP,lm)

      _RL SHZL(ISTRIP,lm),CLRO(ISTRIP,lm)

      _RL CLMO(ISTRIP,lm)

      _RL flx(ISTRIP,lm+1),flxclr(ISTRIP,lm+1)

      _RL cldlz(istrip,lm)

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

      _RL taual(istrip,lm,10)

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

      _RL tmpstrip(istrip,lm)

      _RL tmpimjm(im,jm,lm)

     .                  im,jm,lm,

        do L =1,lm