SUBROUTINE SETORG(X,Y,Z,XT,YT,ZT) LOGICAL JUMP COMMON/CSETOR/JUMP,EX,EY,EZ,AX,AY,AZ,D,R,COSBE,COSAL,COSGA AX=X AY=Y AZ=Z EX=XT EY=YT EZ=ZT DX=AX-EX DY=AY-EY DZ=AZ-EZ D=SQRT(DX*DX+DY*DY+DZ*DZ) COSAL=DX/D COSBE=DY/D COSGA=DZ/D BE=ACOS(COSBE) GA=ACOS(COSGA) SINGA=SIN(GA) IF(SINGA.LT.0.0001) GOTO 10 R=1./SINGA JUMP=.FALSE. RETURN 10 SINBE=SIN(BE) R=1./SINBE JUMP=.TRUE. RETURN END