VMD TCL Position Restraint Generator

# Generate position Restrained files for namd simulations
# In VMD:
# source gen_pr
# gen_pr (# of restraint files to generate)
# The files are named #.pdb
# In the beta(B-factor) field a number is placed
# as a force restraint for each atom
# can be changed to CA by changing "set a.." line


proc gen_pr { num } {
   set a [atomselect top protein]
   set b [atomselect top all]

    for {set i $num } {$i >= 0} {incr i -1} {

        $a set beta $i
        $b writepdb $i.pdb
    }

 }