tomnewport / alchex

Phospholipid Biosynthetic Proteins Project
1 stars 0 forks source link

Simulation container object #7

Closed tomnewport closed 8 years ago

tomnewport commented 8 years ago

Most processes one might want to do involve:

The simulation container will be based around the grompp step. Several GROMACS commands produce files in the working directory.

structure = GROFile("structure.gro")
topology = TOPFile("topol.top")
em_mdp = MDPFile("em.mdp")
md_mdp = mdp.MD_MDP

s = SimulationContainer()
s.add(structure)
s.add(topology)
s.add(em_mdp)

em = s.folder("em")
em.gromacs.grompp(c="/structure.gro", p="/topol.top", f="/em.mdp", o="em.tpr")
em.gromacs.mdrun(deffnm="em")

md = s.folder("md")
s.md.gromacs.grompp(c="/em/em.gro", f="/md.mdp" p="/topol.top")
s.md.bash("echo llama >> llama.txt")
s.md.bash(["echo", "llama"])

In the future this container could also keep a json dictionary to store:

It should also be possible to save and load a container. Later work could allow the container to resume after a crash, or synchronise with a server using SAGA.

Tasks: