
 
function plt(seis,dt,dx,sca,ntskip)
%To generate a shaded plot of your glorified seismic data.
%In the likely case of problems, contact KIM BAK VISUALIZATION INC.
%Otherwise, good luck, and shade them traces!
echo on
%
%dt=.144;   %time discretization
%t=127.;    %total time
%nt=887;    %# of time steps
%dx=.8;     %spatial discretization
%len=98.4;  %profile length
%ntra=124;  %# of traces
%sca=4.;    %scaling factor
%ntskip=1;  %# of timesteps to be skipped
%
hold on;
%matrix of multiplexed seismic data (nt by ntra)
[nt,ntra]=size(seis);t=nt*dt;len=ntra*dx;
dt=dt*ntskip;
nt=nt/ntskip;
nt=floor(nt);
cl=[ 1. 1. 1.];
c=(0.:dt:t+t/10.);
x=[0 nt*dt 0 len];
for i=1:ntra
nn=-ntskip+1;
for k=1:nt
nn=nn+ntskip;
u(k,i)=seis(nn,i)/sca+i*dx;
ux(k,i)=c(k);
uz(k,i)=u(k,i);
if (u(k,i)>i*dx),
uz(k,i)=i*dx;
end;
end;
end;
for i=1:ntra
for k=nt+1:2*nt-1
ux(k,i)=c(2*nt-k);
uz(k,i)=u(2*nt-k,i);
end;
end;
for i=1:ntra
fill(ux(1:2*nt-1,i),uz(1:2*nt-1,i),cl);
end;
plot(c(1:nt),u(1:nt,1:ntra),'-');
xlabel('TIME (SECONDS)');
ylabel('DISTANCE (KM)');
axis(x);
echo off

