%filtro
clear
addpath /Francia/Velocidades
addpath /Francia/annie
addpath /Francia/Estadisticas

load u60
load v60
load u230
load v230
load u800
load v800
load u1190
load v1190
load u1300
load v1300
load v15010w
load u15010w

for i=1:length(u60)-10
    u60m(i)=mean(u60(i:i+10));
       v60m(i)=mean(v60(i:i+10));
   end
   

for i=1:length(u230)-10
    u230m(i)=mean(u230(i:i+10));
       v230m(i)=mean(v230(i:i+10));
       u1190m(i)=mean(u1190(i:i+10));
       v1190m(i)=mean(v1190(i:i+10));
       u15010wm(i)=mean(u15010w(i:i+10));
       v15010wm(i)=mean(v15010w(i:i+10));

   end 
   for i=1:length(u800)-10
       u800m(i)=mean(u800(i:i+10));
       v800m(i)=mean(v800(i:i+10));
       %u1300m(i)=mean(u1300(i:i+15));
       %v1300m(i)=mean(v1300(i:i+15));
       

end

tt23=1:372;


tt=tt23(6:end-5);
x=tt(1):30:tt(end);
a=datestr(x,12);

subplot(2,1,1)
plot(tt,u230m)
hold on
plot(tt(1:length(u800m)),u800m,'k')
hold on
plot(tt,u1190m,'r')
hold on
%plot(tt(1:length(u800m)),u1300m,'g')
hold on
plot(tt(1:length(u60m)),u60m,'c')
hold on
plot(tt,u15010wm,'g')



%set(gca,'xlim',[min(tt1)-3 max(tt1)+3],'ylim',[-20 30],'xtick',x,'xticklabel',[])
%ylabel('\delta u cm/s','fontsize',12)
%xlabel('r=-0.8')
subplot(2,1,2)
plot(tt,v230m)
hold on
plot(tt(1:length(u800m)),v800m,'k')
hold on
plot(tt,v1190m,'r')
hold on
%plot(tt(1:length(u800m)),v1300m,'g')
hold on
plot(tt(1:length(u60m)),v60m,'c')
hold on
plot(tt,v15010wm,'g')


%set(gca,'xlim',[min(tt1)-3 max(tt1)+3],'ylim',[-20 20],'xtick',x,'xticklabel',a,'fontsize',11)
%ylabel('\delta u cm/s','fontsize',12)
load vientochristinescasa
load topexscasa

for i=1:length(topexscasa)-100
    topexscasaf(i)=mean(topexscasa(i:i+100));
       vientochristinescasaf(i)=mean(vientochristinescasa(i:i+100));
   end
   
   load zwstotal
   load ttotal
   for i=1:length(zwstotal)-50
    zwstotalf(i)=mean(zwstotal(i:i+50));
   end
   ttotalf=ttotal(26:end-25);