examine (plot) a few shots
rm sunfil.ras
sioseis << eof
procs segddin filter agc plot end
segddin
#   ipath R000186_1316606048.RAW end   # 1233
#   ipath R000187_1316606061.RAW end   # 1234
#   ipath R000188_1316606074.RAW end   #1235
#   ipath R000300_1316607520.RAW end   # 1347
#   ipath R000400_1316608768.RAW end  # 1447
#   ipath R000500_1316610022.RAW stime 3 end # 1547
   ipath R001000_1316616311.RAW stime 1 secs 4 end # 2047
#   ipath R002000_1316628777.RAW stime 2 secs 3 end # 3047
#   ipath R003000_1316641454.RAW stime 2 secs 3 end # 4063
#   ipath R004000_1316653879.RAW stime 2 secs 3 end # 5063
#   ipath R005000_1316666400.RAW stime 2 secs 3 end # 6063
#   ipath R006000_1316678915.RAW stime 1 secs 4 end  #  7065
#   ipath R007000_1316691381.RAW stime 2 secs 3 end  # 8065
#   ipath R007900_1316702597.RAW stime 1 end # 8965
end
filter
   ftype 0 pass 10 90 end
end
agc
  winlen .250 end
end
plot
   trpin 20 ann sh&tr taginc 2
   vscale 2 def .01
   srpath sunfil.ras
   end
end
end
eof
display -rotate 90 sunfil.ras &


gather (sort by rp) using UKOOA nav and killing bad traces
sioseis << eof
procs segddin weight geom prout gather diskoa end
segddin
   secs 8
   listpath list-new end
end
weight
    fno 0 lno 9999999 twp 30 0 257 0 293 0 301 0 343 0 349 0 377 0 396 0 437 0 end
end
geom
  dbrps 6.25
  type 19 navfil /home/Paul/data/mgl1112/18--MGL1112MCS--04.p190 end
end
gather
   maxtrs 150 end
end
diskoa
   opath gathers.segy end
end
prout
   fno 0 lno 999999 ftr 1 ltr 1 noinc 100 end
end
end
eof


gather (sort by rp) using LDEO realtime GPS and killing bad traces
sioseis << eof
procs segddin weight geom prout gather diskoa end
segddin
   secs 8 
   listpath list-new end
end
weight 
      twp 30 0 169 0 186 0 192 0 294 0 301 0 343 0 349 0 377 0 396 0 
         408 1.25 409 1.25 427 1.5 428 1.5 429 1.5 430 1.5 437 0 end 
end 
geom
  dbrps 6.25 type 9 ggx -12.5 gxp 468 -220 end
end
gather
   end
end
diskoa
   opath gathers-gps.segy end
end
prout
   fno 0 lno 999999 ftr 1 ltr 1 noinc 100 end
end
end
eof


vpick
#! /bin/bash

#  vpick - velocity picker using sioseis and matlab
#  Usage: vpick filename start_rpnum [ rpnum_inc end_rpnum]
#
#  you may want to add the diskin  SET parameter (start and end time)
#  I collect all the cmp gathers into a separate file to keep the disk I/O down.
#  Paul Henkart, May 2006

declare -i start_rpnum=0 rpnum_inc=0 end_rpnum=0 
declare -i rpnum

num_args="$#"
if [ $num_args = "0" ]; then
   echo "Usage: vpick filename start_rpnum [ rpnum_inc end_rpnum]"
   exit
fi
file=$1
if [ ! -f $file ]; then
   echo "***  ERROR  ***  Could not open input file: $file"
   exit
fi
start_rpnum=$2
rpnum_inc=$3
if [ $rpnum_inc == 0 ]; then
     rpnum_inc=1
fi
end_rpnum=$4
if [ $end_rpnum == 0 ]; then
     end_rpnum=$start_rpnum
fi


matfile="vpick.mat"
vpickfile="$file.vpick"

for (( rpnum=start_rpnum; rpnum <= end_rpnum; rpnum=rpnum+rpnum_inc ))
do

   sioseis << eof
   noecho procs diskin filter agc plot velan end
   diskin
      fno $rpnum lno $rpnum allno no
      ipath $file end END
   filter ftype 0 pass 40 150 dbdrop 48 end end
   agc
      winlen .25 center .01  end END
   prout
       fno 0 lno 99999 ftr 0 ltr 99999 end
   end
   velan
       vels 1450 10 2200 nrp 1 type spec winlen .048 opath $matfile
       end
   end
plot
   def .02
!  scalar 3.E-01
  clip .03
  srpath sunfil.ras
  vscale 1 trpin 15
  ann range ann2 sh&tr taginc 1
  end
  end
   end
#  eof must be in column 1 - bash is stupid
eof
   display -rotate 90 sunfil.ras &



