function plotpl
global ephpdir utc_tai
[d,f]=fileparts(which('ephpl')); ephpdir=fullfile(d,'data'); utc_tai=-32; 
figure, box on, axis equal, hold on
t0=caltomjd([2000,1,1]);
for pl=1:9
    for t=1:365*20
        p(:,t)=ephpl(t0+t-1,pl);
    end
    plot3(p(1,:),p(2,:),p(3,:),'.','markersize',3);
end
set(gca,'position',[0.06,0.02,0.88,0.96]);
