plot plate boundary 〜プレート境界を入れる 

  使うコマンド  psxy



bash スクリプト例
#  plot plate boundary

# parameter setting
region=115/155/0/45                   # map region east/west/south/north
proj=M10                                # map projection and scale
ticks=a10f1g10                            # boundary tick info
frame=WSne+tplate_boundary                     # boundary frame info
ridgefile=ridge.d
transfile=fracture.d
trenchfileL=trenchesL
trenchfileR=trenchesR
psfile=plate.ps                  # output postscript file name
#
gmt psxy $transfile -R$region -J$proj -Wthinner,green -K -V > $psfile
gmt psxy $ridgefile -R$region -J$proj -Wthin,red -K -O -V >> $psfile
gmt psxy $trenchfileL -R$region -J$proj -Sf0.3/0.1+l+t -Wthin,blue -Gblue -K -O -V >> $psfile
gmt psxy $trenchfileR -R$region -J$proj -Sf0.3/0.1+r+t -Wthin,blue -Gblue -K -O -V >> $psfile
gmt pscoast -R$region -J$proj -Di -Wthin -Ggray -O -K >> $psfile
gmt psbasemap -R$region -J$proj -B$ticks -B$frame -O >> $psfile
#


Tips