/Applications/MATLAB_R2012a.app/bin/matlab -nodesktop << eof1
fid = fopen('$vpickfile','a');
load $matfile;
n = size(vel);nt = n(1);
nv = n(2);
rpno = vel(1,1);
st = vel(2,1);
dt = vel(3,1);
sv = vel(4,1);
dv = vel(5,1);
vel(1,1) = 0;vel(2,1) = 0;vel(3,1) = 0;vel(4,1) = 0;vel(5,1) = 0;
x = sv:dv:sv+nv*dv-dv;
y = -st:-dt:-(st+nt*dt-dt);
scrsz = get(0,'ScreenSize');
% figure('Position',[left, bottom, width, height])
figure('Position',[1 scrsz(4) scrsz(3)*.75 scrsz(4)])
contour(x,y,vel,20)
title 'rp $rpnum';
xx=[];
yy=[];
n = 0;
but = 1;
while but == 1
   [xi,yi,button] = ginput(1);
   n = n + 1;
   xx(n,1) = xi;
   yy(n,1) = yi;
   if button > 1
      but = 0;
   end
end
pause
n = size(xx);
fprintf(fid,'  fno %.0f vtp ',$rpnum);
for i = 1:n(1)
   fprintf(fid,'%.0f %.3f ',xx(i),-yy(i));
end
fprintf(fid,'end\n');
st = fclose(fid);
quit
eof1

done
echo " "

Stack
#! /bin/bash

   sioseis << eof
   procs diskin weight wbt nmo mute stack diskoa prout end
   diskin
      fno 5200 lno 22200
      ipath /home/Paul/data/mgl1112/segy/gathers-gps.segy end END
   mute  addwb yes xtp 0 -.1 2000 -.1 3000 1 4000 2 5000 3 6000 4 end end
   weight
      fno 0 lno 9999999 xwp 1770 0 2858 0 3720 0 end end
nmo
####  first good
  fno 5200 vtp 1465 2.083 1480 2.306 1520 2.6 1530 2.7 1590 2.8 1640 2.947 1730 3.1 1900 4 2202 8.005 end
  fno 5700 vtp 1465 1.818 1480 2.181 1580 2.557 1840 2.905 1947 3.072 2202 8.005 end
  fno 6200 vtp 1465 1.386 1490 1.72 1562 2.000 1640 2.278 1779 2.613 2199 8 end
  fno 6700 vtp 1465 1.3 1490 1.5  1540 1.7 1660 2.139 1855 2.585 2199 3.685 end
  fno 7200 vtp 1465 1.3  1490 1.5 1540 1.707 1590 1.902 1670 2.125 1765 2.306 2198 3.616 end
  fno 7700 vtp 1465 1.3  1490 1.65 1540 1.7 1560 1.818 1670 2.167   1765 2.5   2100 3 2199 8 end
  fno 8200 vtp 1465 1.317 1500 1.7 1580 1.958 1680 2.306 1700 2.376 2018 2.975 2201 7 end
  fno 8700 vtp 1465 1.4   1510 1.7 1630 2.055 1680 2.195 1860 2.640 2000 2.947 2201 7 end
  fno 9200 vtp 1465 1.3  1490 1.6 1530 1.7 1640 2.000 1830 2.557 1900 3.0 2199 7 end
  fno 9700 vtp 1470 1.4  1500 1.7 1560 1.8 1590 1.9 1720 2.222 1920 2.891 2178 7 end
  fno 10200 vtp 1470 1.219 1480 1.4 1500 1.7           1620 1.9   1700 2.125 1840 2.5 1900 3 2199 7 end
  fno 10700 vtp 1465 1.3 1480 1.5 1500 1.7    1640 2.0   1690 2.260          1850 3 2201 7 end
  fno 11200 vtp 1465 1.5 1480 1.8           1560 2.05             1710 2.4 1850 2.9 2199 7 end
  fno 11700 vtp 1480 1.95  1550 2.2           1570 2.6   1700 3.           2184 7 end
  fno 12200 vtp 1465 2.2 1510 2.571 1590 2.724 1640 4 2199 7 end
  fno 12700 vtp 1465 2.35 1490 2.65 1525 2.808 1600 3.407 1650 4 2202 7 end
  fno 13200 vtp 1470 2.45  1515 2.860 1610 3.3 1850 3.825  1880 4.5 2197 7 end
  fno 13700 vtp 1470 2.5 1480 2.682 1543 3.128 1680 3.5 1830 3.853 1900 4.8 2201 7 end
  fno 14200 vtp 1470 2.53 1480 2.752 1535 3.017 1560 3.128 1700 3.600 1830 4.00 2070 4.480 2199 8 end
  fno 14700 vtp 1470 2.55 1510 2.905 1550 3.184 1770 3.853 1850 4.243 1950 5 2202 8 end
  fno 15200 vtp 1470 2.515 1481 2.738 1542 3.142 1650 3.532 1810 4.020 1950 4.201 2000 5 2200 8 end
  fno 15700 vtp 1470 2.55 1500 2.933 1600 3.449 1740 3.881 1800 4.05 1900 5 2199 8 end
  fno 16200 vtp 1470 2.55 1480 2.836 1490 2.9 1600 3.407 1700 3.769 1850 5 2201 8 end
  fno 16700 vtp 1470 2.55 1475 2.738 1570 3.2 1650 3.4 1590 3.6 1554 3.713 1511 3.894 1700 4.8 2201 7 end
  fno 17200 vtp 1470 2.45 1496 2.640 1550 3.100 1580 3.7 1680 4.0 1700 5 2201 7. end
  fno 17700 vtp 1470 2.7 1500 3.058 1539 3.212 1557 3.337 1650 3.560 1600 5.2 2199 8 end
  fno 18200 vtp 1470 2.7 1500 3.0 1590 3.5 1670 3.825 1870 4.299 2113 4.591 2201 8 end
  fno 18700 vtp 1470 2.7 1510 3.128 1531 3.254 1680 3.894 1850 4.368 1900 5.2 2202 8 end
  fno 19200 vtp 1470 2.7 1480 2.90 1500 3.058 1580 3.463 1700 3.881 1760 4.062 1780 5 2201 8 end
  fno 19700 vtp 1470 2.7 1485 2.975 1530 3.340 1590 3.532 1670 3.783 1800 5. 2198 8 end
  fno 20200 vtp 1470 2.599 1480 2.905 1520 3.212 1590 3.5 1660 3.769 1800 5. 2198 8 end
  fno 20700 vtp 1470 2.696 1480 2.989            1600 3.463 1640 3.6 1740 5 2201 8 end
  fno 21200 vtp 1470 2.780 1480 3.031 1560 3.449 1642 3.825 1690 4.131 1760 5.5 1800 6.1 2199 8 end
  fno 21700 vtp 1470 2.7 1520 3.212 1550 3.4 1610 3.6 1630 3.713 1700 4.5 1800 5.5 2199 8 end
  fno 22200 vtp 1465 2.55 1480 2.919  1560 3.2 1580 3.602 1600 3.727 1750 5 2201 8 end
