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 |
implicit none |
7 |
do 9 i=1,im |
integer im, i, j |
8 |
9 d(i)=vec(i,j)*vec(i,j) |
real vec(im, im), d(im) |
9 |
sum=ssum(im,d,1) |
real sum, ssum |
10 |
sum=sqrt(sum) |
|
11 |
do 10 i=1,im |
DO j = 1, im |
12 |
10 vec(i,j)=vec(i,j)/sum |
DO i = 1, im |
13 |
return |
d(i) = vec(i, j)*vec(i, j) |
14 |
end |
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 |