CPPTRAJ Manual

User Manual:

Open the PDF directly: View PDF PDF.
Page Count: 207

DownloadCPPTRAJ Manual
Open PDF In BrowserView PDF
CPPTRAJ
Daniel R. Roe
February 22, 2019

https://github.com/Amber-MD/cpptraj

Contents
1 Introduction
1.1

6

Manual Syntax Format . . . . . . . . . . . . . . . . . . . . . . . .

2 Running Cpptraj

8

8

2.1

Command Line Syntax . . . . . . . . . . . . . . . . . . . . . . . .

8

2.2

Commands

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

9

2.3

Getting Help

. . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

2.4

Batch mode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

2.5

Interactive mode

. . . . . . . . . . . . . . . . . . . . . . . . . . .

11

2.6

Trajectory Processing Run . . . . . . . . . . . . . . . . . . . . .

11

2.6.1

11

2.7

Actions and multiple topologies . . . . . . . . . . . . . . .

Parallelization . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

12

2.7.1

MPI Trajectory Parallelization

12

2.7.2

OpenMP Parallelization . . . . . . . . . . . . . . . . . . .

12

2.7.3

CUDA Parallelization

13

. . . . . . . . . . . . . . .

. . . . . . . . . . . . . . . . . . . .

3 General Concepts

13

3.1

Units . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.2

Atom Mask Selection Syntax

. . . . . . . . . . . . . . . . . . . .

14

3.3

Ranges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

16

3.4

Parameter/Reference Tagging . . . . . . . . . . . . . . . . . . . .

16

4 Variables and Control Structures

13

17

4.1

for

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

17

4.2

set

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

18

4.3

show . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

19

5 Data Sets and Data Files

19

5.1

Data Set Selection Syntax . . . . . . . . . . . . . . . . . . . . . .

20

5.2

Data Set Math

21

. . . . . . . . . . . . . . . . . . . . . . . . . . . .

1

6 Data File Options

23

6.1

Standard Data File Options . . . . . . . . . . . . . . . . . . . . .

24

6.2

Grace Data File Options . . . . . . . . . . . . . . . . . . . . . . .

25

6.3

Gnuplot Data File Options

. . . . . . . . . . . . . . . . . . . . .

26

6.4

Amber REM Log Options . . . . . . . . . . . . . . . . . . . . . .

26

6.5

Amber MDOUT Options

27

6.6

Evecs File Options . . . . . . . . . . . . . . . . . . . . . . . . . .

27

6.7

Vector psuedo-traj Options

27

6.8

OpenDX le options

. . . . . . . . . . . . . . . . . . . . . . . . .

27

6.9

CCP4 le options . . . . . . . . . . . . . . . . . . . . . . . . . . .

28

. . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . .

6.10 Charmm REPD log options . . . . . . . . . . . . . . . . . . . . .

28

6.11 Amber Constant pH Out options . . . . . . . . . . . . . . . . . .

28

7 Coordinates (COORDS) Data Set Commands
7.1

combinecrd

7.2

28

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

29

crdaction

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

30

7.3

createcrd

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

30

7.4

crdout . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

30

7.5

loadcrd

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

30

7.6

loadtraj

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

7.7

permutedihedrals . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

7.8

reference . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

33

7.9

rotatedihedral . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

33

8 General Commands

33

8.1

activeref . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

34

8.2

calc

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

34

8.3

clear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

35

8.4

create

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

35

8.5

createset . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

35

8.6

datale

35

8.7

datalter

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

36

8.8

dataset . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

36

8.9

debug | prnlev . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

8.10 ensextension . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

8.11 exit | quit

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

8.12 go | run
8.13 help

8.14 list . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

8.15 noexitonerror . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

8.16 noprogress . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

40

8.17 precision . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

40

8.18 readdata . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

40

8.19 readensembledata . . . . . . . . . . . . . . . . . . . . . . . . . . .

41

8.20 readinput

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

8.21 removedata

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2

41
41

8.22 rst

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

8.23 runanalysis

41

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

42

8.24 select . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

43

8.25 selectds

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

43

8.26 sortensembledata . . . . . . . . . . . . . . . . . . . . . . . . . . .

43

8.27 write | writedata

. . . . . . . . . . . . . . . . . . . . . . . . . . .

43

8.28 System Commands . . . . . . . . . . . . . . . . . . . . . . . . . .

44

9 Topology File Commands

44

9.1

angleinfo | angles | printangles . . . . . . . . . . . . . . . . . . . .

45

9.2

atominfo | atoms | printatoms . . . . . . . . . . . . . . . . . . . .

45

9.3

bondinfo | bonds | printbonds . . . . . . . . . . . . . . . . . . . .

46

9.4

change . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

46

9.5

charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

46

9.6

comparetop

47

9.7

dihedralinfo | dihedrals | printdihedrals . . . . . . . . . . . . . . .

