#! /bin/csh
#
#   run_fkrsort.eqdist in Tohoku
#
# Window out the eight vectors for many distances
# Remember to set NT and DT CORRECTLY!
# USE wvint9 (The flipped traces are corrected within this
# code rather than by the external program flip
# set path=( $path ~/FK )
#
#     CHANGE to where you keep the executables etc.
#
set path=( $path /data/a1/fumiko/MT.dir/Doug.dir/MTCODE/BIN )
#         for window
#
# set path=( $path /data/04/dreger/Bin)
#           for mkHelm
#
set dt=0.50
# set npts=512
#        ....
set npts=1024
set dist=30
set depth=$1
set count=0
set j=1
set vshift=0
set i=0
# set loopend=75 <--- number of the distances
set loopend=95
set nvec=0
rehash
#
#
@ nvec=($loopend - $count) * 10
while ($count < $loopend)
@ vshift=$i + 7
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp1$$
@ vshift=$i + 4
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp2$$
@ vshift=$i + 6
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp3$$
@ vshift=$i + 3
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp4$$
@ vshift=$i + 1
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp5$$
@ vshift=$i + 5
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp6$$
@ vshift=$i + 2
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp7$$
@ vshift=$i
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp8$$
cat tmp1$$ tmp2$$ tmp3$$ tmp4$$ tmp5$$ tmp6$$ tmp7$$ tmp8$$ > junk
echo $j
mkHelm format="(6e12.5)" ntr=8 dt=$dt nt=$npts < junk > gil7_{$dist}d{$depth}.disp
compress gil7_{$dist}d{$depth}.disp
#
#          /\ CHANGE THE MODEL NAME
#
\rm tmp*$$
@ i += 10
@ count++
@ j++
@ dist += 5
end