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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】新手求助,写了一段关于空间中的2点的距离满足一点条件的程序老是出错
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 654  |  回复: 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的回帖

jylisace

新虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+1,VIP+0):感谢参与！ 1-4 09:52
如果不满足(sum < 4.0)，肯定无法输出文件。没有你的输入文件，没法帮你验证。
赞一下(11人) 回复此楼
8楼2010-01-04 09:33:01
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学354求调剂 +3 101次希望 2026-03-28 3/150 2026-03-29 00:28 by 544594351
[考研] 295求调剂 +4 wei-5 2026-03-26 4/200 2026-03-28 23:20 by 小木虫tim
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 295材料工程专硕求调剂 +7 1428151015 2026-03-27 7/350 2026-03-28 19:58 by S240
[考研] 321求调剂 +6 璞玉~~ 2026-03-25 7/350 2026-03-28 17:48 by 璞玉~~
[考研] 材料与化工272求调剂 +9 阿斯蒂芬2004 2026-03-28 9/450 2026-03-28 15:21 by 1018329917
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 0703化学求调剂，各位老师看看我！！！ +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 求调剂 +6 林之夕 2026-03-24 6/300 2026-03-27 08:38 by hypershenger
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 机械学硕310分，数一英一，一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 300分，材料，求调剂，英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
信息提示
关闭
请填处理意见
关闭