24小时热门版块排行榜

返回列表

fzuvivi

金虫 (小有名气)

应助: 4 (幼儿园)
金币: 1168
帖子: 224
在线: 330.1小时
虫号: 881678

[交流] 【求助】CRASH: from init：error # 1 q not allowed*

利用ph.x顺利跑完动力学矩阵后，想用q2r.x计算力常数矩阵，可是出现一个CRASH：

from init：error # 1
q not allowed

这该怎么解决。谢谢大家了～

q2r.in:

&input
zasr='simple',
fildyn='CO.dyn',
flfrc='CO.fc',
la2F=.false.
/

回复此楼

» 猜你喜欢

2026年国自然面上资助率已经有20人回复
一个有机合成实验室都需要哪些设备？已经有5人回复
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！已经有22人回复
微信指数没变化，科研之友没阅读已经有18人回复
基础研究怎么拉横向，学校到款任务越来越多，难以完成拉横向，都有哪些途径啊已经有8人回复
HXDI做水性聚氨酯乳液，是不是特别容易出渣已经有3人回复
系统今天又提示维护了，估计离放榜不远了已经有15人回复
你们的时间戳变了吗已经有4人回复
产物和副产物价值比较已经有5人回复
时间戳他又来了已经有17人回复

» 本主题相关价值贴推荐，对您同样有帮助:

解晶体时遇到：Cell contents from UNIT instuction and atom list do not agree 已经有5人回复
fcc格子，执行q2r.x时出错，提示，q not allowed 已经有7人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-02-28 10:34:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

léon_WANG

木虫 (小有名气)

应助: 5 (幼儿园)
金币: 2577.5
帖子: 156
在线: 226.1小时
虫号: 2458360

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

楼主问题解决了吗？我也碰到这问题。这是代码 DO nq = 1,nqs WRITE(6,'(a,3f12.8)') ' q= ',(q(i,nq),i=1,3) lq = .TRUE. DO ipol=1,3 xq = 0.0d0 DO icar=1,3 xq = xq + at(icar,ipol) * q(icar,nq) * nr(ipol) END DO lq = lq .AND. (ABS(NINT(xq) - xq) .LT. eps) iq = NINT(xq) ! m(ipol)= MOD(iq,nr(ipol)) + 1 IF (m(ipol) .LT. 1) m(ipol) = m(ipol) + nr(ipol) END DO IF (.NOT.lq) CALL errore('init','q not allowed',1)，就是读不懂它是怎么来判断的

赞一下