| 查看: 869 | 回复: 24 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
【求助】新手求助,不知道文件中数据有多少行怎么半?
|
|||
|
新手求助!! 我要写一段程序,读取文件中的每一行数据,调用子程序进行计算,可是我并不知道这个文件有多少行,这段程序应该怎么写,请大家帮忙... 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),像我这样写可以吗?? |
» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有14人回复
三甲基碘化亚砜的氧化反应
已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
请问有评职称,把科研教学业绩算分排序的高校吗
已经有5人回复
2025冷门绝学什么时候出结果
已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
AI论文写作工具:是科研加速器还是学术作弊器?
已经有3人回复
论文投稿,期刊推荐
已经有4人回复
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复

|
谢谢你的回复,经过大家的指导我已经能计算出一个文件的行数了,在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 |

19楼2009-11-19 21:45:53
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 441.5小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
2楼2009-11-17 16:19:06

4楼2009-11-17 17:02:53
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 441.5小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
5楼2009-11-17 19:35:32













回复此楼