24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1891  |  回复: 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 popular289 2026-05-18 6/300 2026-05-18 15:17 by gy116024
[基金申请] 国自然上会要求 +3 无名者登山 2026-05-18 5/250 2026-05-18 14:16 by bnullh
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 cjf4bx70cj 2026-05-14 7/350 2026-05-17 18:49 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考博] 现在不知道怎么办,感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 售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 6/300 2026-05-16 11:29 by h3oerqvkv9
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见