24小时热门版块排行榜

返回列表

当前主题已经存档。

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

tpp001

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 749.3
散金: 5
红花: 1
帖子: 1040
在线: 155.9小时
虫号: 423029
注册: 2007-07-21
性别: GG
专业: 化学

[交流] 【求助】新手求助,写了一段关于空间中的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
..................

回复此楼