from phaser import *
i = InputNMA()
i.setROOT("beta_nma")
i.addENSE_PDB_ID("beta","beta.pdb",1.0)
i.setNMAP_MODES([7,10])
i.setNMAP_FORWARD()
r = runNMA(i)
print r.summary()
print "Normal Mode Analysis"
for i in range(0,r.getNum()):
  print "PDB file = ", r.getPdbFile(i)
  displacement = r.getDisplacements(i)
  mode = r.getModes(i)
  for j in range(0,mode.size()):
    print " Mode = " , mode[j], " Displacement = ", displacement[j]