| 查看: 786 | 回复: 21 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】新手请教一段程序的疑问....【已完结】
|
|||
|
我的疑问是,在do i=1,natom, read (1,*) molname,atomname,index,x(i),y(i),z(i),这里do 做 一个循环,我后面的输出要求是输出,x(i-1),y(i-1),z(i-1), x(i),y(i),z(i),x(i+1),y(i+1),z(i+1)这3个坐标,这样输出可以吗?? open(1,file='benz-22.gro',status='unknown') !打开benz-co2.gro文件，并读取 do istep=1, nstep read(1,*) sysname read(1,*) natom do i=1, natom read(1,*) molname,atomname,index,x(i),y(i),z(i) if (atomname .eq.'Cb') then r2=(x(i)-2.623)**2+(y(i)-2.554)**2 if (r2 .LT. 0.01932) then h2=(z(i)-2.608)**2 if (h2 .LT. 0.5776 .AND. z(i) .NE. 2.608) then open(20,file='aromaticco2.txt') write(20,*)'Oa',x(i-1),y(i-1),z(i-1) !输出这个Cb原子的坐标和Cb编号的前面和后面的一个原子的坐标。 write(20,*)'Cb',x(i),y(i),z(i) write(20,*)'Oc',x(i+1),y(i+1),z(i+1) [ Last edited by 波不动 on 2009-11-25 at 23:52 ] |
回复此楼» 猜你喜欢
AI论文写作工具：是科研加速器还是学术作弊器？
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555170
- 散金: 27828
- 红花: 2147
- 沙发: 66666
- 帖子: 1602247
- 在线: 65200.8小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
21楼2009-11-25 22:56:23
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51002.6
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
2楼2009-11-24 12:43:59
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 441.5小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
3楼2009-11-24 16:49:36
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 441.5小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
open(1,file='benz-22.gro',status='unknown') !打开benz-co2.gro文件，并读取 do istep=1, nstep read(1,*) sysname read(1,*) natom I_emuch=1 do i=1, natom read(1,*) molname,atomname,index,x(i),y(i),z(i) if (atomname .eq.'Cb') then r2=(x(i)-2.623)**2+(y(i)-2.554)**2 if (r2 .LT. 0.01932) then h2=(z(i)-2.608)**2 if (h2 .LT. 0.5776 .AND. z(i) .NE. 2.608) then Index_emuch(I_emuch)=i endif endif endif I_emuch=I_emuch+1 Enddo Do J_emuch=1,I_emuch open(20,file='aromaticco2.txt') if(ii.gt.1)write(20,*)'Oa',x(ii-1),y(ii-1),z(ii-1) !输出这个Cb原子的坐标和Cb编号的前面和后面的一个原子的坐标。 write(20,*)'Cb',x(ii),y(ii),z(ii) write(20,*)'Oc',x(ii+1),y(ii+1),z(ii+1) Enddo Enddo |
4楼2009-11-24 17:04:49