Esope source files

   1 : C CDRF      SOURCE    CHAT      05/01/12    21:51:41     5004           
   2 : C CDRF     SOURCE    INSL       24/10/96
   3 :       SUBROUTINE CDRF(S1,PHI,V1)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8 (A-H,O-Z)
   7 :       DIMENSION S1(3),V1(3),TS(3,3)
   8 : C
   9 : C         ************************************
  10 : C         * DETERMINE LE VECTEUR CONTRAINTE  *
  11 : C         * DANS UN REPERE LIE  A LA FISSURE *
  12 : C         ************************************
  13 : C
  13 :       CALL ZERO(V1,3,1)
  14 :       PI=4.D0*DATAN(1.D0)
  15 :       PHIC=(PHI+90.D0)*PI/180.D0
  16 :       A=DCOS(PHIC)
  17 :       B=DSIN(PHIC)
  18 : *
  19 :       TS(1,1)=A*A
  20 :       TS(1,2)=B*B
  21 :       TS(1,3)=2.D0*A*B
  22 :       TS(2,1)=B*B
  23 :       TS(2,2)=A*A
  24 :       TS(2,3)=-2.D0*A*B
  25 :       TS(3,1)=-1.D0*A*B
  26 :       TS(3,2)=B*A
  27 :       TS(3,3)=A*A-B*B
  27 :       CALL BST(TS,S1,3,3,V1)
  28 : C
  29 :       RETURN
  30 :       END
  31 : 
  32 :  
  33 :