#   last checked
end
wbt
5550 1.85
5700 1.8
5850 1.65
5975 1.6
6150 1.4
6250 1.35
6450 1.25
6700 1.2
7200 1.2
7800 1.25
8050 1.3
8800 1.3
9300 1.2
9650 1.2
9700 1.3
9900 1.25
10000 1.2
10200 1.2
10250 1.25
10700 1.25
11200 1.45
11450 1.55
11700 1.9
11850 2.15
12400 2.15
12850 2.35
14300 2.45
15000 2.5
16200 2.5
16400 2.48
16850 2.48
17040 2.25
17150 2.25
17450 2.6
17800 2.65
20000 2.65
20100 2.6
20500 2.6
20700 2.65
20900 2.75
21300 2.75
21550 2.65
21950 2.6
22100 2.5
22200 2.5
end
prout
   fno 0 lno 99999 ftr 1 ltr 1 noinc 100 end
   end
diskoa
   opath stack.segy end
end
end
eof

FDMIGR
#! /bin/bash

   sioseis << eof
procs diskin fdmigr diskoa prout end
   diskin
   set 1 7
      random 0
#      fno 5200 lno 14200
     fno 13200 lno 22200
      ipath stack.segy end END
fdmigr
  dx 6.25 bpad 20 epad 20 tsteps 1 .1 3 .2 nx 10000
