[gmx-developers] issues in calling python with subprocess in git master code
Shirts, Michael R. (mrs5pt)
mrs5pt at eservices.virginia.edu
Fri Nov 28 05:59:47 CET 2014
I'm trying to write a python script using subprocess to call grompp. It
works perfectly fine when called from the command line or when debugging;
this only occurs with python (Anaconda, Python 2.7).
For example, my code looks like:
It produces to stdout:
GROMACS: gmx grompp, VERSION 5.1-dev-20141126-f367cf6 (double precision)
grompp -f r0000.noiter.mdp -c trpo.gro -p trpo.top -o r0000.noiter.tpr
However, the errors that it then gives and mdout.mdp indicate that the mdp
is not actually being read in, and the default values are being used
instead. It even reports that the errors are in the mdp file I fed in,
even though the mdp does not contain those commands.
All files work perfectly well when called from the command line.
The problems also occur the same way when calling as:
So that's not the issue. Same thing happened with os.system.
I've tried a bunch of different ways of calling it in python, and all fail
to have the mdp read in. All other files are read in, and it doesn't
matter what order they are in.
I'm just wondering if other people have experienced this sort of problem
calling gromacs with Python, or if there is something weird I'm doing.
Department of Chemical Engineering
University of Virginia
michael.shirts at virginia.edu
More information about the gromacs.org_gmx-developers