fortran 配列の読み込み

Question

質問

行数をカウントしたのち、以下のデータを読み込みたいのですがうまくいきません。なぜか最後の行だけ読み込まれます。
データをちゃんと読み込むにはどうしたらいいですか？
###読み込みたいデータ
1, 3.88000000e-003, -9.50325916e-019, 3.00000000e-002
2, 3.88000000e-003, -9.50325916e-019, 0.00000000e+000
3, -3.88000000e-003, 4.75162958e-019, 0.00000000e+000
4, -3.88000000e-003, 2.49220638e-016, 3.00000000e-002
5, 3.22000000e-003, -7.88672539e-019, 3.00000000e-002
6, 3.22000000e-003, -7.88672539e-019, 0.00000000e+000

ソースコード

program p4
  implicit none
  integer :: i, n
  real x, y, z, w

  open (17, file='p5.txt', status='old')
  n = 0
  read(17,'()')
  do
    read (17, *, end=100)
    n = n + 1
  end do
100 continue
  rewind (17)
  read(17,'()')

  print *, 'NumRec =', n

  do i = 1, n
    read (17, *) x, y, z, w
  end do
  close (17)

open(20,file='p4.csv',status='replace')

 do i = 1, n
  write (20,*) x,",", y,",", z,",", w
 end do

close(20)
end program p4

結果

6	3.22E-03	-7.89E-19	0
6	3.22E-03	-7.89E-19	0
6	3.22E-03	-7.89E-19	0
6	3.22E-03	-7.89E-19	0
6	3.22E-03	-7.89E-19	0

Answer 1

program p3
implicit none
integer n, i
real,dimension(70000,4) :: a　!十分に大きい行数を設定しました．

open (17, file='p3.txt', status='old')
n = 0

do
read (17, *, end=100)
n = n + 1
end do

100 continue
rewind (17)

print *, 'NumRec =', n
do i = 1, n
read (17, *) a(i,1), a(i,2), a(i,3), a(i,4)
end do
close (17)

open(20,file='p3.csv',status='replace')

do i = 1, n
write (20,*) a(i,1),',', a(i,2),',', a(i,3),',', a(i,4)
end do

close(20)
end program p3

と組んだところ読み込めましたが，
行の最初に空白があるとうまく読み込めません...

Comments

[Y.H.]

ここは質問に対する回答を書くところです。
別途新しい質問としましょう。

Answer 2

読み込まれてますが、毎回、前回読んだものを上書きしていますので、最後に読んだ物だけが残っています。

「配列」というものを聞いたことは無いでしょうか？全部保存するにはそれを使います。

Comments

[MYJ]

アドバイスありがとうございます。配列を意識して組んだところうまく読み込めました．
行の最初に空白があるデータを読み込む場合にエラーが出てしまうのですが，対策方法などありますか？

[otn]

書式 * であればエラーにならないと思いますが、
・どんなプログラムで
・どんなデータを読んだら
・どんなエラーになったのか
を書いてください。

Answer 3

x,y,z,wは配列ではなく単なる変数で、読み込みループで毎回上書きされます。
（読み込みループ終了時点ではx,y,z,wには最後に読み込んだ行の値が格納されています。）

p4.csvに書き出しているのはx,y,z,w（最後に読み込んだ行の値が格納）を読み込んだ行数分となります。

配列の読み込み

とタイトルにありますが、配列に読み込んでないです。
配列に読み込んで、書き出し時に配列の内容を書き出しましょう。

Comments

[MYJ]

アドバイスありがとうございます。配列を意識して組んだところうまく読み込めました．
行の最初に空白があるデータを読み込む場合にエラーが出てしまうのですが，対策方法などありますか？

[Y.H.]

それは、この質問
> 行数をカウントしたのち、以下のデータを読み込みたいのですがうまくいきません。なぜか最後の行だけ読み込まれます。

とは別の問題になるので別の質問とした方が良いでしょう。