create_MIG_panel.sh

#!/bin/sh
# Course: Geophysical Exploration
#
# Sample script : Migration
# Apply migration with increasing velocity
# 
# by W.T.
#

# Edit the name of the input and oudput data
indata=Nzero.su
outdata=cvmpanel.su

# Parameters
vmin=1000 # Minimum velocity (m/s)
vmax=2000 # Maximum velocity (m/s)
dv=100    # Increment of velocity (m/s)
perc=98
#--

nt=`surange < $indata | grep ns | awk '{print $2}'`
dt=`surange < $indata | grep dt | awk '{print $2}'`

cdpmin=900
cdpmax=1300
dxcdp=16.667
tmig=6.0

# Original zero-offset section
cat $indata > $outdata
sunull nt=$nt dt=$dt ntr=20 >>$outdata

# Migrated section
v=$vmin
while [ $v -le $vmax ]; do
    echo $v
    sustolt cdpmin=$cdpmin cdpmax=$cdpmax dxcdp=$dxcdp \
	    vmig=$v tmig=$tmig lstaper=20 lbtaper=100 \
            < $indata >> $outdata
    sunull nt=$nt dt=$dt ntr=5 >>$outdata
    v=`expr $v + $dv`
done

# Plot
suximage < $outdata perc=$perc \
	 x1beg=5.0 x1end=9.0 \
	 label1='Time (s)'&

exit