WANtaro雑記帳（Win_GMT：エラーバー，プロット寸法，プロット色）

rem ***************************************************************************
rem * エラーバー，プロット寸法，プロット色                                    *
rem ***************************************************************************
del .*

set range=0/100/0/100
set scale=10/10
set xga=g10a10
set yga=g10a10
set xlabel="Name of x-axis"
set ylabel="Name of y-axis"
rem
gmtset ANOT_FONT_SIZE 10
gmtset LABEL_FONT_SIZE 10
gmtset TICK_LENGTH 0c



rem ***************************************************************************
rem * 入力用ファイル「inp_wg_31.txt」の作成                                      *
rem * (x座標 y座標 y方向エラーバー長さ)                                            *
rem ***************************************************************************
echo  10 20 5 > inp_wg_31.txt
echo  20 30 5 >> inp_wg_31.txt
echo  30 50 5 >> inp_wg_31.txt
echo  40 80 5 >> inp_wg_31.txt
echo  50 90 5 >> inp_wg_31.txt
echo  60 70 5 >> inp_wg_31.txt
echo  70 60 5 >> inp_wg_31.txt
echo  80 30 5 >> inp_wg_31.txt
echo  90 10 5 >> inp_wg_31.txt
rem
set inp_dat=inp_wg_31.txt
set fig_out=fig_wg_31.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig_out%
psxy %inp_dat% -R -JX -W2 -K -O >> %fig_out%
psxy %inp_dat% -R -JX -SC0.2 -G0 -Ey -N -O >> %fig_out%
rem ***************************************************************************
rem * y方向のエラーバーを描画                                                 *
rem * -SC0.2 -G0 入力された数値群の位置に(C)0.2cmの円を描き(G)黒で塗りつぶす．*
rem * -Ey        y方向にエラーバーを描画．                                    *
rem *            入力ファイルの3列目が片側の長さに対応．                      *
rem ***************************************************************************



rem ***************************************************************************
rem * 入力用ファイル「inp_wg_32.txt」の作成                                      *
rem * (x座標 y座標 x方向エラーバー長さ y方向エラーバー長さ)                   *
rem ***************************************************************************
echo  10 20 5 3> inp_wg_32.txt
echo  20 30 5 3 >> inp_wg_32.txt
echo  30 50 5 3 >> inp_wg_32.txt
echo  40 80 5 3 >> inp_wg_32.txt
echo  50 90 5 3 >> inp_wg_32.txt
echo  60 70 5 3 >> inp_wg_32.txt
echo  70 60 5 3 >> inp_wg_32.txt
echo  80 30 5 3 >> inp_wg_32.txt
echo  90 10 5 3 >> inp_wg_32.txt
rem
set inp_dat=inp_wg_32.txt
set fig_out=fig_wg_32.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig_out%
psxy %inp_dat% -R -JX -W2 -K -O >> %fig_out%
psxy %inp_dat% -R -JX -SC0.2 -G0 -Ex -Ey -N -O >> %fig_out%
rem ***************************************************************************
rem * x,y方向のエラーバーを描画                                               *
rem * -SC0.2 -G0 入力された数値群の位置に(C)0.2cmの円を描き(G)黒で塗りつぶす．*
rem * -Ex        x方向にエラーバーを描画．                                    *
rem *            入力ファイルの3列目が片側の長さに対応．                      *
rem * -Ey        y方向にエラーバーを描画．                                    *
rem *            入力ファイルの4列目が片側の長さに対応．                      *
rem ***************************************************************************



rem ***************************************************************************
rem * 入力用ファイル「inp_wg_33.txt」の作成                                      *
rem * (x座標 y座標 z値)：z値はプロットの大きさ                                *
rem ***************************************************************************
echo  10 20 0.1 > inp_wg_33.txt
echo  20 30 0.3 >> inp_wg_33.txt
echo  30 50 0.5 >> inp_wg_33.txt
echo  40 80 0.7 >> inp_wg_33.txt
echo  50 90 1.0 >> inp_wg_33.txt
echo  60 70 0.8 >> inp_wg_33.txt
echo  70 60 0.6 >> inp_wg_33.txt
echo  80 30 0.4 >> inp_wg_33.txt
echo  90 10 0.2 >> inp_wg_33.txt
rem
set inp_dat=inp_wg_33.txt
set fig_out=fig_wg_33.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig_out%
psxy %inp_dat% -R -JX -W2 -K -O >> %fig_out%
psxy %inp_dat% -R -JX -SC -G255 -W3 -N -O >> %fig_out%
rem ***************************************************************************
rem * プロットの大きさを変える                                                *
rem * -SC -G255 -W3 入力された数値群の位置に(C,G,W)白抜き円を描く．           *
rem *               (C)の次にプロットの大きさは指定しない．                   *
rem *               プロットの大きさは入力ファイルの3列目(z値)に対応．        *
rem *               この事例ではz値に応じて直径の異なる円がプロットされる．   *
rem ***************************************************************************



