%desfase
clear
addpath /Francia/Velocidades
load uvmoy
addpath /Francia/Estadisticas%aqui esta nanmean
vq=uvmoy(4,378:end)';vq=completo_series1(vq,5);
vq=vq(2:end-1);
uq=uvmoy(4,1:377)';uq=completo_series1(uq,5);
uq=uq(2:end-1);
load N3u
load S1u
load S1v
load Y3v
load resnueva
addpath /Francia/vientos
load Pente_vent.mat
topexn=completo_series1(TP_J_J_pente1_v,5);
vienton=completo_series1(uwind23WTP,5);
fet=date_windTP; feti=fet(1):fet(end);
topex=interp1(fet,topexn,feti);
viento=interp1(fet,vienton,feti);
topex=completo_series1(topex,5);
viento=completo_series1(viento,5);
load topex
load vientochristine
load zwstotalf
zwsersquick=zwstotalf(431:end-652);

fety1=feti;
xxx=feti(1):300:feti(end);
fechaA=datestr(xxx,12);
%y2=resnueva(1:294,4);
%y1=vq(1:375);
%y2=uq(1:375);
y1=vientochristine';
%y1=topex';
y2=zwsersquick';
%y1=uq;
%y2=vq;

exp1=2;
exp2=10;
inter=40;
y=y1;
abc
%aqui se escoge el periodo en el que se quiere el desfase
%aa=aa(180:190);
%period=period(180:190);
%aa=aa(1:300);
%period=period(1:300);
%aa=aa(252:256);
%period=period(252:256);
conv_fft_abc

y1=y;
convolution1=wave;
y=y2;
%abc
conv_fft_abc;
y2=y;
convolution2=wave;
integration_essai;
%phase=unwrap(phase);
delay=(t.*phase/(2*pi));
for i=1:n1;
   xx(i)=i+delay(i);
%xx(i)=i-delay(i);
end;


subplot(2,1,1)
plot(fety1,y1)
hold on
plot(fety1,y2,'r')
set(gca,'xlim',[min(feti)-3  max(feti)+3],'xtick',xxx,'xticklabel',[])
title('Normalized meridional component of the north and south signals')

subplot(2,1,2)
plot(fety1,delay,'k')
ylabel('delay (days)')

set(gca,'xlim',[min(feti)-3 max(feti)+3],'xtick',xxx,'xticklabel',fechaA)

figure
aver=1:321;
contourf(feti,aver,abs(cross_wave))
caxis([0 1])
view(0,-90)
%colorbar
hold on
%plot([1 200],[75 75],'w')
%ylim([0 240])
%r=30:30:240;
ylim([0 300])
r=30:30:300;
%set(gca,'YTick',[6 46 70 86 99 109 118 126 133 139] ,'YTickLabel',[10:10:100])
set(gca,'YTick',r,'YTickLabel',floor(period(r)))
set(gca,'xlim',[min(feti) max(feti)],'xtick',xxx,'xticklabel',fechaA)

xlabel('days')
