% plot_vitesses.m
ici=pwd;

%cd e:\Gatti2\golts\Survey\golts_2002
cd e:\Gatti2\golts\Survey\golts_1204
load TGS0416aphi_new.mat

x=(1:1:nEnsembles)';
zc=(-244:4:-8)'; 

% conversion cm/s en m/s
ve=ve/100;
vn=vn/100;
vz=vz/100;
verr=verr/100;

figure

ax1=axes('position',[0.25 0.71 0.45 0.2]);
pcolor(x,zc,flipud(ve'))
shading interp, c1=colorbar;
eval(['title([''GOLTS December 2002 - TGS0416aphi_-new''])'])
ylabel(['ve (m.s-1)'])
set(gca,'XTickLabel',[])

ax2=axes('position',[0.25 0.5 0.45 0.2]);
pcolor(x,zc,flipud(vn'))
shading interp, c2=colorbar;
ylabel(['vn'])
set(gca,'XTickLabel',[])

ax3=axes('position',[0.25 0.29 0.45 0.2]);
pcolor(x,zc,flipud(vz'))
shading interp, c2=colorbar;
ylabel(['vz'])
set(gca,'XTickLabel',[])

ax4=axes('position',[0.25 0.08 0.45 0.2]);
pcolor(x,zc,flipud(verr'))
shading interp, c2=colorbar;
ylabel(['verr'])

figure

eastCurrentComponent=eastCurrentComponent/100;
northCurrentComponent=northCurrentComponent/100;
speedCurrentComponent=speedCurrentComponent/100;

ax1=axes('position',[0.25 0.71 0.45 0.2]);
pcolor(x,zc,flipud(eastCurrentComponent'))
shading interp, c1=colorbar;
eval(['title([''GOLTS December 2002 - TGS0416aphi_-new.mat''])'])
ylabel(['U (m.s-1)'])
set(gca,'XTickLabel',[])

ax2=axes('position',[0.25 0.5 0.45 0.2]);
pcolor(x,zc,flipud(northCurrentComponent'))
shading interp, c2=colorbar;
ylabel(['V'])
set(gca,'XTickLabel',[])

ax3=axes('position',[0.25 0.29 0.45 0.2]);
pcolor(x,zc,flipud(speedCurrentComponent'))
shading interp, c2=colorbar;
ylabel(['module'])
set(gca,'XTickLabel',[])

ax4=axes('position',[0.25 0.08 0.45 0.2]);
pcolor(x,zc,flipud(dirCurrentComponent'))
shading interp, c2=colorbar;
ylabel(['argument'])


figure
cd e:\Gatti2\golts\Survey\golts_1204
load TGS0416skaphi.mat
eastCurrentComponent=eastCurrentComponent/100;
northCurrentComponent=northCurrentComponent/100;
speedCurrentComponent=speedCurrentComponent/100;

ax1=axes('position',[0.25 0.71 0.45 0.2]);
pcolor(x,zc,flipud(eastCurrentComponent'))
shading interp, c1=colorbar;
eval(['title([''GOLTS December 2002 - TGS0416skaphi.mat''])'])
ylabel(['U (m.s-1)'])
set(gca,'XTickLabel',[])

ax2=axes('position',[0.25 0.5 0.45 0.2]);
pcolor(x,zc,flipud(northCurrentComponent'))
shading interp, c2=colorbar;
ylabel(['V'])
set(gca,'XTickLabel',[])

ax3=axes('position',[0.25 0.29 0.45 0.2]);
pcolor(x,zc,flipud(speedCurrentComponent'))
shading interp, c2=colorbar;
ylabel(['module'])
set(gca,'XTickLabel',[])

ax4=axes('position',[0.25 0.08 0.45 0.2]);
pcolor(x,zc,flipud(dirCurrentComponent'))
shading interp, c2=colorbar;
ylabel(['argument'])
cd (ici)