%ondeletes para ver escalas verticales
%13 feb 2003

addpath D:\Francia\Programitas
addpath D:\Francia\Chantal

load eq99_uvctd

s=eq99_uvctd(2201:end-400,5);
t=eq99_uvctd(2201:end-400,4);
z=eq99_uvctd(2201:end-400,2);
u=eq99_uvctd(:,9);
v=eq99_uvctd(2201:end-400,9);
uu=reshape(u,200,18);
b=-30:3:30;
zz=1:10:2000;
x=[-6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.67 -.33 0 .33 .67 1 1.5];
figure
contourf(x,-zz,uu,b)
colorbar
set(gca, 'ylim',[-1600 -600])
load uu6e
y=uu6e(20:end);
%y=restemp(:,3);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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


figure;
%initfig;
r=20:20:200;

subplot(2,1,1)
plot(y,'r')
xlabel('depth (m)')
ylabel('u component Equator E99')
set(gca,'XTick',r,'XTickLabel',r*10)

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 250])
%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)*16))
set(gca,'XTick',r,'XTickLabel',r*16)

xlabel('depth (m)')
ylabel('meters')

wave1=wave;
wave(1:95,:)=0;
wave(105: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)