rem ***************************************************************************
rem * 入力用ファイル「inp_wg_34.txt」の作成                                      *
rem * (x座標 y座標 z値)：z値は色区分に対応                                    *
rem ***************************************************************************
echo  10 20 20 > inp_wg_34.txt
echo  20 30 30 >> inp_wg_34.txt
echo  30 50 50 >> inp_wg_34.txt
echo  40 80 80 >> inp_wg_34.txt
echo  50 90 90 >> inp_wg_34.txt
echo  60 70 70 >> inp_wg_34.txt
echo  70 60 60 >> inp_wg_34.txt
echo  80 30 30 >> inp_wg_34.txt
echo  90 10 10 >> inp_wg_34.txt
rem ***************************************************************************
rem * カラーパレットファイル「inp_wg_34cpt.cpt」の作成                           *
rem * (z1値 R1 G1 B1 z2値 R2 G2 B2)                                           *
rem ***************************************************************************
echo   0   0   0   0  20   0   0   0 > inp_wg_34cpt.cpt
echo  20 255   0   0  40 255   0   0 >> inp_wg_34cpt.cpt
echo  40   0 255   0  60   0 255   0 >> inp_wg_34cpt.cpt
echo  60   0   0 255  80   0   0 255 >> inp_wg_34cpt.cpt
echo  80 255   0 255 100 255   0 255 >> inp_wg_34cpt.cpt
rem
set inp_dat=inp_wg_34.txt
set fig_out=fig_wg_34.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig_out%
psxy %inp_dat% -R -JX -W2 -K -O >> %fig_out%
psxy %inp_dat% -R -JX -SC0.5 -Cinp_wg_34cpt.cpt -N -O >> %fig_out%
rem ***************************************************************************
rem * プロットの色を変える                                                    *
rem * -SC0.5          入力された数値群の位置に(C)0.5cmの円を描く              *
rem * -Cinp_wg_34cpt.cpt カラーパレットファイルinp_wg_34cpt.cptで指定された色       *
rem *                 で円を塗りつぶす．                                      *
rem *                 色の区分は入力ファイル3列目の値(z値)と                  *
rem *                 カラーパレットファイルでの定義に従い行われる            *
rem * RGBでの基本色表示                                                       *
rem *     Black      0   0   0                                                *
rem *     Red      255   0   0                                                *
rem *     Lime       0 255   0                                                *
rem *     Blue       0   0 255                                                *
rem *     Yellow   255 255   0                                                *
rem *     Fuchsia  255   0 255                                                *
rem *     Aqua       0 255 255                                                *
rem *     White    255 255 255                                                *
rem ***************************************************************************



rem ***************************************************************************
rem * 入力用ファイル「inp_wg_35.txt」の作成                                      *
rem * (x座標 y座標 za値 zb値)：za値は色区分に対応，zb値は大きさに対応         *
rem ***************************************************************************
echo  10 20 20 0.1 > inp_wg_35.txt
echo  20 30 30 0.3 >> inp_wg_35.txt
echo  30 50 50 0.5 >> inp_wg_35.txt
echo  40 80 80 0.7 >> inp_wg_35.txt
echo  50 90 90 1.0 >> inp_wg_35.txt
echo  60 70 70 0.8 >> inp_wg_35.txt
echo  70 60 60 0.6 >> inp_wg_35.txt
echo  80 30 30 0.4 >> inp_wg_35.txt
echo  90 10 10 0.2 >> inp_wg_35.txt
rem
set inp_dat=inp_wg_35.txt
set fig_out=fig_wg_35.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig_out%
psxy %inp_dat% -R -JX -W2 -K -O >> %fig_out%
psxy %inp_dat% -R -JX -SC -Cinp_wg_34cpt.cpt -N -O >> %fig_out%
rem ***************************************************************************
rem * プロットの色と大きさを変える                                            *
rem *     色の区分は入力ファイル3列目の値(za値)による                         *
rem *     プロットの大きさは入力ファイル4列目の値(zb値)による                 *
rem * -SC             入力された数値群の位置に(C)円を描く                     *
rem * -Cinp_wg_34cpt.cpt カラーパレットファイルinp_wg_34cpt.cptで指定された色       *
rem *                 で円を塗りつぶす．                                      *
rem ***************************************************************************


rem ***************************************************************************
rem * 変数のクリア                                                            *
rem ***************************************************************************
set range=
set scale=
set xga=
set yga=
set xlabel=
set ylabel=
set fig_out=
set inp_dat=

del .*