| 查看: 1045 | 回复: 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),像我这样写可以吗?? |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
|
谢谢你的回复,经过大家的指导我已经能计算出一个文件的行数了,在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 (幼儿园)
- 金币: 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
