CPPTRAJ: Trajectory Analysis. V16.00b ___ ___ ___ ___ | \/ | \/ | \/ | _|_/\_|_/\_|_/\_|_ | Date/time: 06/21/16 15:48:44 | Available memory: 2.483 GB INPUT: Reading input from 'pt-analysis.chi-pu-py.in' [parm ../../aa/tip3p/ff12sb/run1/build/noWAt.topo.hmr [apa]] Reading '../../aa/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../ag/tip3p/ff12sb/run1/build/noWAt.topo.hmr [apg]] Reading '../../ag/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../ga/tip3p/ff12sb/run1/build/noWAt.topo.hmr [gpa]] Reading '../../ga/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../gg/tip3p/ff12sb/run1/build/noWAt.topo.hmr [gpg]] Reading '../../gg/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../ac/tip3p/ff12sb/run1/build/noWAt.topo.hmr [apc]] Reading '../../ac/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../au/tip3p/ff12sb/run1/build/noWAt.topo.hmr [apu]] Reading '../../au/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../gc/tip3p/ff12sb/run1/build/noWAt.topo.hmr [gpc]] Reading '../../gc/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../gu/tip3p/ff12sb/run1/build/noWAt.topo.hmr [gpu]] Reading '../../gu/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../ca/tip3p/ff12sb/run1/build/noWAt.topo.hmr [cpa]] Reading '../../ca/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../cg/tip3p/ff12sb/run1/build/noWAt.topo.hmr [cpg]] Reading '../../cg/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../ua/tip3p/ff12sb/run1/build/noWAt.topo.hmr [upa]] Reading '../../ua/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../ug/tip3p/ff12sb/run1/build/noWAt.topo.hmr [upg]] Reading '../../ug/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../cc/tip3p/ff12sb/run1/build/noWAt.topo.hmr [cpc]] Reading '../../cc/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../cu/tip3p/ff12sb/run1/build/noWAt.topo.hmr [cpu]] Reading '../../cu/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../uc/tip3p/ff12sb/run1/build/noWAt.topo.hmr [upc]] Reading '../../uc/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [parm ../../uu/tip3p/ff12sb/run1/build/noWAt.topo.hmr [upu]] Reading '../../uu/tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [trajin ../../ac/cluster-all/km-7/ctraj.c1 parm [apc]] Reading '../../ac/cluster-all/km-7/ctraj.c1' as Amber NetCDF [trajin ../../au/cluster-1-1/km-5/ctraj.c1 parm [apu]] Reading '../../au/cluster-1-1/km-5/ctraj.c1' as Amber NetCDF [trajin ../../gc/cluster-1-1/km-6-bbC2/ctraj.c1 parm [gpc]] Reading '../../gc/cluster-1-1/km-6-bbC2/ctraj.c1' as Amber NetCDF [trajin ../../gu/cluster-1-1/km-6/ctraj.c2 parm [gpu]] Reading '../../gu/cluster-1-1/km-6/ctraj.c2' as Amber NetCDF [dihedral c1 :1@O4' :1@C1' :1@N9 :1@C4 type chi] DIHEDRAL: [:1@O4']-[:1@C1']-[:1@N9]-[:1@C4] Output range is -180 to 180 degrees. [dihedral c2 :2@O4' :2@C1' :2@N1 :2@C2 type chi] DIHEDRAL: [:2@O4']-[:2@C1']-[:2@N1]-[:2@C2] Output range is -180 to 180 degrees. [run] ---------- RUN BEGIN ------------------------------------------------- PARAMETER FILES (16 total): 0: [apa] noWAt.topo.hmr, 65 atoms, 2 res, box: Orthogonal, 1 mol 1: [apg] noWAt.topo.hmr, 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [gpa] noWAt.topo.hmr, 66 atoms, 2 res, box: Orthogonal, 1 mol 3: [gpg] noWAt.topo.hmr, 67 atoms, 2 res, box: Orthogonal, 1 mol 4: [apc] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 5: [apu] noWAt.topo.hmr, 62 atoms, 2 res, box: Orthogonal, 1 mol 6: [gpc] noWAt.topo.hmr, 64 atoms, 2 res, box: Orthogonal, 1 mol 7: [gpu] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 8: [cpa] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 9: [cpg] noWAt.topo.hmr, 64 atoms, 2 res, box: Orthogonal, 1 mol 10: [upa] noWAt.topo.hmr, 62 atoms, 2 res, box: Orthogonal, 1 mol 11: [upg] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 12: [cpc] noWAt.topo.hmr, 61 atoms, 2 res, box: Orthogonal, 1 mol 13: [cpu] noWAt.topo.hmr, 60 atoms, 2 res, box: Orthogonal, 1 mol 14: [upc] noWAt.topo.hmr, 60 atoms, 2 res, box: Orthogonal, 1 mol 15: [upu] noWAt.topo.hmr, 59 atoms, 2 res, box: Orthogonal, 1 mol INPUT TRAJECTORIES (4 total): 0: 'ctraj.c1' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 109439 of 109439) 1: 'ctraj.c1' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 19943 of 19943) 2: 'ctraj.c1' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 24389 of 24389) 3: 'ctraj.c2' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 12720 of 12720) Coordinate processing will occur on 166491 frames. BEGIN TRAJECTORY PROCESSING: ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (2 actions): 0: [dihedral c1 :1@O4' :1@C1' :1@N9 :1@C4 type chi] [:1@O4'](1) [:1@C1'](1) [:1@N9](1) [:1@C4](1) 1: [dihedral c2 :2@O4' :2@C1' :2@N1 :2@C2 type chi] [:2@O4'](1) [:2@C1'](1) [:2@N1](1) [:2@C2](1) ----- ctraj.c1 (1-109439, 1) ----- 0% 10% 20% 30% 40% 50% 60% ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (2 actions): 0: [dihedral c1 :1@O4' :1@C1' :1@N9 :1@C4 type chi] [:1@O4'](1) [:1@C1'](1) [:1@N9](1) [:1@C4](1) 1: [dihedral c2 :2@O4' :2@C1' :2@N1 :2@C2 type chi] [:2@O4'](1) [:2@C1'](1) [:2@N1](1) [:2@C2](1) ----- ctraj.c1 (1-19943, 1) ----- 70% ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (2 actions): 0: [dihedral c1 :1@O4' :1@C1' :1@N9 :1@C4 type chi] [:1@O4'](1) [:1@C1'](1) [:1@N9](1) [:1@C4](1) 1: [dihedral c2 :2@O4' :2@C1' :2@N1 :2@C2 type chi] [:2@O4'](1) [:2@C1'](1) [:2@N1](1) [:2@C2](1) ----- ctraj.c1 (1-24389, 1) ----- 80% 90% ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (2 actions): 0: [dihedral c1 :1@O4' :1@C1' :1@N9 :1@C4 type chi] [:1@O4'](1) [:1@C1'](1) [:1@N9](1) [:1@C4](1) 1: [dihedral c2 :2@O4' :2@C1' :2@N1 :2@C2 type chi] [:2@O4'](1) [:2@C1'](1) [:2@N1](1) [:2@C2](1) ----- ctraj.c2 (1-12720, 1) ----- 100% Complete. Read 166491 frames and processed 166491 frames. TIME: Avg. throughput= 409527.6280 frames / second. ACTION OUTPUT: DATASETS (2 total): c1 "c1" (double, torsion(chin)), size is 166491 c2 "c2" (double, torsion(chin)), size is 166491 RUN TIMING: TIME: Init : 0.0001 s ( 0.01%) TIME: Trajectory Process : 0.4065 s ( 99.97%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 0.0000 s ( 0.00%) TIME: Data File Write : 0.0000 s ( 0.00%) TIME: Other : 0.0000 s ( 0.00%) TIME: Run Total 0.4066 s ---------- RUN END --------------------------------------------------- [hist c1,-180,180,0.5 c2,-180,180,0.5 out 2dhist_c1_c2_uy.gnu norm] Hist: 2dhist_c1_c2_uy.gnu: Set up for 2 dimensions using the following datasets: [ c1 c2 ] norm: Sum over bins will be normalized to 1.0. ---------- RUN BEGIN ------------------------------------------------- PARAMETER FILES (16 total): 0: [apa] noWAt.topo.hmr, 65 atoms, 2 res, box: Orthogonal, 1 mol 1: [apg] noWAt.topo.hmr, 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [gpa] noWAt.topo.hmr, 66 atoms, 2 res, box: Orthogonal, 1 mol 3: [gpg] noWAt.topo.hmr, 67 atoms, 2 res, box: Orthogonal, 1 mol 4: [apc] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 5: [apu] noWAt.topo.hmr, 62 atoms, 2 res, box: Orthogonal, 1 mol 6: [gpc] noWAt.topo.hmr, 64 atoms, 2 res, box: Orthogonal, 1 mol 7: [gpu] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 8: [cpa] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 9: [cpg] noWAt.topo.hmr, 64 atoms, 2 res, box: Orthogonal, 1 mol 10: [upa] noWAt.topo.hmr, 62 atoms, 2 res, box: Orthogonal, 1 mol 11: [upg] noWAt.topo.hmr, 63 atoms, 2 res, box: Orthogonal, 1 mol 12: [cpc] noWAt.topo.hmr, 61 atoms, 2 res, box: Orthogonal, 1 mol 13: [cpu] noWAt.topo.hmr, 60 atoms, 2 res, box: Orthogonal, 1 mol 14: [upc] noWAt.topo.hmr, 60 atoms, 2 res, box: Orthogonal, 1 mol 15: [upu] noWAt.topo.hmr, 59 atoms, 2 res, box: Orthogonal, 1 mol INPUT TRAJECTORIES (4 total): 0: 'ctraj.c1' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 109439 of 109439) 1: 'ctraj.c1' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 19943 of 19943) 2: 'ctraj.c1' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 24389 of 24389) 3: 'ctraj.c2' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 12720 of 12720) Coordinate processing will occur on 166491 frames. BEGIN TRAJECTORY PROCESSING: ----- ctraj.c1 (1-109439, 1) ----- 0% 10% 20% 30% 40% 50% 60% ----- ctraj.c1 (1-19943, 1) ----- 70% ----- ctraj.c1 (1-24389, 1) ----- 80% 90% ----- ctraj.c2 (1-12720, 1) ----- 100% Complete. Read 166491 frames and processed 166491 frames. TIME: Avg. throughput= 558880.2992 frames / second. ACTION OUTPUT: ANALYSIS: Performing 1 analyses: 0: [hist c1,-180,180,0.5 c2,-180,180,0.5 out 2dhist_c1_c2_uy.gnu norm] Calculating bins from min=-180 max=180 step=0.5. Dim c1: -180.000000->180.000000, step 0.500000, 720 bins. Calculating bins from min=-180 max=180 step=0.5. Dim c2: -180.000000->180.000000, step 0.500000, 720 bins. Hist: 166491 data points in each dimension. Hist: Allocating histogram, total bins = 518400 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 166491 TIME: Analyses took 0.0068 seconds. DATASETS (3 total): c1 "c1" (double, torsion(chin)), size is 166491 c2 "c2" (double, torsion(chin)), size is 166491 Hist_00018 "Hist_00018" (double matrix), size is 518400 DATAFILES (1 total): 2dhist_c1_c2_uy.gnu (Gnuplot File): Hist_00018 RUN TIMING: TIME: Init : 0.0000 s ( 0.00%) TIME: Trajectory Process : 0.2979 s ( 26.31%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 0.0068 s ( 0.60%) TIME: Data File Write : 0.8275 s ( 73.07%) TIME: Other : 0.0002 s ( 0.00%) TIME: Run Total 1.1324 s ---------- RUN END --------------------------------------------------- TIME: Total execution time: 1.6042 seconds. -------------------------------------------------------------------------------- To cite CPPTRAJ use: Daniel R. Roe and Thomas E. Cheatham, III, "PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data". J. Chem. Theory Comput., 2013, 9 (7), pp 3084-3095.