PRO PoCarot,az,za,x4,y4,z4

;rdcols,'PoCa.fov',h,az1,za1
;az = [az,az1]
;za = [za,za1]
r=1.  
x2=0. & x3=0.& x4=0.& y2=0.& y3=0.& y4=0.
az = az*!dtor
za = za*!dtor
a=-13.00*!dtor
b=-12.00*!dtor 
c=+06.00*!dtor

x1 = r*cos(90-az)*sin(za)
y1 = r*sin(90-az)*sin(za)
z1 = r*cos(za)

rot_3d,3,x1,y1,z1,a,x2,y2,z2
rot_3d,2,x2,y2,z2,b,x3,y3,z3
rot_3d,1,x3,y3,z3,c,x4,y4,z4
az = az * !radeg
za = za * !radeg
RETURN
END