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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】如何实现数字紧贴字母输入
51/1返回列表
查看: 2006  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +7 人间唯你是清欢 2026-02-28 10/500 2026-02-28 20:16 by 公瑾逍遥
[考研] 284求调剂 +3 天下熯 2026-02-28 3/150 2026-02-28 20:14 by 公瑾逍遥
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考博] 博士推荐 +5 花儿笑? 2026-02-21 6/300 2026-02-28 18:53 by nxgogo
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[考研] 311求调剂 +7 南迦720 2026-02-28 7/350 2026-02-28 18:28 by leonnulll
[考研] 材料学调剂 +4 提神豆沙包 2026-02-28 4/200 2026-02-28 18:26 by houyaoxu
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 5/250 2026-02-28 18:20 by liqiongjy
[考研] 材料类求调剂 +4 wana_kiko 2026-02-28 4/200 2026-02-28 18:08 by djennjx
[考研] 材料调剂 +3 爱擦汗的可乐冰 2026-02-28 3/150 2026-02-28 18:06 by houyaoxu
[考研] 求调剂 +3 repeatt?t 2026-02-28 3/150 2026-02-28 18:00 by houyaoxu
[考研] 化工专硕348，一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 290求调剂 +4 材料专硕调剂； 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