EW9914 SIOSEIS "realtime" stack script

run on "heezen" in directory:    /ldata/realtime/brute_stack

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

sioseis << eof

procs segdin prout geom header diskob gather nmo mute diskod stack diskoe filter diskoa agc plot end

segdin
# take all shots 
   ffilen 99999   
#   ftr 1 ltr 480 trinc 2 fcset 1 lcset 1
   ftr 241 ltr 480 fcset 1 lcset 1
   stime 1 secs 5.0 
   offline yes  # eject after the rewind after EOT
   iunit 43 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 -196.6   # RESET the closest group only.
  ggx -12.5         # Used to extrapolate gxp!
 dfls 37.5 dbrps 6.25 smear 6.25
#  dfls 37.5 dbrps 12.5 smear 12.5
  rpadd 1000 end
end

mute
   fno 1 lno 999999
   addwb yes xtp 200 -.1 1500 -.1 3000 1 6200 2 end
end
  
diskoa   # Write out disk file
  opath LAUBASIN.$LINENO end
end
diskob 
# write every 50th shot to a "circular" file 
    fno 1 lno 999999 noinc 50 rewind 1 
    opath /ldata/realtime/shots/latest.shot end 
end
diskoc 
    set 3 6 opath /ldata/realtime/cmps/cmps.$LINENO end 
end
diskod 
    fno 1 lno 999999 noinc 50 rewind 1 
    opath /ldata/realtime/shots/latest.mute end 
end
diskoe   # Write out disk file
  opath stack.$LINENO 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
   end
end

gather
#   maxtrs 90 maxrps 500 end
   maxtrs 50 maxrps 250 end
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 30 
  def 0.01 trpin 80 wiggle 0 
  stime 2 nsecs 4 vscale 5 
  opath siopltfil.$LINENO end 
 end 


end
eof