
本来は、プログラムをくんで自動化しているのでしょうがここではGMTの描き方と言うことで了承してください。
「地図とグラフもこんな風に一緒に書けるんだよ」っていうこと。
シェルスクリプトはこちら
#/bin/csh!
gmtset HEADER_FONT 4
###########################################
set file=rensyu2a.dat
set file1=rensyu2b.dat
set file2=rensyu2c.dat
set fileps=sample6.ps
set resol=1
set range=138/145/35/46
set range1=-100/0/35/46
set range2=138/145/-100/0
set scale=m0.5
set fontsize=25
set fonttipe=7
set l=1/30
set g=2/2
#############################################
pstext -JX5/5 -R-5/10/-5/10 -K -Y2.8 -X1 -P << END >! $fileps
-5 -3 40 0 $fonttipe 5 Sample1
END
####################
pscoast -J$scale -R$range -Ba1g1/a1g1:."":neWS -G155/220/0 -S220/220/255 -W1 -Dl -O -K >> $fileps
psxy -J$scale -R -W -G255/0/0 -Sc0.1 -L -O -K -P $file2 >> $fileps
####################
psxy -Jx-0.02/0.65 -R$range1 -Ba20f20g20/a1f1g1neWS -W10 -O -K -P -X4.5 -Y0 $file >> $fileps
psxy -Jx-0.02/0.65 -R -W -G255/0/0 -Sc0.1 -L -O -K -P $file >> $fileps
pstext -JX -R-5/10/-5/10 -X-3.3 -Y2 -O -K -P << END >> $fileps
3.2 -3 $fontsize 90 $fonttipe 5 Latitude(Degrees)
END
pstext -JX -R-5/10/-5/10 -X3.3 -Y-4.3 -O -P -K << END >> $fileps
-4 0 $fontsize 0 $fonttipe 5 Depth(km)
END
####################
psxy -Jx0.5/0.015 -R$range2 -Ba1f1g1/a20f20g20neWs -W10 -O -K -P -Y0.5 -X-4.5 $file1 >> $fileps
psxy -Jx0.5/0.015 -R -W -G255/0/0 -Sc0.1 -L -O -K -P -Y-0.5 $file1 >> $fileps
pstext -JX -R-5/10/-5/10 -X-1 -O -K -P << END >> $fileps
-4 -4 $fontsize 90 $fonttipe 5 Depth(km)
END
pstext -JX -R-5/10/-10/10 -X1 -O -P << END >> $fileps
0 -9 $fontsize 0 $fonttipe 5 Longitude(Degrees)
END
gs $fileps
#
モドリマスカ?