--- trunk/filtrez/acc.f 2013/11/15 18:45:49 76 +++ trunk/filtrez/acc.f90 2014/03/05 14:38:41 81 @@ -1,14 +1,19 @@ -! -! $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) - dimension vec(im,im),d(im) - do 10 j=1,im - do 9 i=1,im - 9 d(i)=vec(i,j)*vec(i,j) - sum=ssum(im,d,1) - sum=sqrt(sum) - do 10 i=1,im - 10 vec(i,j)=vec(i,j)/sum - return - end + +! $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) + DIMENSION vec(im, im), d(im) + + 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