回答編集履歴
1
誤字修正 プログラム改
answer
CHANGED
|
@@ -2,7 +2,7 @@
|
|
|
2
2
|
|
|
3
3
|
csv は list-directed (* format のこと）で読み込めばいい。ただし、色々落とし穴があるので要注意だが。（ブービートラップだらけのべトコンの村を捜索するつもりで行けば問題ない。）
|
|
4
4
|
|
|
5
|
-
とりあえずワーク配列に一
|
|
5
|
+
とりあえずワーク配列に一気にに読み込んで、その後で変形するのが楽。
|
|
6
6
|
Fortran は列優先なので、これは頭に叩き込んでおく必要がある。
|
|
7
7
|
wk に一気に読み込むと、入力ファイルの一列目は wk(1, :) に入り、二列目は wk(2,:) に入る。
|
|
8
8
|
|
|
@@ -22,10 +22,12 @@
|
|
|
22
22
|
real(8) :: wk(2, (nc+2)*(nd+2))
|
|
23
23
|
open(10, file = 'input.dat')
|
|
24
24
|
read(10, *) wk
|
|
25
|
+
|
|
26
|
+
a = reshape(wk(1, :), shape(a))
|
|
27
|
+
b = reshape(wk(2, :), shape(b))
|
|
28
|
+
! サイズを手動で与えてもいい a = reshape(wk(1, :), [nc+2, nd+2])
|
|
29
|
+
|
|
25
30
|
|
|
26
|
-
a = reshape(wk(1, :), [nc+2, nd+2])
|
|
27
|
-
b = reshape(wk(2, :), [nc+2, nd+2])
|
|
28
|
-
|
|
29
31
|
!check
|
|
30
32
|
block
|
|
31
33
|
integer :: i
|