_RL   offset, alpha, rasmax

      offset = 0.10

     2   tmfrc,cp,gravity,alhl,gamfac,cldlz,RHCRIT,offset,alpha)

      _RL tmscl,tmfrc,cp,gravity,alhl,gamfac,offset,alpha

     1                                                rhc,offset,alpha)

      _RL offset,alpha

      rhcrit = rhc(i,L) - offset

      ratio  = alpha*(rh-rhcrit)/offset

     .  CLSBTH,TMSTP,one,cp,grav,alhl,gamfac,cldlz,rhcrit,offset,alpha)

     .                                              rhcrit,offset,alpha)

      integer   i,j,Lena,offset

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      offset = offset+1

        b(i,j) = a(i+offset-1,j)

        b(i,j) = a(i+offset-1,j)

        b(Lena+i,j) = a(Lena+offset-1,j)

      integer i,j,Lena,offset

      offset = ib*(k-1)

      Lena = min(ib,irun-offset)

      offset = offset+1

        b(i,j) = a(i+offset-1,j)

        b(i,j) = a(i+offset-1,j)

        b(Lena+i,j) = a(Lena+offset-1,j)

      integer ib,ia,L,k,irun,Lena,offset

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      offset = offset+1

      a(i+offset-1,j) = b(i,j)

      integer ib,ia,L,k,irun,Lena,offset

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      offset = offset+1

      a(i+offset-1,j) = a(i+offset-1,j) + float(b(i,j))

      integer ib,ia,L,k,irun,Lena,offset

      offset = ib*(k-1)

      Lena    = min(ib,irun-offset)

      offset = offset+1

      a(i+offset-1,j) = a(i+offset-1,j) + b(i,j)

      integer i,k,Lena,offset

      offset = ib*(npeice-1)

      Lena    = min(ib,irun-offset)

      offset = offset+1

        b(i,k) = a(indx(i+offset-1),k)

        b(i,k) = a(indx(i+offset-1),k)

        b(Lena+i,k) = a(indx(Lena+offset-1),k)

      integer i,L,offset,Lena

      offset = ib*(npeice-1)

      Lena    = min(ib,numpts-offset)

      offset = offset+1

       a(indx(i+offset-1),L) = a(indx(i+offset-1),L) + b(i,L)*chfr(i)

      integer i,L,offset,Lena

      offset = ib*(npeice-1)

      Lena    = min(ib,numpts-offset)

      offset = offset+1

        if( a(indx(i+offset-1),L).eq.undef .or. b(i,L).eq.undef ) then

         a(indx(i+offset-1),L) = undef

         a(indx(i+offset-1),L)=a(indx(i+offset-1),L) + b(i,L)*chfr(i)

        a(indx(i+offset-1),L)=a(indx(i+offset-1),L) + b(i,L)*chfr(i)

          offset = ib*(npeice-1)

          Lena   = MIN(ib,numpts-offset)

             ij = indx(ivt+offset) - 1

              ij = indx(ivt+offset) - 1

           offset = ib*(npeice-1)

           Lena    = MIN(ib,numpts-offset)

             ij = indx(ivt+offset)

              ij = indx(ivt+offset)

              ij = indx(ivt+offset)

      INTEGER offset, Lena

      offset     = 0. _d 0

                offset  = offset + RA(i,j,bi,bj)*bpdifanom(i,j,bi,bj)

      _GLOBAL_SUM_RL( offset     , mythid )

     &                                 - offset/offset_sum

      _RL offset

      offset     = 0. _d 0

                  offset=offset+RA(i,j,bi,bj)*mean_psMssh_all(i,j,bi,bj)

                  offset=offset+RA(i,j,bi,bj)

      _GLOBAL_SUM_RL( offset     , mythid )

     &      offset,offset_sum

      offset = offset / offset_sum

     &                 mean_psMssh_all(i,j,bi,bj) - offset

     &             psmean(i,j,bi,bj) -mdtob(i,j,bi,bj) - offset

     &            psmean(i,j,bi,bj)=psmean(i,j,bi,bj)-offset

      _RL offset

      offset     = 0. _d 0

              offset     = offset +

      _GLOBAL_SUM_RL( offset     , mythid )

      offset = offset / offset_sum

              prof_etan_mean(i,j,bi,bj)=offset+psmean(i,j,bi,bj)

              junk = psmean(i,j,bi,bj) - mdt(i,j,bi,bj) + offset

     &                 (tpobs(i,j,bi,bj)+mdt(i,j,bi,bj)-offset) )

     &                 (ersobs(i,j,bi,bj)+mdt(i,j,bi,bj)-offset) )

     &                 (gfoobs(i,j,bi,bj)+mdt(i,j,bi,bj)-offset) )

      _RL offset