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