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 eofReturn to EW0008 example. Return to SIOSEIS examples.
Go to the list of seismic processes. Go to SIOSEIS introduction.