3-3 波形比较

RR1=R-1/1/0/0.5
RR2=R-1/1/0/0.5
JJ=JX1.5c/-15c
PS=wave_compare.ps
TXT1=point.txt
TXT2=4_VZ_t2.txt
cp ../$TXT1 .
cp ../$TXT2 .
gmt psxy -JX10c/-15c -BWN -Bya0.1+l"time(s)" -K -$RR1 -P >$PS <<EOF
1.0 1.0
EOF
for((k=1;k<7;k++)); do
    #kk=$[$k-1]
    gmt psxy $TXT1 -D$[k]c/0c -W1.3p -i$k,0 -$JJ -$RR1 -K -O  >>$PS;
    gmt psxy $TXT2 -D$[k]c/0c -W1.3p  -A -i$k,0 -$JJ -$RR2 -Wred -K -O >>$PS
done
 gmt psxy $TXT1 -D7c/0c -i7,0 -$JJ -W1.3p  -$RR1  -K -O  >>$PS;
 gmt psxy $TXT2 -D7c/0c -i7,0 -$JJ -W1.3p -$RR2 -Wred -K -O >>$PS
# is fault
gmt pslegend -DjTL+w10c/2c/+o0.5c/-0.8c -JX10c/15c -$RR1 -O  <<EOF >>$PS 
N 2
S 0.2c - 0.5c - 1p,black 0.5c original wave 
S 0.2c - 0.5c - 1p,red  0.5c extracted wave 
EOF
gmt psconvert -A -Teg $PS

wave_compare.png
更多资料
GMT目录