xyz2tms

 [xp,yp,zp]=xyz2tms(x,y,z,l)

 x-east, y-north, z-up seems the natural way
 to define array coords and baselines to me.
 Convert these for an array at latitude l to the coord
 frame defined by TMS fig 4.1 - x,y in equatorial plane
 with x to meridian, y east and z to north celestial pole

0001 function [xp,yp,zp]=xyz2tms(x,y,z,l)
0002 % [xp,yp,zp]=xyz2tms(x,y,z,l)
0003 %
0004 % x-east, y-north, z-up seems the natural way
0005 % to define array coords and baselines to me.
0006 % Convert these for an array at latitude l to the coord
0007 % frame defined by TMS fig 4.1 - x,y in equatorial plane
0008 % with x to meridian, y east and z to north celestial pole
0009 
0010 xp=-y; yp=x;
0011 
0012 x=xp; y=yp;
0013 t=pi/2-l;
0014 
0015 sint=sin(t); cost=cos(t);
0016 
0017 xp=+cost.*x+sint.*z;
0018 yp=y;
0019 zp=-sint.*x+cost.*z;
0020 
0021 return