####  first good
  fno 5200 vtp 1465 2.083 1480 2.306 1520 2.6 1530 2.7 1590 2.8 1640 2.947 1730 3.1 1900 4 2202 8.005 end
  fno 5700 vtp 1465 1.818 1480 2.181 1580 2.557 1840 2.905 1947 3.072 2202 8.005 end
  fno 6200 vtp 1465 1.386 1490 1.72 1562 2.000 1640 2.278 1779 2.613 2199 8 end
  fno 6700 vtp 1465 1.3 1490 1.5  1540 1.7 1660 2.139 1855 2.585 2199 3.685 end
  fno 7200 vtp 1465 1.3  1490 1.5 1540 1.707 1590 1.902 1670 2.125 1765 2.306 2198 3.616 end
  fno 7700 vtp 1465 1.3  1490 1.65 1540 1.7 1560 1.818 1670 2.167   1765 2.5   2100 3 2199 8 end
  fno 8200 vtp 1465 1.317 1500 1.7 1580 1.958 1680 2.306 1700 2.376 2018 2.975 2201 7 end
  fno 8700 vtp 1465 1.4   1510 1.7 1630 2.055 1680 2.195 1860 2.640 2000 2.947 2201 7 end
  fno 9200 vtp 1465 1.3  1490 1.6 1530 1.7 1640 2.000 1830 2.557 1900 3.0 2199 7 end
  fno 9700 vtp 1470 1.4  1500 1.7 1560 1.8 1590 1.9 1720 2.222 1920 2.891 2178 7 end
  fno 10200 vtp 1470 1.219 1480 1.4 1500 1.7           1620 1.9   1700 2.125 1840 2.5 1900 3 2199 7 end
  fno 10700 vtp 1465 1.3 1480 1.5 1500 1.7    1640 2.0   1690 2.260          1850 3 2201 7 end
  fno 11200 vtp 1465 1.5 1480 1.8           1560 2.05             1710 2.4 1850 2.9 2199 7 end
  fno 11700 vtp 1480 1.95  1550 2.2           1570 2.6   1700 3.           2184 7 end
  fno 12200 vtp 1465 2.2 1510 2.571 1590 2.724 1640 4 2199 7 end
  fno 12700 vtp 1465 2.35 1490 2.65 1525 2.808 1600 3.407 1650 4 2202 7 end
  fno 13200 vtp 1470 2.45  1515 2.860 1610 3.3 1850 3.825  1880 4.5 2197 7 end
  fno 13700 vtp 1470 2.5 1480 2.682 1543 3.128 1680 3.5 1830 3.853 1900 4.8 2201 7 end
  fno 14200 vtp 1470 2.53 1480 2.752 1535 3.017 1560 3.128 1700 3.600 1830 4.00 2070 4.480 2199 8 end
  fno 14700 vtp 1470 2.55 1510 2.905 1550 3.184 1770 3.853 1850 4.243 1950 5 2202 8 end
  fno 15200 vtp 1470 2.515 1481 2.738 1542 3.142 1650 3.532 1810 4.020 1950 4.201 2000 5 2200 8 end
  fno 15700 vtp 1470 2.55 1500 2.933 1600 3.449 1740 3.881 1800 4.05 1900 5 2199 8 end
  fno 16200 vtp 1470 2.55 1480 2.836 1490 2.9 1600 3.407 1700 3.769 1850 5 2201 8 end
  fno 16700 vtp 1470 2.55 1475 2.738 1570 3.2 1650 3.4 1590 3.6 1554 3.713 1511 3.894 1700 4.8 2201 7 end
  fno 17200 vtp 1470 2.45 1496 2.640 1550 3.100 1580 3.7 1680 4.0 1700 5 2201 7. end
  fno 17700 vtp 1470 2.7 1500 3.058 1539 3.212 1557 3.337 1650 3.560 1600 5.2 2199 8 end
  fno 18200 vtp 1470 2.7 1500 3.0 1590 3.5 1670 3.825 1870 4.299 2113 4.591 2201 8 end
  fno 18700 vtp 1470 2.7 1510 3.128 1531 3.254 1680 3.894 1850 4.368 1900 5.2 2202 8 end
  fno 19200 vtp 1470 2.7 1480 2.90 1500 3.058 1580 3.463 1700 3.881 1760 4.062 1780 5 2201 8 end
  fno 19700 vtp 1470 2.7 1485 2.975 1530 3.340 1590 3.532 1670 3.783 1800 5. 2198 8 end
  fno 20200 vtp 1470 2.599 1480 2.905 1520 3.212 1590 3.5 1660 3.769 1800 5. 2198 8 end
  fno 20700 vtp 1470 2.696 1480 2.989            1600 3.463 1640 3.6 1740 5 2201 8 end
  fno 21200 vtp 1470 2.780 1480 3.031 1560 3.449 1642 3.825 1690 4.131 1760 5.5 1800 6.1 2199 8 end
  fno 21700 vtp 1470 2.7 1520 3.212 1550 3.4 1610 3.6 1630 3.713 1700 4.5 1800 5.5 2199 8 end
  fno 22200 vtp 1465 2.55 1480 2.919  1560 3.2 1580 3.602 1600 3.727 1750 5 2201 8 end
#   last checked
end
prout
   fno 0 lno 99999 ftr 1 ltr 1 noinc 100 end
   end
diskoa
   opath fdmigr2.segy end
end
end
eof

filter gains plot
sioseis << eof
procs diskin filter gains plot end
diskin
#   ntodo 1000
   ipath fdmigr2.segy end
end
filter
   ftype 0 pass 20 120 dbdrop 48 end
end
gains
   subwb yes type 5 alpha .6 end
end
plot
   scalar 5.5E-04
   stime 1 nsecs 6
   srpath sunfil.ras
   trpin 200 wiggle 0 vscale 1
   def .02 clip .015
    ann gmtint anninc 5 ann2 rpno
   end
end
end
eof
display -rotate 90 sunfil.ras &