24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1208  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

shen614

新虫 (初入文坛)

[求助] fortran 从文件读入问题 很基础

没有FORTRAN基础
想从文件里面读取数据

程序
WRITE(*,*) 'The Input Filename=?  '
          READ(*,999) FILE1
         
999     format(a12)
          OPEN(1,FILE=FILE1)
          read(1,*) npp,nppp
            read(1,*) (zong(i),i=1,2)
          read(1,*) (const(i),i=1,npp)
          read(1,*) (para(i),i=1,nppp)
          read(1,*) ((p(i,j),j=1,3),i=1,4)
          read(1,*) err


          ENDFILE 1

数据文件 a.txt
9,14
4.52,7.39
1,18,273,6000,1337,40500,0.25,0.1,0.2
3.9724,1.2322,-1.8977,-1.8977,0.5056,1.0387,-0.5280,-0.5280,2.02,-0.4962,0,0,298.15,10**-3
6.46,3.62,2.04,6.56,3.72,2.14,6.66,3.82,2.24,6.34,3.52,1.94
10**-6

但是老是出现 run-time error F6501:READ(a)-end of file encountered
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

你数据文件在执行程序在同一文件夹吗?
6楼2013-06-07 11:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
dbb627: 金币+2, 谢谢指导 2013-06-07 19:54:40
READ(*,999) FILE1
999     format(a12)
可以不要。 然后
OPEN(1,FILE=‘a.txt’)
          read(1,*) npp,nppp
            read(1,*) (zong(i),i=1,2)
          read(1,*) (const(i),i=1,npp)
          read(1,*) (para(i),i=1,nppp)
          read(1,*) ((p(i,j),j=1,3),i=1,4)
          read(1,*) err
close(1)
如果数据个数对,就应该可以
2楼2013-06-07 10:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shen614

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by pippi6 at 2013-06-07 10:44:12
READ(*,999) FILE1
999     format(a12)
可以不要。 然后
OPEN(1,FILE=‘a.txt’)
          read(1,*) npp,nppp
            read(1,*) (zong(i),i=1,2)
          read(1,*) (const(i),i=1,npp)
          ...

还是不行啊
数据个数是对的啊,贴在上面了呀
和输入格式有关嘛?
而且有不同的数据的,不能每次改程序,会很烦的
3楼2013-06-07 10:53:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
3楼: Originally posted by shen614 at 2013-06-07 10:53:12
还是不行啊
数据个数是对的啊,贴在上面了呀
和输入格式有关嘛?
而且有不同的数据的,不能每次改程序,会很烦的...

10**-3 这个不能这么写。要写成 1e-3, 10**-6 为 1e-6
数据个数你要自己检查

» 本帖已获得的红花(最新10朵)

4楼2013-06-07 11:08:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见