psgain

su 的源代码库

#! /bin/sh
# Try some display gains on the muted, filtered data
demo=gain

perc=90

f=0,15,25,30
amps=0.,0.,1.,1.

xmute=1,49,97
tmute=1.2,.2,1.2
ntaper=5

tpow=1
sufilter <data.su f=$f amps=$amps |
sumute key=tracl xmute=$xmute tmute=$tmute ntaper=$ntaper |
sugain tpow=$tpow |
supsimage  perc=$perc \
	title="filter, mute & tpow=$tpow"\
	label1="Time (sec)" label2="Offset (km)" >${demo}.1.eps

tpow=2
sufilter <data.su f=$f amps=$amps |
sumute key=tracl xmute=$xmute tmute=$tmute ntaper=$ntaper |
sugain tpow=$tpow |
supsimage  perc=$perc \
	title="filter, mute & tpow=$tpow"\
	label1="Time (sec)" label2="Offset (km)" >${demo}.2.eps

tpow=3
sufilter <data.su f=$f amps=$amps |
sumute key=tracl xmute=$xmute tmute=$tmute ntaper=$ntaper |
sugain tpow=$tpow |
supsimage  perc=$perc \
	title="filter, mute & tpow=$tpow"\
	label1="Time (sec)" label2="Offset (km)" >${demo}.3.eps

tpow=2.2
sufilter <data.su f=$f amps=$amps |
sumute key=tracl xmute=$xmute tmute=$tmute ntaper=$ntaper |
sugain tpow=$tpow |
supsimage  perc=$perc \
	title="filter, mute & tpow=$tpow"\
	label1="Time (sec)" label2="Offset (km)" >${demo}.4.eps

echo 
echo "Output PostScript files:"

echo -n "     "
for i in ${demo}.*.eps
do
	 echo -n "$i "
done
echo
echo "You may view these with your PostScript Previewer"

exit 0