47

9.8

mass . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

47

9.9

molinfo

48

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

9.10 parm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.10.1 PDB format:

. . . . . . . . . . . . . . . . . . . . . . . . .

48
49

9.11 parmbox . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

49

9.12 parminfo

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

50

9.13 parmstrip

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

50

9.14 parmwrite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

50

9.15 resinfo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

51

9.16 scaledihedralk . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

51

9.17 solvent . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

52

10 Trajectory File Commands

52

10.1 ensemble . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

53

10.2 ensemblesize . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

54

10.3 reference . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

54

10.4 trajin

55

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10.4.1 Options for Amber NetCDF, Amber NC Restart, Amber

Restart: . . . . . . . . . . . . . . . . . . . . . . . . . . . .

57

10.4.2 Options for CHARMM DCD: . . . . . . . . . . . . . . . .

57

10.5 trajout . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

58

10.5.1 Options for pdb format:

. . . . . . . . . . . . . . . . . . .

10.5.2 Options for Amber ASCII format:
10.5.3 Options for Amber NetCDF format:

59

. . . . . . . . . . . . .

60

. . . . . . . . . . . .

60

10.5.4 Options for Amber Restart/NetCDF Restart format:

. . .

60

10.5.5 Options for CHARMM DCD: . . . . . . . . . . . . . . . .

61

10.5.6 Options for GROMACS TRX/XTC format: . . . . . . . .

61

10.5.7 Options for mol2 format:

. . . . . . . . . . . . . . . . . .

61

10.5.8 Options for SQM input format: . . . . . . . . . . . . . . .

61

3

11 Action Commands

62

11.1 angle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

65

11.2 areapermol

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

65

11.3 atomiccorr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

66

11.4 atomicuct | rmsf . . . . . . . . . . . . . . . . . . . . . . . . . . .

66

11.5 atommap

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

68

11.6 autoimage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

69

11.7 average

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

69

11.8 avgcoord . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

70

11.9 bounds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

71

11.10box . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

71

11.11center

72

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

11.12check | checkoverlap | checkstructure . . . . . . . . . . . . . . . .

72

11.13checkchirality . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

73

11.14closest | closestwaters . . . . . . . . . . . . . . . . . . . . . . . . .

74

11.15cluster . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

75

11.16clusterdihedral

. . . . . . . . . . . . . . . . . . . . . . . . . . . .

75

11.17contacts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

76

11.18createcrd

77

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

11.19createreservoir

. . . . . . . . . . . . . . . . . . . . . . . . . . . .

77

11.20density . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

78

11.21diusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

79

11.22dihedral

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

81

11.23dihedralscan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

81

11.24dipole

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

11.25distance

81

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

81

11.26drms | drmsd (distance RMSD) . . . . . . . . . . . . . . . . . . .

82

11.27dssp

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

84

11.28energy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

84

11.29esander

86

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

11.30lter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

87

11.31xatomorder

. . . . . . . . . . . . . . . . . . . . . . . . . . . . .

87

11.32ximagedbonds . . . . . . . . . . . . . . . . . . . . . . . . . . . .

88

11.33gist (Grid Inhomogeneous Solvation Theory)
11.34grid

. . . . . . . . . . .

88

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

97

11.35hbond

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

11.36image

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103

99

11.37jcoupling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
11.38lessplit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
11.39lie

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106

11.40lipidorder
11.41lipidscd

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108

11.42makestructure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
11.43mask . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
11.44matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
11.45mindist

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113

4

11.46minimage
11.47molsurf

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113

11.48multidihedral
11.49multivector

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115

11.50nastruct . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
11.51nativecontacts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
11.52outtraj . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
11.53pairdist

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125

11.54pairwise . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125
11.55principal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126
11.56projection . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 127
11.57pucker . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128
11.58radgyr | rog . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
11.59radial | rdf . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
11.60randomizeions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
11.61replicatecell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
11.62rms | rmsd . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
11.63rms2d | 2drms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
11.64rmsavgcorr

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135

11.65rmsf | atomicuct . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
11.66rotate

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 136

11.67rotdif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
11.68runavg | runningaverage . . . . . . . . . . . . . . . . . . . . . . . 137
11.69scale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
11.70secstruct . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
11.71spam . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
11.72setvelocity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
11.73stfcdiusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
11.74strip . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
11.75surf

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143

11.76symmrmsd

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144

11.77temperature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
11.78trans | translate . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
11.79unstrip . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
11.80unwrap
11.81vector

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147

11.82velocityautocorr . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
11.83volmap . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 150
11.84volume . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151
11.85watershell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151

12 Analysis Commands

152

12.1 autocorr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
12.2 avg . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
12.3 calcstate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
12.4 cluster . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157

5

12.5 cphstats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
12.6 corr | correlationcoe
12.7 crank | crankshaft
12.8 crduct

