subroutine scopy(n,sx,incx,sy,incy)
c
      IMPLICIT NONE
c
      integer n,incx,incy,ix,iy,i
      !real sx((n-1)*incx+1),sy((n-1)*incy+1)     
      real sx(*),sy(*)   ! TB18
c
      iy=1
      ix=1
      do 10 i=1,n
         sy(iy)=sx(ix)
         ix=ix+incx
         iy=iy+incy
         !print*, 'TB18 scop',i,sx(ix),sy(iy)
10    continue
c
      return
      end