gmt gmtset MAP_ORIGIN_Y 15c
RR=R0/0.013/0/0.013
JJ=JX8c/8c
ps=merge.ps
FI=slowness2
TXT=summation_slowness_t2.txt
cp ~/matlabfile/graduate/jieguo3/TXT -FI.grd
gmt grd2cpt FI.cpt
gmt grdimage -K RR -FI.cpt -BWSne -Bxa0.004+l"U@-s@-(s/m)" -Bya0.004+l"u@-r@-(s/m)" -P >JJ -ps
0.0065 0.0120 t=0.285s
EOF

second pic

rm gmt.*
FI2=slowness3
TXT2=summation_slowness_t3.txt
cp ~/matlabfile/graduate/jieguo3/TXT2 -FI2.grd
gmt grd2cpt --COLOR_BACKGROUND=0/0/127 FI2.cpt
gmt grdimage -X9c -K -O RR -FI2.cpt -BWSne -Bxa0.004+l"U@-s@-(s/m)" >>RR -ps
0.0065 0.0120 t=0.366s
EOF

the third pic

FI3=slowness4
TXT3=summation_slowness_t4.txt
cp ~/matlabfile/graduate/jieguo3/TXT3 -FI3.grd
gmt grd2cpt FI3.cpt
gmt grdimage -X-9c -Y-10c -K -O RR -FI3.cpt -BWSne -Bxa0.004+l"U@-s@-(s/m)" -Bya0.002+l"u@-r@-(s/m)" >>JJ -ps
0.0065 0.0120 t=0.410s
EOF

gmt psscale -O -B2+l"relative power(dB)" -D12c/0c+w8c/0.5c -Cps
gmt psconvert -A -Teg $ps
rm gmt.* *.cpt *.grd

8.png
更多资料
GMT目录