. . . . . . . . . . . . . . . . . . . . . . . . . 167
. . . . . . . . . . . . . . . . . . . . . . . . . . 168

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168

12.9 crosscorr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
12.10curvet

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 169

12.11diagmatrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
12.12divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
12.13t

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172

12.14hist | histogram . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
12.15integrate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 174
12.16ired

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175

12.17kde . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 176
12.18lifetime

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 177

12.19lowestcurve
12.20meltcurve
12.21modes

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180

12.22multicurve . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 184
12.23multihist . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
12.24phipsi

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185

12.25regress . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
12.26remlog . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 187
12.27rms2d | 2drms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 188
12.28rmsavgcorr

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189

12.29rotdif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 190
12.30runningavg
12.31spline

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 194

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 194

12.32statistics | stat

. . . . . . . . . . . . . . . . . . . . . . . . . . . . 195

12.32.1 Torsion Analysis
12.32.2 Distance Analysis

. . . . . . . . . . . . . . . . . . . . . . . 196
. . . . . . . . . . . . . . . . . . . . . . 196

12.32.3 Pucker Analysis . . . . . . . . . . . . . . . . . . . . . . . . 197
12.33ti . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
12.34timecorr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 198
12.35vectormath
12.36wavelet

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 199

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 200

13 Analysis Examples

203

13.1 Cartesian covariance matrix calculation and projection (PCA) . . 203
13.2 Dihedral covariance matrix calculation and projection for backbone phi/psi (PCA)

1

. . . . . . . . . . . . . . . . . . . . . . . . . 204

Introduction

Cpptraj [1] (the successor to ptraj ) is the main program in Amber for processing
coordinate trajectories and data les. Cpptraj has a wide range of functionality,

6

and makes use of OpenMP/MPI to speed up many calculations, including processing ensembles of trajectories and/or conducting multiple analyses in parallel
with MPI.[2]
Here are several notable features of cpptraj :
1. Trajectories with dierent topologies can be processed in the same run.
2. Several actions/analyses in cpptraj are OpenMP parallelized; see section
2.7.2 for more details.
3. Trajectory and ensemble reads can be MPI parallelized.
4. Almost any le read or written by cpptraj can be compressed (with the exception of the NetCDF trajectory format). So for example gzipped/bzipped
topology les can be read, and data les can be written out as gzip/bzip2
les. Compression is detected automatically when reading, and is determined by the lename extension (.gz and .bz2 respectively) on writing.
5. The format of output data les can be specied by extension. For example,
data les can be written in xmgrace format if the lename given has a
'.agr' extension. A trajectory can be written in DCD format if the '.dcd'
extension is used.
6. Multiple output trajectories can be specied, and can be written during
action processing (as opposed to only after) via the

outtraj

command.

In addition, output les can be directed to write only specic frames from
the input trajectories.
7. Multiple reference structures can be specied. Specic frames from trajectories may be used as a reference structure.
8. The

rmsd action allows specication of a separate mask for the reference

structure. In addition, per-residue RMSD can be calculated easily.
9. Actions that modify coordinates and topology such as the

strip /closest

actions can often write an accompanying fully-functional stripped topology
le.
10. Users usually are able to ne-tune the output format of data les declared
in actions using the 

out

keyword (for example, the precision of the

numbers can be changed). In addition, users can control which data sets
are written to which les (e.g.

if two actions specify the same data le

with the 'out' keyword, data from both actions will be written to that
data le).
11. Users can manipulate data sets using mathematical expressions (with some
limitations), see 5.2 on page 21 for details.
12. There is some support for creating internal loops over e.g. mask expressions and setting internal variables (see

7

for , set , and show commands).

1.1

Manual Syntax Format

The syntax presented in this manual uses the following conventions:

<>
[]

Denotes a variable.

Denotes something is optional.

{|}

Denotes several choices separated by the '|' character; one of the choices
must be specied.

...

Denotes the preceding option can be repeated.

Everything else is as printed.

2

Running Cpptraj

Cpptraj can be run in either interactive mode or in batch mode.

2.1

Command Line Syntax

