% seis1ddata.m inputs the data for seis1d
% change for various models
nl=5; % number of layers (number of primary reflections will be nl-1
h=[0.5 1.0 0.5 0.75 1.5];  % thicknesses of layers
%Change densities to change reflection coefficients
rho=[1.0 1.0 1.0 0.6 0.6]; % densities of layers
v=[1.5 2.0 2.5 3.0 4.0];  % velocities of layers
%
% input other parameters
%
npts=400;	% number of points in seismogram
dt=0.01;	% sample interval
maxt=4.0;	% maximum time for seismogram plot
maxz=4.5;	% maximum depth for model plot
vmin=0.0;	% minimum velocity for model plot
vmax=4.2;	% maximum velocity for model plot
nwavelet=1; % Select wavelet; nwavelet = 1, Ricker wavelet; otherwise, Signal wavelet
% Signal wavelet input (not used if nwavelet = 1).
nw1=100; % number of points in Signal wavelet
next=9; % number of extrema
ts=0.91; % length of wavelet in seconds
thalf=0.2; % if thalf > 0, the wavelet is weighted by an exponential such that at 
%         thalf s the amplitude is multipied by 0.5
%         thalf is usually chosen to be ~t or t/2
1
% Ricker wavelet input (not used if nwavelet ~= 1).
tshift=0.25;	% time shift for center of wavelet = arrival time
% The wavelet is shifted by half the wavelet length so that the travel
% times for reflections are located at the center of the wavelet.
nw=50;	% number of points in ricker wavelet
% Change the following line:
freq=8;		% frequency of wavelet
nsw=10;	% length of cosine taper on each end of wavelet
% Change the following line for multiples:
nmult=0; % 1 = include multiples; 0 = primaries only
%
mtitle='"Standard Model"';	% title for model plot
ftitle=['Standard Model, 4 layers,',int2str(freq),' Hz ricker wavelet']; % title for figure
% end of input
