  % Decomposition en ondelettes de la NAO. 

clear all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Chargement des donnees a etudier

%load /usr/home/lbulod/Francia/Velocidades/S4v

load /Francia/annie/u225pr

addpath /Francia/Estadisticas%aqui esta nanmean
addpath /Francia/Velocidades
load /Francia/imagenest/ll

y=u225pr';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

exp1=2;
exp2=7;
inter=40;
abc
conv_fft_abc


figure;
%initfig;

subplot(2,1,1)
plot(y,'r')
xlabel('days')
ylabel('u component S1 840m')
hold on

subplot(2,1,2)
contourf(abs(wave),[0:0.05:1])
caxis([0 1])
view(0,-90)
%colorbar
hold on
%plot([1 200],[75 75],'w')
ylim([0 200])
r=30:30:200;
%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)))
xlabel('days')
ylabel('period(days)')

wave1=wave;
wave(1:100,:)=0;
wave(110:end,:)=0;
maxmax
plot(dilat_number,'w','LineWidth',2)


subplot(2,1,1)
plot(ridge,'k','LineWidth',2)
%NAO_filtree=real(wave(75,:));
%plot(NAO_filtree,'r')
A=corrcoef(y,ridge);
pourvar=round(A(1,2)^2*1000)/10
tt1=['Signal a 33 ans extrait (' num2str(pourvar) '% de variance expliquee)']

%legend('NAO moyenne annuelle',tt1)



