sioseis << eof
procs diskin filter diskoa end
diskin
   fno 1 lno 1 ftr 1 ltr 1 set .1 .2
   ipath source.sig.test.sgy end
end
filter
   ftype 0 pass 20 125 end
end
diskoa
  opath data end
end
end
eof

matlab << eof
a = zeros(401,3);
a(:,1) = (.1:.00025:.2)';
rdsegy data 1 1;
a(:,2) = ans;
plot(a(:,1),a(:,2),a(:,1),a(:,3));
xlabel('time');
ylabel('counts');
title('Watergun with 20x125 filter');
print water6.ps
rdsegy data 1 1;
a = fft(ans,2048);
b = zeros(1025,2);
b(:,1) = (0:1.953125:2000)';
b(:,2) = abs(a(1:1025));
plot(b(:,1),b(:,2));
print water7.ps
plot(b(1:257,1),b(1:257,2));
print water8.ps             
plot(b(1:65,1),b(1:65,2));  
print water9.ps
quit
eof
convert water6.ps water6.gif
convert water7.ps water7.gif
convert water8.ps water8.gif
convert water9.ps water9.gif