スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

gnuplotのmultiplotのメモ

Windows版のgnuplot(wgnuplot)を使って、一つの図の中に複数のグラフを描画する際のmultiplotを利用するためのTIPについて書きます。
  • multiplotの中でterminalの変更はできない
  • replotは使えないのでplotですべてプロット。改行は,\
  • 横軸が共通のときはmarginで合わせる
  • xrangeをデフォルトへ戻すのはset xrange [*:*]
  • wgnuplot.exeのプロセスが残ってしまうなら unset multiplot の後に何でもいいからplot

001_2015101506141594f.png

Fig.1: AkaiKKRでγ-Mnの反強磁性の結果をmultiplotで一つのグラフにまとめたもの



gnuplotのmultiplot


gnuplotではmultiplotを利用して一枚の図の中に複数のグラフを描くことが可能です。この際にいくつか注意しなければならない点があります。今回は、それらのポイントを(最小限のものだけ)まとめました。multiplotの詳しい使い方は、米澤 進吾さんのホームページなどを参照されるのが良いと思います。

set terminal pngcairo enhanced color size 520,600
set output "fccMn.png"

## *** マルチプロットの開始 ***
set multiplot layout 2,1
set lmargin 8 # 左マージン
set rmargin 2 # 右マージン

## *** 磁気モーメントのプロット ***
set key left top # 凡例の位置
set yrange [0:4]
set xlabel "Lattice constant a (Bohr)"
set ylabel "Local spin moment ({/Symbol m}_B)"
plot 'fccMn-Moment.txt' u 1:4 w lp title "Non-mag",\
'fccMn-Moment.txt' u 1:3 w lp title "Ferro",\
'fccMn-Moment.txt' u 1:2 w lp title "Antiferro"

## *** 全エネルギーのプロット ***
set key right top
unset ytics
set yrange [*:*] # yの範囲をデフォルトに戻す
set xlabel "Lattice constant a (Bohr)"
set ylabel "Total Energy (arb. unit)"
plot 'fccMn-NMG.txt' u 1:2 w lp title "Non-mag",\
'fccMn-FMG.txt' u 1:2 w lp title "Ferro",\
'fccMn-AFM.txt' u 1:($2/2.0) w lp title "Antiferro"

## *** マルチプロットの終了 ***
unset multiplot

## *** 後処理 ***
# 何かプロットしておくとwgnuplot.exeのプロセスが残らない
set terminal windows
plot sin(x)


terminalの切り替え不可


multiplot環境下ではterminalの切り替えができません。
私はこれまでmultiplotを使わない状況では、windowsのterminalで画面上で描画を確認した後、terminalをpngcairo等に変更して画像出力をしていましたがmultiplotを使う場合は、最初からpngcairoやpostscriptなど画像出力用のterminalを指定して一発勝負にします。

一つのグラフにつき一回のplot


multiplot環境の中では、一回のplotに付きひとつのパネルが作成されます。
一つのパネルの中に複数の線を描きたい場合、一回のplotで行う必要があります。
plot 'fccMn-NMG.txt' u 1:2 w lp title "Non-mag",\
'fccMn-FMG.txt' u 1:2 w lp title "Ferro",\
'fccMn-AFM.txt' u 1:($2/2.0) w lp title "Antiferro"

この際に、改行を表す ,\ を利用します。

グラフの横位置を合わせる


複数のパネルが共通のx軸を持っている場合、位置がずれて見栄えが悪くなることがあります。
そういう場合は lmarginrmargin でマージンを設定します。
set lmargin 8
set rmargin 2


描画範囲のリセット


描画範囲を設定するときには set yrange [0:4] のように指定します。
これをデフォルトの設定に戻すには set yrange [*:*]のように指定します。
set xrange [*:*]
set yrange [*:*]

なんとなく unset を使いたくなりますが、こちらは間違いです。

wgnuplot.exeのプロセスが残ってしまう場合


私の環境(のうちのひとつ)では wgnuplot.exe のプロセスが残ってしまうことがあるようです。
multiplot環境を抜けた後、terminalをwindowsに切り替えて、何でも良いのでプロットするとプロセスが残りません。
set terminal windows
plot sin(x)

