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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 求助计算机高人，编写个Fortran程序
141/2返回列表12下一页
查看: 1989  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

q1036210267

木虫 (正式写手)

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`
1楼 2011-09-06 16:05:16
已阅   回复此楼   关注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
赞一下(3人) 回复此楼
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楼还是很厉害
赞一下(1人) 回复此楼
充满挑战的生活只有主动出击才能赢得胜算！！
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
你确定你运行的是我的程序和你给出的那个示例文件吗？我这里最后一列没有变化啊……

需要注意的是，这个程序对输入文件的格式有着非常严格的要求，所以一旦输入文件格式发生变化，结果就会不同，甚至得不到正确的结果……
赞一下(1人) 回复此楼
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 列了……
赞一下(1人) 回复此楼
7楼2011-09-07 16:03:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
余泽成(金币+1): 谢谢参与应助！ 2011-09-07 22:01:12
你把源码命令为 file.f90，然后再编译下看看……记得用 f90 或更新版本的编译器……
赞一下(1人) 回复此楼
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 语法的书的看一下吧
赞一下(1人) 回复此楼
10楼2011-09-07 16:29:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员94tPwG 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +6 墨墨漠 2026-03-18 6/300 2026-03-20 00:02 by 23Postgrad
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 281求调剂（0805） +9 烟汐忆海 2026-03-16 19/950 2026-03-19 11:42 by laoshidan
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
信息提示
关闭
请填处理意见
关闭