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

24小时热门版块排行榜    

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

tjyl

金虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2,VIP+0):谢谢应助！ 11-18 14:36
intel的编译器应该这样可以的，不过gfortran不行，可能eof返回的类型不对的原因。
另外可以使用iostat来判断。
open(15, file='foo.txt')
do
   read(fh, iostat=stat) buf
   if (stat /= 0) exit
   ! process buf
end do
close(15)

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

tjyl

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+1,VIP+0):谢谢评论！ 11-18 14:36
以前是看到有的人这样用，只是觉得这种方式不好。
引用回帖:
Originally posted by maomao1210 at 2009-11-18 13:19:

你试了嘛？呵呵，咋知道不行呢？

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

jjdg

版主 (知名作家)

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
反正程序这种东西不上机实测是不知道好用不好用的！
赞一下(1人) 回复此楼
努力学习！以正当途径！获得需要的知识！
13楼2009-11-18 23:30:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by tjyl at 2009-11-18 13:08:
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 ...

你好,我用compad visual fortran 运行了这段,可是还是出现错误..错误信息是:
D:\rdf\gro\benz\ljinteraction.f90(14) : Error: Constants and expressions are invalid in read-only I/O lists.   [X]
read(1,*) atomname,x(i),y(i),z(i)
-------------------^
D:\rdf\gro\benz\ljinteraction.f90(14) : Error: Constants and expressions are invalid in read-only I/O lists.   [Y]
read(1,*) atomname,x(i),y(i),z(i)
------------------------^
D:\rdf\gro\benz\ljinteraction.f90(14) : Error: Constants and expressions are invalid in read-only I/O lists.   [Z]
read(1,*) atomname,x(i),y(i),z(i)
我不太清楚这个数组x(i),y(i),z(i),怎么定义好,我在program main的开头这样定义的
program main
real x, y, z
我觉得是不是这里有问题呀...就是这个数组到底应该怎么半呢?
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
14楼2009-11-19 15:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (知名作家)
本帖仅楼主可见
15楼2009-11-19 16:07:57
已阅   申请程序强帖   回复此楼   编辑   查看我的主页

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

comments on L14

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1,VIP+0):辛苦了！ 11-19 18:12
real x,y,z 只是定义了实数而非数组。
子程序里x,y,z是数组吗？数组是必须的吗？
也许可以用动态数组，但我长时间不用了，基本都忘了。看看CVF的user manual吧。
赞一下(11人) 回复此楼
16楼2009-11-19 16:16:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
子程序里面的x,y,z,就是一个变量......
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
17楼2009-11-19 16:49:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1,VIP+0):辛苦了！ 11-19 21:52
不好意思，没有注意到数组x,y,z还没有申明。
可以用一个链表来做。
或者笨一点的方法是用先那样读一次，计算出行号i。
然后再allocate(x(i),y(i),z(i))
引用回帖:
Originally posted by tpp001 at 2009-11-19 15:58:

你好,我用compad visual fortran 运行了这段,可是还是出现错误..错误信息是:
D:\rdf\gro\benz\ljinteraction.f90(14) : Error: Constants and expressions are invalid in read-only I/O lists.   [X]
read(1 ...

赞一下(11人) 回复此楼
18楼2009-11-19 21:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
引用回帖:
Originally posted by tjyl at 2009-11-19 21:26:
不好意思，没有注意到数组x,y,z还没有申明。
可以用一个链表来做。
或者笨一点的方法是用先那样读一次，计算出行号i。
然后再allocate(x(i),y(i),z(i))


谢谢你的回复,经过大家的指导我已经能计算出一个文件的行数了,在linux用cat -n 就可以.
可是我的程序又有其他问题了.望帮忙谢谢
我下面这个程序 compile没有错误,可是execute就不行,
我程序的主题就是,调用一个函数,计算函数值,然后累加函数数值.,输出函数数值的总和!
program main
implicit none
integer i
double precision,external :: cljinteraction !申明函数
character*5 atomname
double precision x(33972),y(33972),z(33972),line(33972)
double precision ljcc !函数输出数值
double precision vicc
    vicc=0.0
     open(1,file='co2ben.txt',status='unknown')
    do i=1, 33972
   read(1,*) line(i),atomname,x(i),y(i),z(i)
if (atomname .eq. 'Cb') then !如果原子是Cb,就调用下面函数
    ljcc=cljinteraction(x(i),y(i),z(i))
     vicc=vicc+ljcc !累加函数数值
endif
enddo
write(6,*) vicc !输出函数数值的总和
end

function cljinteraction(x,y,z)
implicit none
double precision cljinteraction
double precision x,y,z !输入变量
double precision ri2c, sr2c, sr6c, sr12c !中间使用量
    ri2c=(x-3.08)**2+(y-3.37)**2+(z-3.129)**2
   sr2c=0.32535/ri2c
   sr6c=sr2c*sr2c*sr2c
   sr12c=sr6c**2
   cljinteraction=sr12c-sr6c
return
end
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
19楼2009-11-19 21:45:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)

这是我的输入文件 'co2ben.txt'

这是我的输入文件 'co2ben.txt' 部分数据
Oa   3.106000       3.385000       3.639000   
Cb   3.119000       3.271000       3.627000   
Oc   4.060000      -2.784000       1.788000   
Oa   2.935000       3.083000       3.510000   
Cb   2.845000       3.101000       3.577000   
Oc   2.701000       3.042000       3.654000   
Oa   3.074000       3.273000       3.460000   
Cb   3.038000       3.260000       3.569000   
Oc   2.751000       3.116000       3.642000   
Oa   3.347000       3.438000       3.404000   
Cb   3.322000       3.381000       3.501000   
Oc   3.004000       3.250000       3.678000   
Oa   3.433000       3.481000       3.513000   
Cb   3.347000       3.406000       3.515000   
Oc   3.295000       3.326000       3.598000   
Oa   3.288000       3.359000       3.679000   
Cb   3.254000       3.374000       3.571000   
Oc   3.261000       3.329000       3.513000
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
20楼2009-11-19 21:47:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tpp001 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +10 Micky11223 2026-03-25 14/700 2026-03-28 15:48 by Micky11223
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 299求调剂 +7 嗯嗯嗯嗯2 2026-03-27 7/350 2026-03-28 13:09 by 唐沐儿
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 292求调剂 +14 鹅鹅鹅额额额额� 2026-03-25 15/750 2026-03-28 08:45 by WYUMater
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +8 曼111 2026-03-24 9/450 2026-03-28 07:58 by YYYYX1234
[考研] 272求调剂 +7 脚滑的守法公民 2026-03-27 7/350 2026-03-27 17:23 by laoshidan
[考研] 279 分 求调剂 +4 睡个好觉_16 2026-03-24 4/200 2026-03-27 15:05 by 醉在风里
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 一志愿陕师大生物学071000，298分，求调剂 +5 SYA！ 2026-03-23 5/250 2026-03-27 09:29 by 不吃魚的貓
[考研] 343求调剂 +4 赠我一本书 2026-03-23 4/200 2026-03-27 00:40 by wxiongid
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 340求调剂 +3 Amber00 2026-03-26 3/150 2026-03-26 18:57 by 不吃魚的貓
[考研] 297求调剂 +6 田洪有 2026-03-26 6/300 2026-03-26 15:55 by 不吃魚的貓
[考研] 材料与化工328分调剂 +6 。，。，。，。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
信息提示
关闭
请填处理意见
关闭