Esope source files

   1 : C INVAT2    SOURCE    CHAT      05/01/13    00:42:30     5004           
   2 : 
   3 :       SUBROUTINE INVARI_T2 (SIG,N,T) !!-30:30
   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=asin(t)/3.D0
  13 :       endif
  14 :       return
  15 :       end
  16 :  
  17 :  
  18 :