acc.f90 Source File

!
! $Header$
!
  subroutine acc(vec,d,im)
    implicit none
    integer :: im
    real :: vec(im,im),d(im)
    integer :: i,j
    real ::sum
    real,external :: ssum
    do j=1,im
      do i=1,im
        d(i)=vec(i,j)*vec(i,j)
      enddo
      sum=ssum(im,d,1)
      sum=sqrt(sum)
      do i=1,im
        vec(i,j)=vec(i,j)/sum
      enddo
    enddo
    return
end subroutine acc