Sources Esope | Cast3M

invaj2
C INVAJ2    SOURCE    CB215821  22/11/16    21:15:02     11500          
C INVIJ2    SOURCE    CHAT      05/01/13    00:43:06     5004
 
      SUBROUTINE INVARI_J2 (SIG,N,J2)
      IMPLICIT INTEGER(I-N)
      integer n,i
      real*8 sig(n),J2,desv(6)
      call Desviador(sig,desv,n)
      J2=0.D0
      do i=1,3
         J2=J2+desv(i)*desv(i)
      enddo
      J2=J2/2.D0
      if (n.lt.4) return
      do i=4,n
         J2=J2+desv(i)*desv(i)
      enddo
      return
      end