% plot_cap.m
ici=pwd;

cd e:\Gatti2\golts\Survey\golts_2002
ses='212';
eval(['load TGS0',ses,'skaphi.mat'])

x=(1:1:nEnsembles)';
zc=(-244:4:-8)'; 
intmoy=(int1+int2+int3+int4)/4;
cormoy=(cor1+cor2+cor3+cor4)/4;

% avoir les courants en m.s-1
eastCurrentComponent=eastCurrentComponent/100;
northCurrentComponent=northCurrentComponent/100;
speedCurrentComponent=speedCurrentComponent/100;

cap=sensors(:,1);

figure
subplot(311)
pcolor(x,zc,flipud(speedCurrentComponent'))
shading interp, c1=colorbar;
Pc=get(c1,'position');
Pc(3)=0.03;
Pp1=get(gca,'position');
set(c1,'position',Pc,'fontsize',12)
title(['GOLTS Dec 2002 - TGS0',ses,''],'fontsize',12)
ylabel(['amplitude (m.s-1)'])
set(gca,'xlim',[x(1) x(end)],'Xtick',[x(1):50:x(end)],'layer','top','fontsize',12)

subplot(312)
plot(x,cap)
hold on
plot(x,vesselRefSpeed(:,1),'r')
plot(x,speedCurrentComponent(:,19)*1000,'g')
Pp=get(gca,'position');
Pp(3)=Pp1(3);
set(gca,'position',Pp)
set(gca,'xlim',[x(1) x(end)],'Xtick',[x(1):50:x(end)],'layer','top','fontsize',12)

subplot(313)
plot(x,-btSpeed(:,2))
hold on
plot(x,bottomRange(:,1),'+k')
plot(x,vtgSpeed(:,2),'r')
plot(x,ggaSpeed(:,2),'g')
Pp=get(gca,'position');
Pp(3)=Pp1(3);
set(gca,'position',Pp)
set(gca,'xlim',[x(1) x(end)],'Xtick',[x(1):50:x(end)],'layer','top','fontsize',12)

cd (ici)



