%
% Overlays ECMWF ensemble tracks on WRT/AER events sets. 
% "Alltrmap.m" must be run first.
% Copyright WindRiskTech L.L.C., October, 2010
%
%--------------------------
file='tracks.mat';
load(file)
load('sorted')
[nens,nt]=size(la);
hold on
for i=1:nens,
         clear x y
         y=ay-ay.*(la(i,1:nn(i))+90)./180;
         x=ax.*(lo(i,1:nn(i))+360-displacement)/360;
         for j=1:nn(i),
             if x(j) > ax
                 x(j)=x(j)-ax;
             end    
         end   
    h=plot(x,y,'w','linewidth',1);
end
hold off