%Vientos a 23 west
%Series sin ciclo anual correspondientes a la fecha de los mouillages
clear
load v23windcompletosca %viento 23 W sin ciclo anual pero la serie completa, tres anos.
load win23w.dat %todos los datos, las dos componentes, fecha... aqui lo tengo por lo de las fechas
load u23windcompletoscasa
load tf
in1=find(datenum(2001,12,18)==tf); %fecha de la serie cortada
in2=find(datenum(2002,12,17)==tf); %fecha final de la serie cortada
time=tf(in1:in2);
 
w23cortadoscaf=v23windcompletosca(in1:in2);
wu23cortadoscasaf=u23windcompletoscasa(in1:in2);
figure
plot(time,w23cortadoscaf)
hold on
plot(time,wu23cortadoscasaf,'c')
xx=time(1):30:time(end);
aa=datestr(xx,12);
set(gca,'xtick',xx,'xticklabel',aa)