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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 求助计算机高人，编写个Fortran程序
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 2038  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700，311，求调剂 +13 冬十三 2026-04-04 14/700 2026-04-05 22:23 by 醉翁wl
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 一志愿青科085500，初试295分，公共课213分 +3 遇到的人愿望都� 2026-04-05 3/150 2026-04-05 18:45 by 蓝云思雨
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 22408，264求调剂 +3 ywh729 2026-04-03 4/200 2026-04-04 11:04 by ywh729
[考研] 求调剂，一志愿南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-04-03 10/500 2026-04-04 09:01 by T可可西里T
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 初试成绩337找调剂 +3 ??? ?. ? 2026-04-03 3/150 2026-04-03 11:43 by 土木硕士招生
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 一志愿北交大材料工程，总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 358求调剂 +3 王向阳花 2026-03-31 3/150 2026-04-01 09:56 by zzchen2000
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135，现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
信息提示
关闭
请填处理意见
关闭