求助，tcl脚本使用foreach遍历出现未统计完就断掉的情况

liucw

最近使用vmd编写脚本时发现这样的一个问题，统计时出现这样的情况，前面统计出现未完全统计就意外结束，只有最后一次循环完全统计结束。请看图，脚本我贴在后面，请大家帮我看看，具体也是根据网上的脚本自己改了一下，加了循环。

1. #---------------------------------------------------
   
2. proc lbond {sel fp1 step} {
   
3. set select1 [atomselect top "$sel"]
   
4. set inde [$select1 get index]
   
5. #---------------------------------------------------
   
6. foreach k $inde {
   
7. set sel1 [atomselect top "index [expr $k-1]"]
   
8. set sel2 [atomselect top "index $k"]
   
9. if {[expr $k % 2] == 0} {
   
10. set outfile [open lbond/[expr $k-1]-$k.dat w]
    
11. set nf [molinfo top get numframes]
    
12. for { set i $fp1 } { $i <= 4700 } { incr i $step} {  
    
13.         $sel1 frame $i
    
14.         set ct1 [measure center "$sel1"]
    
15.         $sel2 frame $i
    
16.         set ct2 [measure center "$sel2"]
    
17.         set vsub [vecsub $ct1 $ct2]
    
18.         set vlen [veclength $vsub]
    
19.         if {$vlen < 4} {
    
20.         puts $outfile "$i $vlen"
    
21.         } else {puts $outfile "$i 0"}
    
22. }
    
23. }
    
24. puts "done!"
    
25. }
    
26. close $outfile
    
27. puts "All Done!"
    
28. }

复制代码