#!  /bin/csh -fx
#  Script to run the basic 3D inverse Q estimation/filter

#  qest (estimation)
#  pipe Qs into scale to invert Q's (to absorption). This is analogous to
#  smoothing slownesses rather than velocities. We then smooth temporally
#  with an 11-point boxcar and then pipe directly into the data slicer to
#  get to the time slice domain.

qest -Nclassdata | scale -exp-1 | mixn -m11 -w1 |\
qdslice -Otmp1

#  In TI domain do an 11x11 spatial smoother

fkkstrip -Ntmp1 -Otmp2 -l11 -pass -dup -xtaper0 -ytaper0

#  Slice TI's back into LI/DI domain

qdslice -Ntmp2 -R | scale -exp-1 | hdrswap -N2classdata -Oclassdata_qsm

#  qifl (inverse Q application)

qifl -Nclassdata -Q -qtclassdata_qsm -Oclassdata_iqf -fc70