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.
/

回复此楼

» 猜你喜欢

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

解晶体时遇到：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)，就是读不懂它是怎么来判断的

赞一下