Script Library

This page is designed to serve as a repository for files that may be useful to group members. For information on uploading files, see notes below.

Queue submission scripts.

The apoa.sh model script is for a simple 8-processor job in NAMD. The multiple_images.sh script is to run 8 single-processor jobs in NAMD simultaneously, while the sequential-jobs-separately-submitted.sh script is for submitting a series of jobs in series, as would be done in aimless shooting, for example.

The qstatus script gives a summary of the state of all nodes in group0.

MD Analysis scripts

Basic analysis for NAMD trajectories

basic_analysis.sh A simple script for analyzing NAMD output files and trajectories. This script extracts energy, temperature, and pressure values from a simulation log file and plots them as a timeseries, along with running average of energy, RMSD of system and of a selection (specified at command line), and density. Users should feel free to adapt this script to their own needs.

Note that this requires runavg.awk (assumed to be in /home/$(whoami)/tools), and note that the formula for calculating density uses only simulation cell lengths: (volume) = a*b*c. Run with options "--usage" or "--help" for more details:

Usage: basic_analysis.sh --name SIMNAME [ --logfile FILE ] [ --psf PSFFILE ] [ --outputdir ODR ] [ --dcdfreq N ] [ --rmsdsel SELTEXT ]

Performs simple analysis tasks for simulation without output SIMNAME.
OPTIONS:
	--logfile FILE : Use FILE as NAMD log file.  Default value is SIMNAME.log
 	--psf PSFFILE  : Use PSFFILE (either CHARMM or AMBER format) as system topology file.
			 Should have .psf or .parm7 extension for filetype recognition by VMD.
			 Default value: SIMNAME.psf
	--outputdir ODR: Write analysis output to directory ODR. Figures will go to ODR/figures/.
			 Default value: $(dirname SIMNAME)/analysis_output/
	--dcdfreq N    : Frequency at which frames were written to DCD file SIMNAME.dcd. Default 100.
	--rmsdsel SELTXT: VMD selection text for RMSD calculation. Must be enclosed in quotes.
			 Default value: "all".

Sample scripts for basic tasks in CHARMM

The scripts.tar.gz tar ball contains following scripts, along with a small set of auxiliary (data) files. These scripts are also available on this link, and questions about how to use these scripts are addressed in this thread on the CHARMM forums.

VMD TCL script to calculate density profile normal to z-axis

denprofile_analysis.tcl TCL script for VMD to calculate density profile normal to z-axis. Run the script with argument "usage" for more info, or see following example:

vmd -dispdev text -psf system.psf -dcd sample_run.dcd -args seltext "not water" outfile zinfo_timeseries.txt histfile denprofile.txt first 100000 last 200000 freq 100 

Note that the -psf PSFfile topology description could be replaced with other filetype/filename options, e.g. -parm AMBERparmfile. Bin width and number are specified in the script, although they could be adjusted to be set at command line.

Basic and miscellaneous scripts

autocorr.awk AWK script to calculate autocorrelation of a time series.
histogram_dec.awk AWK script to create histogram.
histogram_discrete.awk AWK script to create histogram of discrete values (e.g. a variable taking integral values).
bootstrap_mean.awk AWK script to re-sample mean values in the bootstrap method.
runavg.awk AWK script to calculate running average of energy from NAMD output. Requires setting variables at command line.

Notes on uploading files

When uploading scripts, please remember that comments are essential to helping others understand your work. All uploaded files should include author information, a date or version number, a description of the file's purpose and methodology, a description of the expected input (if applicable), and a description of the output format.

To upload: Choose "Attachment" under the "Add" menu at right. Then, link to the attached file in this page using the following wiki syntax: [^filename.sh], followed by a description of the file.

To use: download the file. Note that some browsers will append unwanted extensions to downloaded files. Remember that all every user is responsible for understanding at a detailed level all the tools that he/she uses!