WANtaro雑記帳（Win_GMT：片対数グラフにプロットする）

rem ***************************************************************************
rem * 片対数グラフにプロットする                                              *
rem ***************************************************************************

rem ***************************************************************************
rem * 入力用ファイル「inp_wg_11.txt」の作成                                   *
rem ***************************************************************************
echo   0.01  -8 > inp_wg_11.txt
echo   0.10  -4 >> inp_wg_11.txt
echo   1      0 >> inp_wg_11.txt
echo   5      8 >> inp_wg_11.txt
echo  10      7 >> inp_wg_11.txt
echo  20      6 >> inp_wg_11.txt
echo  50      5 >> inp_wg_11.txt
echo  75      4 >> inp_wg_11.txt
echo 100      3 >> inp_wg_11.txt

rem ***************************************************************************
rem * 入力用ファイル「inp_wg_12.txt」の作成                                   *
rem ***************************************************************************
echo   0.01  3 > inp_wg_12.txt
echo   0.10  4 >> inp_wg_12.txt
echo   1     5 >> inp_wg_12.txt
echo   5     6 >> inp_wg_12.txt
echo  10     7 >> inp_wg_12.txt
echo  20     8 >> inp_wg_12.txt
echo  50     0 >> inp_wg_12.txt
echo  75    -4 >> inp_wg_12.txt
echo 100    -8 >> inp_wg_12.txt


del .*

set range=0.01/100/-10/10
set scale=8l/4
set xga=g3a1
set yga=g2a4
set xlabel="Name of x-axis"
set ylabel="Name of y-axis"
rem
set inp_dat=inp_wg_11.txt
rem
gmtset ANOT_FONT_SIZE 12
gmtset LABEL_FONT_SIZE 8
gmtset TICK_LENGTH 0c
rem
set fig_out=fig_wg_11.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 -O >> %fig_out%
rem ***************************************************************************
rem * psbasemapコマンドで作成した枠線内に，psxyコマンドを用いてデータプロット *
rem * を重ね書きする．                                                        *
rem * 入力データファイルはinp_datに格納された「inp_wg_11.txt」．              *
rem * 出力ファイルはfig_outに格納された「fig_wg_11.eps」                      *
rem * -W2    入力された数値群を2dpiの線で結ぶ                                 *
rem * -SC0.2 入力された数値群の位置に(C)0.2cmの円を描く                       *
rem * -K,-O,>,>>について                                                      *
rem *     バッチコマンドの1行目には-Kをいれ，>で出力ファイルを指定．          *
rem *     バッチコマンドの2行目以降には-K，-Oをいれ，>>で出力ファイルを指定． *
rem *     バッチコマンドの最終行には-Oをいれ，>>で出力ファイルを指定．        *
rem * 記号の線が細く，記号の下に折れ線が見えていて，カッコいいグラフでは      *
rem * ありません．                                                            *
rem * またx=0.01とx=100上の丸が半分しか描画されていません．                   *
rem ***************************************************************************
rem
set fig_out=fig_wg_12.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 -N -O >> %fig_out%
rem ***************************************************************************
rem * -SC0.2 -G0 入力された数値群の位置に(C)0.2cmの円を描き(G)黒で塗りつぶす．*
rem * -N         -Nオプションによりグラフ枠の外にも描画可能にする．           *
rem *            これによりx=0.01とx=100上の黒丸が完全に描画されます．        *
rem ***************************************************************************



set inp_dat1=inp_wg_11.txt
set inp_dat2=inp_wg_12.txt
set fig_out=fig_wg_13.eps
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig_out%
psxy %inp_dat1% -R -JX -W3 -K -O >> %fig_out%
psxy %inp_dat1% -R -JX -SC0.2 -G0 -N -K -O >> %fig_out%
psxy %inp_dat2% -R -JX -W3t10_10:0 -K -O >> %fig_out%
psxy %inp_dat2% -R -JX -SS0.2 -G255/255/255 -W3 -N -K -O >> %fig_out%
echo 0.02 0.95  8 0 0 TL %inp_dat1% | pstext -R0/1/0/1 -JX8/4 -N -K -O >> %fig_out%
echo 0.02 0.85  8 0 0 TL %inp_dat2% | pstext -R0/1/0/1 -JX8/4 -N -K -O >> %fig_out%
echo 0.99 1.05 10 0 5 BR %fig_out%  | pstext -R0/1/0/1 -JX8/4 -N -O >> %fig_out%
rem ***************************************************************************
rem * 入力データファイルとして，inp_dat1に「inp_wg_11.txt」，                 *
rem *                           inp_dat2に「inp_wg_12.txt」を格納．           *
rem * 出力ファイルとして，fig_outに「fig_test13.eps」を格納．                 *
rem * psxy %inp_dat1% -R -JX -W3 -K -O >> %fig_out%                           *
rem *     「inp_wg_11.txt」より数値群を読み込み，(-W3)3dpiの実線で結ぶ．      *
rem * psxy %inp_dat1% -R -JX -SC0.2 -G0 -N -K -O >> %fig_out%                 *
rem *     「inp_wg_11.txt」より数値群を読み込み，(-SC0.2 -G0)黒丸を表示する． *
rem * psxy %inp_dat2% -R -JX -W3t10_10:0 -K -O >> %fig_out%                   *
rem *     「inp_wg_12.txt」より数値群を読み込み，                             *
rem *       (-W3t10_10:0)3dpiの点線で結ぶ．                                   *
rem * psxy %inp_dat2% -R -JX -SS0.2 -G255/255/255 -W3 -N -K -O >> %fig_out%   *
rem *     「inp_wg_12.txt」より数値群を読み込み，                             *
rem *       (-SS0.2 -G255/255/255 -W3)白抜きの正方形を表示する．              *
rem *      -SS0.2        (S)0.2cmの正方形表示                                 *
rem *      -G255/255/255 (G)正方形を白で塗りつぶし                            *
rem *      -W3           (W)正方形の外形を3dpiの実線で囲む                    *
rem * テキストの表示                                                          *
rem * pstext テキスト表示コマンド                                             *
rem * pstextは以下の7個のデータ組を読み取りテキストを描画します．             *
rem * x座標 y座標 フォントサイズ 回転角(度) フォントNo 位置関係 テキスト      *
rem * 位置関係は以下のアルファベットの組み合わせを指定します．                *
rem *     TL(top-left)    TC(top-center)    TR(top-right)                     *
rem *     ML(middle-left) MC(middle-center) MR(middle-right)                  *
rem *     BL(bottom-left) BC(bottom-center) BR(bottom-right)                  *
rem * echo 0.02 0.95  8 0 0 TL %inp_dat1%                                     *
rem *     コマンドプロンプト画面でpstextコマンドに送る7個のパラメータを       *
rem *     指定します．                                                        *
rem *     TLは座標(0.02,0.95)に文字列の上端(top)と左端(left)を一致させる指定．*
rem *     inp_dat1は文字列「inp_text1.txt」を示す．                           *
rem * pstext -R0/1/0/1 -JX8/4 -N -K -O >> %fig_out%                           *
rem *     -R0/1/0/1  描画範囲をx軸を0から1，y軸を0から1に再定義．             *
rem *     -JX8/4     x軸長さを普通軸8cm，y軸長さを普通軸4cmに再定義．         *
rem *     上記再定義された軸の座標(0.02,0.95)を基準に文字列を描画する．       *
rem ***************************************************************************

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

del .*