%set up the model (global variables)

global  time dt
global  line_length depth
global  sources receivers spacing
global  velocity1 layer 
global  tref wref

load wavelet.mat
% load model1.mat

dt = 2;          % time sampling interval
time = 0:dt:800;  % time sampling
amplitude = 50;

sources = [250 500 750];     % positions of three sources
receivers = 0:25:1000;    % positions of 41 receivers
depth = 500;
velocity1 = 2;

line_length = receivers(end);                 % total length of seismic line
spacing = line_length/(length(receivers)-1);  % receiver spacing

% draw and pick a reflecting layer
 layer = picklayer(line_length,depth);

% print out the picked points
layer
%S - source positions (can be an array)
%   R - array of receiver positions
section1 = reflection(sources(1),receivers,layer,velocity1);
section2 = reflection(sources(2),receivers,layer,velocity1);
section3 = reflection(sources(3),receivers,layer,velocity1);
plot_section(receivers,section1,'Distance(m)','b-')
plot_section(receivers,section2,'Distance(m)','r-')
plot_section(receivers,section3,'Distance(m)','g-')
plot(sources,[0 0 0],'ko')
% save the model and resulting section

save model1.mat layer*   section*