24小时热门版块排行榜    

查看: 2622  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
请教前辈,如何让后面的数字紧贴着字母输入呢后面还有几百行
回复此楼
已阅   回复此楼   关注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 ]
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楼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时不要加参数
4楼2011-01-12 13:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[教师之家] 课题组招2027级博士 上海工程技术大学 激光智能制造方向 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:37 by oxidpl
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +4 13108017953 2026-09-21 4/200 2026-09-24 10:06 by kudofaye
信息提示
请填处理意见