cpptraj [-p ] [-i ] [-y ] [-x ]
[-ya ] [-xa ] []
[-c ] [-d ] [-w ] [-o ]
[-h | --help] [-V | --version] [--defines] [-debug <#>]
[--interactive] [--log ] [-tl]
[-ms ] [-mr ] [--mask ] [--resmask ]
* denotes a flag may be specified multiple times.

-p * Load  as a topology file.
-i * Read input from .
-y * Read from trajectory file ; same

as input 'trajin '.
-x * Write trajectory file ; same as
input 'trajout '.
-ya * Input trajectory file arguments.
-xa * Output trajectory file arguments.
<le>* A topology, input trajectory, or file
containing cpptraj input.
-c * Read  as reference
coordinates; same as input 'reference '.
-d * Read data in from file  ('readdata
').
-w  Write data from  as file
 ('writedata ).
8

-o  Write CPPTRAJ STDOUT output to file

.
-h | help Print command line help and exit.
-V | version Print version and exit.
denes Print compiler defines and exit.
-debug <#> Set global debug level to <#>; same as
input 'debug <#>'.
interactive Force interactive mode.
log  Record commands to  (interactive
mode only). Default is 'cpptraj.log'.
-tl Print length of trajectories specified with '-y' to
STDOUT. The total number of frames is written out as
'Frames: '
-ms  Print selected atom numbers to STDOUT.
Selected atoms are written out as 'Selected= 1 2 3
...'
-mr  : Print selected residue numbers to STDOUT.
Selected residues are written out as 'Selected= 1 2
3 ...'
mask  Print detailed atom selection to STDOUT.
resmask  : Print detailed residue selection to
STDOUT.
Note that unlike ptraj, in cpptraj it is not required that a topology le be
specied on the command line as long as one is specied in the input le with
the 'parm' keyword.

Multiple topology/input les can be specied by use of

multiple '-p' and '-i' ags. All topology and coordinate ags will be processed
before any input ags.

2.2

Commands

Input to cpptraj is in the form of commands, which can be categorized in to 2
types: immediate and queued. Immediate commands are executed as soon as
they are encountered. Queued commands are initialized when they are encoun-

run or go command.
reference ) are queued

tered, but are not executed until a Run is executed via a
Actions, Analyses, and Trajectory commands (except

commands; however, they can also be run immediately via commands such as

crdaction , runanalysis , loadcrd , etc. See 7 on page 28 for more details.
Commands fall into seven categories:

General

(Immediate) These commands are executed immediately when en-

tered.

System

(Immediate) These are unix system commands (e.g. 'ls', 'pwd', etc).

9

Coords

(Immediate) These commands are used to manipulate COORDS data

sets; see 7 on page 28 for more details.

Trajectory

(Queued) These commands prepare cpptraj for reading or writing

trajectories during a Run.

Topology

(Immediate) These commands are used to read, write, and modify

topology information.

Action

(Queued) These commands specify actions that will be performed on

coordinate frames read in from trajectories during a Run.

Analysis

(Queued) These commands specify analyses that will be performed

on data that has been either generated from a Run or read in from an
external source.

Control

(Immediate) These commands set up control blocks that can be used

to e.g. loop over a set of commands.
In addition to normal commands, cpptraj now has the ability to perform certain
basic math operations, even on data sets. See 5.2 on page 21 for more details.
Commands in cpptraj can be read in from an input le or from the interactive
command prompt. A '#' anywhere on a line denotes a comment; anything after
'#' will be ignored no matter where it occurs. A '\' allows the continuation of
one line to another. For example, the input:

# Sample input
trajin mdcrd # This is a trajectory
rms first out rmsd.dat \
:1-10
Translates to:

trajin mdcrd
rms first out rmsd.dat :1-10
2.3

Getting Help

If in interactive mode, 'help ' can be used to get the associated
keywords as well as an abbreviated description of the command.

Most com-

mands have a corresponding test which also serves as an example of how to use
the command.

See $AMBERHOME/AmberTools/test/cpptraj/README for

more details.

2.4

Batch mode

In batch mode, cpptraj is executed from the command line with one or more
input les containing commands to be processed or STDIN. The syntax of  is similar to that of ptraj. Keywords specifying dierent commands

10

are given one per line. Lines beginning with '#' are ignored as comments. Lines
can also be continued through use of the '\' character. This is the only allowed
mode for cpptraj.MPI.

2.5

Interactive mode

In interactive mode users can enter commands in a UNIX-like shell. Interactive
mode is useful for running short and simple analyses or for trying out new kinds
of analyses. If cpptraj is run with 'interactive', no arguments, or no specied
input le:

cpptraj
cpptraj --interactive
cpptraj 
cpptraj -p 
this brings up the interactive interface. This interface supports command history (via the up and down arrows) and tab completion for commands and le
names. If no log le name has been given (with 'log '), all commands
used in interactive mode will be logged to a le named 'cpptraj.log', which can
subsequently be used as input if desired.

When starting cpptraj, command

histories will be read from any existing logs.

2.6

Trajectory Processing Run

Like ptraj, a trajectory processing Run is one of the main ways to run cpp-

traj. First the Run is set up via commands read in from an input le or the
interactive prompt.

Trajectories are then read in one frame at a time (or in

the case of ensemble processing all frames from a given step are read). Actions
are performed on the coordinates stored in the frame, after which any output
coordinates are written.

At the end of the run, any data sets generated are

written, and any queued Analyses are performed.

2.6.1

Actions and multiple topologies

Since cpptraj supports multiple topology les, during a Run actions are set up
every time the topology changes in order to recalculate things like what atoms
are in a mask etc. Actions that are not valid for the current topology are skipped
for that topology. So for example given two topology les with 100 residues, if
the rst topology le processed includes a ligand named MOL and the second
one does not, the action:

distance :80 :MOL out D_80-to-MOL.dat
will be valid for the rst topology but not for the second, so it will be skipped
as long as the second topology is active.

11

2.7

Parallelization

Cpptraj has many levels of parallelization that can be enabled via the '-mpi',
'-openmp', and/or '-cuda' congure ags for MPI, OpenMP, and CUDA parallelization respectively. At the highest level, trajectory and ensemble reads are
parallelized with MPI. In addition, certain time consuming actions have been
parallelized with OpenMP and/or CUDA.
Note that any combination of the '-openmp', '-cuda', and '-mpi' ags may
be used to generate a hybrid MPI/OpenMP/CUDA binary; however this may
require additional runtime setup (e.g.

setting OMP_NUM_THREADS for

OpenMP) to work properly and not oversubscribe cores.

2.7.1

MPI Trajectory Parallelization

Cpptraj has two levels of MPI parallelization for reading input trajectories.

trajin ' trajectory input, where the trajectory read is divided as

The rst is for '

evenly as possible among all input frames. For example, if given two trajectories
of 1000 frames each and 4 MPI threads, thread 0 reads frames 1-500 of trajectory
1, thread 1 reads frames 501-1000 of trajectory 1, thread 2 reads frames 1-500
of trajectory 2, and thread 3 reads frames 501-1000 of trajectory 2.
The second is for

'ensemble' trajectory input, where the reading/processing/writing

of each member of the ensemble is divided up among MPI threads. The number of MPI threads must be a multiple of the ensemble size. If the number of
threads is greater than the ensemble size then the processing of each ensemble
member will be divided among MPI threads. For example, given an ensemble
of 4 trajectories and 8 threads, threads 0 and 1 are assigned to the rst ensemble trajectory, threads 2 and 3 are assigned to the second ensemble trajectory,
and so on. When using ensemble mode in parallel it is recommended that the

ensemblesize

command be used prior to any ensemble command as this will

make set up far more ecient.
In order to use the MPI version, Amber/cpptraj should be congured with
the '-mpi' ag. You can tell if cpptraj has been compiled with MPI as it will
print 'MPI' in the title, and/or by calling 'cpptraj denes' and looking for
'-DMPI'.

2.7.2 OpenMP Parallelization
Some of the more time-consuming actions/analyses in cpptraj have been parallelized with OpenMP to take advantage of machines with multiple cores. In order to use OpenMP parallelization Amber/cpptraj should be congured with the
'-openmp' ag. You can easily tell if cpptraj has been compiled with OpenMP
as it will print 'OpenMP' in the title, and/or by calling 'cpptraj denes' and
looking for '-D_OPENMP'. The following actions/analyses have been OpenMP
parallelized:

2drms/rms2d
atomiccorr
12

checkstructure
closest
cluster (pair-wise distance calculation and sieved frame restore only)
dssp/secstruct
energy
gist (non-bonded calculation)
hbond
kde
lipidscd
mask (distance-based masks only)
matrix (coordinate covariance matrices only)
minimage
radial
replicatecell
rmsavgcorr
spam
surf
velocityautocorr
volmap
watershell
wavelet
By default OpenMP cpptraj will use all available cores. The number of OpenMP
threads can be controlled by setting the OMP_NUM_THREADS environment
variable.

2.7.3

CUDA Parallelization

Some time-consuming actions in cpptraj have been parallelized with CUDA
to take advantage of machines with NVIDIA GPUs.

In order to use CUDA

parallelization Amber/cpptraj should be congured with the '-cuda' ag. You
can easily tell if cpptraj has been compiled with CUDA as it will print 'CUDA'
and details on the current graphics device in the title, and/or by calling 'cpptraj
denes' and looking for '-DCUDA'. The following actions have been CUDA
parallelized:

closest
watershell

3
3.1

General Concepts
Units

Cpptraj uses the AKMA system of units.

The execption is time, which is

typically expressed in ps (except where noted).

13

3.2

Variable

Unit

Length

Angstrom

Energy

kcal/mol

Mass

AMU

Charge

electron

Time

ps (typically)

Force

kcal/mol*Angstrom

Atom Mask Selection Syntax

The mask syntax is similar to ptraj. Note that the characters ':', '@', and '*'
are reserved for masks and should not be used in output le or data set names.
All masks are case-sensitive.

Either names or numbers can be used.

Masks

can contain ranges (denoted with '-') and comma separated lists. The logical
operands '&' (and), '|' (or), and '!' (not) are also supported.
The syntax for elementary selections is the following:

@{atom numlist}

e.g. '@12,17', '@54-85', '@12,54-85,90'

@{atom namelist}

e.g. '@CA', '@CA,C,O,N,H'

@%{atom type name}

e.g. '@%CT'

@/{atom_element_name}
:{residue numlist}
:{residue namelist}
:/{chain id}

e.g. '@/N'

e.g. ':1-10', ':1,3,5', ':1-3,5,7-9'
e.g. ':LYS', ':ARG,ALA,GLY'

e.g. ':/B', ':/A,D'. Requires chain ID information be present in

the topology.

:;{pdb residue number}

e.g. ':;2-4,8'. Requires a PDB loaded as topology,

or Amber topology with embedded PDB information (see

^{molecule numlist}

?? on page ??).

e.g. '^1-10', ':23,84,111'



Selection by distance, see below.

Several wildcard characters are supported:

'*'

Zero or more characters.

'='
' ?'

Same as '*'
One character.

The wildcards can also be used with numbers or other mask characters, e.g.
':?0' means :10,20,30,40,50,60,70,80,90, ':*' means all residues and '@*' means
all atoms.
Compound expressions of the following type are allowed:

14

:{residue numlist | namelist}@{atom namelist | numlist}
and are processed as:

:{residue numlist | namelist} & @{atom namelist | numlist}
e.g. ':1-10@CA' is equivalent to :1-10 & @CA.
More examples:

:ALA,TRP

All alanine and tryptophan residues.

:5,10@CA

CA carbon in residues 5 and 10.

:*&!@H=

All non-hydrogen atoms (equivalent to "!@H=").

@CA,C,O,N,H
!@CA,C,O,N,H

All backbone atoms.
All non-backbone atoms (=sidechains for proteins only).

:1-500@O&!(:WAT|:LYS,ARG)

All backbone oxygens in residues 1-500 but

not in water, lysine or arginine residues.

^1-2:ASP

All residues named 'ASP' in the rst two molecules.

:/A,D@CA

All atoms named 'CA' in chains A and D.

Distance-based Masks
There are two very important things to keep in mind when using distance based
masks:
1. Distance-based masks that update each frame are currently only supported
by the

mask action.

mask action requires dening
reference ; distances are then calculated using the

2. Selection by distance for everything but the
a reference frame with

specied reference frame only. This reference frame can be changed using
the

activeref

command.

 expression followed by a
 followed by a  (which is in Angstroms).
The  consists of 2 characters: '<' (within) or '>' (without) followed by either '^' (molecules), ':' (residues), or '@' (atoms). For
example, '<:3.0' means residues within 3.0 Angstroms etc. For residue- and
The syntax for selection by distance is a

molecule-based distance selection, if any atom in that residue/molecule matches
the given distance criterion, the entire residue/molecule is selected.
In plain language, the entire distance mask can be read as Select

operator>  of .
:11-17<@2.4
15

', 'nastruct resrange ', 'rmsd perres range ');

ify a range and ',' to separate dierent ranges. For example 1-2,4-6,9 species
1 to 2, 4 to 6, and 9, i.e. '1 2 4 5 6 9'.

3.4

Parameter/Reference Tagging

Parameter and reference les may be 'tagged' (i.e. given a nickname); these tags
can then be used in place of the le name itself. A tag in cpptraj is recognized
by being bounded by brackets ('[' and ']'). This can be particularly useful when
reading in many parameter or reference les.

For example, when reading in

multiple reference structures:

trajin Test1.crd
reference 1LE1.NoWater.Xray.rst7 [xray]
reference Test1.crd lastframe [last]
reference Test2.crd 225 [open]
rms Xray ref [xray] :2-12@CA out rmsd.dat
rms Last ref [last] :2-12@CA out rmsd.dat
rms Open ref [open] :2-12@CA out rmsd.dat
This denes three reference structures and gives them tags [xray], [last], and
[open]. These reference structures can then be referred to by their tags instead
of their lenames by any action that uses reference structures (in this case the
RMSD action).
Similarly, this can be useful when reading in multiple parameter les:

parm tz2.ff99sb.tip3p.truncoct.parm7 [tz2-water]
parm tz2.ff99sb.mbondi2.parm7 [tz2-nowater]
trajin tz2.run1.explicit.nc parm [tz2-water]
reference tz2.dry.rst7 parm [tz2-nowater] [tz2]
rms ref [tz2] !(:WAT) out rmsd.dat
This denes two parm les and gives them tags [tz2-water] and [tz2-nowater],
then reads in a trajectory associated with one, and a reference structure associated with the other. Note that in the 'reference' command there are two tags;

16

the rst goes along with the 'parm' keyword and species what parameter le
the reference should use, the second is the tag given to the reference itself (as
in the previous example) and is referred to in the subsequent RMSD action.

4

Variables and Control Structures

As of version 18, CPPTRAJ has limited support for script variables and 'for'
loops. Script variables are referred to by a dollar sign ('$') prex and are replaced
when they are processed. Note that to use script variables in CPPTRAJ input
that is also inside e.g. a BASH script, they can be escaped with the '\' character,
e.g.

#!/bin/bash
TOP=MyTop.parm7
cpptraj < inmask  [parm  | parmindex <#> | <#>] ... |
=;[;][] ... }
END KEYWORD: 'done'
Available 'end OP'
: '<' '>'
Available 'increment OP' : '++', '--', '+=', '-='

atoms|residues|molecules|molrstres|mollastres  inmask 
Loop over atoms/residues/molecules/first residue in
molecules/last residue in molecules selected by the
given mask expression, set as script variable .

parm  | parmindex <#> <#> Select

topology that  should be based on (default
first topology).

=;[;][]
Loop over integer script variable  starting

17

from , optionally ending at , increment
by .
Create a for loop using one or more mask expressions or integers. Loops can be
nested, but currently inner loops cannot refer to output loop variables. Integer
loops may be used without an end condition, but in that case at least one
descriptor in the loop should have an end condition or refer to a mask. Loops
are ended by the

done keyword.

Note that unlike C-style for loops the variable

is not incremented on the nal execution, so the nal value of

 (or the last selection in the mask).



will be

For example:

for atoms A0 inmask :1-3@CA i=1;i++
distance d$i :TCS $A0 out $i.dat
done
This loops over all atoms in the mask expression ':1-3@CA' (all atoms named
CA in residues 1 to 3) and creates a variable named 'i' that starts from 1 and is
incremented by 1 each iteration. Inside the loop, the mask selection is referred to
by

$A0 and the integer by $i.

This is equivalent to doing 3 distance commands

like so:

distance d1 :TCS :1@CA out 1.dat
distance d2 :TCS :2@CA out 2.dat
distance d3 :TCS :3@CA out 3.dat
4.2

set

set {    |
  {atoms|residues|molecules} inmask 
[parm  | parmindex <#> | <#>]
  trajinframes }
Available  : '=', '+='

   Set or append a script
variable.

  {atoms|residues|molecules} inmask 
Set/append a script variable to/by the total number
of atoms/residues/molecules selected by given mask
expression.

parm  | parmindex <#> | <#> Topology to
which mask should correspond (default first).

  trajinframes Set/append a script

variable to/by the total number of frames in
trajectories currently loaded by trajin commands.

18

Set ( = '=') or append ( = '+=') a script variable. Script variables
are referred to in CPPTRAJ input by using a dollar sign '$' prex.
For example, the following input will print info for the last 10 atoms in a
topology to 'last10.dat':

set Natom = atoms inmask *
last10 = $Natom - 10
show
atoms "@$last10 - $Natom" out last10.dat
4.3

show

show
Show all current script variables and their values.

5

Data Sets and Data Files

In cpptraj, Actions and Analyses can generate one or more data sets which are
available for further processing. For example, the
a data set containing distances vs time.

distance

command creates

The data set can be named by the

user simply by specifying a non-keyword string as an additional argument. If
no name is given, a default one will be generated based on the action name and
data set number. For example:

distance d1-2 :1 :2 out d1-2.dat
will create a data set named d1-2. If a name is not specied, e.g.:

distance :1 :2 out d1-2.dat
the data set will be named Dis_00000.
Data les are created automatically by most commands, usually via the



out keyword.

and

create

Data les can also be explicitly created with the

write /writedata
readdata

commands. Data can also be read in from les via the

command. Cpptraj currently recognizes the formats listed in 1, although it cannot write in all formats. In addition, a data set must be valid for the data le
format. For example, 3D data (such as a grid) can be written to an OpenDX
format le but not a Grace format le.
The default le format is called 'Standard', which simply has data in columns,
like ptraj, although multiple data sets can be directed to the same output le.
The format of a le can be changed either by specifying a recognized keyword
(either on the command line itself or later via a 'datale' command) or by giving
the le an extension corresponding to te format, so 'lename.agr' will output
in Grace format, and 'lename.gnu' will output in Gnuplot contour, and so on.
The xmgrace/gnuplot output is particularly nice for the secstruct sumout and

19

Format

Filename Extensions

Keyword

Valid Dimensions
1D, 2D, 3D

Standard

.dat

dat

Grace

.agr, .xmgr

grace

1D

Gnuplot

.gnu

gnu

1D, 2D

Xplor

.xplor, .grid

xplor

3D

OpenDX

.dx

opendx

3D

Amber REM log

.log

remlog

-

R

Amber MDOUT

.mdout

mdout

-

Energy infor

Modes data set only

Amber Evecs

.evecs

evecs

Amber Constant pH output

.cpout

cpout

pH data only

Vector pseudo-traj

.vectraj

vectraj

Vector data set only.

W

Gromacs XVG

.xvg

xvg

-

R

CCP4

.ccp4

ccp4

3D

Charmm REPD log

.exch

charmmrepd

-

R

Charmm Output

.charmmout

charmmout

-

Energy infor

Table 1: DataFile formats recognized by cpptraj. 'Valid Dimensions' shows what
dimensions the format is valid for (e.g. you cannot write a 1D data set with
OpenDX format).

rmsd perresout les.

Additional options for data les can be found in 6 on

page 23.
Any action using the out keyword will allow data sets from separate commands to be written into the same le. For example, the commands:

dihedral phi :1@C :2@N :2@CA :2@C out phipsi.dat
dihedral psi :2@N :2@CA :2@C :3@N out phipsi.dat
will assign the phi and psi data sets generated from each action to the standard data output le phipsi.dat:

#Frame
5.1

phi

psi

Data Set Selection Syntax

Many analysis commands can be used to analyze multiple data sets. The general
format for selecting data sets is:

[]:
The '*' character can be used as a wild-card for entire names (no partial
matches).

ˆ :

The data set name, usually specied in the action (e.g.

'distance d0 @1 @2' the data set name is d0).

20

in

ˆ :

Optional; this is set for certain data sets internally in or-

der to easily select subsets of data.

The brackets are required.

For

example, when using 'hbond series', both solute-solute and solute-solvent
hydrogen bond time series may be generated. To select all solute-solute
hydrogen bonds one would use the aspect  [solutehb]; to select solutesolvent hydrogen bonds the aspect  [solventhb] would be used. Aspects
are hard-coded and are listed in the commands that use them.

ˆ :

Optional; for actions that generate many data sets (such as

'rmsd perres') an index is used. Depending on the action, the index may
correspond to atom #s, residue #s, etc. A number range (comma and/or
dash separated) may be used.
For example: to select all data sets with aspect  [shear] named NA_00000:

NA_00000[shear]
To select all data sets with aspect  [stagger] with any name, indices 1 and 3:

*[stagger]:1,3
In ensemble mode, data set selection has additional syntax:

[]:%
Where  is the ensemble member number starting from 0.

5.2

Data Set Math

As of version 15, cpptraj can perform basic math operations, even on data sets
(with some limitations). Currently recognized operations are:

Operation

Symbol

Minus

-

Plus

+

Divide

/

Multiply

*

Power

^

Negate

-

Assign

=

Several functions are also supported:

21

Function

Form

Square Root

sqrt()

Exponential

exp()

Natural Logarithm

ln()

Absolute Value

abs()

Sine

sin()

Cosine

cos()

Tangent

tan()

Summation

sum()

Average

avg()

Standard Deviation

stdev()

Minimum

min()

Maximum

max()

Numbers can be expressed in scientic notation using E notation, e.g. 1E-5
= 0.00001. The parser also recognizes PI as the number pi. Expressions can also
be enclosed in parentheses. So for example, the following expression is valid:

> 1 - ln(sin(PI/4) * 2)^2
Result: 0.879887
Results of numerical calculations like the above can be assigned to a variable
(essentially a data set of size 1) for use in subsequent calculations, e.g.

> R = 1 - ln(sin(PI/4) * 2)^2
Result stored in 'R'
> R + 1 Result: 1.879887
Data sets can be specied in expressions as well.

Currently data sets in an

expression must be of the same type and only 1D, 2D, and 3D data sets are
supported. Functions are applied to each member of the data set. So for example, given two 1D data sets of the same size named D0 and D1, the following
expression:

> D2 = sqrt( D0 ) + D1
would take the square root of each member of D0, add it to the corresponding
member of D1, and assign the result to D2.

The following table lists which

operations are valid for data set types. If a type is not listed it is not supported:

Data Set Type

Supported Ops

Supported Funcs

1D (integer, double, oat)

All

All

1D (vector)

+, -, *, /, =

None

2D (matrices)

+, -, /, *, =

sum, avg, stdev, min, max

3D (grids)

+, -, /, *, =

sum, avg, stdev, min, max

22

Notes
'*' is dot product

6

Data File Options

Data le output can be handled multiple ways in cpptraj. Output data les can
be created by Actions/Analyses/Commands, or can be explicitly created with

writedata ( 8.27 on page 43) or create ( 8.4 on page 35) commands. Reading
data from les is only done via the readdata command ( 8.18 on page 40).
In general, data les which have been declared with an

'out'

recognize data le write keywords on the same command line.
the

'time'

keyword will
For example,

argument can be passed directly to the output from a

distance

command:

distance d0 :1 :2 out d0.agr time 0.001
The data le format can be changed from standard implicitly by using specic
lename extensions or keywords. If the extension is not recognized or no keyword
is give the default format is 'Standard'. Keywords and extensions for data le
formats recognized by cpptraj are shown in 1.

Note that the use of certain

options may be restricted for certain data le formats. These options can also
be passed to data les via the

datale command ( 8.6 on page 35).

[]
[{xlabel|ylabel|zlabel}