C CDRF      SOURCE    CHAT      05/01/12    21:51:41     5004           
C CDRF     SOURCE    INSL       24/10/96
      SUBROUTINE CDRF(S1,PHI,V1)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION S1(3),V1(3),TS(3,3)
C
C         ************************************
C         * DETERMINE LE VECTEUR CONTRAINTE  *
C         * DANS UN REPERE LIE  A LA FISSURE *
C         ************************************
C
      CALL ZERO(V1,3,1)
      PI=4.D0*DATAN(1.D0)
      PHIC=(PHI+90.D0)*PI/180.D0
      A=DCOS(PHIC)
      B=DSIN(PHIC)
*
      TS(1,1)=A*A
      TS(1,2)=B*B
      TS(1,3)=2.D0*A*B
      TS(2,1)=B*B
      TS(2,2)=A*A
      TS(2,3)=-2.D0*A*B
      TS(3,1)=-1.D0*A*B
      TS(3,2)=B*A
      TS(3,3)=A*A-B*B
      CALL BST(TS,S1,3,3,V1)
C
      RETURN
      END