  % Decomposition en ondelettes de la NAO. 

clear all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Chargement des donnees a etudier
%load /usr/lodyc/decennal/jmignot/BERGEN/ATMOS/EOF/E75_2075-2274_eof_slp_25N80N_ANNUALmean_d3.mat
load /Francia/Velocidades/N2v
load /Francia/Velocidades/N2u
load /Francia/Velocidades/uvmoy

load /Francia/Temperatura/te80
load /Francia/Temperatura/te100
addpath /Francia/Programas
%y=PC(:,1)';
%y=reshape(y,12,200);
%y=squeeze(mean(y,1));
%y=res(:,6);
indice=isnan(te100); indi2=find(indice==0);te100=te100(indi2);
indice=isnan(te80); indi2=find(indice==0);te80=te80(indi2);
u=uvmoy(4,1:377)';in1=find(isnan(u)==0);u=u(in1);
v=uvmoy(4,378:754)';v=v(in1);
%mimo=[te80,te100];
%mimo=[N2u,N2v];
mimo=[u,v];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:2
    y=mimo(:,i);
exp1=0;
exp2=7;
inter=40;
abc
conv_fft_abc


figure(2)
%initfig;

%subplot(2,1,1)
%plot(y,'k')
%xlabel('jours')
%hold on

subplot(1,2,i)
contourf(abs(wave),[0:0.05:1])
%caxis([0 1])
view(0,-90)
%colorbar
hold on
%plot([1 200],[75 75],'w')
%ylim([0 140])
set(gca,'YTick',[35:34:280] ,'YTickLabel',[2 4 7 12 22 39 71 128])
xlabel('jours')
ylabel('periodes(jours)')
end
suptitle('Mouillage sud 1525 m u et v')
%wave1=wave;
%wave(1:230,:)=0;
%wave(245: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 62-75 extrait (' num2str(pourvar) '% de variance expliquee)']

