应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】新手求助,写了一段关于空间中的2点的距离满足一点条件的程序老是出错
51/1返回列表
查看: 842  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tpp001

金虫 (著名写手)

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
迷茫在知识的海洋里，需要你的指导。thankyou
1楼 2009-11-04 16:55:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haixiang

金虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+2,VIP+0):感谢应助！ 12-17 18:57
将sum值打印出来看看是否满足条件
另外 sum=sqrt(n)
if (sum < 4.0) then
这行写错了
赞一下(11人) 回复此楼
7楼2009-12-17 16:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1,VIP+0):谢谢应助！ 11-4 17:53
是编译没通过吗？请把4个error信息也给出来看看。
引用回帖:
j=abs(x(i)**2-2.615**2)
k=abs(y(i)**2-2.572**2)
m=abs(z(i)**2-2.601**2)

这里是istep还是i?
赞一下(7人) 回复此楼
2楼2009-11-04 17:03:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
nono2009(金币+0,VIP+0):但是你的这些i都不在循环内，其实都等于natom+1. 11-4 19:17
是i ,这里的 i 其实就是原子的编号，因为我一个构象里面有4500个原子。
istep是步数。。分子动力学每跑一步，就出来一个构象。。。。。
赞一下(10人) 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
3楼2009-11-04 17:44:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tpp001

金虫 (著名写手)
出现的错误信息是:
D:\fortran\compile\getdistance.f90(20) : Error: This name does not have a type, and must have an explicit type.   [Z]
read(1,*) molname,atomname,index,x(i),y(i),z(i)
我觉得好象是没有定义x(i),y(i),z(i)把,不知道应该怎么定义他们,请大家指教
赞一下 回复此楼
迷茫在知识的海洋里，需要你的指导。thankyou
4楼2009-11-04 21:35:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[基金申请] 申请科研项目纯靠实力能成功吗？ 200+4 David1982 2026-10-13 7/350 2026-10-14 14:24 by 13108017953
[硕博家园] 极限拉扯，如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