Go to the list of seismic processes.      Go to SIOSEIS introduction.
               A later stack script
Changes:
1) stacking 120 of 480 traces (2km of 6km streamer)
2) A new velocity function with a low velocity zone.
3) New geom option type 9 - binning by realtime lat/long
4) New gains parameter tadd used with subwb yes


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

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

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

procs segdin geom prout gather nmo mute stack
      diskoe filter gains plot end

segdin
   ffilen 99999 
   fcset 1 lcset 1
   ftr 361 ltr 480 ! drop long ranges (1-360 ~= 2km to 4km)
   stime 3 secs 4
   offline yes  # eject after the rewind after EOT
   iunit 4 end
end

tredit
  limits -2000 0 kill outside lhdr 10  # kill 2km to 6km
  fno 0 lno 999999 end
end

prout
  fno 0 lno 999999 ftr 480 ltr 480 noinc 10 end
end

weight
   fno 0 lno 9999999 twp 51 0 end 
end

geom
  type 9  # The new dfls from the realtime nav block
#  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 $DX dbrps 6.25 smear 6.25
  dfls 37.5 dbrps 6.25 smear 6.25
  gxp 480 -184.25 ggx -12.5 dfls 37.5 dbrps 6.25 smear 6.25
  bgp 1 6 2 18 3 30 4 42 5 66 6 90 7 114 8 138 9 162 10 186 11 210
     12 234 13 258 14 282 15 306 16 330 17 354 18 378 19 402 20 426
     21 438 22 459 23 462 24 474
  rpadd 1000 end
end

mute
   fno 1 lno 999999
     addwb yes xtp 2000 -.03 2500 .6 6000 .7 end
end

plot
   nsecs 4
   scalar 2.E-04
   tlines 0.1 0.5 1 nibs 7224
   def 0.150 clip .1 trpin 75 wiggle -40 vscale 5
   ann gmtint anninc 5
   opath /export/home/EW0008/Realtimestack/Line$LINENO/siopltfil.Line$LINENO
  end
end
  
diskoe   # Write out disk file
  opath /export/home/EW0008/Realtimestack/Line$LINENO/stack.line$LINENO 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 subwb yes tadd 3 end
end

nmo
 addwb yes stretc 1.
  vtp 1510 0 1530 .3 1550 .4 1510 .5 1560 .6 1580 .9 1650 1.5 end
end

end

eof
Return to EW0008 example.      Return to SIOSEIS examples.

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