A bug in replica.namd?

From: Ada Zhan (zyqfrog10_at_msn.com)
Date: Wed Nov 07 2012 - 16:37:13 CST

Hi folks,

I am testing the replica exchange with abf in namd 2.9.
It works for many steps initially until the following message:
====
TCL: expected integer but got ""
FATAL ERROR: expected integer but got ""
    while executing
"format "$restart_root.tcl" """
    ("while" body line 69)
    invoked from within
"while {$i_run < $num_runs} {

  run $steps_per_run
  save_array
  incr i_step $steps_per_run
  set TEMP $saved_array(TEMP)
  set POTENTIAL [expr $save..."
    (file "~/NAMD_2.9_Source/lib/replica/replica.namd" line 98)
    invoked from within
"source ~/NAMD_2.9_Source/lib/replica/replica.namd "
    invoked from within
"if { ! [catch numPes] } { source ~/NAMD_2.9_Source/lib/replica/replica.namd }"

=====
I checked the replica.namd script. The error happens at the line :

set rfile [open [format "$restart_root.tcl" ""] "w"]

====
There is indeed a "" after "format ...". But what integer should be expected here ? Is that a bug in this script?

Thanks, all!
-Ada

This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:22:44 CST