[gmx-users] Velocity as a function of distance Z (Mark Abraham)

Kamps, M. m.kamps at student.rug.nl
Wed Mar 1 14:11:00 CET 2017


Dear Mark,

Thanks again for your reply. I'm sorry for asking these probably
stupid questions, but I'm not able to figure it out.

The objects I am interested in are small polymers, consisting of three
residue groups; Eth EthE and EthB (corresponding to both the end
groups of a polymer, and the middle section). I total these groups
consist of around 25000 atoms.

With gmx select I will enter the following command, at an arbitrarily
chosen time:
gmx select -f input.trr -s input.tpr -on output.ndx -b 7.5 -e 7.5
In the following selection screen my preferred molecules are listed as
three different groups (Eth, EthE and EthB), therefore I manually
select these groups with the command: resname Eth EthE EthB, press
enter and end with Ctrl-D. It will then proceed to process the frames,
where GROMACS tells me it has analyzed 1 frame, at the 30th timestep,
which is 7.5.

I will then switch to gmx traj where I enter the following:
gmx traj -f input.trr -s input.tpr -n output.ndx -ov output.xvg
It will then process all frames, and not only the frame I selected
during gmx select. Do I need to specify -b and -e again during gmx
traj?

Then, the output could be read with xmgrace, however looking at the
files it will try to plot 25000 different lines, which it probably
cant. Opening the .xvg file plots the velocity/time plot of only 1
atom, while opening the  .xvg file for the time selected output (via
gmx traj -b -e) shows an empty plot for only one atom.

I'm getting lost in all the options, any help would be appreciated.

Mark


More information about the gromacs.org_gmx-users mailing list