| 查看: 652 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】新手求助,写了一段关于空间中的2点的距离满足一点条件的程序老是出错
|
|||
|
program get-atom-coordinate implicit none integer istep,nstep !构象的个数 integer natom !每个构象中原子的总数 character sysname,molname,atomname !体系的名字，分子的名字，原子的名字 integer i integer index !原子序号 dimension x(4500),y(4500),z(4500) real xbox,ybox,zbox !盒子的大小 real j,k,m,n,sum nstep=250 !共有250个构象 open(1,file='co2.gro',status='unknown') !打开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) enddo read(1,*) xbox,ybox,zbox if (atomname.eq.'Cb') then !原子的名字atomname 为Cb的原子的坐标满足与空间中一点（2.615，2.572，2.601）的距离小于4 j=abs(x(i)**2-2.615**2) k=abs(y(i)**2-2.572**2) m=abs(z(i)**2-2.601**2) n=j+k+m sum=sqrt(n) if (sum < 4.0) then open(2,file='co2ben.txt') write(2,*)x(i-1),y(i-1),z(i-1),x(i),y(i),z(i) !输出这个Cb原子的坐标和Cb编号的前面的一个原子的坐标。 end if end if enddo close(1) end 上面的程序老是出错4个errors,请大家对应上面的注释和下面的输入文件帮 帮忙谢谢......下面是输入文件中的其中一个构象 Generated by trjconv : cinna-co2mix t= 1000.00000 4500 4DRG Oa 1 0.775 1.784 0.721 4DRG Cb 2 0.667 1.750 0.699 4DRG Oc 3 0.560 1.714 0.676 5DRG Oa 4 4.774 0.611 0.616 5DRG Cb 5 4.674 0.667 0.603 5DRG Oc 6 4.575 0.723 0.591 6DRG Oa 7 1.307 2.386 1.515 6DRG Cb 8 1.208 2.405 1.460 6DRG Oc 9 1.109 2.428 1.406 7DRG Oa 10 0.716 1.682 0.398 7DRG Cb 11 0.665 1.603 0.332 7DRG Oc 12 0.611 1.525 0.268 8DRG Oa 13 0.346 3.926 3.246 8DRG Cb 14 0.400 4.023 3.218 8DRG Oc 15 0.456 4.120 3.190 9DRG Oa 16 4.243 2.495 3.733 9DRG Cb 17 4.325 2.444 3.670 9DRG Oc 18 4.404 2.392 3.605 10DRG Oa 19 1.790 1.759 2.631 10DRG Cb 20 1.867 1.742 2.547 10DRG Oc 21 1.941 1.729 2.460 11DRG Oa 22 1.073 4.489 1.693 11DRG Cb 23 1.167 4.549 1.721 11DRG Oc 24 1.259 4.612 1.747 12DRG Oa 25 0.407 4.290 3.671 12DRG Cb 26 0.318 4.339 3.727 12DRG Oc 27 0.230 4.388 3.782 .................. |
回复此楼» 猜你喜欢
0703化学调剂，求导师收
已经有10人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有5人回复
-
一志愿太原理工安全工程300分，求调剂
已经有6人回复
-
0856求调剂
已经有4人回复
-
11408软件工程求调剂
已经有3人回复
-
一志愿211 初试270分 求调剂
已经有8人回复
-
300求调剂，材料科学英一数二
已经有6人回复
-
学硕274求调剂
已经有9人回复
-
本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以
已经有7人回复
-
317分 一志愿南理工材料工程 本科湖工大 求调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
5楼2009-11-04 21:55:50
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2009-11-04 17:03:27
3楼2009-11-04 17:44:58
4楼2009-11-04 21:35:34
