% seis1ddata.m inputs the data for seis1d
% change for various models
nl=12; % number of layers (number of primary reflections will be nl-1
h=[1.0 1.5 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 3.0];  % thicknesses of layers
v=[2.0 4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0];  % velocities of layers
%
% input other parameters
%
mtitle='"Gradient Model -- Model 1"';	% title for model plot
ftitle='Gradient model, 3 layers, 5 Hz ricker wavelet'; % title for figure
npts=400;	% number of points in seismogram
dt=0.01;	% sample interval
maxt=4.0;	% maximum time for seismogram plot
maxz=4.0;	% maximum depth for model plot
vmin=0.0;	% minimum velocity for model plot
vmax=8.0;	% 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
%
% 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
freq=5;		% frequency of wavelet
nsw=10;	% length of cosine taper on each end of wavelet
nmult=0; % 1 = include multiples; 0 = primaries only
% end of input