Backup MD data processing related scripts.
These scripts will do following tasks in order:
- convert gromacs .trr files to .xtc files, and merge into 'merge.xtc'.
- extract Atom Group ‘protein’ from original .gro file to a new file 'npt_ab_protein.gro' for subsequent processing.
- periodic boundary condition treatment by methods:
- protocol1: whole->atom
- protocol2: whole->atom->nojumps
- center & align and generate .xtc files
- protocol1: atom_rottrans.xtc
- protocol2: atom_rottrans_nojump.xtc
- plot information in log files: temperature/pressure/kinetic/potential/viral/lbox
- automate visualization by vmd tcl scripts.
Main script. Copy post_process_mv.sh/trr2xtc.py/atom_only_pbc.sh/draw_log_plots.py/useme.tcl to work directory and execute processing protocol.
Move trr files to ./coor_nc, log files to ./log_txt.
convert trr to xtc and merge all xtc files.
Deal with PBC in two protocols. Check atom index for PBC treatment and centering manually and interact with scripts. Output with Atom Group 'protein' by default, update this script for other types of entities.
Visualization by vmd:
- Change display method to orthographic
- save visualization state to 'useme.vmd'
- load trajectory 'atom_rottrans.xtc'
Read log files in directory 'log_txt', plot temperature/pressure/kinetic/potential/viral/lbox
- updata variable SCRIPTS_DIR in md_process.sh to access to all neccesery scripts.
- copy md_process.sh to data's directory, for example 'md_post_output', source it.
- input initial gro file
- input atom index for PBC treatment and centering.
- update automatic center atom selection, for example by centroid.
- pbc nojump is not robust.