はじめに
メインで使っているのはMacbookなのですが、どうしても会社のWindows7機を使わざるを得ず、必要なソフト類をインストールしました。
これらの中で、GMT5は、Windowsコマンドプロンプトでバッチファイルを実行するため、シェルスクリプトで書いたものは使えず、Windowsコマンドプロンプト用に修正しなければなりません。
特に最初はインストールしたソフトが期待通りの動作をするかを確認するためのテストランを行うわけですが、このテストラン用バッチファイルを作るのが面倒。いつもその場その場でマシンが変わるたびにMac用スクリプトを書き換えていたのですが、この際、今回作成したWindows用データをアップしておくことにしました。
動作確認したGMTのバージョンは5.2.1です。GMT4は対象外ですので注意してください。
一般的なグラフ作成
a_gra.bat
echo 0.01 -6 > _inp_a.txt
echo 0.1 -4 >> _inp_a.txt
echo 1 0 >> _inp_a.txt
echo 20 8 >> _inp_a.txt
echo 50 4 >> _inp_a.txt
echo 100 3 >> _inp_a.txt
echo 0.01 3 > _inp_b.txt
echo 0.1 4 >> _inp_b.txt
echo 1 5 >> _inp_b.txt
echo 7 2 >> _inp_b.txt
echo 50 -4 >> _inp_b.txt
echo 100 -8 >> _inp_b.txt
set range=0.01/100/-10/10
set scale=5l/3
set xga=a1g3
set yga=g2a4
set xlabel="x-label"
set ylabel="y-label"
set fig=fig_gra_slog.eps
gmt set FONT_ANNOT_PRIMARY 12
gmt set FONT_LABEL 8
gmt set MAP_TICK_LENGTH_PRIMARY 0c
gmt psbasemap -R%range% -JX%scale% -Bx%xga%+l"%xlabel%" -By%yga%+l"%ylabel%" -BWSen -P -K > %fig%
gmt psxy _inp_a.txt -R -J -W1.5 -K -O >> %fig%
gmt psxy _inp_a.txt -R -J -SC0.2 -G#000000 -N -K -O >> %fig%
gmt psxy _inp_b.txt -R -J -W1.5,3_3:0 -K -O >> %fig%
gmt psxy _inp_b.txt -R -J -SC0.2 -W1,#0000ff -G#ffffff -N -K -O >> %fig%
echo 0.01 -10 | gmt psxy -R -J -Sp -O >> %fig%
echo 2.5 1 > _inp_a.txt
echo 7.5 5 >> _inp_a.txt
echo 12.5 8 >> _inp_a.txt
echo 17.5 3 >> _inp_a.txt
echo -75 40 > _inp_b.txt
echo -25 20 >> _inp_b.txt
echo 25 -30 >> _inp_b.txt
echo 70 30 >> _inp_b.txt
set rangeSW=0/20/0/10
set rangeNE=-100/100/-50/50
set scale=5/3
set afgS=a5f1g5
set afgW=a2f1g2
set afgN=a50f10
set afgE=a20f10
set labelS="Label-S"
set labelW="Label-W"
set labelN="Label-N"
set labelE="Label-E"
set fig=fig_gra_2axes.eps
gmt set FONT_ANNOT_PRIMARY 8
gmt set FONT_LABEL 12
gmt set MAP_TICK_LENGTH_PRIMARY 0.2c
gmt psbasemap -R%rangeSW% -JX%scale% -Bx%afgS%+l"%labelS%" -By%afgW%+l"%labelW%" -BWS -P -K > %fig%
gmt psxy _inp_a.txt -R -J -W2 -K -O >> %fig%
echo 17.5 3 10p 0 LM SW | gmt pstext -R -J -F+f+a+j -N -K -O >> %fig%
gmt psbasemap -R%rangeNE% -JX%scale% -Bx%afgN%+l"%labelN%" -By%afgE%+l"%labelE%" -BNE -K -O >> %fig%
gmt psxy _inp_b.txt -R -J -W2,5_5:0 -K -O >> %fig%
echo 70 30 10p 0 LM NE | gmt pstext -R -J -F+f+a+j -N -K -O >> %fig%
echo -100 -50 | gmt psxy -R -J -Sp -O >> %fig%
何種類かの世界地図を描画
a_world.bat
rem Hemisphere map
gmt pscoast -Rg -JA75/30/4.5i -Bg -Dc -A1000 -Gnavy -P > fig_map_lambert_az_hemi.eps
rem Orthographic projection (-Jg -JG)
gmt pscoast -Rg -JG75/30/4.5i -Bg -Dc -A5000 -Gpink -Sthistle -P > fig_map_orthographic.eps
rem Mercator projection (-Jm -JM)
gmt set MAP_FRAME_TYPE fancy
gmt pscoast -R-90/270/-75/75 -Jm1.2e-2i -Bxa30f15 -Bya30f15 -Dc -A5000 -Gred ^
-P > fig_map_mercator.eps
rem Cylindrical equidistant projection (-Jq -JQ)
gmt pscoast -Rg -JQ90/4.5i -B30f30g30 -Dc -A5000 -Gtan4 -Slightcyan -P > fig_map_equi_cyl.eps
rem Cylindrical equal-area projections (-Jy -JY)
gmt pscoast -R-90/270/-90/90 -JY90/30/4.5i -B30g30 -Dc -A10000 -Sdodgerblue ^
-Wthinnest -P > fig_map_general_cyl.eps
rem Miller Cylindrical projection (-Jj -JJ)
gmt pscoast -R-90/270/-80/90 -Jj1:400000000 -Bx30g30 -By30g30 -Dc -A10000 ^
-Gkhaki -Wthinnest -P -Sazure > fig_map_miller.eps
rem Cylindrical stereographic projections (-Jcyl_stere -JCyl_stere)
gmt set FORMAT_GEO_MAP dddA
gmt pscoast -R-90/270/-60/80 -JCyl_stere/90/45/4.5i -Bxa30f30g30 -Bya30g30 -Dc -A5000 ^
-Wblack -Gseashell4 -Santiquewhite1 -P > fig_map_gall_stereo.eps
rem Hammer projection (-Jh -JH)
gmt pscoast -Rg -JH75/4.5i -Bg -Dc -A10000 -Gblack -Scornsilk -P > fig_map_hammer.eps
rem Mollweide projection (-Jw -JW)
gmt pscoast -Rd -JW75/4.5i -Bg -Dc -A10000 -Gtomato1 -Sskyblue -P > fig_map_mollweide.eps
rem Winkel Tripel projection (-Jr -JR)
gmt pscoast -Rd -JR75/4.5i -Bg -Dc -A10000 -Gburlywood4 -Swheat1 -P > fig_map_winkel.eps
rem Robinson projection (-Jn -JN)
gmt pscoast -Rd -JN75/4.5i -Bg -Dc -A10000 -Ggoldenrod -Ssnow2 -P > fig_map_robinson.eps
rem Eckert IV and VI projection (-Jk -JK)
gmt pscoast -Rg -JKf75/4.5i -Bg -Dc -A10000 -Wthinnest -Givory -Sbisque3 -P > fig_map_eckert4.eps
gmt pscoast -Rg -JKs75/4.5i -Bg -Dc -A10000 -Wthinnest -Givory -Sbisque3 -P > fig_map_eckert6.eps
以上