シェルはたまにしか使わないので、すぐ忘れてしまう。ので、ここにメモすることにする。[Unix|Linux]コマンドの使い方は[[linux_command]]へ．

====== スクリプトを保存する時，自動的に chmod +x を行う ======
.emacs に以下を記述するとファイルの先頭に#!が含まれているとき、自動的に chmod +x を行ってくれる。
<code emacs>
(add-hook 'after-save-hook
          'executable-make-buffer-file-executable-if-script-p)
</code>



====== シェルからgnuplotで図を作成(epsやpsなどにする) =====
たくさん図を作るときなどに。
<code bash>
make && ./a.out
gnuplot <<EOF
set term post eps
set out "hoge.eps"
plot 'fort.88'
EOF
</code>

====== シェルからgnuplotで図を作成(Xで表示) =====
=====決定版？=====
次のようにやると、./do.shとやると数値解と解析解を比べてくれる。
<code bash>
#do.shの中身
#!/bin/sh
make && ./a.out
gnuplot -persist <<EOF
plot 'suuchikai.dat','kaisekikai.dat'
EOF
</code>

=====↑に到達するまでのやりとり=====
以下のdo.shとplot.gpを用意すると、./do.shとやるだけで、数値解と解析解を比べてくれる。

<code bash>
#do.shの中身
make && ./a.out
./plot.gp
</code>
<code gnuplot>
#plot.gpの中身
#!/thin/local/bin/gnuplot -persist
plot 'suuchikai.dat','kaisekikai.dat'
</code>

-persistとやると、Xの画面が勝手に消えなくなる。
#!/thin/local/bin/gnuplotの部分は各自のgnuplotが置いてある場所を指定すること。
<code bash>
which gnuplot
</code>
とやって出てきたものを書き込めばよろしい。
バグ取り時に大活躍(のはず)。
ただし、スクリプトが二つのファイルになってしまう。
もっとsmartなやり方があるのかもしれない。FIXME


gnuplotのパスに依存しないやり方はこんな感じかしら．
<code bash>
#!/bin/sh
make && ./a.out
gnuplot -persist plot.gp
</code>

====== readする値をshellで食わす =====
  * fortranなんかのread文に値を食わせるときに使うshell．
  * ファイルを作ってリダイレクトするのが面倒なときに．

  * 変数１つの場合
<code bash>
echo 5 | ./a.out
</code>
  * 変数が2つ以上
<code bash>
echo "5\n6" | ./a.out
</code>
あるいは
<code bash>
./a.out << EOF
5
6
EOF
</code>

====== 文字を置換する(sed) ======
inputの中のすべてのappleをorangeに置換して、outputに放り込む。
<code bash>
sed -e 's/apple/orange/g' input > output
</code>

====== シェルで実数計算 ======
bcと、scaleを使う。scaleには欲しい桁数を入力する。
`は'ではないので注意が必要．
<code bash>
sq3=`echo "scale=17; sqrt(3)" |bc`
a=`echo "scale=17; $sq3/10.0" |bc`
</code>

bcについて補足． 
-lオプションで20桁まで計算＋関数の使用が可能になる．
<code bash>
echo "e(1)" | bc -l
</code>
e(x)はexp(1)のこと．

====== ファイルから何かを読み取って変数に代入 ======
hoge.datが一行からなる情報を持っているとき。
<code bash>
a=`cat hoge.dat`
</code>

複数行からなるpiyo.datから一行ごとに読み取りたければreadを使う。
<code>
cat piyo.dat | while read i
do
a=`echo $i`
(aを使ってなにか作業)
done
</code>

====== data file から必要な列だけを取り出す ======
awkを使う。n列目とm列目が欲しいなら、
<code bash>
awk '{print $n,$m}' hoge.dat
</code>


====== 0.1刻みで変数iを増加させながら何か実行 ======
<code fortran>
do i=1,100
   di=1.d-1*dble(i)
   write(*,*) di
   何か作業
enddo
</code>
のようなことをやりたければ、
<code bash>
i=1
while [ `echo "$i <= 10" |bc`  -eq 1 ]
  do
  echo $i
  何かコマンド
  i=`echo "$i + 0.1" |bc`
done
</code>
とすればよい。

====== 1刻みで変数aを増加させながら何か実行 ======
aを0から99まで動かすとき、
<code bash>
a=0
while test $a -ne 100
do 
  echo $a
　何か実行
  a=`expr $a + 1`
done
</code>

====== 変数を適当に変えてケーススタディー ======
周波数omegaを適当に振って計算した例。
<code bash>
make
for omega in 3.125 3.135 3.145 3.155 6.265 6.275 6.285 6.295
  do
  rm -fr output_linc_beta0_w${omega}  
  echo ${omega} | ./a.out
  cp -r output output_linc_beta0_w${omega}
done
</code>

====== 変数を端末から読み込みたい ======
readを使う。
<code bash>
echo "what's your name"
read name
echo your name is $name.
</code>

====== 3分待ってからプログラムを実行 ======
sleepを使う。
<code bash>
sleep 180 && ./a.out
</code>