sioseis << eof
procs diskin diskoa filter diskob diskoc end
diskin
   fno 331 lno 331 set 3.5 7.5 ftr 25 ltr 25
   ipath nbp9702tape38.segy end 
end
diskoa
    opath data end
end
filter
   ftype 0 dbdrop 48 pass 5 80 end
end
diskob
   opath dataf end
end
diskoc
   decimf 2 opath datafr end
end
end
eof



matlab5.1 << eof1
a = (0:.4882812:250)';
rdsegy data 331 25;
x = fft(ans,1024);
b = abs(x(1:513,1));
rdsegy dataf 331 25;
x = fft(ans,1024);
c = abs(x(1:513,1));
 plot(a,b,'b',a,c,'r');
rdsegy datafr 331 25;
x = fft(ans,1024);
d = abs(x(1:513,1));
d = d * 2.;
e = (0:.2441406:125);
plot(a,b,'b',a,c,'r',e,d,'g');
xlabel ('frequency (Hz.)');
title('Decimation to 4ms ok?');        
gtext('blue = raw trace');
gtext('red = 5x80 filter');
gtext('green = filter and decimation');
print -dpsc palmer5.ps;
eof1