[gmx-users] writing script to analyze gmx data

[Jian Zou]

Hi all,

I have a question about writing a script to analyze gmx data..

some analysis utilities of gromacs need keyboard input after running the 
command.

how could I incorporate this keyboard input into the script?


for example, I want use "g_dist" to calculate the distance between Group 2
and Group 5, I write the following script:

----start of the script----
#!/bin/sh
g_dist -f md.trr -s b4md.tpr -n index.ndx -o dist.xvg
----end of the script----

the input file "md.trr", "b4md.tpr", "index.ndx"
and output the distance between Group 2 and Group 5 into the file
"dist.xvg".

after running the script, it will prompt to select the groups, I will input
2 and 5 for this case.

I want to write sth. to put these selection into the script, then what
should I do?


Thanks a lot in advance.


Regards,

Jian Zou