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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】新手求助,不知道文件中数据有多少行怎么半?
252/3返回列表上一页123下一页
查看: 1169  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-16 22:25 by lindoudou
[考博] 现在不知道怎么办，感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 19:10 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-16 18:45 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 3/150 2026-05-16 17:50 by j6b2pdz07o
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 3/150 2026-05-16 17:30 by j6b2pdz07o
[考博] 光量子物理方向 博士招生 1人（2026.09） +3 sandyworld 2026-05-15 3/150 2026-05-16 17:11 by zznnnj
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 11:36 by h3oerqvkv9
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 cjf4bx70cj 2026-05-14 6/300 2026-05-16 11:16 by h3oerqvkv9
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +11 dd921ww 2026-05-12 13/650 2026-05-16 09:16 by Howard28
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 04:29 by k37jurhrau
[考博] 西南大学考核制博士 +4 lijunjie84 2026-05-11 7/350 2026-05-15 23:20 by 同仁堂教主
[基金申请] 青C资助名额大幅增加！ +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
关闭
请填处理意见
关闭