昨日の続き

gnuplotの自動化
この間勉強したshellスクリプトを使うとグラフがばしばし作れる.
move.plt
set terminal png
set output "calcname.png"
set xlabel "position"
set ylabel "probability"
set yrange [0:0.2]
plot "p202.dat" every :::col::col u 2:3 title "T=col with lines"

move.sh
#!/bin/bash
max_line=`cat p202.dat | wc -l`;
max_time=`expr $max_line / 100`;
for *1
do
j=`expr $i + 10000`;
sed "s/name/$j/g" move.plt | sed "s/col/$i/g" | gnuplot;
done
ちょっと行数がずれるけどまあ良いや.
そしてanimate.自分のマシンでやると遅くなるのでnatureに送り込む.頑張れ〜

http://nature.tsukuba.ac.jp/~n0310255/calc2.gif
http://nature.tsukuba.ac.jp/~n0310255/calc3.gif
http://nature.tsukuba.ac.jp/~n0310255/calc4.gif
http://nature.tsukuba.ac.jp/~n0310255/calc5.gif
成果物

*1:i=0;i<${max_time};i++