24小时热门版块排行榜    

查看: 272  |  回复: 0

dreamtiger

铜虫 (初入文坛)

[求助] FATAL ERROR: force not a vector

用namd模拟计算时,出现下列错误
Sample Text
    while executing
"addforce "0.0 0.0 $f""
    (procedure "calcforces" line 8)
    invoked from within
"calcforces 296855 1 18.0 -18.0 -0.001"
运用的脚本

    * tclBC  on
    * tclBCScript{
    * proc calcforces{step unique zp zn} {
    *   if { $step % 20 == 0 } {
    *     cleardrops
    *   }
    *   while {[nextatom]} {
    *     set x [lindex [getcoord] 0]
    *     set y [lindex [getcoord] 1]
    *     set z [lindex [getcoord] 2]
    *     set r2 [expr "$x * $x + $y * $y"]
    *     set m [getmass]
    *     if { ($z>$zp || $z<$zn) && $r2 < $RADII2 && $m<16.0 && $m>15.9} {
    *       addforce { 0 0 -0.01}
    *       continue
    *     }
    *     dropatom
    *   }
    * }
    * }
    * tclBCArgs{   }


bay_gulf68 告诉我让我在add 前加print ,看一看,
我将脚本改成
  if { ($z>$zp || $z<$zn) && $r2 < $RADII2 && $m<16.0 && $m>15.9} {
     print  $f
      addforce { 0 0 -0.01}
再次续跑,结果时间变化很快,一会功夫完成了 未命名.jpg
下载 (2.14 KB)
前天 14:19

不知道运用print是否正确,还有如果论坛中有高手知道怎么解决FATAL ERROR: force not a vector这个问题,希望能够多多指教。
Sample Text
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 dreamtiger 的主题更新
信息提示
请填处理意见