例1：GMTに最初から含まれている地図データで書いたもの

ーシェルスクリプトはこちらー



#/bin/csh!
######################
set psfile=map.ps
set range=0/360/40/90
set scale=A180/90/7.5
set fontsize=25
set fontsize1=15
set fonttipe=1
set fonttipe1=7
set w=20/255/0/0
########################
pstext -JX6/6 -R-5/10/-5/10 -Y8 -X1.5 -P -K << END >! $psfile
0 0 40 0 $fonttipe 5 MAP
END
pscoast -J$scale -R$range -Ba20g10/a20g10:."":newS -G155/220/0 -S220/220/255 -W1 -Dc -A1000 -X-1.1 -Y-6.5 -O -K>> $psfile
psxy bound-ba.dat -J$scale -R$range -O -W$w -K >> $psfile
psxy bound-be.dat -J$scale -R$range -O -W$w -K >> $psfile
psxy bound-la.dat -J$scale -R$range -O -W$w -K >> $psfile
psxy bound-gr.dat -J$scale -R$range -O -W$w -K >> $psfile
psxy bound-ok.dat -J$scale -R$range -O -W$w -K >> $psfile
pstext -JX10/10 -R-5/10/-5/10 -X-3 -Y-3 -P -O << END >> $psfile
3 6.5 $fontsize 0 $fonttipe 5 Barents
5 7 $fontsize 0 $fonttipe 5 Greenland
7 6 $fontsize 0 $fonttipe 5 Labrador
4.5 1.5 $fontsize 0 $fonttipe 5 Bering
2 1.8 $fontsize 0 $fonttipe 5 Okhotsk
########
2.1 0.5 $fontsize1 0 $fonttipe1 5 Hokkaido
3.5 0 $fontsize1 65 $fonttipe1 5 Kamchatka Pen.
5.7 2.7 $fontsize1 0 $fonttipe1 5 ALASKA
9 3 $fontsize1 65 $fonttipe1 5 CANADA
3 6 $fontsize1 0 $fonttipe1 5 Novaya Zemlya
4.5 0 $fontsize1 0 $fonttipe1 5 PACIFIC OCEAN
4.2 5 $fontsize1 0 $fonttipe1 5 ARCTIC OCEAN
7.7 5.5 $fontsize1 0 $fonttipe1 5 Hudoson Bay
2 5 $fontsize1 -70 $fonttipe1 5 RUSSIAN
6 9 $fontsize1 -30 $fonttipe1 5 ATLANTIC OCEAN
END
gs $psfile


例２：特別に地図データを持ってきてそれを図化したもの
#! /bin/tcsh -f
gmtset BASEMAP_TYPE FANCY
gmtset DEGREE_FORMAT 3
gmtset HEADER_FONT 7
gmtset LABEL_FONT 9
gmtset LABEL_FONT_SIZE 10
#####################
set file=map_gmt
set gfile=$file.grd
set psfile=$file.ps3
set resol=0.01
set range=135/145/42/50
set scale=m1.5
#set range=140/142/42/44
#set scale=m6
set color=cpt.dat
###########################
xyz2grd $file -G$gfile -I$resol -R$range -V
makecpt -T-100/2500/100 -Ctopo -V > $color
grdimage $gfile -J$scale -R$range -C$color -Y8 -P -K >! $psfile
psbasemap -R$range -J$scale -Ba1f1:"Longitude":/a1f1:"Latitude"::."Gtopo $range":neWS -O -K >> &$psfile
psscale -C$color -P -D15/5/10/0.5 -O >> $psfile
rm $gfile $color
gs $psfile



例３：元のデータが粗いとモザイク状になってしまう(シェルスクリプトは省略) モドリマスカ？