版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

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

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

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

回复此楼

» 猜你喜欢

真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有3人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

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

gibbs程序包，fortran版，CPC库程序，没有错误，可直接编译运行已经有9人回复
求一个简单的读写文件的小程序。。。已经有4人回复
用Fortran编写程序输出的数据一般用什么程序作图比较好？已经有28人回复
求助：用Fortran写一个关于晶粒个数统计和编号的程序。卡了我好久了。已经有4人回复
求助：intel visual fortran的格式问题已经有14人回复
fortran输出求助已经有6人回复
一个主程序，一个subroutine程序，要如何将两个程序组合在fortran里运行已经有4人回复
求fortran95编写的数值程序已经有4人回复
fortran中到底该如何比较两个实数是否相等已经有8人回复
求高人帮忙看一下我这个Fortran程序的错误出在哪里？已经有12人回复
大家帮我看一个fortran的程序，我总是计算不出正确的结果已经有20人回复
我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？已经有4人回复
求助：编写程序已经有5人回复
有没有人关注fortran语言的GPU并行呢？已经有12人回复
【求助】如何写FORTRAN程序实现求平均最近邻距离已经有34人回复
【求助】FORTRAN写文件时每行开头都多了一个空格？已经有6人回复
【求助】求Fortran90编译器？已经有11人回复
【求助】初学fortran，已经有15人回复
【求助】如何实现编写的fortran程序让别人使用，但不让别人知道源代码已经有6人回复
【分享】fortran程序设计视频教程【已搜索无重复】已经有15人回复

生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`

1楼2011-09-06 16:05:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

引用回帖:

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的回帖

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

引用回帖:

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的回帖

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

引用回帖:

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的回帖

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

引用回帖:

10楼: Originally posted by snoopyzhao at 2011-09-07 16:29:46:
这样吧，你把
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))
再 ...

恩，好的。谢谢你了~~

赞一下

回复此楼

生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`

12楼2011-09-07 16:39:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