% plottrlat_superp_1204.m
% plot transects de haut en bas: amplitude
% X=lon
% cree par JG le 26/09
% xx numero de station depart
% yy numero de station fin
% zz lettre pour indiquer si c'est le a: 1er passage, b 2e etc...

ici=pwd;
cd g:\golts\Survey\Survey291102\ADCPTethys\tronca
xx1=6;
yy1=7;
zz1='a';
eval(['load g1202tr',int2str(xx1),'t',int2str(yy1),zz1,'t.mat;'])
xyt1=xyt;
cd (ici)
[ampl1,angle1]=umod(uv(1:2:end-1,:),uv(2:2:end,:));
clear xyt zc uv
min1=min(xyt1(:,1));
max1=max(xyt1(:,1));

cd g:\golts\Survey\Survey291102\ADCPTethys\tronca
xx2=4;
yy2=9;
zz2='a';
eval(['load g1202tr',int2str(xx2),'t',int2str(yy2),zz2,'t.mat;'])
xyt2=xyt;
cd (ici)
[ampl2,angle2]=umod(uv(1:2:end-1,:),uv(2:2:end,:));
clear xyt zc uv
min2=min(xyt2(:,1));
max2=max(xyt2(:,1));

cd g:\golts\Survey\Survey291102\ADCPTethys\tronca
xx3=7;
yy3=6;
zz3='a';
eval(['load g1202tr',int2str(xx3),'t',int2str(yy3),zz3,'t.mat;'])
xyt3=xyt;
cd (ici)
[ampl3,angle3]=umod(uv(1:2:end-1,:),uv(2:2:end,:));
clear xyt uv
min3=min(xyt3(:,1));
max3=max(xyt3(:,1));

% cd g:\golts\Survey\Survey291102\ADCPTethys\tronca
% xx4=4;
% yy4=9;
% zz4='a';
% eval(['load g1202tr',int2str(xx4),'t',int2str(yy4),zz4,'t.mat;'])
% xyt4=xyt;
% cd (ici)
% [ampl4,angle4]=umod(uv(1:2:end-1,:),uv(2:2:end,:));
% clear xyt uv
% min4=min(xyt4(:,1));
% max4=max(xyt4(:,1));

mint=min([min1,min2,min3]);
maxt=max([max1,max2,max3]);

figure

ax1=axes('position',[0.25 0.71 0.45 0.2]);
pcolor(xyt1(:,1),zc,flipud(ampl1))
hold on
contour(xyt1(:,1),zc,flipud(ampl1),[0.2 0.2],'w');
shading interp, caxis([0 0.7]), c1=colorbar;
%set(c3,'ytick',[0:0.5:1.5],'yticklabel',[0 0.5 1 1.5])
set(gca,'XTickLabel',[])
set(gca,'layer','top')
ylabel(['',int2str(xx1),' to ',int2str(yy1),' ',zz1,''])
title('GOLTS Decembre 2002  intensité des courants (m.s-1)')
axis([mint maxt -200 0])
hold off

ax2=axes('position',[0.25 0.5 0.45 0.2]);
pcolor(xyt2(:,1),zc,flipud(ampl2))
hold on
contour(xyt2(:,1),zc,flipud(ampl2),[0.2 0.2],'w');
shading interp, caxis([0 0.7]), c2=colorbar;
set(gca,'XTickLabel',[])
set(gca,'layer','top')
ylabel(['',int2str(xx2),' to ',int2str(yy2),' ',zz2,''])
axis([mint maxt -200 0])
hold off

ax3=axes('position',[0.25 0.29 0.45 0.2]);
pcolor(xyt3(:,1),zc,flipud(ampl3))
hold on
contour(xyt3(:,1),zc,flipud(ampl3),[0.2 0.2],'w');
shading interp,caxis([0 0.7]), c3=colorbar;
%set(gca,'XTickLabel',[])
set(gca,'layer','top')
ylabel(['',int2str(xx3),' to ',int2str(yy3),' ',zz3,''])
axis([mint maxt -200 0])
hold off

% ax4=axes('position',[0.25 0.08 0.45 0.2]);
% pcolor(xyt4(:,1),zc,flipud(ampl4))
% hold on
% contour(xyt4(:,1),zc,flipud(ampl4),[0.2 0.2],'w');
% shading interp,caxis([0 0.7]), c4=colorbar;
% ylabel(['',int2str(xx4),' to ',int2str(yy4),' ',zz4,''])
% axis([mint maxt -200 0])
% hold off

set(gcf,'unit','centimeter',...
		'papertype','A4',...
		'PaperPosition',[-3.51 3.34 30 23])
 

%print -dpsc2 -P\\00SAMBACOM\laser4550bc