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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】新手求助,写了一段关于空间中的2点的距离满足一点条件的程序老是出错
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 644  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 学硕274求调剂 +8 Li李鱼 2026-03-26 8/400 2026-03-28 16:53 by xuxiang
[考研] 调剂310 +6 温柔的晚安 2026-03-25 7/350 2026-03-28 16:04 by qq2000
[考研] 本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以 +5 芝士雪baoo 2026-03-28 6/300 2026-03-28 15:25 by 1018329917
[考研] 0703化学 +10 妮妮ninicgb 2026-03-27 10/500 2026-03-28 14:07 by 哒哒哒呱呱呱
[考研] 一志愿太原理工安全工程300分，求调剂 +4 0857求调剂. 2026-03-24 5/250 2026-03-28 13:45 by 求知若渴lz
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 305求调剂 +5 哇卢卡库 2026-03-26 5/250 2026-03-27 14:01 by laoshidan
[考研] 316求调剂 +5 Pigcasso 2026-03-24 5/250 2026-03-27 12:10 by zhshch
[考研] 一志愿郑大085600，310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 材料学硕，求调剂 6+5 糖葫芦888ll 2026-03-22 10/500 2026-03-27 08:18 by hypershenger
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 一志愿河工大 081700 276求调剂 +4 地球绕着太阳转 2026-03-23 4/200 2026-03-26 14:27 by zzll406
[考研] 263求调剂 +6 yqdszhdap－ 2026-03-22 10/500 2026-03-26 13:11 by 公瑾逍遥
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 机械学硕总分317求调剂！！！！ +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 0854AI CV方向招收调剂 +4 章小鱼567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 材料/农业专业，07/08开头均可，过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪！？！
[考研] 315分，诚求调剂，材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
信息提示
关闭
请填处理意见
关闭