% close all
% clear all
% 
% load ./Profile_4/vfinal_sch_8.mat
figure;imagesc(vfinal_sch_8,[1200 2000]);
hold on;
contour(vfinal_sch_8,30,'w');
colorbar
xlabel('X')
ylabel('Depth')
title('Traveltime Velocity Tomogram');
axis equal

[gvx,gvz]=gradient(vfinal_sch_8,1,1);
%gv=gvz(find(abs(gvz)>0.0001))=1;
gvz(find(abs(gvz)>0.0001))=1;
gvz(find(abs(gvz)<0.0001))=0;
gv=vfinal_sch_8.*gvz;
figure(2);imagesc(gv,[1200 2000]);hold on;contour(gv,60,'w');
colorbar
xlabel('X')
ylabel('Depth')
title('Traveltime Velocity Tomogram');
% axis equal

[gvx,gvz]=gradient(vfinal_sch_8,1,1);
gvz(find(abs(gvz)>0))=1;
gvz(find(abs(gvz)==0))=0;
gv=vfinal_sch_8.*gvz;
figure;imagesc(gv,[1200 2000]);hold on;contour(gv,30,'w');
colorbar
xlabel('X')
ylabel('Depth')
title('Traveltime Velocity Tomogram');
axis equal