tms2xyz

 [xp,yp,zp]=tms2xyz(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]=tms2xyz(x,y,z,l)
0002 % [xp,yp,zp]=tms2xyz(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 sinl=sin(l); cosl= cos(l);
0017 
0018 xp=+cost.*x-sint.*z;
0019 yp=y;
0020 zp=+sint.*x+cost.*z;
0021 
0022 return