PRO trsk2,x,y,z,olt,oln,oht,alpha,xp,yp,zp

; Modelled on subroutine TRSK2.FOR in GEOLB.FOR
; Angles input in degrees, converted to radians for computations

trsk1,olt,oln,oht,xo,yo,zo
;x3=x-xo		; These lines commented out 950302
;y3=y-yo		; to save memory, at the expense of
;z3=z-zo		; redundant computations below
colt=COS(olt*!DTOR)
solt=SIN(olt*!DTOR)
coln=COS(oln*!DTOR)
soln=SIN(oln*!DTOR)
cang=COS(alpha*!DTOR)
sang=SIN(alpha*!DTOR)

tmp1=(x-xo)*coln+(y-yo)*soln
tmp2=-(x-xo)*soln+(y-yo)*coln
xp=-cang*((tmp1*solt)-(z-zo)*colt)+sang*tmp2
yp=-sang*((tmp1*solt)-(z-zo)*colt)-cang*tmp2
zp=tmp1*colt+(z-zo)*solt
RETURN
END