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 250 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 20x250 zero phase filter');
print water10.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 water11.ps
plot(b(1:257,1),b(1:257,2));
print water12.ps             
plot(b(1:65,1),b(1:65,2));  
print water13.ps
quit
eof
convert water10.ps water10.gif
convert water11.ps water11.gif
convert water12.ps water12.gif
convert water13.ps water13.gif