DtransuにはTECPLOT用の変換ツール(trans_tec.f)がついています。
今日もこれを使って変換し、濃度分布のチェックをしていたのですが、どうも薄い濃度の表現ができてないように感じました。変換されたTXTデータをチェックすると、小数点以下4桁しか書き出されていません。fortranソースを見ても倍精度になっていましたので、おかしいなあと思いながら読んでました。
もしかしてDtransu側が4桁しか吐き出していないのか?そっちを見ないといけないのか?と不安になりかけた頃、後ろのほうでformatされていました。そりゃそうですよね。
200 format(f10.3,4(a,f10.3),a,f10.4,a,e12.4,a,e12.4,a,e12.4)
これを、8桁にしとけば今回は十分な精度になるので、変えました。
200 format(f10.3,4(a,f10.3),a,f10.8,a,e12.4,a,e12.4,a,e12.4)
指数表示のほうがよかったかな?
人の書いたソースって、読めないんですよね。ま、自分の書いたソースも後でわからくなるときありますが。めんどくさくてもコメント行を書かないといけないですね。
0 件のコメント:
コメントを投稿