Esope source files

   1 : C INVA      SOURCE    CHAT      05/01/13    00:42:09     5004           
   2 : C INVA     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE INVA(S1,DJ2,DI1,CO3T,DJ3,NSTRS)
   4 : C             ****************************************************
   5 : C             *  CALCUL DES IVARIANT : I1, J2, J3 et COS(3TETA)  *
   6 : C             *           (  MODEL D'OTTOSEN   )                 *
   7 : C             ****************************************************
   8 :       IMPLICIT INTEGER(I-N)
   9 :       IMPLICIT REAL*8(A-H,O-Z)
  10 :       DIMENSION S1(NSTRS)
  11 : C
  11 :       CALL PARINV(S1,SX,SY,SZ,SXY,SXZ,SYZ,NSTRS)
  12 : C
  13 :       DI1=SX+SY+SZ
  14 :       SSX=SX-DI1/3.D0
  15 :       SSY=SY-DI1/3.D0
  16 :       SSZ=SZ-DI1/3.D0
  17 : C
  18 :       DJ2=(SSX*SSX+SSY*SSY+SSZ*SSZ)/2.D0+SXY*SXY+SXZ*SXZ+SYZ*SYZ
  19 :       DJ=SSX*SSY*SSZ+2.D0*SXY*SXZ*SYZ-SSX*SYZ*SYZ-SSY*SXZ*SXZ
  20 :       DJ3=DJ-SSZ*SXY*SXY
  21 :       CO3T=1.5D0*DSQRT(3.D0)*DJ3/DSQRT(DJ2*DJ2*DJ2)
  22 : C
  23 :       IF(CO3T .LT. -1.D0) CO3T=-1.D0
  24 :       IF(CO3T .GT. 1.D0) CO3T=1.D0
  25 : C---------------------------------------------------------------------
  26 :       IECR=0
  27 :       IF(IECR.EQ.3) THEN
  28 :       WRITE(*,100) SX,SY,SZ,SXY,SXZ,SYZ,DI1,DJ2,DJ3,CO3T
  29 :   100 FORMAT('SX=',E9.3,1X,'SY=',E9.3,1X,'SZ=',E9.3,1X,'SXY=',E9.3,1X
  30 :      *,'SXZ=',E9.3,1X,'SYZ=',E9.3,2X,/,'DI1=',E9.3,2X,'DJ2=',E9.3,2X,
  31 :      * 'DJ3=',E9.3,1X,'CO3T=',E10.4)
  32 :       ENDIF
  33 : C---------------------------------------------------------------------
  34 :       RETURN
  35 :       END
  36 :  
  37 :