Fortran 左揃えで桁数が違う場合の読み込み方

Question

＜実現したいこと＞
fortranで下の形式のデータを読み込みたいです．

0001 0.00000e+000 4/PART_4
0158 0.00000e+000 4/PART_4
1179 0.00000e+000 4/PART_4
1180 0.00000e+000 4/PART_4
1181 0.00000e+000 4/PART_4
1449 -2.32840e+000 5/PART_5
1450 -3.48224e+000 5/PART_5
1451 -3.19070e+000 5/PART_5
1452 -2.89924e+000 5/PART_5
1453 -3.58439e+000 5/PART_5
1580 0.00000e+000 5/PART_5
1681 0.00000e+000 5/PART_5

＜問題＞
真ん中の列が左揃えになっており，下の命令で読み込んだところ，
負号のところでエラーが出ます．

read(*,'(i4,4X,e13.5,a11)')

＜質問＞
左揃えのデータのデータ長が，負号によって１桁違う場合の読み込み方を教えていただけないでしょうか．
よろしくお願いします．

Answer 1

スマートなやり方は入力ファイルの桁数をそろえることです。
おそらくファイル作成時にそのような融通が利かないと思われるので、プログラム実行前に変換してしまいましょう。マイナス記号の手前の空白を1つ削ります。
sedが一番楽でしょう。

$ sed -e 's/ -/-/g' < sample.dat 
0001    0.00000e+000    4/PART_4
0158    0.00000e+000    4/PART_4
1179    0.00000e+000    4/PART_4
1180    0.00000e+000    4/PART_4
1181    0.00000e+000    4/PART_4
1449   -2.32840e+000    5/PART_5
1450   -3.48224e+000    5/PART_5
1451   -3.19070e+000    5/PART_5
1452   -2.89924e+000    5/PART_5
1453   -3.58439e+000    5/PART_5
1580    0.00000e+000    5/PART_5
1681    0.00000e+000    5/PART_5

これで、どの行も'(i4,3X,e13.5,a12)'のフォーマット指定で読めるようになります。
e指定子を読む場合、空白が混ざっていても大丈夫です。

入力ファイルをいじれない場合、マイナスがあるかどうかで書式指定子を変えるしかないでしょうか。

program hello
   character(len=33) linebuf
   integer a
   real b
   character(len=12) c

   do
     read (*, "(A)", end=999 ) linebuf
        if (linebuf(9:9) == "-") then
            read(linebuf,'(i4,4X,e13.5,a12)') a,b,c
        else
            read(linebuf,'(i4,4X,e12.5,a12)') a,b,c
        end if
        print *, a, b, c
   end do
   999 print *, "finish."
end program Hello

一旦linebufに入力データを保存し、9文字目がハイフンかどうかでe指定子の長さを変えています。
この方法だとほかのケースが混ざってきたら都度フォーマットを作らないといけないのでスマートじゃない。
私は前者の方法を強く推します。

Comments

[Sayuki]

非常に丁寧で分かりやすい説明、ありがとうございました。
このプログラムの使用者はwindows環境で使うそうなので、とりあえずは後者で対応させていただきます。
自分でも、もう少しフォーマットの指定の仕方など研究してみるつもりですので、よりよい方法を見つけた場合はそちらをベストアンサーにさせていただくかもしれませんが、ご了承ください。