| 查看: 946 | 回复: 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 ] |
回复此楼» 猜你喜欢
0703 化学 求调剂，一志愿山东大学 342 分
已经有3人回复
-
0703化学调剂，求导师收
已经有9人回复
-
0856求调剂
已经有3人回复
-
085701求调剂初试286分
已经有4人回复
-
求调剂
已经有8人回复
-
316求调剂
已经有6人回复
-
315求调剂
已经有5人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有11人回复
-
332求调剂
已经有4人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有11人回复
我的输入文件的部分内容是这样的
|
我的输入文件中共有2500个构象,每一个构象由一个溶质分子(包含18个原子)和1500个CO2分子组成,我想让构象中只要CO2分子中的C原子满足某一个条件,就输出这个CO2分子中的3个原子的坐标..... cinna-co2mix 4518 1CIN C1 1 3.374 3.315 3.175 1CIN C2 2 3.468 3.213 3.177 1CIN C3 3 3.422 3.082 3.177 1CIN C4 4 3.285 3.054 3.176 1CIN C5 5 3.187 3.155 3.177 1CIN C6 6 3.237 3.288 3.177 1CIN H7 7 3.410 3.418 3.174 1CIN H8 8 3.574 3.236 3.178 1CIN H9 9 3.495 3.002 3.177 1CIN H10 10 3.257 2.949 3.176 1CIN H11 11 3.171 3.373 3.180 2CIN C12 12 3.042 3.118 3.176 2CIN H13 13 3.024 3.011 3.176 2CIN C14 14 2.929 3.193 3.176 2CIN H15 15 2.928 3.339 3.177 3CIN C16 16 2.826 3.152 3.176 3CIN H17 17 2.809 3.042 3.176 3CIN O18 18 2.726 3.221 3.176 4DRG Oa 19 1.388 0.106 0.427 4DRG Cb 20 1.274 0.106 0.427 4DRG Oc 21 1.159 0.105 0.426 5DRG Oa 22 1.388 0.106 0.856 5DRG Cb 23 1.274 0.106 0.855 5DRG Oc 24 1.159 0.105 0.854 6DRG Oa 25 1.388 0.106 1.284 6DRG Cb 26 1.274 0.106 1.283 6DRG Oc 27 1.159 0.105 1.283 |
7楼2009-11-24 21:16:21
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51005.2
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
2楼2009-11-24 12:43:59
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
3楼2009-11-24 16:49:36
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 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
