%cdf1='44013s-1a.cdf';
cdf1='3412-a1h.cdf';
u1='east';
v1='north';
files=[3401];
u2='east';
v2='north';
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');
