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 .*
|