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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 求助计算机高人，编写个Fortran程序
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 2039  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336材料与化工085600求调剂 +8 水星记infp 2026-04-05 8/400 2026-04-06 06:09 by houyaoxu
[考研] 277求调剂 数一104分 +6 瓶子PZ 2026-04-05 6/300 2026-04-05 20:38 by 啵啵啵0119
[考研] 材料调剂 +6 一样YWY 2026-04-05 6/300 2026-04-05 20:30 by 南航~万老师
[考研] 301求调剂 +5 细胞相关蛋白 2026-04-02 9/450 2026-04-05 20:14 by lys0704
[考研] 301求调剂 +3 XYPLR 2026-04-05 4/200 2026-04-05 19:07 by XYPLR
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 17:59 by 猪会飞
[考研] 一志愿南航，数一英一学硕317求调剂！！ +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] 求调剂，一志愿郑州大学材料与化工专硕，英二数二342分，求老师收留 +18 v12abo 2026-04-02 20/1000 2026-04-05 11:37 by a8144223
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 311分 22408 求调剂 +3 bing_bot 2026-04-03 3/150 2026-04-05 00:43 by chongya
[考研] 求调剂：085600材料与化工，考材科基，总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 材料调剂 +11 吴棂颖！ 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 265求调剂 +17 林深温澜 2026-04-01 20/1000 2026-04-04 01:09 by userper
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 一志愿9初试366 本双非求调剂 +4 运气来得若有似� 2026-04-02 4/200 2026-04-02 09:56 by guanxin1001
[考研] 【调剂】一志愿厦大生物与医药调剂 +3 Echo虾米 2026-03-31 3/150 2026-04-01 08:40 by JourneyLucky
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
信息提示
关闭
请填处理意见
关闭