应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】如何实现数字紧贴字母输入
51/1返回列表
查看: 2465  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xk6891

至尊木虫 (著名写手)

[交流] 【求助】如何实现数字紧贴字母输入

H  1
H  2
H  3
H  4
H  5
H  6
H  7
H  8
H  9
H 10
H 11
H 12
请教前辈，如何让后面的数字紧贴着字母输入呢后面还有几百行
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-01-12 00:16:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
★ ★
jjdg(金币+2): 感谢参与 2011-06-17 22:09:03
引用回帖:
Originally posted by snoopyzhao at 2011-01-13 19:59:40:


我觉得 90 的话，入门可以看 “fortran 90 编程指南”，袁驷，叶康生译，高等教育出版社出版的。但有些术语可能翻译的比较别扭……

好的，找来看一看

fortran对于格式的要求很严格，刚开始读取数据的时候出现的是
At line 10 of file tiH.f90 (unit = 17, file = 'wa_xyz.dat')
Fortran runtime error: End of file
这个错误修改之后，对于读取和写入数据采取不同的格式之后，看出了问题所在，下面是
[code:]
program main
implicit none
integer i,j
character B
real A
dimension B(768),A(768,3)

open(17,file='wa_xyz',status='old')
do i=1,768
read (17,*) B(i),A(i,1),A(i,2),A(i,3)
end do
close(17)

open(18,file='nem_wa_xyz',status='unknown')
do i=1,768
if(mod(i,3)/=0)then
write(18,*)(B(i),A(i,j),j=1,3)
end if
end do
close(18)

stop
end program
[code:]
读取这样一段数据：
O         3.668000       10.082000       15.904000
H         3.224000       10.451000       15.101000
H         3.092000       10.379000       16.627001
O         5.186000       16.695999       12.072000
H         6.083000       16.979000       11.816000
H         5.337000       15.832000       12.490000
O         3.786000       16.318001        1.310000
H         3.781000       15.388000        1.196000
H         4.528000       16.580000        0.781000
O         6.205000        6.991000       18.441999
H         6.782000        7.347000       19.152000
H         5.309000        7.062000       18.799999
O        16.125999       14.701000        4.577000
H        16.070000       14.272000        3.726000
………………………………………………………………………………
得到的如下
O   3.6680000     O   10.082000     O   15.904000   
H   3.2240000     H   10.451000     H   15.101000   
O   5.1859999     O   16.695999     O   12.072000   
H   6.0830002     H   16.979000     H   11.816000   
O   3.7860000     O   16.318001     O   1.3099999   
H   3.7809999     H   15.388000     H   1.1960000   
O   6.2049999     O   6.9910002     O   18.441999   
H   6.7820001     H   7.3470001     H   19.152000   
O   16.125999     O   14.701000     O   4.5770001   
H   16.070000     H   14.272000     H   3.7260001   
O   2.7500000     O   8.8439999     O   9.6289997   
H   2.1099999     H   8.2900000     H   9.9980001   
O   10.674000     O   8.9090004     O   4.9860001   
H   10.998000     H   8.8710003     H   5.9320002   
O  2.00000009E-03 O   5.6669998     O   14.591000   
H  0.48300001     H   5.9899998     H   13.830000   
………………………………………………………………………………
这就可以明显的看出问题是出在哪里了，本来读取数据的时候采用的是和写入一样的方式(B(i),A(i,j),j=1,3),一直出现上面的错误，记录一下，以备日后学习汲取教训

[ Last edited by xk6891 on 2011-3-20 at 13:20 ]

进一步修改成这样就顺多了，首先要能读，然后才是写入数据
code:
program main
implicit none
integer i,j
character B
real A
dimension B(768),A(768,3)

open(17,file='wa_xyz',status='old')
do i=1,768
read(17,*)B(i),A(i,1),A(i,2),A(i,3)
end do
close(17)

open(18,file='nem_wa_xyz',status='unknown')
do i=1,768
if(mod(i,3)/=0)then
write(18,"(A3,3(F16.6))"B(i),A(i,1),A(i,2),A(i,3)
end if
end do
close(18)

stop
end program
code:
得到的结果
  O        3.668000       10.082000       15.904000
  H        3.224000       10.451000       15.101000
  O        5.186000       16.695999       12.072000
  H        6.083000       16.979000       11.816000
  O        3.786000       16.318001        1.310000
  H        3.781000       15.388000        1.196000
  O        6.205000        6.991000       18.441999
  H        6.782000        7.347000       19.152000
  O       16.125999       14.701000        4.577000
  H       16.070000       14.272000        3.726000
  O        2.750000        8.844000        9.629000

[ Last edited by xk6891 on 2011-3-20 at 13:25 ]
赞一下(1人) 回复此楼
22楼2011-01-13 23:56:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

uiui

新虫 (小有名气)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xk6891(金币+5):谢谢前辈，问题解决了。用了第二种方法 2011-01-12 10:32:06
余泽成(金币+3):谢谢参与应助！ 2011-01-13 19:47:34
program align_left_output_1
  
   do i=1,200
      L=int(log10(i*1.0))+1
      write(*,100) i
   enddo
   100 format('H'I

end program


program align_left_output_2
   character*8 tmpchar

   do i=1,200
      write(tmpchar,'(I8)') i
      write(*,*) 'H'//adjustl(tmpchar)
   enddo

end program

第一种方法形式简单，容易理解。但变宽度输出依赖于编译器，可移植性不好；而且用到了对数运算，在数据量很大时效率较低。
第二种方法借用了临时字符串做内部文件，效率可能更高一些。
赞一下(2人) 回复此楼
2楼2011-01-12 09:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
引用回帖:
Originally posted by uiui at 2011-01-12 09:17:49:
program align_left_output_1
  
   do i=1,200
      L=int(log10(i*1.0))+1
      write(*,100) i
   enddo
   100 format('H'I

end program


program align_left_output_2
   charact ...

请教前辈一个问题，call random_seed产生随机数的时候，网上说的会根据系统的日期和时间提供种子，但是我每次编译都没有改变，求教一下这是为何，如何使它能生成不同的随机数列。
赞一下 回复此楼
3楼2011-01-12 10:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uiui

新虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
在call random_seed时不要加参数
赞一下(1人) 回复此楼
4楼2011-01-12 13:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +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
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 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
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