	SUBROUTINE ROTA(R,PSI,TETA,PHI)
C
C	SOUS PROGRAMME DE CALCUL DE NOUVELLES COORDONNEES
C	APRES UNE ROTATION DES 3 ANGLES D'EULER(PSI,TETA,PHI)
C	CES ANGLES SONT EXPRIMES EN RADIANS.
C
C	SOUS PROGRAMMES REQUIS POUR LE TKB: NEANT.
C
	IMPLICIT REAL*8(A-H,O-Z)
	DIMENSION R(3)
	X1=R(1)
	Y1=R(2)
	Z1=R(3)
	C1=DCOS(PSI)
	C2=DCOS(TETA)
	C3=DCOS(PHI)
	S1=DSIN(PSI)
	S2=DSIN(TETA)
	S3=DSIN(PHI)
	A=C3*C1-S3*S1*C2
	B=C3*S1+S3*C1*C2
	C=S3*S2
	X2=A*X1+B*Y1+C*Z1
	A=-S3*C1-C3*S1*C2
	B=-S3*S1+C3*C1*C2
	C=C3*S2
	Y2=A*X1+B*Y1+C*Z1
	A=S1*S2
	B=-S2*C1
	C=C2
	Z2=A*X1+B*Y1+C*Z1
	R(1)=X2
	R(2)=Y2
	R(3)=Z2
	RETURN
	END
