24小时热门版块排行榜    

查看: 601  |  回复: 2

xk6891

至尊木虫 (著名写手)

[求助] fortran如何消去相同行,使其只显示其中一行

有下列这样的一组数,希望能将多余,相同内容的行删去(只用fortran,不玩shell),譬如(下面1,3,5行相同)
  C   7.273878   8.481164   4.787954
  H   6.812623   8.515214   5.770640
  C   7.273878   8.481164   4.787954
  H   7.588283   9.573033   4.554239
  C   7.273878   8.481164   4.787954
  H   7.102666   8.004816   3.811218
看了网上的这种解法
[code]
   character aa(1000000)*256,bb*256
    dimension idd(1000000)
    open(1,file='text.txt',status='old')
    icc=0
1    continue
    read(1,'(a)'.,end=2,err=1) bb*256
    icc=icc+1
    aa(icc)(1:256)=bb(1:256)
2   continue
    idd(1)=0
    do 10 i=1,icc-1
    do 20 j=i+1,icc
    idd(j)=0
    if(aa(i)(1:256).eq.aa(j)(1:256)) idd(j)=1
20  continue
10  continue
    do 30 i=1,icc
    if(idd(i).eq.1) goto 30
    write(*,'(A)') aa(i)
30  continue
    end
[code]
鄙人愚钝,有几处表示不理解
aa(1000000)*256表示什么意思
aa(icc)(1:256)=bb(1:256) 不理解
最关键的是我认为if(aa(i)(1:256).eq.aa(j)(1:256)) idd(j)=1这句话,是否在比较的时候会产生歧义,因为嵌套循环如果只取一个标号(j)作比较标志的话,肯定重复了,而且do 30 i=1,icc,因为已经有消去的行了,如果在用icc的循环输出,那肯定多了
还请大侠们指点一下,给个详细些的思路或者解法,晚辈拜谢。

[ Last edited by xk6891 on 2011-6-16 at 21:33 ]
回复此楼
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)

明白了,是我自己想的方向有错误,现在实验成功了,只是还没有想明白
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
2楼2011-06-16 21:45:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)

shell的话还真简单...
3楼2011-06-16 22:11:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xk6891 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见