function  [dir0,k]=finddir(aa,bb,g,f,dir,si,ni)
%判斷兩點之間是否有往返方向觀測值
m=length(dir);
k=0;
dir0=0;
for i=1:m
        if(g(i)==aa & f(i)==bb)
                dir0=dir(i);
                k=1;
                return
        end
end
return