rm sunfil.ras
sioseis << eof
procs diskin header geom weight gather nmo stack filter plot end
diskin
   ftr 1 ltr 30  ! only use first 30 traces
   ipath /home/phenk/data/Tj37/raw_segy.dir/day203-1413z.segy end
   ftr 1 ltr 30  ! only use first 30 traces
   ipath /home/phenk/data/Tj37/raw_segy.dir/day203-1440z.segy end
   ftr 1 ltr 30  ! only use first 30 traces
   ipath /home/phenk/data/Tj37/raw_segy.dir/day203-1505z.segy end
end
header
    fno 0 lno 999999 ftr 0 ltr 99999
    i79 = i79 + 20    ! correct the year
    end
end
geom
   rpadd 10000
   type 6 navfil sio_line34
   ntrcs 48 offset 45
    gxp 1 -12.5 ggx 12.5 dbrps 6.25 end
end
gather   ! must be given so default parameters can be set
   frp 9952 ! 10000-48 see FRP and NTRCS in GATHER doc
   end
end
weight
    fno 0 lno 999999
    twp 15 0 16 0 17 0 18 0 20 0 25 2000 29 0 30 0 end
end
nmo  
    fno 9969 lno 9969  vtp  1500  0.173  1550  0.243  1600  0.417 end
    fno 10189 lno 10189  vtp  1500  0.162  1550  0.268  1600  0.360 end
    fno 10392 lno 10392  vtp  1500  0.229  1700  0.394  1750  0.527 end
    fno 10592 lno 10592  vtp  1500  0.256  1600  0.451  1850  0.588 end
    fno 10784 lno 10784  vtp  1500  0.268  1600  0.360  1700  0.521  1800  0.631 end
    fno 10945 lno 10945  vtp  1500  0.259  1550  0.326  1700  0.455 end
    fno 11133 lno 11133  vtp  1500  0.256  1600  0.394 end
    fno 11325 lno 11325  vtp  1500  0.261  1600  0.385 end
    fno 11496 lno 11496  vtp  1500  0.268  1600  0.392  1950  0.489 end
    fno 11655 lno 11655  vtp  1500  0.261  1550  0.331  1750  0.437 end
end
prout
 fno 1 lno 99999 noinc 100 ftr 1 ltr 1 end
end
filter
   ftype 0 pass 30 120 dbdrop 48 end
end
plot
   ann gmtint anninc 1 ann2 rpno
   trpin 300
   nibs 2858 srpath sunfil.ras taginc 2
   def .01 clip .005 wiggle 0  nsecs 1 vscale 2
   end
end
end
eof
display -rotate 90 sunfil.ras &
convert -rotate 90 sunfil.ras short_range_stk.png