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. & z1=0.
az = az*!DTOR
za = za*!DTOR
a=-13.20*!DTOR
b=04.00*!DTOR 
c=+00.00*!DTOR

polrec3d,r,za,((90*!dtor)-az),x1,y1,z1 ; from jhu library

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