Sources Esope | Cast3M

tuyfrm
C TUYFRM    SOURCE    CHAT      05/01/13    03:57:13     5004
      SUBROUTINE TUYFRM(IWRK,KERRE,VEC1,VEC2)
C-----------------------------------------------------------------------
C  ROUTINE DE REMISE A JOUR DES VECTEURS DEFINISSANT LE REPERE LOCAL
C           ET L'AXE DE LA FISSURE POUR LE TUYAU FISSURE
C                HYPOTHESE DES PETITES ROTATIONS
C-----------------------------------------------------------------------
C  ENTREE
C          IWRK    POINTEUR SUR SEGMENT DE TRAVAIL
C          VEC1    LE VECTEUR DIRECTION A REACTUALISER
C          VEC2    LE VECTEUR ORIENTATION A REACTUALISER
C  SORTIE
C          KERRE   INDICE D'ERREUR  ( 0 SI TOUT EST OK )
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
      SEGMENT/IWRK/(XDDL(LRE)*D,XE(3,NBNN)*D,WORK(LW)*D)
      DIMENSION VEC1(*),VEC2(*),OMEG(3)
C
      KERRE=0
      IF (IDIM.EQ.3) THEN
      DO 1 I=1,IDIM
  1   OMEG(I)=0.5D00*(XDDL(3+I)+XDDL(9+I))
      CALL ROTVEC(VEC1,OMEG,IDIM,KERRE)
      IF(KERRE.NE.0) GO TO 666
      CALL ROTVEC(VEC2,OMEG,IDIM,KERRE)
      ENDIF
C
 666  RETURN
      END