[gmx-users] The energy output code

Mark Abraham Mark.Abraham at anu.edu.au
Wed Jun 21 14:43:24 CEST 2006


Uwe Richter wrote:
> Hi,
> 
> I'd rather resort to the following.
> 
> Extract each single frame from your trajectory.
> Check whether your criteria for including certain atoms/indices are
> fulfilled.
> Create a separate index file for each frame.
> Recalculate the energy.
> 
> With some scripting ability this approach could be automatized easily.

Indeed. If you extra the single frames to a .gro or .pdb format to make
life easy for the script, remember to take enough precision from the
binary trajectory to the text gro/pdb so that your subsequent energy
calculations are not meaningfully dependent on your text precision.
Check the manpage of tpbconv for details about this.

Mark



More information about the gromacs.org_gmx-users mailing list