function [xi,yi,zi]=pdrog(xd,yd,zd,jd,s1,s2);
 jd1=julian(s1);
 ind=tind(jd,s1,s2);
 xi=xd(ind,:);
 yi=yd(ind,:);
 zi=zd(ind,:);
 n=length(ind);
 line(xi,yi);
 line(xi(1,:),yi(1,:),'linestyle','x','markersize',5,'color',[1 1 0]);
 line(xi(n,:),yi(n,:),'linestyle','o','markersize',5,'color',[1 1 .99]);


