% plot_int.m
ici=pwd;


cd D:\golts\Survey\golts_0604\
%cd D:\golts\Valor\Valor080303\golts_aphi2\
%cd D:\golts\Survey\golts_2002
%cd E:\Gatti2\Theveny\golts_valor_mars
%cd E:\Gatti2\golts\Valor\Valor040904\golts
ses='406';
eval(['load TGS0',ses,'skaphi.mat'])

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

% intmin=60;
% [ind1]=find(int1<intmin);
% [ind2]=find(int2<intmin);
% [ind3]=find(int3<intmin);
% [ind4]=find(int4<intmin);
% int1(ind1)=NaN;
% int2(ind2)=NaN;
% int3(ind3)=NaN;
% int4(ind4)=NaN;

figure
ax1=axes('position',[0.25 0.71 0.45 0.2]);
pcolor(x,zc,flipud(int1'))
 %cmin=min(min(int1))
 %cmax=max(max(int1))
cmin=35;
cmax=100;
caxis([cmin cmax])
shading interp, c1=colorbar;
title(['GOLTS Juin 2004 - TGS0',ses,''])
ylabel(['int1'])
set(gca,'XTickLabel',[])

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

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

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

nomfigjpg=['int_',ses,'.jpg'];
nomfig=['int_',ses,'.fig'];
saveas(gcf,nomfigjpg,'jpg')
saveas(gcf,nomfig,'fig')

figure
plot(x,int1(:,37))

cd (ici)