24小时热门版块排行榜    

查看: 2370  |  回复: 13

q1036210267

木虫 (正式写手)

[求助] 求助计算机高人,编写个Fortran程序

求助计算机高人
编写Fortran程序,将waterbox.xyz文件中所有氧原子的X+1.0, Y+2.0, Z+3.0,而氢原子的保持不变。改变后的结果保存为一样的格式,并且命名为waterbox.are.

1500为水分子数目 ,然后第一竖行为第一个,第二个。。第N个水分子(这里我只写到40)第二竖行为氧原子和氢原子,接下来的三个竖行是需要改变的(第三四五竖行的氧原子坐标改变)。以后其他竖行的都不需要改变。
希望计算机高手给我编下,不胜感激
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

生活中坚持快乐的原则不改变,学习工作中踏实上进!`~`
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
xzhdty: 欢迎常来程序语言看看 2011-09-06 22:23:55
ben_ladeng: 2011-09-07 00:23:25
微尘、梦想(金币+5): 感谢应助~ 2011-09-07 12:14:46
q1036210267(金币+20): 辛苦了 2011-09-07 16:40:38
很简单的嘛……
CODE:
program waterbox
   implicit none
   integer :: num, no, i, i1, i2, i3
   character(len=3) :: ele
   character(len=80) :: line
   real :: x, y, z

   open(unit=11, file='waterbox.xyz', status='old')
   open(unit=12, file='waterbox.are', status='unknown')

   read(11,'(i6)') num
   write(12,'(i6)') num

   do i=1,num
      if (mod(i,3) == 1) then
         read(11, '(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x,y,z,i1,i2,i3
         write(12,'(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x+1.,y+2.,z+3.,i1,i2,i3
      else
         read(11,'(a)') line
         write(12,'(a)') trim(line)
      end if
   end do

end program waterbox

2楼2011-09-06 17:41:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

865504913

木虫 (著名写手)

xzhdty: 欢迎常来 2011-09-06 22:24:19
引用回帖:
2楼: Originally posted by snoopyzhao at 2011-09-06 17:41:00:
很简单的嘛……
[code]
program waterbox
   implicit none
   integer :: num, no, i, i1, i2, i3
   character(len=3) :: ele
   character(len=80) :: line
   real :: x, y, z

   open(unit=11, fi ...

高人无处不在啊  本科学的FORTRAN  现在有点淡忘了 不过 2楼还是很厉害
充满挑战的生活只有主动出击才能赢得胜算!!
3楼2011-09-06 18:33:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

q1036210267

木虫 (正式写手)

引用回帖:
2楼: Originally posted by snoopyzhao at 2011-09-06 17:41:00:
很简单的嘛……
[code]
program waterbox
   implicit none
   integer :: num, no, i, i1, i2, i3
   character(len=3) :: ele
   character(len=80) :: line
   real :: x, y, z

   open(unit=11, fi ...

不好意思现在才回复你,最后那几列数字应该不变,具体的参数怎么设啊。
我是刚开始学,所以很多都不懂。麻烦你了



生活中坚持快乐的原则不改变,学习工作中踏实上进!`~`
4楼2011-09-07 09:29:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
xzhdty(金币+1): 欢迎常来 2011-09-07 13:46:53
你确定你运行的是我的程序和你给出的那个示例文件吗?我这里最后一列没有变化啊……

需要注意的是,这个程序对输入文件的格式有着非常严格的要求,所以一旦输入文件格式发生变化,结果就会不同,甚至得不到正确的结果……
5楼2011-09-07 10:46:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

q1036210267

木虫 (正式写手)

引用回帖:
5楼: Originally posted by snoopyzhao at 2011-09-07 10:46:04:
你确定你运行的是我的程序和你给出的那个示例文件吗?我这里最后一列没有变化啊……

需要注意的是,这个程序对输入文件的格式有着非常严格的要求,所以一旦输入文件格式发生变化,结果就会不同,甚至得不到正确 ...

我就是按照这个程序输入的啊。。结果还是刚才显示的那样,为什么呢? 是不是哪里   出错拉


      program waterbox
      implicit none
      integer :: num, no, i, i1, i2, i3
      character(len=3) :: ele
      character(len=80) :: line
      real :: x, y, z

      open(unit=11, file='waterbox.xyz', status='old')
      open(unit=12, file='waterbox.are', status='unknown')

      read(11,'(i6)') num
      write(12,'(i6)') num

      do i=1,num
      if (mod(i,3) == 1) then
         read(11, '(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x,y,z,i1,i2,i3
         write(12,'(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x+1.,y+2.,z+3.,i1,i2,i3
      else
         read(11,'(a)') line
         write(12,'(a)') trim(line)
      end if
      end do

      end program waterbox
生活中坚持快乐的原则不改变,学习工作中踏实上进!`~`
6楼2011-09-07 12:56:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
余泽成(金币+1): 谢谢参与应助! 2011-09-07 22:00:58
引用回帖:
6楼: Originally posted by q1036210267 at 2011-09-07 12:56:13:
我就是按照这个程序输入的啊。。结果还是刚才显示的那样,为什么呢? 是不是哪里   出错拉


      program waterbox
      implicit none
      integer :: num, no, i, i1, i2, i3
      character(len= ...

你的源程序是不是用的固定格式?或者你的编译器默认使用固定格式?

你看一下
         write(12,'(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x+1.,y+2.,z+3.,i1,i2,i3
这一行是不是超出 72 列了……
7楼2011-09-07 16:03:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
余泽成(金币+1): 谢谢参与应助! 2011-09-07 22:01:12
你把源码命令为 file.f90,然后再编译下看看……记得用 f90 或更新版本的编译器……
8楼2011-09-07 16:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

q1036210267

木虫 (正式写手)

引用回帖:
8楼: Originally posted by snoopyzhao at 2011-09-07 16:04:44:
你把源码命令为 file.f90,然后再编译下看看……记得用 f90 或更新版本的编译器……

write(12,'(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x+1.,y+2.,z+3.,i1,i2,i3

中 i3 显示为绿色,我截个图给你看看。这是为什么啊?

生活中坚持快乐的原则不改变,学习工作中踏实上进!`~`
9楼2011-09-07 16:22:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★
余泽成(金币+1): 谢谢参与应助! 2011-09-07 22:01:27
这样吧,你把
         write(12,'(i6,a3,tr2,3(f12.6),3(i6))') no,ele,x+1.,y+2.,z+3.,i1,i2,i3
改为
         write(12,100) no,ele,x+1.,y+2.,z+3.,i1,i2,i3
100   format(i6,a3,tr2,3(f12.6),3(i6))
再试试吧

我觉得你最好还是找本基础 Fortran 语法的书的看一下吧
10楼2011-09-07 16:29:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员94tPwG 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见