24小时热门版块排行榜    

查看: 2063  |  回复: 8

carrolyb

新虫 (初入文坛)

[求助] fortran如何获取文件中指定列的数据 已有2人参与

我想读取给定第二列的数据(文件数据如下):
  4.902380       24.09600       19.19362   
   5.117340       24.64800       19.53066   
   5.376030       24.58470       19.20867   
   1.905530       7.827580       5.922050   
   1.605610       8.201700       6.596090   
   1.946720       7.271230       5.324510   
   1.742370       10.20920       8.466830   
   1.901910       9.991350       8.089440   
   1.525720       10.93020       9.404480   
   1.705040       12.12240       10.41736   
   1.564290       11.48820       9.923910   
   1.638400       12.97950       11.34110   
   1.404010       20.71810       19.31409   
   1.571210       21.79690       20.22569   
   1.633620       23.00130       21.36768   
   1.478230       11.46750       9.989270   
   1.470280       14.52380       13.05352   
   1.670590       15.37190       13.70131   
编写的程序为:
program main
implicit none
integer i
real ::a(10,3),b(10)
open(1,file='C:\Users\user\Desktop\11\zongie.txt')
do i=1,10
read(1,*) a(i,2)
b(i)=a(i,2)
enddo

close(1)
write(*,*) b
end
但是运行后总是输出的是第一列的数据,不知道程序哪里有问题,请高手帮忙看看,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浔荆

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
read(1,*) a(i,2)
这个代码会读取编号1文件中某行的所有数据,然后逐个赋值给后面的变量列表。但由于变量列表只有一个变量,所以只会将该行数据中的第一个赋值给变量,也就是你的第一列数据。
改成这样应该就行了,其他不用变:
read(1,*) a(i,

另:文件编号不要用个位数,可能会与系统默认编号冲突
2楼2015-10-26 11:45:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浔荆

木虫 (正式写手)

【答案】应助回帖

read(1,*) a(i,1:3)
上面的冒号被系统吃了。。。
3楼2015-10-26 11:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

www.cctv

银虫 (小有名气)

你这个有规律,可以采用代码。

发自小木虫Android客户端
有幸相遇 不负相见
4楼2015-10-26 19:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liupw789k

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
Do j=1,3
   Do i=1,10
      Read(1,*)a(i,j)
   Enddo
Enddo
B(=a(:,2)
Fortran是按列读取数据的!
手机打字真麻烦!!

[ 发自小木虫客户端 ]
Purescienceissomethingworththinkingabout!
5楼2015-10-27 00:19:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liupw789k

金虫 (小有名气)

引用回帖:
5楼: Originally posted by liupw789k at 2015-10-27 00:19:20
Do j=1,3
   Do i=1,10
      Read(1,*)a(i,j)
   Enddo
Enddo
B(=a(:,2)
Fortran是按列读取数据的!
手机打字真麻烦!!

2楼是对的!我这个输出了好使!读数据我就没试过了!!嘿嘿~

[ 发自小木虫客户端 ]
Purescienceissomethingworththinkingabout!
6楼2015-10-27 00:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

庐山活水

铜虫 (小有名气)

自信人生二百年,会当击水三千里
7楼2015-10-27 00:25:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carrolyb

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 浔荆 at 2015-10-26 11:46:06
read(1,*) a(i,1:3)
上面的冒号被系统吃了。。。

明白了,谢谢
8楼2015-10-29 18:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carrolyb

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 浔荆 at 2015-10-26 11:46:06
read(1,*) a(i,1:3)
上面的冒号被系统吃了。。。

9楼2015-10-29 18:29:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 carrolyb 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 F4zcb7xRCEJA 2026-08-11 3/150 2026-08-11 21:51 by pMlf9MbdMHyc
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 27/1350 2026-08-11 21:12 by klxdzh
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 yPv6ndDtiRxQ 2026-08-11 4/200 2026-08-11 21:03 by thgx1vqMqrLw
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 8/400 2026-08-11 19:56 by thgx1vqMqrLw
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GpvwGNVshCOp 2026-08-10 7/350 2026-08-11 19:39 by thgx1vqMqrLw
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
请填处理意见