j.vstack

来自「seismic software,very useful」· VSTACK 代码 · 共 23 行

#! /bin/ksh
# Velocity analyses for the cmp gathers
#set -x

SCRATCH=/scratch/jkc/Nmodata

# Determine velocity sampling.
nv=100 dv=50.0 fv=4000.0
integer cdpmin=1500 cdpmax=3000 dcdp=100 fold=64 cdp

# Get a velocity stack.
>vstack
cdp=cdpmin
while ((cdp<=cdpmax))
do
	suwind <$SCRATCH/traces.su key=cdp min=$cdp max=$cdp count=$fold |
	suvelan nv=$nv dv=$dv fv=$fv |
	sustack >>vstack
	cdp=cdp+dcdp
done

suxwigb <vstack