Esope source files

   1 : C CLVECA    SOURCE    CHAT      05/01/12    22:08:58     5004           
   2 :       SUBROUTINE CLVECA (SIGMA,VECZ,ANGLE,IANGLE)                               
   3 : C                                                                               
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   6 : C                                                                               
   7 :       DIMENSION SIGMA(*),SIG(3),VECZ(*)                                         
   8 : C                                                                               
   9 : C INITIALISATION                                                                
  10 : C                                                                               
  11 :       UNIT=0.0174532925199432957692D0                                           
  12 :       SIG(1)=SIGMA(1)                                                           
  13 :       SIG(2)=SIGMA(2)                                                           
  14 :       SIG(3)=SIGMA(4)                                                           
  15 : C                                                                               
  16 :       IF(IANGLE.EQ.0) THEN                                                      
  16 :         CALL VPMAX3 (SIGMA,SIGMAX,VECZ)                                         
  17 :         RETURN                                                                  
  18 :       ENDIF                                                                     
  19 : C                                                                               
  20 :       IF(IANGLE.GT.0) THEN                                                      
  20 :         CALL DIAGOD (SIG)                                                       
  21 :         ANGLE=SIG(3)                                                            
  22 :         IF(IANGLE.EQ.2) ANGLE=ANGLE*UNIT                                        
  23 :         RETURN                                                                  
  24 :       ENDIF                                                                     
  25 : C                                                                               
  26 :       END                                                                       
  27 :