#! /bin/csh
# 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 ../BIN)
#
set dt=1.00
set npts=512
##
set dist=(604 768 958)
set loopend=3
set depth=10
##
set count=0
set j=1
set vshift=0
set i=0
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$$
@ vshift=$i + 8
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp10$$
@ vshift=$i + 9
window v0=$vshift e0=0 nt=$npts nx=$nvec nv=1 < vec > tmp9$$
cat tmp1$$ tmp2$$ tmp3$$ tmp4$$ tmp5$$ tmp6$$ tmp7$$ tmp8$$ tmp9$$ tmp10$$ > junk
echo $j
mkHelm format="(6e13.5)" ntr=10 dt=$dt nt=$npts < junk > socal_{$dist[$j]}d{$depth}.disp
\rm tmp*$$
@ i += 10
@ count++
@ j++
end
#putmech in=test_604d10.disp out=test1 azimuth=247.7 strike=93. rake=5. dip=83. moment=1.2e25
#putmech in=test_768d10.disp out=test2 azimuth=99.6 strike=93. rake=5. dip=83. moment=1.2e25
#putmech in=test_958d10.disp out=test3 azimuth=7.5 strike=93. rake=5. dip=83. moment=1.2e25
##PAYG
#sac2bin in=test1.tan out=junk
#bin2sac par=b2s.par stime=60.4 < junk > test1.tan
#sac2bin in=test1.rad out=junk
#bin2sac par=b2s.par stime=60.4 < junk > test1.rad
#sac2bin in=test1.ver out=junk
#bin2sac par=b2s.par stime=60.4 < junk > test1.ver
##OTAV
#sac2bin in=test2.tan out=junk
#bin2sac par=b2s.par stime=76.8 < junk > test2.tan
#sac2bin in=test2.rad out=junk
#bin2sac par=b2s.par stime=76.8 < junk > test2.rad
#sac2bin in=test2.ver out=junk
#bin2sac par=b2s.par stime=76.8 < junk > test2.ver
##HDC
#sac2bin in=test3.tan out=junk
#bin2sac par=b2s.par stime=95.8 < junk > test3.tan
#sac2bin in=test3.rad out=junk
#bin2sac par=b2s.par stime=95.8 < junk > test3.rad
#sac2bin in=test3.ver out=junk
#bin2sac par=b2s.par stime=95.8 < junk > test3.ver
