clear all
close all
wl_SS=[407,439,485,507,527,594,651,715,878];
list=dir('*.mat');
gamma_all=ones(59,26)*NaN;
bb527_all=ones(59,26)*NaN;
k=1
for i=1:length(list)
    load(list(i).name)
    [m n]=size(datT);
    if n>6
        x(k)=cast_time;
        gamma_all(uint8(floor(prrT)/2)+1,k)=gamma;
        bb527_all(uint8(floor(prrT)/2)+1,k)=datT(:,11);
        k=k+1;
    end
end
I=find(gamma_all(:)>2);
gamma_all(I)=NaN;

I=find(isnan(gamma_all)==0);
y=0.5:2:116.5;
xx=x'*ones(1,59);
xx=xx';
yy=y'*ones(1,26);
scatter(xx(I),-yy(I),40,gamma_all(I),'filled')
colorbar
axis tight;
datetick
xlabel('date')
ylabel('depth')

I=find(bb527_all(:)>4e-3);
bb527_all(I)=NaN;

figure
I=find(isnan(bb527_all)==0);
y=0.5:2:116.5;
xx=x'*ones(1,59);
xx=xx';
yy=y'*ones(1,26);
scatter(xx(I),-yy(I),40,bb527_all(I),'filled')
colorbar
axis tight;
datetick
xlabel('date')
ylabel('depth')
axis tight;

figure
plot(bb527_all,-y,'.')