| 查看: 949 | 回复: 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 ] |
回复此楼» 猜你喜欢
11408软件工程求调剂
已经有3人回复
-
一志愿211 初试270分 求调剂
已经有8人回复
-
300求调剂，材料科学英一数二
已经有6人回复
-
学硕274求调剂
已经有9人回复
-
本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以
已经有7人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有3人回复
-
317分 一志愿南理工材料工程 本科湖工大 求调剂
已经有7人回复
-
081200-11408-276学硕求调剂
已经有5人回复
-
316求调剂
已经有7人回复
-
322求调剂
已经有7人回复
|
哈哈,的确是"read(1,*) xbox,ybox,zbox !读入盒子的大小" 这个的位置问题.. 但是我还是很奇怪,为什么这个位置要放在下面程序的红色这个位置才行,而放在兰色这个位置就不行呢??? 全程序如下... program getdistance implicit none integer istep,nstep !构象的个数 integer natom !每个构象中原子的总数 character*80 sysname,molname,atomname !体系的名字，分子的名字，原子的名字 integer i,p,q,r,atomindex(500) !定义了一个存编号的数组 integer index !原子序号 real x(4514),y(4514),z(4514) !原子坐标 real xbox,ybox,zbox !盒子的大小 real j,k,m,n,sum nstep=2500 !共有2500个构象 open(20,file='co2ben.txt',status='unknown') !打开一个输出文件 open(1,file='benz-7.gro',status='unknown') !打开co2.gro文件，并读取 do istep=1, nstep p=0 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 !原子的名字atomname 为Cb的原子的坐标满足与空间中一点（2.615，2.572，2.601）的距离小于4 j=abs(x(i)-3.019)**2 k=abs(y(i)-3.426)**2 m=abs(z(i)-3.129)**2 n=j+k+m sum=sqrt(n) if (sum < 0.56) then p=p+1 atomindex(p)=i end if end if enddo read(1,*) xbox,ybox,zbox !读入盒子的大小do q=1,p r=atomindex(p) write(20,*)'Oa',x(r-1),y(r-1),z(r-1) !输出这个Cb原子的坐标和Cb编号的前面的一个原子的坐标。 write(20,*)'Cb',x(r),y(r),z(r) write(20,*)'Oc',x(r+1),y(r+1),z(r+1) read(1,*) xbox,ybox,zbox !读入盒子的大小 enddo enddo close(1) end |
19楼2009-11-25 22:53:06
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
