# load multi-frame pdb file, storing B factors from each frame in user.


proc pdbbfactormodif { fname } {
  mol new $fname waitfor all
  mol delrep 0 top
  mol representation CPK
  mol modmaterial 0 0 Glossy

  mol addrep top

  set all [atomselect top all]
  set frame 0
  set numat 0
  set listnumat {}
  set in [open $fname r]
  set beta {}
  while { [gets $in line] != -1 } {
    switch -- [string range $line 0 3] {
      END {
        $all frame $frame
        $all set user $beta
        set beta {}
        incr frame
      }
      ATOM -
      HETA {
        incr numat
        
        lappend beta [string range $line 61 66]
        if { [expr [string range $line 61 66]] > 0.5 } {

          set sel [atomselect 0 {index $numat}]

          $sel set radius 10
        }
      }
    }
  }
}