1 : C TUYFRM SOURCE CHAT 05/01/13 03:57:13 5004 2 : SUBROUTINE TUYFRM(IWRK,KERRE,VEC1,VEC2) 3 : C----------------------------------------------------------------------- 4 : C ROUTINE DE REMISE A JOUR DES VECTEURS DEFINISSANT LE REPERE LOCAL 5 : C ET L'AXE DE LA FISSURE POUR LE TUYAU FISSURE 6 : C HYPOTHESE DES PETITES ROTATIONS 7 : C----------------------------------------------------------------------- 8 : C ENTREE 9 : C IWRK POINTEUR SUR SEGMENT DE TRAVAIL 10 : C VEC1 LE VECTEUR DIRECTION A REACTUALISER 11 : C VEC2 LE VECTEUR ORIENTATION A REACTUALISER 12 : C SORTIE 13 : C KERRE INDICE D'ERREUR ( 0 SI TOUT EST OK ) 14 : C----------------------------------------------------------------------- 15 : IMPLICIT INTEGER(I-N) 16 : IMPLICIT REAL*8(A-H,O-Z) 17 : -INC CCOPTIO 18 : SEGMENT/IWRK/(XDDL(LRE)*D,XE(3,NBNN)*D,WORK(LW)*D) 19 : DIMENSION VEC1(*),VEC2(*),OMEG(3) 20 : C 21 : KERRE=0 22 : IF (IDIM.EQ.3) THEN 23 : DO 1 I=1,IDIM 24 : 1 OMEG(I)=0.5D00*(XDDL(3+I)+XDDL(9+I)) 24 : CALL ROTVEC(VEC1,OMEG,IDIM,KERRE) 25 : IF(KERRE.NE.0) GO TO 666 25 : CALL ROTVEC(VEC2,OMEG,IDIM,KERRE) 26 : ENDIF 27 : C 28 : 666 RETURN 29 : END 30 :
© Cast3M 2003 - All rights reserved.
Disclaimer