求助tcl，循环出问题

chw985550192

for {set j 96} {$j <= 117} {incr j 1} {

set outfile [open num-$j.dat w]

set O5 [[atomselect top {not water and name "O.*" and within 3 of index $j}] get index]
set  O5num [llength $O5]
puts  $outfile [format "%5d" $O5num]

}

close $outfile

我想统计index 96 to 117 附近的氧原子

我直接在tk console中输入 set O5 [[atomselect top {not water and name "O.*" and within 3 of index 96}] get index] ,输出的结果为120 125 136 158 169

可是做上述循环的时候输出 1 2 10，输出结果竟然和上述不一样，且输出是错误的。


请问上述的循环中那里出问题了呢，谢谢。

enthalpy

{set j 96} {$j <= 117} {incr j 1} {

set outfile [open num-$j.dat w]

set O5 [[atomselect top “not water and name "O.*" and within 3 of index $j"] get index]

set  O5num [llength $O5]

puts  $outfile [format "%5d" $O5num]

}

你要把花括号替换成双引号，里面的变量j 才能替换。

chw985550192

enthalpy 发表于 2020-7-10 15:09
{set j 96} {$j

我给的例子是对当前帧数进行输出，并没有涉及到具体那一帧。我只是想问上述的循环那里出问题。

enthalpy

chw985550192 发表于 2020-7-10 15:18
我给的例子是对当前帧数进行输出，并没有涉及到具体那一帧。我只是想问上述的循环那里出问题。

我刚才看错了，你看上面修改的回复。

chw985550192

enthalpy 发表于 2020-7-10 15:19
我刚才看错了，你看上面修改的回复。

谢谢你的回答，你说的会有一个问题，我定义的 name "O.*",已经有双引号了，如果再有一个双引号就会报错，请问这个问题要怎么解决呢。

wyb

每次读取一帧的时候需要清空之前的选择