sioseis << eof
procs diskin  diskoa end
diskin
   fno 1 lno 1 ftr 1 ltr 1 set .1 .2
   ipath source.sig.test.sgy end
end
diskoa
  opath data end
end
end
eof

matlab << eof
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 water3.ps
plot(b(1:257,1),b(1:257,2));
print water4.ps             
plot(b(1:65,1),b(1:65,2));  
print water5.ps
quit
eof
convert water3.ps water3.gif
convert water4.ps water4.gif
convert water5.ps water5.gif