Esope source files

   1 : C SIJ       SOURCE    CHAT      05/01/13    03:18:33     5004           
   2 :       SUBROUTINE SIJ(TENS,SOR,VAL)                                              
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION TENS(3),SOR(3)                                                  
   6 : C                                                                               
   7 : C------------------------------------------                                     
   8 : C    CALCUL DU TENSEUR SIJ ET DE SIGMA*                                         
   9 : C------------------------------------------                                     
  10 : C                                                                               
  11 :       SOR(1)=(2.D0*TENS(1)-TENS(2))/3.D0                                        
  12 :       SOR(2)=(2.D0*TENS(2)-TENS(1))/3.D0                                        
  13 :       SOR(3)=TENS(3)                                                            
  14 :       VAL=SQRT(TENS(1)*TENS(1)+TENS(2)*TENS(2)-TENS(1)*TENS(2)                  
  15 :      1 +3.D0*TENS(3)*TENS(3))                                                   
  16 :       RETURN                                                                    
  17 :       END                                                                       
  18 :