Question

您好，遇到一个比较急的问题，关于Fortran90的直接文件读取，希望向您请教！谢谢！

csv数据格式如下：

a,b,year,month,day
0.4324756287533950 ,0.9,1980,1,1
0.3428334118681960 ,0,1980,1,2
0.7745397378446070 ,0,1980,1,3
0.8918107602794460 ,0,1980,1,4
。。。
包括表头共1828行，希望能从这个csv文件中直接读取某一行某一列的值赋给某个参数用于计算。请问应怎么做呢？
另外，为了测试，自己写了个程序来读取这个文件，但输出结果都是0，请问问题出在哪里？应如何解决呢？
代码如下
program test
implicit none
integer::i
real::a(1827)
open(1,file='aa.csv',access='direct',recl=18,action='readwrite',form='formatted')
read(1,*)
do i=1,3
write(*,*) a(i),b(i)
enddo
end program
谢谢您！
抱歉，程序real一行少打了。应为

real::a(1827),b(1827)
在下本意是用两个函数a(i)，b(i)直接获取文件中对应位置的a、b值用于计算。望不吝赐教！谢谢！

Answer 1

	program test
c	implicit none
	integer i
	real b(1827),c(1827),d(1827),e(1827)
	double precision a(1827)
c	open(1,file='aa.csv',access='direct',recl=99,form='formatted')
	open(1,file='aa.csv',status='old')
	read(1,*)
	do 5 i=1,3
	read(1,*)a(i),b(i),c(i),d(i),e(i)
5	continue
	close(1)
	do 10 i=1,3
	write(*,100)a(i),b(i),c(i),d(i),e(i)
10	continue
100	format(1x,1f18.15,1f10.7,3f10.0)
	end
	
	以上为fortran77代码。调试通过，供参考。
	
	program test
c	implicit none
	integer i
	real b(1827),c(1827),d(1827),e(1827)
	double precision a(1827)
	open(1,file='aa.csv',access='direct',recl=99,form='formatted')
c	open(1,file='aa.csv',status='old')
	read(1,*)
	do 5 i=1,3
	read(1,*)a(i),b(i),c(i),d(i),e(i)
5	continue
	close(1)
	do 10 i=1,3
	write(*,100)a(i),b(i),c(i),d(i),e(i)
10	continue
100	format(1x,1f18.15,1f10.7,3f10.0)
	end
	以上一样调试通过。

追答

这样的数据文件，还是以顺序文件处理为宜！

追问

您好！我用自由格式（只会这个）改写后编译ok，但跑不出。forrtl:severe(59):list-direct I/O syntax error。文件本身是excel转的，应该没问题的。若顺序文件，如何用函数子程序a(i)获取文件中对应位置的a呢？谢谢！

追答

要取第i行，先空读文件i次，再读取的就是！

追问

好的，谢谢您！还有一些小问题想请教一下：请问为什么只有a用双精度？您的程序是怎么写成这种自带行数和背景的样式的呢？

追答

这个𣎴是

自带的，是百度提供的功能！你选择某种语言提交成果就行，我选的是delphi。

追问

求教为何结果为0呢？字数限制，抱歉只能传图了。

追答

因为你的csv文件的第一个数位数较多，需要双精度

追问

谢谢，修改了下程序，但运行结果为0，还请指教>。<！

追答

为什么要用子程序呢？

追问

因为需要进行迭代计算，简单来说例如x(n)=f[P(n-1)]这样，用子程序比较方便。不知您能否帮改正一下这个程序和子程序呢？或者您不用子程序如何实现呢？

追答

program pn
integer i,n
real b(1827),c(1827),d(1827),e(1827),x
double precision a(1827)
open(1,file='et.csv',status='old')
read(1,*)
do 5 i=1,1827
read(1,*)a(i),b(i),c(i),d(i),e(i)
5 continue
close(1)

read(*,*)n
c n=1
x=10*b(n)
write(*,*)x
stop
end