% plot_cor.m
ici=pwd;

%cd e:\Gatti2\golts\Survey\golts_0604\
cd e:\Gatti2\golts\Survey\golts_2002
%cd E:\Gatti2\golts\Valor\Valor040904\golts
load TGS0202skaphi.mat

% [ind1]=find(cor1<40);
% [ind2]=find(cor2<40);
% [ind3]=find(cor3<40);
% [ind4]=find(cor4<40);
% cor1(ind1)=NaN;
% cor2(ind2)=NaN;
% cor3(ind3)=NaN;
% cor4(ind4)=NaN;

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

figure
ax1=axes('position',[0.25 0.71 0.45 0.2]);
pcolor(x,zc,flipud(cor1'))
shading interp, c1=colorbar;
eval(['title([''GOLTS June 2002 - TGS0202''])'])
ylabel(['cor1'])
set(gca,'XTickLabel',[])

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

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

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

cd (ici)