#! /bin/csh if( $#argv < 1 ) then echo "Usage: stack line-number" exit 1 endif set LINENO = $1 # This should be run in /export/home2/scratch/ew0207/brute so the # sioseis tmp files do not conflict with other sioseis scripts /export/home2/scratch/ew0207/bin/sioseis << eof procs segddin prout weight geom header gather nmo mute diskod stack diskoe filter agc plot end segddin ftr 241 ltr 480 fcset 1 lcset 1 stime 1 secs 6.0 format seisnet stack /export/home2/scratch/ew0207/logs/latest end end weight fno 0 lno 999999 twp 236 0 end end geom type 2 # increment the shot loaction based on the shot number fs 1 ls 999999 # all shot have the same parameters (preset) gxp 480 -180 # RESET the closest group only. ggx -12.5 # Used to extrapolate gxp! dfls 37.5 dbrps 6.25 smear 6.25 # use this if using trinc 1 # dfls 37.5 dbrps 12.5 smear 12.5 # use this if using trinc 2 rpadd 1000 end end mute fno 1 lno 999999 ! addwb yes xtp 200 -.1 1500 -.1 3000 1 6200 2 end addwb yes xtp 200 -.1 2500 -.1 4000 1 end end diskoa # Write out the filtered stack file disk file opath /export/home2/scratch/ew0207/line$LINENO.stack-filter end end diskob # write every 50th shot to a "circular" file # remember that segdin limited the traces read (ftr/ltr/trinc) fno 1 lno 999999 noinc 50 rewind 1 opath /export/home2/scratch/ew0207/shots/latest.shot end end diskod # write out the muted gather fno 1 lno 999999 noinc 50 rewind 1 opath /export/home2/scratch/ew0207/shots/latest.mute1 end end diskoe # Write out disk file opath /export/home2/scratch/ew0207/stacked_lines/line$LINENO.stack end end prout fno 0 lno 99999 ftr 479 ltr 479 noinc 10 end end header fno 0 lno 9999999 ftr 1 ltr 9999 r50 r54 / 750. # convert water depth to water time for addwb (mute) end end gather # maxtrs 90 maxrps 500 end maxtrs 50 maxrps 250 end # half the streamer end nmo # real time nmo, replace interpolation by RP to WB depth in Meters. # If water depth changes by > 500 m, use previous value. Water-depth # velocity functions derived from ESP5, interpolation by iso-velocity layering vtrkwb 500 stretc 0.50 fno 1000 lno 1000 vtp 1500 1.333 1557 1.414 1607 1.443 1789 1.492 2346 1.645 2638 1.746 2900 1.846 2971 1.872 3150 1.983 3141 2.102 3264 2.362 4228 3.742 4343 3.892 4898 4.393 end fno 1500 lno 1500 vtp 1500 2.0 1539 2.081 1574 2.110 1705 2.159 2137 2.312 2379 2.413 2603 2.513 2665 2.539 2827 2.650 2834 2.769 2967 3.029 3939 4.409 4053 4.559 4596 5.060 end fno 2000 lno 2000 vtp 1500 2.667 1529 2.748 1557 2.777 1659 2.826 2012 2.979 2218 3.080 2414 3.180 2468 3.206 2614 3.317 2629 3.436 2761 3.696 3711 5.076 3823 5.226 4351 5.727 end fno 2500 lno 2500 vtp 1500 3.333 1524 3.414 1546 3.443 1629 3.492 1928 3.645 2108 3.746 2282 3.846 2330 3.872 2463 3.983 2481 4.102 2608 4.362 3526 5.742 3636 5.892 4146 6.393 end fno 3000 lno 3000 vtp 1500 4.0 1520 4.080 1538 4.110 1609 4.159 1868 4.312 2028 4.413 2184 4.513 2228 4.539 2350 4.650 2368 4.769 2489 5.029 3373 6.409 3479 6.559 3972 7.060 end fno 3500 lno 3500 vtp 1500 4.667 1517 4.748 1533 4.777 1595 4.826 1823 4.979 1967 5.080 2108 5.180 2148 5.206 2260 5.317 2279 5.436 2395 5.696 3243 7.076 3346 7.226 3822 7.727 end fno 4000 lno 4000 vtp 1500 5.333 1515 5.414 1529 5.443 1583 5.492 1788 5.645 1919 5.746 2048 5.846 2085 5.872 2189 5.983 2208 6.102 2317 6.362 3131 7.742 3231 7.892 3692 8.393 end fno 4500 lno 4500 vtp 1500 6.0 1513 6.081 1526 6.110 1574 6.159 1760 6.312 1879 6.413 1999 6.513 2033 6.539 2130 6.650 2148 6.769 2252 7.029 3034 8.409 3131 8.559 3577 9.060 end fno 5000 lno 5000 vtp 1500 6.667 1512 6.748 1523 6.777 1567 6.826 1737 6.979 1847 7.080 1958 7.180 1990 7.206 2080 7.317 2098 7.436 2197 7.696 2948 9.076 3042 9.226 3474 9.727 end end filter pass 5 40 ftype 0 dbdrop 48 end end agc winlen .5 center .1 end end plot scalar 1.e-07 tlines 0.5 1 nibs 7224 ann gmtint anninc 5 def 0.01 trpin 80 wiggle 0 stime 3.0 nsecs 4 vscale 5 opath /export/home2/scratch/ew0207/stack_plots/line$LINENO.atlantek srpath /export/home2/scratch/ew0207/stack_plots/line$LINENO.ras end end end eofReturn to SIOSEIS examples.
Go to the list of seismic processes. Go to SIOSEIS introduction.