#! /bin/sh
set -x

fx=0 nx=721 dx=12.5
fxo=0 nxo=185 dxo=50
ft=0 nt=501 dt=0.005
fzo=0 nzo=61 dzo=50
fz=0 nz=601 dz=5

veldata=model.200
veldata2=model.200o
veldata3=junk2
suaddhead <$veldata ns=$nzo |
sushw key=dt,tracl,cdp,tracr a=4000,1,1,1 b=0,1,1,1 >junk
suinterp <junk linear=0 ninterp=3 |
suwind count=721 |
sustrip >$veldata3
#ximage <$veldata3 n1=$nzo
convtz <$veldata3 type=1 vtype=1 vfile=$veldata3 \
        nx=$nx nz=$nzo dz=$dzo nt=$nt dt=$dt >junk
convtz <junk type=0 vtype=0 vfile=junk \
        nx=$nx nz=$nz dz=$dz nt=$nt dt=$dt >$veldata2
#ximage <$veldata2 n1=$nz

sugain <kzmig_all.stack.taui_ti pbal=1 |
sustrip |
conv_tz type=0 vtype=0 vfile=junk int_type=0 \
	nx=$nx nz=$nz dz=$dz nt=$nt dt=$dt |
ximage perc=95 wbox=700 hbox=400 f2=0 d2=0.0125 x2beg=3.0 \
title="$1 Pre-stack depth migration $1" titlesize=25 \
label1="Depth (km)" label2="Midpoint (km)" \
f1=0 d1=$dz n1=$nz d2num=1 n2=$nx &
#supsimage >kzmig.stack.ps