% plottrlat.m
% trace en contours la coupe verticale des courants mesurés par l'ADCP de coque
% pour un transect choisi
% de haut en bas: UU (haut), VV, amplitude, angle [0,360] (bas) avec echelles independantes
% x = lat
% cree par AP le 10/09/02
% modifie par JG le 16/01/04

% xx numero de station depart
% yy numero de station fin
% zz lettre pour indiquer si c'est le a: 1er passage, b 2e etc...

clear
xx=1;
yy=5;
zz='a';

ici=pwd;
cd e:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
eval(['load g1203tr',int2str(xx),'t',int2str(yy),zz,'t.mat;'])

ax1=axes('position',[0.25 0.71 0.45 0.2]);
contourf(xyt(:,2),zc,flipud(uv(1:2:end-1,:)),[-1:0.2:1]);
c1=colorbar;
%set(c1,'ytick',[-1:0.5:2],'yticklabel',[-1 -0.5 0 0.5 1 1.5 2])
set(gca,'XTickLabel',[])
eval(['title([''GOLTS December 2003 - Transect ',int2str(xx),' to ',int2str(yy),' ',zz,'.''])'])
text(min(xyt(:,2))+0.01,-200,'U (Est>0)');

ax2=axes('position',[0.25 0.5 0.45 0.2]);
contourf(xyt(:,2),zc,flipud(uv(2:2:end,:)),[-1:0.2:1])
c2=colorbar;
%set(c2,'ytick',[0:0.5:3],'yticklabel',[0 0.5 1 1.5 2 2.5 3])
set(gca,'XTickLabel',[])
text(min(xyt(:,2))+0.01,-200,'V (Nord>0)');

cd (ici);
[ampl,angle]=umod(uv(1:2:end-1,:),uv(2:2:end,:));

ax3=axes('position',[0.25 0.29 0.45 0.2]);
contourf(xyt(:,2),zc,flipud(ampl),[0:0.2:1]);
c3=colorbar;
%set(c3,'ytick',[0:0.5:3],'yticklabel',[0 0.5 1 1.5 2 2.5 3])
set(gca,'XTickLabel',[])
text(min(xyt(:,2))+0.01,-200,'amplitude');

ax4=axes('position',[0.25 0.08 0.45 0.2]);
contourf(xyt(:,2),zc,flipud(angle),[0:90:360])
c4=colorbar;
set(c4,'ytick',[0:90:360],'yticklabel',[0 90 180 270 360])
xlabel(['latitude']);
text(min(xyt(:,2))+0.01,-200,'angle');

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