function eigs=mode_traj2(fname)
% Call eigs=mode_traj2('fname');
format compact
eval(['load ',fname])
eval(['eigs=',fname,';'])
[nrow,ncol]=size(eigs);
%ii0=find(eigs(:,3)>=0 | eigs(:,3)==-2);
ii0=find(eigs(:,3)>=0);
iiall=find(eigs(:,3)>=-2 & eigs(:,3)<=1);
iimiss=find(eigs(:,3)==-1);
iieig=find(eigs(:,3)==1);
iitraj=find(eigs(:,3)==0);
figure
hold on
% Plot branch points:
iibp=find(eigs(:,3)==-5);
%plot(eigs(ii,1),eigs(ii,2),'g', ...%eigs(ii0,1),eigs(ii0,2),'r-', ...
%   eigs(ii1,1),eigs(ii1,2),'g+')
plot_cfill2(eigs(iibp,1),eigs(iibp,2),'s','m',6);
text(eigs(iibp,1),eigs(iibp,2),' BP','Color','m')
plot(eigs(iiall,1),eigs(iiall,2),'g-');
plot_cfill2(eigs(iimiss,1),eigs(iimiss,2),'s','g',4);
plot_cfill2(eigs(iitraj,1),eigs(iitraj,2),'o','b',4);
plot_cfill2(eigs(iieig,1),eigs(iieig,2),'o','r',6);
text(eigs(iieig,1),eigs(iieig,2),int2str([1:length(iieig)]'))
%for j=1:length(iieig)
%   text(eigs(iieig(j),1),eigs(iieig(j),2),int2str(j))
%end
xlabel('COS(THref)'); ylabel('ATTEN (dB/km)'); grid; 
mytitle(fname);
xlabel('COS(THref)'); ylabel('ATTEN (dB/km)'); mytitle(fname);
zoom
%
% COPYRIGHT (c) 2008
% 
% Applied Research Laboratories, The University of Texas at Austin
% (ARL:UT)
% All rights reserved.
% 
% This material may be reproduced by or for the U.S. Government
% only. This notice must appear in all copies of this file and its
% derivatives. Further distribution outside the U.S. Government is
% prohibited without the explicit written authorization from the
% Director of the Environmental Sciences Laboratory, ARL:UT.
% 
% THIS SOFTWARE IS PROVIDED "AS IS" WITHOUT WARRANTY OF ANY KIND,
% EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO THE
% IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR
% PURPOSE, OR NONINFRINGEMENT.
% 
% IN NO EVENT WILL ARL:UT BE LIABLE FOR ANY DAMAGES WHATSOEVER ARISING
% OUT OF THE USE OR INABILITY TO USE THIS SOFTWARE, EVEN IF ARL:UT HAS
% BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
