fortran [Applied Mechanics Lab. WIKI]

この文書は読取専用です。文書のソースを閲覧することは可能ですが、変更はできません。もし変更したい場合は管理者に連絡してください。

====== Fortranのtips ======

===== 連番のファイル名をつける。 =====
ファイル名をcharacter変数に代入して使う。ファイル名の変更にはwrite文を使う。
逐次計算に使える。
<code fortran>
character(len=12) filename
...
do i=1,100
  ...
  write(filename,'(a,i4.4,a)') 'hoge',i,'.dat'
  open(10,file=filename)
  ...
    write(10,*) ...
  ...
  close(10)
end do
</code>

このファイルをperlに食わせて処理すると便利かも．
以下，サンプル．

<code perl>
#!/usr/bin/perl
$gpfile .= sprintf("%s\n","set terminal gif");
$gpfile .= sprintf("%s\n","set pm3d map");
for($i=1;$i < 30;$i++){
    $file = sprintf("hoge%04d",$i*10);
    $gpfile .= sprintf("%s\n","set output \"$file.gif\"");
    $gpfile .= sprintf("%s\n","splot \"$file.dat\" using 1:2:5");
};
open (GP, "| gnuplot" ) or die "no gnuplot";
print GP "$gpfile";
close GP;
system("gifsicle --delay 50 hoge*gif > anime.gif");
</code>

gifsicle はgifアニメを作成するツール．アニメを見るにはgifview．
<code bash>
gifview -a anim.gif
</code>