Skip to content

Instantly share code, notes, and snippets.

@hainm
Created September 13, 2016 05:46
Show Gist options
  • Save hainm/7614bfdd85094b3670de66b687c4f169 to your computer and use it in GitHub Desktop.
Save hainm/7614bfdd85094b3670de66b687c4f169 to your computer and use it in GitHub Desktop.
import parmed as pmd
import sander
# make default pme input
inp = sander.pme_input()
parm_file = '../ubiquitin/prmtop'
rst7_file = '../ubiquitin/restrt.save'
parm = pmd.load_file(parm_file, rst7_file)
with sander.setup(parm_file, parm.coordinates, parm.box, mm_options=inp):
e, f = sander.energy_forces()
print(e.elec)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment