stack.WtoE


#!/bin/csh -f
if( $#argv != 1 ) then
    echo "Usage: stack line-number "
    exit 1
endif
set LINENO = $1

mkdir /export/home/EW0008/Realtimestack/Line$LINENO
chmod a+rw /export/home/EW0008/Realtimestack/Line$LINENO

/home/henkart/bin/sioseis-fast << eof

procs segdin geom prout header gather diskoa tredit nmo stack
      diskoe filter gains plot end

segdin
   ffilen 99999 
   fcset 1 lcset 1
   ftr 1 ltr 324
   stime 3 secs 4
   offline yes  # eject after the rewind after EOT
   iunit 4 end
end

tredit
  limits -1000 0 kill outside lhdr 10    # kill ranges outside 0 - 1000
  fno 0 lno 999999 end
end

prout
  fno 0 lno 99999 ftr 324 ltr 324 noinc 10 end
end

weight
  fno 0 lno 999999 twp 187 0 247 0 295 0 307 0 end
end

geom
  type 2 # increment the shot location based on the shot number
  fs 1 ls 999999  # all shot have the same parameters (preset)
  ggx -12.5         # Used to extrapolate gxp!
  dfls 37.5 dbrps 6.25 smear 6.25
  gxp 324 -54.25 ggx -12.5 dfls 37.5 dbrps 6.25 smear 6.25
  bgp 1 6 2 7 3 18 4 30 5 42 6 66 7 90 8 102 9 126 10 138
       11 162 12 174 13 186 14 198 15 210 16 222 17 234 18 246
      19 258 20 270 21 282 22 294 23 306 24 318
  rpadd 1000 end
end

mute
   fno 1 lno 999999
   addwb yes xtp 200 -.1 1500 -.1 3000 1 6200 2 end
end

plot
   nsecs 4  dir ltr
   scalar 2.E-05
   tlines 0.1 0.5 1 nibs 7224
   def 0.150 clip .2 trpin 75 wiggle -50 vscale 5
   ann gmtint anninc 5
   opath /export/home/EW0008/Realtimestack/Line$LINENO/siopltfil.Line$LINENO
  end
end
  
diskoa
# write every 50th cmp gather into a file for velocity analysis
   opath /export/home/EW0008/Realtimestack/Line$LINENO/line$LINENO.cmp.inc50.segy
   fno 1 lno 9999999 noinc 50 end
end

diskob 
# write every 20th shot to a "circular" file 
    fno 1 lno 999999 noinc 20 rewind 1 
    opath /export/home/EW0008/Realtimestack/latest.shot.segy end 
end

diskoe   # Write out disk file
  opath /export/home/EW0008/Realtimestack/Line$LINENO/stack.line$LINENO end
end

header
   fno 0 lno 9999999 ftr 1 ltr 9999 
   r50 r54 / 750.  # convert water depth to water time
   bhdr 28 1
   end
end

gather
   maxtrs 90 maxrps 500 end
#   maxtrs 60 maxrps 500 end
end

filter
#   ftype 0 dbdrop 48 pass 20 120 end
   ftype 0 dbdrop 48 pass 40 220 end
end

gains
   type 3 alpha 4.2 end
end

agc 
  winlen 1 pctagc 100 end 
end 

nmo
 addwb yes
 stretc 0.7
vtp
1515 0
1502 0.109
1520 0.145
1520 0.188
1520 0.243
1520 0.352
1520 0.407
1540 0.491
1541 0.57
1540 0.726
1559 0.808
1579 0.912
1626 1.927
end

end

end
eof

Go to the list of seismic processes.      Go to SIOSEIS introduction.