%shot gathers

timelen = zeros(length(time),length(receivers)); 
for (k = 1:length(sources))
    traces = timelen;
    S = [sources(k) 0];
    figure(k);
    hold on;
    subplot(2,1,1), plot(layer(:,1),layer(:,2),'b-');
    axis([receivers(1) receivers(length(receivers)) -1 depth]);
    hold on;
    set(gca,'YDir','reverse');
    
reflections = diprefl(S,receivers,time,layer,velocity1);
for (refnum = 1:size(reflections,1) )
    wav1 = wavelet(amplitude,reflections(refnum,:),time);
    traces = traces + wav1;
end;


for j = 1:length(receivers) %seperate traces by channel - add offset
     for i = 1:size(traces,1)
         traces(i,j)=traces(i,j) + receivers(j);
     end
end


hold on;
subplot(2,1,2) , plot(traces,time);
axis([receivers(1) receivers(length(receivers)) time(1) time(length(time)) ]);
set(gca,'YDir','reverse');
end;