1 : C KNORM2 SOURCE CHAT 05/01/13 00:58:18 5004 2 : SUBROUTINE KNORM2(NES,CM,NCM,U,S) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : DIMENSION CM(NES,NCM),U(*),V(2) 6 : C 7 : DO 10 K=1,NES 8 : V(K)=CM(K,2)-CM(K,1) 9 : 10 CONTINUE 10 : S=0 11 : DO 20 K=1,NES 12 : S=S+V(K)*V(K) 13 : 20 CONTINUE 14 : S=DSQRT(S) 15 : U(1) = -V(2)/S 16 : U(2) = V(1)/S 17 : U(3) = -(U(1)*CM(1,1)+U(2)*CM(2,1)) 18 : RETURN 19 : END 20 :
© Cast3M 2003 - All rights reserved.
Disclaimer