%
% Plots the evolution of maximum wind speed for each ensemble member.
%
%  Copyright WindRiskTech, L.L.C., October 2010
%-----------------------------------------------------------------------
transfac=0.7;  %  Fraction of translation speed to add to circular wind
%-----------------------------------------------------------------------
clf
%
load('temp.mat')
load('sorted')
[ut,vt,jmax] = utrans(latstore,longstore);
vnet=vstore+transfac.*sqrt(ut.^2+vt.^2);
[q,dum]=size(vstore);
hold on
for i=1:q,
    clear time
    time=[0:(jmax(i)-1)]./12;
    h=plot(time,vnet(i,1:jmax(i)),'linewidth',1,'color','b');
end    
hold off
bold
xlabel('Days Since Initialization')
ylabel('Maximum Wind Speed (Knots)')
tit1=[storm(1:2) ' ' storm(3:4) ' ' datea(1:4) ' ' datea(5:6) ' ' datea(7:8)...
    ' ' datea(9:10) ' GMT'];
title(tit1,'fontsize',12,'fontweight','bold')
set(gca,'color',[0.8 0.8 0.8])
