%cdf1='44013s-1a.cdf';
cdf1='44013tau-a.cdf';
u1='east_str';
v1='north_str';
files=[3391 3392 3411 3412 3481 3482 3591 3581];
u2='east';
v2='north';
axis([-8 8 -8 8]);
n=length(files);
for i=1:n,
  cdf2=[int2str(files(i)) '-a1h.cdf']
  [theta(i),maj(i),min(i)]=doellipse(cdf1,u1,v1,cdf2,u2,v2);
end
%cdf1='logan-a.cdf';
%cdf2='3751-a1h.cdf';
%[A,B,hits]=dobinrot(cdf1,'east_str','north_str',cdf2,'east','north');
%cdf2='3752-a1h.cdf';
%[A,B,hits]=dobinrot(cdf1,'east_str','north_str',cdf2,'east','north');
