EW0210 SIOSEIS constant velocity stack script




sioseis << eof
procs diskin header velan header2 stack filter agc plot end
diskin
   fno 1 lno 50 allno no
   ipath /export/home2/scratch/ew0210/seismic/line10/data/dkgather
   end
end
filter
    pass 5 60 ftype 0 dbdrop 48 end
end
agc
   winlen .5 center .1 end
end
header
#  save the rp number because velan renumbers the rp to be the velocity 
number
   fno 0 lno 999999 l5 = l6 end
end
header2
#   restore the original rp number even though there will be lots with t
he same number
   fno 0 lno 999999 l6 = l5 end
end
velan
   nrp 50 vels 1450 25  1600 50 2200 type cvel end
end
plot
   stime 3
   nibs 200 vscale 1  nsecs 8
   trpin 100 fspace 50 nspace 3 spacei 50
   def .01 clip .005
   ann header ihdr 46 ftag 1 taginc 50
   srpath sunfil.ras
   opath siofil hpath headers
   end
end
end
eof