subroutine STRIPIT(a,b,irun,ia,ib,l,k)

      _RL a(ia,l), b(ib,l)

        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)

      subroutine STRIPITINT(a,b,irun,ia,ib,l,k)

      integer ia,ib,irun,l,k,a(ia,l),b(ib,l)

        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)

      subroutine PASTIT(b,a,ib,ia,irun,L,k)

      _RL a(ia,l), b(ib,l)

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

      subroutine PSTBITINT(b,a,ib,ia,irun,l,k)

      integer b(ib,l)

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

      subroutine PSTBMPIT(b,a,ib,ia,irun,l,k)

      _RL a(ia,l), b(ib,l)

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

      subroutine STRIP2TILE(a,indx,b,irun,ia,ib,levs,npeice)

      _RL a(ia,levs), b(ib,levs)

        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)

      subroutine PASTE2GRD_OLD(b,indx,chfr,ib,numpts,a,ia,levs,npeice)

      _RL a(ia,levs), b(ib,levs), chfr(ib)

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

      subroutine PASTE2GRD (b,indx,chfr,ib,numpts,a,ia,levs,npeice,

      _RL a(ia,levs), b(ib,levs), chfr(ib)

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

         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)

      _RL     b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)

          b(i,j,k) = 1. _d 0 - c(i,j,k) - a(i,j,k)

     I                        a, b, c,

            eps4 = b*pLoc + c0

      _RL a0, b, c0

      b  =  -7.61    _d -4

        PtrISt(g)som_P(a,b,c,d,e,f)

     &           + b*pLoc(I,J) + c0

            eps4 = b*pLoc(I,J) + c0

      _RL a0, a1, a2, b, c0

      b  =  -7.61    _d -4

       b  = -7.53     _d -4

          b(i,j,k) = 1. _d 0 - c(i,j,k) - a(i,j,k)

         IF (b(i,j,1).NE.0.) bet(i,j,1) = 1. _d 0 / b(i,j,1)

          IF ( ( b(i,j,k) - a(i,j,k)*gam(i,j,k) ) .NE. 0.) 

     &        bet(i,j,k) = 1. _d 0 / ( b(i,j,k) - a(i,j,k)*gam(i,j,k) )

      _RL b(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)

          b(i,j,k) = 1. _d 0 - c(i,j,k) - a(i,j,k)

         IF (b(i,j,1).NE.0.) bet(i,j,1) = 1. _d 0 / b(i,j,1)

          IF ( ( b(i,j,k) - a(i,j,k)*gam(i,j,k) ) .NE. 0.) 

     &        bet(i,j,k) = 1. _d 0 / ( b(i,j,k) - a(i,j,k)*gam(i,j,k) )

      _RL b(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)

      _RS function smoothMin_R4( a, b )

      _RS a, b

      smoothMin_R4 = .5*( a+b - smoothAbs_R4(a-b) )

      _RL function smoothMin_R8( a, b )

      _RL a, b

      smoothMin_R8 = .5*( a+b - smoothAbs_R8(a-b) )

  and b) to comment out these lines.

  and b) to comment out these lines.

  and b) to comment out these lines.

      double precision z1, b, discri

      b  = z1 + dble(fp)

      if (sngl(b)*sign .gt. 0.0) then

         t    = t + fp*(ta-t)/sngl(b+discri)

         den  = sngl(z1+b+dble(fpa))

         anum = sngl(b-discri)

            b = x2 + k1local*x + k12

            b2=b*b

            fn = k1local*x*diclocal/b +

     &        2.0*diclocal*k12/b +

            df = ((k1local*diclocal*b) - k1local*x*diclocal*db)/b2 -

              b = x2 + k1local*x + k12

              b2=b*b

              fn = k1local*x*diclocal/b +

     &          2.0*diclocal*k12/b +

        _RL  b