SUBROUTINE add_vec(n,x,y,z)
      INTEGER n,i
      REAL x(n),y(n),z(n)
      DO i=1,n
         z(i)=x(i)+y(i)
      ENDDO
      RETURN
      END