C PRINC     SOURCE    CHAT      05/01/13    02:27:24     5004           
C PRINC     SOURCE    INSL       24/10/96    
      SUBROUTINE PRINC(S,T,NSTRS)
C
C   CALCUL DES CONTRAINTES ET DIRECTIONS PRINCIPALES
C    L'ANGLE PHI=T(4) : EST EN DEGRE
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION S(4),T(4)
      PI=4.D0*DATAN(1.D0)
      A1=S(1)+S(2)
      A2=S(1)-S(2)
      A3=S(3)
      A4=S(3)
      IF(NSTRS.EQ.4.OR.NSTRS.EQ.6) A4=S(4)
      IF(DABS(A4) .GE. 1.D-8) GOTO 11
      IF(A2 .LT. 1.D-8) GOTO 121
      T(1)=S(1)
      T(2)=S(2)
      T(3)=S(3)
      T(4)=0.0D+00
      GO TO 10
  121 T(1)=S(2)
      T(2)=S(1)
      T(3)=S(3)
      T(4)=90.0D+00
      GO TO 10
   11 IF(DABS(A2) .GE. 1.D-8) GOTO 111
      T(1)=S(1)+DABS(A4)
      T(2)=S(2)-DABS(A4)
      T(3)=S(3)
      T(4)=45.0D+00*DABS(A4)/A4
      GO TO 10
  111 A5=DSQRT(A2*A2+4.0D+00*A4*A4)
      T(1)=(A1+A5)/2.0D+00
      T(2)=(A1-A5)/2.0D+00
      T(3)=S(3)
      T(4)=DATAN(2.0D+00*A4/A2)*90.0D+00/PI
      IF(A2 .GE. 1.D-8) GOTO 10
      T(4)=T(4)+90.0D+00*DABS(A4)/A4
  10  CONTINUE
      IF(NSTRS.EQ.3.OR.NSTRS.EQ.5) T(3)=0.D0
      RETURN
      END