Esope source files

   1 : C INVART    SOURCE    CHAT      05/01/13    00:42:27     5004           
   2 : 
   3 :       SUBROUTINE INVARI_T (SIG,N,T)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer n
   6 :       real*8 sig(n),t,J2,J3
   7 :       t=0.D0
   7 :       call Invari_J2(sig,n,J2)
   8 :       if (J2.gt.0.D0) then
   8 :          call Invari_J3(sig,n,J3)
   9 :          t=3.D0*dsqrt(3.D0)*J3/(2.D0*J2**1.5D0)
  10 :          if (t.lt.-1.D0) t=-1.D0
  11 :          if (t.gt.1.D0)  t=1.D0
  12 :          t=acos(t)/3.D0
  13 :       endif
  14 :       return
  15 :       end
  16 :  
  17 :  
  18 :