Question

fortran中如何从文件中读取数据并赋值给二维数组

Answer 1

这是一个基本操作。如果你这么问，我建议你系统的阅读一本教科书。

以下是通用的格式：

Open( 12 , File = 文件名 )

Read( 12 , * ) 数组
Close( 12 )

但具体你的文件是什么样子，是文本文件，还是二进制文件？文本文件是成列的，还是成块的？数据是什么类型的？文本字符串？还是实型数字？整型？

因为这些，具体的读取方法也会有差异。

追问

1 0.179854005575 45.4653015137 69.3728027344
1 0.892018020153 45.9355010986 68.210899353
1 2.33839988708 46.7653007507 65.7986984253
.
..
.
.
1 18.1438007355 55.8218994141 39.4707984924
这是txt文件中的数据，不要第一列，把后边三列赋值给二维数组。

追答

Real a(3,行数) , b
integer i
Open(12,File=文件名)
Do i = 1 , 行数
read(12,*) b , a(:,i)
End Do
Close(12)

追问

高手,能不能再帮个忙呢.下边怎么输不出的d，能帮忙看看呢。
program main
implicit none
Real a(3,50640) , b
real x,y,z
integer i,j
Open(12,File='test1.xyz')
open(102,file='d')
Do i=1,50640
do j=1,50640
read(12,*) b , a(:,i)
if (i/=j)then
x=a(1,i)-a(1,j)
y=a(2,i)-a(2,j)
z=a(3,i)-a(3,j)
d=x**2+y**2+z**2
write(102,*)d
end if
end do
end do
stop
end

追答

1.为什么用双重循环？你的文件一共多少行？是 50640 行，还是 2564409600 （二十五亿）行？
2.d 这个变量没有定义，你确定你的编译器能正常编译？

追问

不好意思高手写错了。txt的数据中三列是原子的x,y,z坐标，然后我想计算每两原子间的距离。请问高手肿么办。x,y,z,分别是第一列，第二列，第三列互减。d=x^2+y^2+z^2;帮帮忙高手。

追答

你确定要计算每一个原子和原子间的距离么？这样有 25 亿个结果。如果一个数存 8 位，需要19GB硬盘来存储。

追问

对的,我是要计算任意两个原子之间的距离,高手能否给编一个通用的程序,跪求!!!!!!!!!!!

追答

是任意，还是所有？
这可是有差别的。

Program Baidu_Thlws
Implicit none
Real a(3,50640) , b , x , y , z
integer i , j
Open(12,File='test1.xyz')
open(102,file='d')
Do i = 1 , 50640
read(12,*) b , a(:,i)
End Do
Close(12)
Do i = 1 , 50640
Do j = 1 , 50640
if ( i == j ) Cycle
x=a(1,i)-a(1,j)
y=a(2,i)-a(2,j)
z=a(3,i)-a(3,j)
write(102,*) x**2+y**2+z**2
End Do
End Do
Close(102)
End Program Baidu_Thlws

另外，计算 25亿个结果，一般的计算机可能需要几十小时。请耐心等待。