CPPTRAJ: Trajectory Analysis. V16.00b ___ ___ ___ ___ | \/ | \/ | \/ | _|_/\_|_/\_|_/\_|_ | Date/time: 06/21/16 12:26:20 | Available memory: 3.438 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.c5 parm [apc]] Reading '../../ac/cluster-all/km-7/ctraj.c5' as Amber NetCDF [trajin ../../au/cluster-1-1/km-5/ctraj.c4 parm [apu]] Reading '../../au/cluster-1-1/km-5/ctraj.c4' as Amber NetCDF [trajin ../../gc/cluster-1-1/km-6-bbC2/ctraj.c4 parm [gpc]] Reading '../../gc/cluster-1-1/km-6-bbC2/ctraj.c4' as Amber NetCDF [trajin ../../gu/cluster-1-1/km-6/ctraj.c4 parm [gpu]] Reading '../../gu/cluster-1-1/km-6/ctraj.c4' 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. [dihedral d1 :1@C5' :1@C4' :1@C3' :1@O3' type delta] DIHEDRAL: [:1@C5']-[:1@C4']-[:1@C3']-[:1@O3'] Output range is -180 to 180 degrees. [dihedral e1 :1@C4' :1@C3' :1@O3' :2@P type epsilon] DIHEDRAL: [:1@C4']-[:1@C3']-[:1@O3']-[:2@P] Output range is -180 to 180 degrees. [dihedral z1 :1@C3' :1@O3' :2@P :2@O5' type zeta] DIHEDRAL: [:1@C3']-[:1@O3']-[:2@P]-[:2@O5'] Output range is -180 to 180 degrees. [dihedral a2 :1@O3' :2@P :2@O5' :2@C5' type alpha] DIHEDRAL: [:1@O3']-[:2@P]-[:2@O5']-[:2@C5'] 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.c5' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 43245 of 43245) 1: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 9975 of 9975) 2: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 4050 of 4050) 3: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 4722 of 4722) Coordinate processing will occur on 61992 frames. BEGIN TRAJECTORY PROCESSING: ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (6 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) 2: [dihedral d1 :1@C5' :1@C4' :1@C3' :1@O3' type delta] [:1@C5'](1) [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) 3: [dihedral e1 :1@C4' :1@C3' :1@O3' :2@P type epsilon] [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) [:2@P](1) 4: [dihedral z1 :1@C3' :1@O3' :2@P :2@O5' type zeta] [:1@C3'](1) [:1@O3'](1) [:2@P](1) [:2@O5'](1) 5: [dihedral a2 :1@O3' :2@P :2@O5' :2@C5' type alpha] [:1@O3'](1) [:2@P](1) [:2@O5'](1) [:2@C5'](1) ----- ctraj.c5 (1-43245, 1) ----- 0% 10% 20% 30% 40% 50% 60% ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (6 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) 2: [dihedral d1 :1@C5' :1@C4' :1@C3' :1@O3' type delta] [:1@C5'](1) [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) 3: [dihedral e1 :1@C4' :1@C3' :1@O3' :2@P type epsilon] [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) [:2@P](1) 4: [dihedral z1 :1@C3' :1@O3' :2@P :2@O5' type zeta] [:1@C3'](1) [:1@O3'](1) [:2@P](1) [:2@O5'](1) 5: [dihedral a2 :1@O3' :2@P :2@O5' :2@C5' type alpha] [:1@O3'](1) [:2@P](1) [:2@O5'](1) [:2@C5'](1) ----- ctraj.c4 (1-9975, 1) ----- 70% 80% ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (6 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) 2: [dihedral d1 :1@C5' :1@C4' :1@C3' :1@O3' type delta] [:1@C5'](1) [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) 3: [dihedral e1 :1@C4' :1@C3' :1@O3' :2@P type epsilon] [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) [:2@P](1) 4: [dihedral z1 :1@C3' :1@O3' :2@P :2@O5' type zeta] [:1@C3'](1) [:1@O3'](1) [:2@P](1) [:2@O5'](1) 5: [dihedral a2 :1@O3' :2@P :2@O5' :2@C5' type alpha] [:1@O3'](1) [:2@P](1) [:2@O5'](1) [:2@C5'](1) ----- ctraj.c4 (1-4050, 1) ----- 90% ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (6 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) 2: [dihedral d1 :1@C5' :1@C4' :1@C3' :1@O3' type delta] [:1@C5'](1) [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) 3: [dihedral e1 :1@C4' :1@C3' :1@O3' :2@P type epsilon] [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) [:2@P](1) 4: [dihedral z1 :1@C3' :1@O3' :2@P :2@O5' type zeta] [:1@C3'](1) [:1@O3'](1) [:2@P](1) [:2@O5'](1) 5: [dihedral a2 :1@O3' :2@P :2@O5' :2@C5' type alpha] [:1@O3'](1) [:2@P](1) [:2@O5'](1) [:2@C5'](1) ----- ctraj.c4 (1-4722, 1) ----- 100% Complete. Read 61992 frames and processed 61992 frames. TIME: Avg. throughput= 295563.1204 frames / second. ACTION OUTPUT: DATASETS (6 total): c1 "c1" (double, torsion(chin)), size is 61992 c2 "c2" (double, torsion(chin)), size is 61992 d1 "d1" (double, torsion(delta)), size is 61992 e1 "e1" (double, torsion(epsilon)), size is 61992 z1 "z1" (double, torsion(zeta)), size is 61992 a2 "a2" (double, torsion(alpha)), size is 61992 RUN TIMING: TIME: Init : 0.0001 s ( 0.03%) TIME: Trajectory Process : 0.2097 s ( 99.94%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 0.0000 s ( 0.00%) TIME: Data File Write : 0.0000 s ( 0.01%) TIME: Other : 0.0000 s ( 0.00%) TIME: Run Total 0.2099 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. [hist d1,-180,180,0.5 c2,-180,180,0.5 out 2dhist_d1_c2_uy.gnu norm] Hist: 2dhist_d1_c2_uy.gnu: Set up for 2 dimensions using the following datasets: [ d1 c2 ] norm: Sum over bins will be normalized to 1.0. [hist z1,-180,180,0.5 c2,-180,180,0.5 out 2dhist_z1_c2_uy.gnu norm] Hist: 2dhist_z1_c2_uy.gnu: Set up for 2 dimensions using the following datasets: [ z1 c2 ] norm: Sum over bins will be normalized to 1.0. [hist a2,-180,180,0.5 c2,-180,180,0.5 out 2dhist_a2_c2_uy.gnu norm] Hist: 2dhist_a2_c2_uy.gnu: Set up for 2 dimensions using the following datasets: [ a2 c2 ] norm: Sum over bins will be normalized to 1.0. [hist z1,-180,180,0.5 a2,-180,180,0.5 out 2dhist_z1_a2_uy.gnu norm] Hist: 2dhist_z1_a2_uy.gnu: Set up for 2 dimensions using the following datasets: [ z1 a2 ] 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.c5' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 43245 of 43245) 1: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 9975 of 9975) 2: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 4050 of 4050) 3: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 4722 of 4722) Coordinate processing will occur on 61992 frames. BEGIN TRAJECTORY PROCESSING: ----- ctraj.c5 (1-43245, 1) ----- 0% 10% 20% 30% 40% 50% 60% ----- ctraj.c4 (1-9975, 1) ----- 70% 80% ----- ctraj.c4 (1-4050, 1) ----- 90% ----- ctraj.c4 (1-4722, 1) ----- 100% Complete. Read 61992 frames and processed 61992 frames. TIME: Avg. throughput= 538279.1946 frames / second. ACTION OUTPUT: ANALYSIS: Performing 5 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: 61992 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 61992 1: [hist d1,-180,180,0.5 c2,-180,180,0.5 out 2dhist_d1_c2_uy.gnu norm] Calculating bins from min=-180 max=180 step=0.5. Dim d1: -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: 61992 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 61992 2: [hist z1,-180,180,0.5 c2,-180,180,0.5 out 2dhist_z1_c2_uy.gnu norm] Calculating bins from min=-180 max=180 step=0.5. Dim z1: -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: 61992 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 61992 3: [hist a2,-180,180,0.5 c2,-180,180,0.5 out 2dhist_a2_c2_uy.gnu norm] Calculating bins from min=-180 max=180 step=0.5. Dim a2: -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: 61992 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 61992 4: [hist z1,-180,180,0.5 a2,-180,180,0.5 out 2dhist_z1_a2_uy.gnu norm] Calculating bins from min=-180 max=180 step=0.5. Dim z1: -180.000000->180.000000, step 0.500000, 720 bins. Calculating bins from min=-180 max=180 step=0.5. Dim a2: -180.000000->180.000000, step 0.500000, 720 bins. Hist: 61992 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 61992 TIME: Analyses took 0.0371 seconds. DATASETS (11 total): c1 "c1" (double, torsion(chin)), size is 61992 c2 "c2" (double, torsion(chin)), size is 61992 d1 "d1" (double, torsion(delta)), size is 61992 e1 "e1" (double, torsion(epsilon)), size is 61992 z1 "z1" (double, torsion(zeta)), size is 61992 a2 "a2" (double, torsion(alpha)), size is 61992 Hist_00022 "Hist_00022" (double matrix), size is 518400 Hist_00023 "Hist_00023" (double matrix), size is 518400 Hist_00024 "Hist_00024" (double matrix), size is 518400 Hist_00025 "Hist_00025" (double matrix), size is 518400 Hist_00026 "Hist_00026" (double matrix), size is 518400 DATAFILES (5 total): 2dhist_c1_c2_uy.gnu (Gnuplot File): Hist_00022 2dhist_d1_c2_uy.gnu (Gnuplot File): Hist_00023 2dhist_z1_c2_uy.gnu (Gnuplot File): Hist_00024 2dhist_a2_c2_uy.gnu (Gnuplot File): Hist_00025 2dhist_z1_a2_uy.gnu (Gnuplot File): Hist_00026 RUN TIMING: TIME: Init : 0.0000 s ( 0.00%) TIME: Trajectory Process : 0.1152 s ( 2.19%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 0.0371 s ( 0.70%) TIME: Data File Write : 5.1173 s ( 97.09%) TIME: Other : 0.0011 s ( 0.00%) TIME: Run Total 5.2707 s ---------- RUN END --------------------------------------------------- TIME: Total execution time: 5.5456 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.