trunk/filtrez/acc.f


! $Header: /home/cvsroot/LMDZ4/libf/filtrez/acc.F,v 1.1.1.1 2004/05/19
! 12:53:09 lmdzadmin Exp $

SUBROUTINE acc(vec, d, im)
  implicit none
  integer im, i, j
  real vec(im, im), d(im)
  real sum, ssum

  DO j = 1, im
    DO i = 1, im
      d(i) = vec(i, j)*vec(i, j)
    END DO
    sum = ssum(im, d, 1)
    sum = sqrt(sum)
    DO i = 1, im
      vec(i, j) = vec(i, j)/sum
    END DO
  END DO
  RETURN
END SUBROUTINE acc
1
2	! $Header: /home/cvsroot/LMDZ4/libf/filtrez/acc.F,v 1.1.1.1 2004/05/19
3	! 12:53:09 lmdzadmin Exp $
4
5	SUBROUTINE acc(vec, d, im)
6	implicit none
7	integer im, i, j
8	real vec(im, im), d(im)
9	real sum, ssum
10
11	DO j = 1, im
12	DO i = 1, im
13	d(i) = vec(i, j)*vec(i, j)
14	END DO
15	sum = ssum(im, d, 1)
16	sum = sqrt(sum)
17	DO i = 1, im
18	vec(i, j) = vec(i, j)/sum
19	END DO
20	END DO
21	RETURN
22	END SUBROUTINE acc