| 查看: 1039 | 回复: 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),像我这样写可以吗?? |
回复此楼» 猜你喜欢
数一英一271专硕（085401）求调剂，可跨
已经有8人回复
-
295求调剂
已经有4人回复
-
0856求调剂
已经有5人回复
-
材料学硕333求调剂
已经有10人回复
-
299求调剂
已经有6人回复
-
一志愿郑州大学，080500学硕，总分317分求调剂
已经有6人回复
-
289求调剂
已经有12人回复
-
070305高分子化学与物理 304分求调剂
已经有6人回复
-
0703化学调剂，求导师收
已经有10人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有5人回复
tjyl
金虫 (正式写手)
- 程序强帖: 2
- 应助: 0 (幼儿园)
- 金币: 3218.1
- 红花: 2
- 帖子: 576
- 在线: 156.6小时
- 虫号: 765184
- 注册: 2009-05-07
- 专业: 生物无机化学
- 管辖: 程序语言
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+5,VIP+0):向你学习！ 11-19 23:14
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+5,VIP+0):向你学习！ 11-19 23:14
|
按照你的格式改写好了 用你发的数据测试了一下是可以的。 其实本来根本就不用去操心有多少行的，因为你只是为了计算总和。 直接按照前面的那样只判断到文件末没有就可以了。 这样改写只是为了随便试试而已。如果只是关心总和也不用去计算到底多少行。 使用了eof函数，intel编译器可以。CVF不大清楚，没有安装。 如果不行的话就使用前面说的read的IOSTAT来判断到行末没有。 ================== 这个是会计算行数，然后再分配空间的 ================= program read_data implicit none real*8 ljcc,vicc real*8,external :: cljinteraction real*8,allocatable :: x(: ),y(: ),z(: ) character*5,allocatable :: atomname(: ) integer*4 row,i row=0 open(32,file='co2ben.txt') do while (.not. eof(32)) row=row+1 read(32,*) end do rewind(32) allocate(x(row),y(row),z(row),atomname(row)) vicc=0.0 do i=1, row read(32,*) atomname(i),x(i),y(i),z(i) if (atomname(i) .eq. 'Cb') then !如果原子是Cb,就调用下面函数 ljcc=cljinteraction(x(i),y(i),z(i)) vicc=vicc+ljcc !累加函数数值 endif enddo write(6,*) vicc !输出函数数值的总和 deallocate(x,y,z,atomname) end program function cljinteraction(x,y,z) implicit none real*8 cljinteraction real*8 x,y,z !输入变量 real*8 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 function ＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ 这个就是只计算这个结果的版本 ＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ program read_data implicit none real*8 ljcc,vicc real*8,external :: cljinteraction real*8 x,y,z character*5 atomname vicc=0.0 open(32,file='co2ben.txt') do while (.not. eof(32)) read(32,*) atomname,x,y,z if (atomname .eq. 'Cb') then !如果原子是Cb,就调用下面函数 ljcc=cljinteraction(x,y,z) vicc=vicc+ljcc !累加函数数值 endif enddo write(6,*) vicc !输出函数数值的总和 end program function cljinteraction(x,y,z) implicit none real*8 cljinteraction real*8 x,y,z !输入变量 real*8 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 function [ Last edited by tjyl on 2009-11-19 at 22:37 ] |
24楼2009-11-19 22:25:34
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
2楼2009-11-17 16:19:06
4楼2009-11-17 17:02:53
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
5楼2009-11-17 19:35:32
