%ondeletes para ver escalas verticales usando datos equalant 2000
%13 feb 2003
clear
addpath D:\Francia\Programitas
addpath D:\Francia\Hidrologia
mu=[];
tint=1:16:141*16;
ti=1:10:141*16;ti=ti(1:end-1);
for i=23:27
        %eval(['load ','eq210w',num2str(i)])
                eval(['load ','eq',num2str(i),'.dat'])

u=['eq',num2str(i),'(:,2)'];uu=eval(u);uu=uu(1:141);

z=['eq',num2str(i),'(:,1)'];z=eval(z);z=z(1:141);
uuu=interp1(tint,uu,ti);uuu=uuu(1:200);
mu=[mu;uuu];
end
b=-30:5:30;
zi=1:10:2000;
x=[.67 .33 0 -.33 -.67];
mu=mu';
figure
contourf(x,-zi,mu,b)
hold on
ja=ones(size(zi)); ja=ja'*x;
plot(ja,-zi,'*')
colorbar


y=mu(:,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('meters')
ylabel('u component Equator E00')
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 200])
r=20:20: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)*10))
set(gca,'XTick',r,'XTickLabel',r*10)

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

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

