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