uen2tms

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

 
 
 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
 To local topocentric Up East and North components.

0001 function [u,e,n]=tms2uen(x,y,z,l)
0002 % [xp,yp,zp]=tms2xyz(x,y,z,l)
0003 %
0004 %
0005 %
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 % To local topocentric Up East and North components.
0010 
0011 
0012 %xp=-y; yp=x;
0013 
0014 %x=xp; y=yp;
0015 %t=pi/2-l;
0016 
0017 %sint=sin(t); cost=cos(t);
0018 sinl=sin(l); cosl= cos(l);
0019 
0020 u=+cosl.*x-sinl.*z;
0021 e=y;
0022 n=sinl.*x+cosl.*z;
0023 
0024 return