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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 求助计算机高人，编写个Fortran程序
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 2037  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物与医药求调剂 +5 heguanhua 2026-04-05 6/300 2026-04-05 22:58 by Hdyxbekcb
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 085602调剂 初试总分335 +8 19123253302 2026-04-05 8/400 2026-04-05 20:44 by 南航~万老师
[考研] 0854电子信息319求调剂（接受跨专业调剂） +3 星星不眨眼喽 2026-04-05 3/150 2026-04-05 20:20 by 啵啵啵0119
[考研] 298求调剂 +3 manman511 2026-04-05 3/150 2026-04-05 18:09 by kk112233
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 0703总分331求调剂 +10 ZY-05 2026-04-04 13/650 2026-04-05 11:03 by xiayan13521
[考研] 359求调剂22408 +3 123456789qw 2026-03-31 3/150 2026-04-05 10:09 by zhq0425
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 315求调剂 +13 小羊小羊_ 2026-04-02 14/700 2026-04-04 20:30 by 蓝云思雨
[考研] 309分085801求调剂 +11 MY_angel 2026-03-31 11/550 2026-04-04 19:11 by 蓝云思雨
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 366求调剂一志愿东北大学 +8 运气来得若有似� 2026-04-02 8/400 2026-04-02 21:39 by dongzh2009
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
信息提示
关闭
请填处理意见
关闭