应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】新手求助,不知道文件中数据有多少行怎么半?
251/3返回列表123下一页
查看: 1290  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

tpp001

金虫 (著名写手)

[交流] 【求助】新手求助,不知道文件中数据有多少行怎么半?

新手求助!!
我要写一段程序,读取文件中的每一行数据,调用子程序进行计算,可是我并不知道这个文件有多少行,这段程序应该怎么写,请大家帮忙...
program main
implicit none
integer i, nrow ! nrow是文件的行数
character*10 atomname
open(1,file='co2.gro',status='unknown')
do i=1, nrow !读取文件中的每一行数据,可是我并不知道这个nrow的具体数值.
read(1,*) atomname,x(i),y(i),z(i)
if (atomname .eq. 'cb') then
call sub(x(i),y(i),z(i))
请大家帮忙,另外我想读取每一行数据后,调用子程序进行运算,我想知道把文件中的每一行数值赋值给子程序sub(x,y,z),像我这样写可以吗??
回复此楼

» 猜你喜欢
迷茫在知识的海洋里，需要你的指导。thankyou
1楼 2009-11-17 16:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1,VIP+0):谢谢参与应助！ 11-17 16:55
可以这么写：
L=1
10  CONTINUE
READ(UNIT, *,END=120) A(L),B(L)
CALL SUB(A(L),B(L))
L=L+1
GOTO 10
120  CONTINUE
赞一下(7人) 回复此楼
2楼2009-11-17 16:19:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

读

至尊木虫 (文坛精英)

会冒泡之超级潜水员

文献杰出贡献
senlia(金币+0,VIP+0):偶尔用用goto 其实满好的 比用其它语句代替反而要明白些 11-17 21:18
少用GOTO
回复此楼
3楼2009-11-17 16:51:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by maomao1210 at 2009-11-17 16:19:
可以这么写：
L=1
10  CONTINUE
READ(UNIT, *,END=120) A(L),B(L)
CALL SUB(A(L),B(L))
L=L+1
GOTO 10
120  CONTINUE

我是新手 能不能详细的说一下,另外这里的END=120是什么意思,有点晕,看了这段
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
4楼2009-11-17 17:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by 读 at 2009-11-17 16:51:
少用GOTO

少用不是不用，不然咋不直接把这个功能去掉呢？哈哈。
阁下给个妙招呀。
赞一下(1人) 回复此楼
5楼2009-11-17 19:35:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

senlia

木虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+1,VIP+0):辛苦了！ 11-17 23:57
引用回帖:
Originally posted by tpp001 at 2009-11-17 17:02:

我是新手 能不能详细的说一下,另外这里的END=120是什么意思,有点晕,看了这段

如果不用 end=120 这样的参数设置 那么如果读文件读到数据没有了还往下读的话就会出错 但是如果用了 那就表示 如果读文件出错了 就跳到120指定的语句
赞一下(11人) 回复此楼
书到用时方恨少, 天涯何处无芳草.
6楼2009-11-17 20:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你这样不行的
并没有判断是否到达文件末尾。这样写和直接用死循环没有区别。
引用回帖:
Originally posted by maomao1210 at 2009-11-17 16:19:
可以这么写：
L=1
10  CONTINUE
READ(UNIT, *,END=120) A(L),B(L)
CALL SUB(A(L),B(L))
L=L+1
GOTO 10
120  CONTINUE

赞一下(1人) 回复此楼
7楼2009-11-18 13:03:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★ ★
nono2009(金币+2,VIP+0):谢谢应助！ 11-18 13:21
使用函数EOF。 eof(file_id)如果为真则说明到达文件末尾了，就不读了。
下面是编译器手册上的说明和范例:

EOF

Inquiry Intrinsic Function (Generic): Checks whether a file is at or beyond the end-of-file record.
Syntax

result = EOF (unit)

unit
       

(Input) Must be of type integer. It represents a unit specifier corresponding to an open file. It cannot be zero unless you have reconnected unit zero to a unit other than the screen or keyboard.
Results

The result type is default logical. The value of the result is .TRUE. if the file connected to unit is at or beyond the end-of-file record; otherwise, .FALSE..
Example

! Creates a file of random numbers, reads them back

REAL x, total

INTEGER count

OPEN (1, FILE = 'TEST.DAT')

DO I = 1, 20

CALL RANDOM_NUMBER(x)

WRITE (1, '(F6.3)') x * 100.0

END DO

CLOSE(1)

OPEN (1, FILE = 'TEST.DAT')

DO WHILE (.NOT. EOF(1))

count = count + 1

READ (1, *) value

total = total + value

END DO

100 IF ( count .GT. 0) THEN

WRITE (*,*) 'Average is: ', total / count

ELSE

WRITE (*,*) 'Input file is empty '

END IF

STOP

END
See Also
ENDFILE
BACKSPACE
REWIND
引用回帖:
Originally posted by tpp001 at 2009-11-17 16:08:
新手求助!!
我要写一段程序,读取文件中的每一行数据,调用子程序进行计算,可是我并不知道这个文件有多少行,这段程序应该怎么写,请大家帮忙...
program main
implicit none
integer i, nrow ! nrow是文件的行数 ...

赞一下(10人) 回复此楼
8楼2009-11-18 13:05:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★ ★
nono2009(金币+2,VIP+0):谢谢应助！ 11-18 13:21
program main
implicit none
integer i, nrow ! nrow是文件的行数
character*10 atomname
open(1,file='co2.gro',status='unknown')
i=1
do while (.not. eof(1))
read(1,*) atomname,x(i),y(i),z(i)
if (atomname .eq. 'cb') then
call sub(x(i),y(i),z(i))
i=i+1
end do


end program
没有测试过，自己试试吧。
引用回帖:
Originally posted by tpp001 at 2009-11-17 16:08:
新手求助!!
我要写一段程序,读取文件中的每一行数据,调用子程序进行计算,可是我并不知道这个文件有多少行,这段程序应该怎么写,请大家帮忙...
program main
implicit none
integer i, nrow ! nrow是文件的行数 ...

[ Last edited by tjyl on 2009-11-18 at 13:10 ]
赞一下(10人) 回复此楼
9楼2009-11-18 13:08:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by tjyl at 2009-11-18 13:03:
你这样不行的
并没有判断是否到达文件末尾。这样写和直接用死循环没有区别。


你试了嘛？呵呵，咋知道不行呢？
赞一下(1人) 回复此楼
10楼2009-11-18 13:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tpp001 的主题更新
251/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了，不知道这次命运如何 +8 雨打竹帘 2026-08-30 12/600 2026-09-03 06:24 by sincosx
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 6/300 2026-09-02 23:59 by rM1TE0WVDIIY
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +14 Kittylucky 2026-08-27 15/750 2026-09-02 19:35 by liumei0601
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗？ +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
信息提示
关闭
请填处理意见
关闭