プロセスが残らないなら必要ないと思います。

関連エントリ




参考URL




付録


このエントリで使用したファイルを添付します。ファイル名末尾の".txt"を削除して、"_"を"."に変更すれば使えるはずです。(参考:ねがてぃぶろぐの付録)


フィードバック



にほんブログ村 その他趣味ブログ 電子工作へ

 ↑ 電子工作ブログランキング参加中です。1クリックお願いします。


コメント・トラックバックも歓迎です。 ↓      


 ↓ この記事が面白かった方は「拍手」をお願いします。
スポンサーサイト

tag: gnuplot multiplot AkaiKKR machikaneyama 

FC2カウンター
カテゴリ
ユーザータグ

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoC強磁性CPAPICOPアンプecalj状態密度モンテカルロ解析常微分方程式トランジスタodeDOSインターフェースPDS5022定電流スイッチング回路分散関係半導体シェルスクリプトレベルシフト乱数HP6632A可変抵抗トランジスタ技術R6452AI2C温度解析ブレッドボードバンドギャップ確率論反強磁性セミナーバンド構造数値積分偏微分方程式非線形方程式ソルバ熱設計絶縁三端子レギュレータISO-I2CA/DコンバータシュミットトリガフォトカプラカオスPWscfGW近似LM358LEDマフィンティン半径発振回路USB数値微分TL431PC817Cサーボアナログスイッチ直流動作点解析補間カレントミラー74HC4053bzqltyチョッパアンプFFT2ちゃんねる開発環境量子力学単振り子電子負荷VESTAQuantumESPRESSO標準ロジックパラメトリック解析ブラべ格子イジング模型アセンブラLDA基本並進ベクトルBSchSMPTLP621失敗談六方最密充填構造コバルト位相図QSGWGGAスイッチト・キャパシタewidth状態方程式VCAキュリー温度繰り返し最適化仮想結晶近似不規則合金熱伝導gfortran相対論抵抗FETMaximaQuantum_ESPRESSOcygwinランダムウォークラプラス方程式スピン軌道相互作用スレーターポーリング曲線マントルシュレディンガー方程式ZnO自動計測QNAP固有値問題ダイヤモンドデータロガー井戸型ポテンシャルTLP552CIFxcrysdenゼーベック係数熱力学条件分岐MCU最小値UPS格子比熱最大値ガイガー管平均場近似過渡解析Writer509スーパーセルFXA-7020ZR差し込みグラフ第一原理計算テスタ起電力OpenMP三角波ubuntuLM555NE555ブラウン運動詰め回路ハーフメタルawkfsolveUbuntuフェルミ面TLP521トランスMAS830LPGACK1026OPA2277フィルタトレーナーバトルEAGLEノコギリ波負帰還安定性ナイキスト線図MBEP-10LMC6622SC1815CapSenseAACircuitPIC16F785入出力固定スピンモーメントFSMTeX結晶磁気異方性全エネルギーc/a合金multiplotgnuplot非線型方程式ソルバL10構造正規分布等高線ジバニャン方程式ヒストグラム確率論初期値interp1fcc面心立方構造ウィグナーザイツ胞半金属デバイ模型磁気モーメント電荷密度重積分SIC不純物問題擬ポテンシャル状態図cif2cellPWgui二相共存ウルツ鉱構造edeltquantumESPRESSOフォノンリジッドバンド模型スワップ領域BaO岩塩構造ルチル構造マテリアルデザインspecx.fフラクタルマンデルブロ集合キーボードRealforceクーロン散乱三次元疎行列縮退化学反応関数フィッティング最小二乗法Excel直流解析PCTS-110TS-112日本語パラメータ・モデル等価回路モデル文字列ポケモンGO熱拡散方程式HiLAPW両対数グラフ片対数グラフ陰解法Crank-Nicolson法ifort境界条件連立一次方程式グラフの分割軸ラベルヒストグラム不規則局所モーメントスーパーリーグ円周率Gimp凡例線種シンボルトラックボール

最新コメント
リンク

にほんブログ村 その他趣味ブログ 電子工作へ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。