Esope source files

   1 : C ENDOCA    SOURCE    CHAT      05/01/12    23:38:57     5004           
   2 :       SUBROUTINE ENDOCA(EPOUSI,A,IMOD)                                          
   3 : C_______________________________________________________________________        
   4 : C                                                                               
   5 : C     CHARGEMENT DE A POUR LA DIAGONALISATION JACOB3                            
   6 : C                                                                               
   7 : C     PP 27/8/91                                                                
   8 : C                                                                               
   9 : C     18/10/91 se souvenir que GAMMA a encore frappe !                          
  10 : C     =====> diviser par 2 les termes de cisaillement                           
  11 : C     ET DISCOSSIER EPS ET STR: IMOD=1--->STR IMOD=2--->EPS                     
  12 : C_______________________________________________________________________        
  13 : C                                                                               
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  16 :       DIMENSION EPOUSI(*),A(3,3)                                                
  17 : *                                                                               
  18 : -INC CCOPTIO                                                                    
  19 : C                                                                               
  19 :       CALL ZERO(A,3,3)                                                          
  20 :       A(1,1)=EPOUSI(1)                                                          
  21 :       A(2,2)=EPOUSI(2)                                                          
  22 :       A(3,3)=EPOUSI(3)                                                          
  23 :       A(1,2)=EPOUSI(4)/IMOD                                                     
  24 :       A(2,1)=A(1,2)                                                             
  25 :       IF(IFOUR.GE.1.OR.IFOUR.LE.-3) THEN                                        
  26 :         A(3,1)=EPOUSI(5)/IMOD                                                   
  27 :         A(1,3)=A(3,1)                                                           
  28 :         A(3,2)=EPOUSI(6)/IMOD                                                   
  29 :         A(2,3)=A(3,2)                                                           
  30 :       ENDIF                                                                     
  31 : C                                                                               
  32 :       RETURN                                                                    
  33 :       END                                                                       
  34 :