1 : C CINVJ2 SOURCE CHAT 05/01/12 22:02:56 5004 2 : 3 : SUBROUTINE C_INVARI_J2 (SIG,N,J2) 4 : IMPLICIT INTEGER(I-N) 5 : integer n,i 6 : complex*16 sig(n),J2,desv(6) 6 : call C_Desviador(sig,desv,n) 7 : J2=0.D0 8 : do i=1,3 9 : J2=J2+desv(i)*desv(i) 10 : enddo 11 : J2=J2/2.D0 12 : if (n.lt.4) return 13 : do i=4,n 14 : J2=J2+desv(i)*desv(i) 15 : enddo 16 : return 17 : end 18 : 19 : 20 :
© Cast3M 2003 - All rights reserved.
Disclaimer