CPPTRAJ: Trajectory Analysis. V16.00b ___ ___ ___ ___ | \/ | \/ | \/ | _|_/\_|_/\_|_/\_|_ | Date/time: 06/02/16 18:13:46 | Available memory: 3.756 GB INPUT: Reading input from 'pt-analysis1.in' [parm ../../../tip3p/ff12sb/run1/build/noWAt.topo.hmr] Reading '../../../tip3p/ff12sb/run1/build/noWAt.topo.hmr' as Amber Topology [trajin ../ctraj.c4] Reading '../ctraj.c4' as Amber NetCDF [reference ../../../A-form.pdb [a-form]] Reading '../../../A-form.pdb' as PDB Setting active reference for distance-based masks: 'A-form.pdb' [rmsd rms_2_aform :1,2&!@H= ref [a-form] mass] RMSD: (:1,2&!@H*), reference is "[a-form]:1" (:1,2&!@H*), mass-weighted. Best-fit RMSD will be calculated, coords will be rotated and translated. [distance dist_sugar_ox :1@O2' :2@O4'] DISTANCE: :1@O2' to :2@O4', center of mass. [distance dist_N3-O5p :1@N3 :1@O5'] DISTANCE: :1@N3 to :1@O5', center of mass. [distance COM :1@N1,C2,H2,N3,C4,C5,C6,N2,H21,H22,N7,C8,H8,N9,O6 :2@N1,C6,H6,C5,H5,C4,N4,H41,H42,N3,C2,O2] DISTANCE: :1@N1,C2,H2,N3,C4,C5,C6,N2,H21,H22,N7,C8,H8,N9,O6 to :2@N1,C6,H6,C5,H5,C4,N4,H41,H42,N3,C2,O2, center of mass. [dihedral g1 :1@O5' :1@C5' :1@C4' :1@C3' type gamma] DIHEDRAL: [:1@O5']-[:1@C5']-[:1@C4']-[:1@C3'] 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 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 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. [dihedral b2 :2@P :2@O5' :2@C5' :2@C4' type beta] DIHEDRAL: [:2@P]-[:2@O5']-[:2@C5']-[:2@C4'] Output range is -180 to 180 degrees. [dihedral g2 :2@O5' :2@C5' :2@C4' :2@C3' type gamma] DIHEDRAL: [:2@O5']-[:2@C5']-[:2@C4']-[:2@C3'] Output range is -180 to 180 degrees. [dihedral d2 :2@C5' :2@C4' :2@C3' :2@O3' type delta] DIHEDRAL: [:2@C5']-[:2@C4']-[:2@C3']-[:2@O3'] 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 (1 total): 0: noWAt.topo.hmr, 64 atoms, 2 res, box: Orthogonal, 1 mol INPUT TRAJECTORIES (1 total): 0: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 4050 of 4050) Coordinate processing will occur on 4050 frames. REFERENCE FRAMES (1 total): 0: [a-form]:1 Active reference frame for distance-based masks is '[a-form]:1' BEGIN TRAJECTORY PROCESSING: ..................................................... ACTION SETUP FOR PARM 'noWAt.topo.hmr' (14 actions): 0: [rmsd rms_2_aform :1,2&!@H= ref [a-form] mass] Target mask: [:1,2&!@H*](40) Reference mask: [:1,2&!@H*](40) Warning: Coordinates are being rotated and box coordinates are present. Warning: Unit cell vectors are NOT rotated; imaging will not be possible Warning: after the RMS-fit is performed. 1: [distance dist_sugar_ox :1@O2' :2@O4'] :1@O2' (1 atoms) to :2@O4' (1 atoms), imaged. 2: [distance dist_N3-O5p :1@N3 :1@O5'] :1@N3 (1 atoms) to :1@O5' (1 atoms), imaged. 3: [distance COM :1@N1,C2,H2,N3,C4,C5,C6,N2,H21,H22,N7,C8,H8,N9,O6 :2@N1,C6,H6,C5,H5,C4,N4,H41,H42,N3,C2,O2] :1@N1,C2,H2,N3,C4,C5,C6,N2,H21,H22,N7,C8,H8,N9,O6 (14 atoms) to :2@N1,C6,H6,C5,H5,C4,N4,H41,H42,N3,C2,O2 (12 atoms), imaged. 4: [dihedral g1 :1@O5' :1@C5' :1@C4' :1@C3' type gamma] [:1@O5'](1) [:1@C5'](1) [:1@C4'](1) [:1@C3'](1) 5: [dihedral d1 :1@C5' :1@C4' :1@C3' :1@O3' type delta] [:1@C5'](1) [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) 6: [dihedral e1 :1@C4' :1@C3' :1@O3' :2@P type epsilon] [:1@C4'](1) [:1@C3'](1) [:1@O3'](1) [:2@P](1) 7: [dihedral z1 :1@C3' :1@O3' :2@P :2@O5' type zeta] [:1@C3'](1) [:1@O3'](1) [:2@P](1) [:2@O5'](1) 8: [dihedral c1 :1@O4' :1@C1' :1@N9 :1@C4 type chi] [:1@O4'](1) [:1@C1'](1) [:1@N9](1) [:1@C4](1) 9: [dihedral a2 :1@O3' :2@P :2@O5' :2@C5' type alpha] [:1@O3'](1) [:2@P](1) [:2@O5'](1) [:2@C5'](1) 10: [dihedral b2 :2@P :2@O5' :2@C5' :2@C4' type beta] [:2@P](1) [:2@O5'](1) [:2@C5'](1) [:2@C4'](1) 11: [dihedral g2 :2@O5' :2@C5' :2@C4' :2@C3' type gamma] [:2@O5'](1) [:2@C5'](1) [:2@C4'](1) [:2@C3'](1) 12: [dihedral d2 :2@C5' :2@C4' :2@C3' :2@O3' type delta] [:2@C5'](1) [:2@C4'](1) [:2@C3'](1) [:2@O3'](1) 13: [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.c4 (1-4050, 1) ----- 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. Read 4050 frames and processed 4050 frames. TIME: Avg. throughput= 49422.7906 frames / second. ACTION OUTPUT: DATASETS (14 total): rms_2_aform "rms_2_aform" (double, rms), size is 4050 dist_sugar_ox "dist_sugar_ox" (double, distance), size is 4050 dist_N3-O5p "dist_N3-O5p" (double, distance), size is 4050 COM "COM" (double, distance), size is 4050 g1 "g1" (double, torsion(gamma)), size is 4050 d1 "d1" (double, torsion(delta)), size is 4050 e1 "e1" (double, torsion(epsilon)), size is 4050 z1 "z1" (double, torsion(zeta)), size is 4050 c1 "c1" (double, torsion(chin)), size is 4050 a2 "a2" (double, torsion(alpha)), size is 4050 b2 "b2" (double, torsion(beta)), size is 4050 g2 "g2" (double, torsion(gamma)), size is 4050 d2 "d2" (double, torsion(delta)), size is 4050 c2 "c2" (double, torsion(chin)), size is 4050 RUN TIMING: TIME: Init : 0.0000 s ( 0.06%) TIME: Trajectory Process : 0.0819 s ( 99.86%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 0.0000 s ( 0.00%) TIME: Data File Write : 0.0000 s ( 0.04%) TIME: Other : 0.0000 s ( 0.00%) TIME: Run Total 0.0821 s ---------- RUN END --------------------------------------------------- [avg rms_2_aform out avg-rmsd-all.dat] AVERAGE: Calculating average of 1 data sets. Data set base name 'AVERAGE_00016' Output to to 'avg-rmsd-all.dat' [hist rms_2_aform,0,5,0.1 out rmsd-hist.dat name RMSD norm] Hist: rmsd-hist.dat: Set up for 1 dimensions using the following datasets: [ rms_2_aform ] norm: Sum over bins will be normalized to 1.0. [hist COM,0,10,0.1 out COM-hist.dat name CoM norm] Hist: COM-hist.dat: Set up for 1 dimensions using the following datasets: [ COM ] norm: Sum over bins will be normalized to 1.0. [hist dist_sugar_ox,0,10,0.1 out distO2O4-hist.dat norm] Hist: distO2O4-hist.dat: Set up for 1 dimensions using the following datasets: [ dist_sugar_ox ] norm: Sum over bins will be normalized to 1.0. [hist dist_N3-O5p,0,10,0.1 out distN3-O5p-hist.dat norm] Hist: distN3-O5p-hist.dat: Set up for 1 dimensions using the following datasets: [ dist_N3-O5p ] norm: Sum over bins will be normalized to 1.0. [hist rms_2_aform,0,5,0.1 dist_sugar_ox,0,10,0.1 out 2dhist_rms-dist.gnu norm] Hist: 2dhist_rms-dist.gnu: Set up for 2 dimensions using the following datasets: [ rms_2_aform dist_sugar_ox ] norm: Sum over bins will be normalized to 1.0. [hist g1,-180,180,1 out dihed_g1-hist.dat name gamma_1 norm] Hist: dihed_g1-hist.dat: Set up for 1 dimensions using the following datasets: [ g1 ] norm: Sum over bins will be normalized to 1.0. [hist d1,-180,180,1 out dihed_d1-hist.dat name delta_1 norm] Hist: dihed_d1-hist.dat: Set up for 1 dimensions using the following datasets: [ d1 ] norm: Sum over bins will be normalized to 1.0. [hist e1,-180,180,1 out dihed_e1-hist.dat name epsilon_1 norm] Hist: dihed_e1-hist.dat: Set up for 1 dimensions using the following datasets: [ e1 ] norm: Sum over bins will be normalized to 1.0. [hist z1,-180,180,1 out dihed_z1-hist.dat name zeta_1 norm] Hist: dihed_z1-hist.dat: Set up for 1 dimensions using the following datasets: [ z1 ] norm: Sum over bins will be normalized to 1.0. [hist c1,-180,180,1 out dihed_c1-hist.dat name chi_1 norm] Hist: dihed_c1-hist.dat: Set up for 1 dimensions using the following datasets: [ c1 ] norm: Sum over bins will be normalized to 1.0. [hist a2,-180,180,1 out dihed_a2-hist.dat name alpha_2 norm] Hist: dihed_a2-hist.dat: Set up for 1 dimensions using the following datasets: [ a2 ] norm: Sum over bins will be normalized to 1.0. [hist b2,-180,180,1 out dihed_b2-hist.dat name beta_2 norm] Hist: dihed_b2-hist.dat: Set up for 1 dimensions using the following datasets: [ b2 ] norm: Sum over bins will be normalized to 1.0. [hist g2,-180,180,1 out dihed_g2-hist.dat name gamma_2 norm] Hist: dihed_g2-hist.dat: Set up for 1 dimensions using the following datasets: [ g2 ] norm: Sum over bins will be normalized to 1.0. [hist d2,-180,180,1 out dihed_d2-hist.dat name delta_2 norm] Hist: dihed_d2-hist.dat: Set up for 1 dimensions using the following datasets: [ d2 ] norm: Sum over bins will be normalized to 1.0. [hist c2,-180,180,1 out dihed_c2-hist.dat name chi_2 norm] Hist: dihed_c2-hist.dat: Set up for 1 dimensions using the following datasets: [ c2 ] norm: Sum over bins will be normalized to 1.0. ---------- RUN BEGIN ------------------------------------------------- PARAMETER FILES (1 total): 0: noWAt.topo.hmr, 64 atoms, 2 res, box: Orthogonal, 1 mol INPUT TRAJECTORIES (1 total): 0: 'ctraj.c4' is a NetCDF AMBER trajectory, Parm noWAt.topo.hmr (Orthogonal box) (reading 4050 of 4050) Coordinate processing will occur on 4050 frames. REFERENCE FRAMES (1 total): 0: [a-form]:1 Active reference frame for distance-based masks is '[a-form]:1' BEGIN TRAJECTORY PROCESSING: ----- ctraj.c4 (1-4050, 1) ----- 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. Read 4050 frames and processed 4050 frames. TIME: Avg. throughput= 53374.4514 frames / second. ACTION OUTPUT: ANALYSIS: Performing 16 analyses: 0: [avg rms_2_aform out avg-rmsd-all.dat] 1: rms_2_aform 1: [hist rms_2_aform,0,5,0.1 out rmsd-hist.dat name RMSD norm] Calculating bins from min=0 max=5 step=0.1. Dim rms_2_aform: 0.000000->5.000000, step 0.100000, 50 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 50 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 2: [hist COM,0,10,0.1 out COM-hist.dat name CoM norm] Calculating bins from min=0 max=10 step=0.1. Dim COM: 0.000000->10.000000, step 0.100000, 100 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 100 Warning: Frame 12 Coordinates out of bounds (-1) Warning: Frame 13 Coordinates out of bounds (-1) Warning: Frame 15 Coordinates out of bounds (-1) Warning: Frame 67 Coordinates out of bounds (-1) Warning: Frame 162 Coordinates out of bounds (-1) Warning: Frame 185 Coordinates out of bounds (-1) Warning: Frame 217 Coordinates out of bounds (-1) Warning: Frame 218 Coordinates out of bounds (-1) Warning: Frame 263 Coordinates out of bounds (-1) Warning: Frame 313 Coordinates out of bounds (-1) Warning: Frame 314 Coordinates out of bounds (-1) Warning: Frame 335 Coordinates out of bounds (-1) Warning: Frame 352 Coordinates out of bounds (-1) Warning: Frame 384 Coordinates out of bounds (-1) Warning: Frame 469 Coordinates out of bounds (-1) Warning: Frame 489 Coordinates out of bounds (-1) Warning: Frame 500 Coordinates out of bounds (-1) Warning: Frame 501 Coordinates out of bounds (-1) Warning: Frame 551 Coordinates out of bounds (-1) Warning: Frame 562 Coordinates out of bounds (-1) Warning: Frame 563 Coordinates out of bounds (-1) Warning: Frame 570 Coordinates out of bounds (-1) Warning: Frame 593 Coordinates out of bounds (-1) Warning: Frame 602 Coordinates out of bounds (-1) Warning: Frame 719 Coordinates out of bounds (-1) Warning: Frame 720 Coordinates out of bounds (-1) Warning: Frame 763 Coordinates out of bounds (-1) Warning: Frame 777 Coordinates out of bounds (-1) Warning: Frame 785 Coordinates out of bounds (-1) Warning: Frame 786 Coordinates out of bounds (-1) Warning: Frame 815 Coordinates out of bounds (-1) Warning: Frame 816 Coordinates out of bounds (-1) Warning: Frame 915 Coordinates out of bounds (-1) Warning: Frame 916 Coordinates out of bounds (-1) Warning: Frame 922 Coordinates out of bounds (-1) Warning: Frame 933 Coordinates out of bounds (-1) Warning: Frame 952 Coordinates out of bounds (-1) Warning: Frame 953 Coordinates out of bounds (-1) Warning: Frame 1063 Coordinates out of bounds (-1) Warning: Frame 1070 Coordinates out of bounds (-1) Warning: Frame 1080 Coordinates out of bounds (-1) Warning: Frame 1207 Coordinates out of bounds (-1) Warning: Frame 1262 Coordinates out of bounds (-1) Warning: Frame 1320 Coordinates out of bounds (-1) Warning: Frame 1358 Coordinates out of bounds (-1) Warning: Frame 1374 Coordinates out of bounds (-1) Warning: Frame 1375 Coordinates out of bounds (-1) Warning: Frame 1376 Coordinates out of bounds (-1) Warning: Frame 1395 Coordinates out of bounds (-1) Warning: Frame 1455 Coordinates out of bounds (-1) Warning: Frame 1466 Coordinates out of bounds (-1) Warning: Frame 1486 Coordinates out of bounds (-1) Warning: Frame 1488 Coordinates out of bounds (-1) Warning: Frame 1589 Coordinates out of bounds (-1) Warning: Frame 1649 Coordinates out of bounds (-1) Warning: Frame 1740 Coordinates out of bounds (-1) Warning: Frame 1749 Coordinates out of bounds (-1) Warning: Frame 1754 Coordinates out of bounds (-1) Warning: Frame 1758 Coordinates out of bounds (-1) Warning: Frame 1786 Coordinates out of bounds (-1) Warning: Frame 1795 Coordinates out of bounds (-1) Warning: Frame 1804 Coordinates out of bounds (-1) Warning: Frame 1926 Coordinates out of bounds (-1) Warning: Frame 2004 Coordinates out of bounds (-1) Warning: Frame 2041 Coordinates out of bounds (-1) Warning: Frame 2090 Coordinates out of bounds (-1) Warning: Frame 2094 Coordinates out of bounds (-1) Warning: Frame 2095 Coordinates out of bounds (-1) Warning: Frame 2096 Coordinates out of bounds (-1) Warning: Frame 2181 Coordinates out of bounds (-1) Warning: Frame 2182 Coordinates out of bounds (-1) Warning: Frame 2256 Coordinates out of bounds (-1) Warning: Frame 2276 Coordinates out of bounds (-1) Warning: Frame 2277 Coordinates out of bounds (-1) Warning: Frame 2278 Coordinates out of bounds (-1) Warning: Frame 2304 Coordinates out of bounds (-1) Warning: Frame 2305 Coordinates out of bounds (-1) Warning: Frame 2324 Coordinates out of bounds (-1) Warning: Frame 2358 Coordinates out of bounds (-1) Warning: Frame 2384 Coordinates out of bounds (-1) Warning: Frame 2418 Coordinates out of bounds (-1) Warning: Frame 2419 Coordinates out of bounds (-1) Warning: Frame 2483 Coordinates out of bounds (-1) Warning: Frame 2511 Coordinates out of bounds (-1) Warning: Frame 2512 Coordinates out of bounds (-1) Warning: Frame 2600 Coordinates out of bounds (-1) Warning: Frame 2653 Coordinates out of bounds (-1) Warning: Frame 2725 Coordinates out of bounds (-1) Warning: Frame 2903 Coordinates out of bounds (-1) Warning: Frame 2925 Coordinates out of bounds (-1) Warning: Frame 2942 Coordinates out of bounds (-1) Warning: Frame 2978 Coordinates out of bounds (-1) Warning: Frame 2981 Coordinates out of bounds (-1) Warning: Frame 2989 Coordinates out of bounds (-1) Warning: Frame 3120 Coordinates out of bounds (-1) Warning: Frame 3135 Coordinates out of bounds (-1) Warning: Frame 3188 Coordinates out of bounds (-1) Warning: Frame 3189 Coordinates out of bounds (-1) Warning: Frame 3191 Coordinates out of bounds (-1) Warning: Frame 3192 Coordinates out of bounds (-1) Warning: Frame 3259 Coordinates out of bounds (-1) Warning: Frame 3277 Coordinates out of bounds (-1) Warning: Frame 3294 Coordinates out of bounds (-1) Warning: Frame 3333 Coordinates out of bounds (-1) Warning: Frame 3335 Coordinates out of bounds (-1) Warning: Frame 3405 Coordinates out of bounds (-1) Warning: Frame 3511 Coordinates out of bounds (-1) Warning: Frame 3514 Coordinates out of bounds (-1) Warning: Frame 3537 Coordinates out of bounds (-1) Warning: Frame 3547 Coordinates out of bounds (-1) Warning: Frame 3608 Coordinates out of bounds (-1) Warning: Frame 3680 Coordinates out of bounds (-1) Warning: Frame 3743 Coordinates out of bounds (-1) Warning: Frame 3747 Coordinates out of bounds (-1) Warning: Frame 3754 Coordinates out of bounds (-1) Warning: Frame 3815 Coordinates out of bounds (-1) Warning: Frame 3849 Coordinates out of bounds (-1) Warning: Frame 4008 Coordinates out of bounds (-1) Warning: Frame 4014 Coordinates out of bounds (-1) Warning: Frame 4016 Coordinates out of bounds (-1) Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 3930 3: [hist dist_sugar_ox,0,10,0.1 out distO2O4-hist.dat norm] Calculating bins from min=0 max=10 step=0.1. Dim dist_sugar_ox: 0.000000->10.000000, step 0.100000, 100 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 100 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 4: [hist dist_N3-O5p,0,10,0.1 out distN3-O5p-hist.dat norm] Calculating bins from min=0 max=10 step=0.1. Dim dist_N3-O5p: 0.000000->10.000000, step 0.100000, 100 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 100 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 5: [hist rms_2_aform,0,5,0.1 dist_sugar_ox,0,10,0.1 out 2dhist_rms-dist.gnu norm] Calculating bins from min=0 max=5 step=0.1. Dim rms_2_aform: 0.000000->5.000000, step 0.100000, 50 bins. Calculating bins from min=0 max=10 step=0.1. Dim dist_sugar_ox: 0.000000->10.000000, step 0.100000, 100 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 5000 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 6: [hist g1,-180,180,1 out dihed_g1-hist.dat name gamma_1 norm] Calculating bins from min=-180 max=180 step=1. Dim g1: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 7: [hist d1,-180,180,1 out dihed_d1-hist.dat name delta_1 norm] Calculating bins from min=-180 max=180 step=1. Dim d1: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 8: [hist e1,-180,180,1 out dihed_e1-hist.dat name epsilon_1 norm] Calculating bins from min=-180 max=180 step=1. Dim e1: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 9: [hist z1,-180,180,1 out dihed_z1-hist.dat name zeta_1 norm] Calculating bins from min=-180 max=180 step=1. Dim z1: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 10: [hist c1,-180,180,1 out dihed_c1-hist.dat name chi_1 norm] Calculating bins from min=-180 max=180 step=1. Dim c1: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 11: [hist a2,-180,180,1 out dihed_a2-hist.dat name alpha_2 norm] Calculating bins from min=-180 max=180 step=1. Dim a2: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 12: [hist b2,-180,180,1 out dihed_b2-hist.dat name beta_2 norm] Calculating bins from min=-180 max=180 step=1. Dim b2: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 13: [hist g2,-180,180,1 out dihed_g2-hist.dat name gamma_2 norm] Calculating bins from min=-180 max=180 step=1. Dim g2: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 14: [hist d2,-180,180,1 out dihed_d2-hist.dat name delta_2 norm] Calculating bins from min=-180 max=180 step=1. Dim d2: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 15: [hist c2,-180,180,1 out dihed_c2-hist.dat name chi_2 norm] Calculating bins from min=-180 max=180 step=1. Dim c2: -180.000000->180.000000, step 1.000000, 360 bins. Hist: 4050 data points in each dimension. Hist: Allocating histogram, total bins = 360 Histogram: Normalizing sum of bin populations to 1.0 Sum over all bins is 4050 TIME: Analyses took 0.0017 seconds. DATASETS (36 total): rms_2_aform "rms_2_aform" (double, rms), size is 4050 dist_sugar_ox "dist_sugar_ox" (double, distance), size is 4050 dist_N3-O5p "dist_N3-O5p" (double, distance), size is 4050 COM "COM" (double, distance), size is 4050 g1 "g1" (double, torsion(gamma)), size is 4050 d1 "d1" (double, torsion(delta)), size is 4050 e1 "e1" (double, torsion(epsilon)), size is 4050 z1 "z1" (double, torsion(zeta)), size is 4050 c1 "c1" (double, torsion(chin)), size is 4050 a2 "a2" (double, torsion(alpha)), size is 4050 b2 "b2" (double, torsion(beta)), size is 4050 g2 "g2" (double, torsion(gamma)), size is 4050 d2 "d2" (double, torsion(delta)), size is 4050 c2 "c2" (double, torsion(chin)), size is 4050 AVERAGE_00016[avg] "AVERAGE_00016[avg]" (double), size is 1 AVERAGE_00016[sd] "AVERAGE_00016[sd]" (double), size is 1 AVERAGE_00016[ymin] "AVERAGE_00016[ymin]" (double), size is 1 AVERAGE_00016[ymax] "AVERAGE_00016[ymax]" (double), size is 1 AVERAGE_00016[yminidx] "AVERAGE_00016[yminidx]" (integer), size is 1 AVERAGE_00016[ymaxidx] "AVERAGE_00016[ymaxidx]" (integer), size is 1 AVERAGE_00016[names] "AVERAGE_00016[names]" (string), size is 1 RMSD "RMSD" (double), size is 50 CoM "CoM" (double), size is 100 Hist_00025 "Hist_00025" (double), size is 100 Hist_00026 "Hist_00026" (double), size is 100 Hist_00027 "Hist_00027" (double matrix), size is 5000 gamma_1 "gamma_1" (double), size is 360 delta_1 "delta_1" (double), size is 360 epsilon_1 "epsilon_1" (double), size is 360 zeta_1 "zeta_1" (double), size is 360 chi_1 "chi_1" (double), size is 360 alpha_2 "alpha_2" (double), size is 360 beta_2 "beta_2" (double), size is 360 gamma_2 "gamma_2" (double), size is 360 delta_2 "delta_2" (double), size is 360 chi_2 "chi_2" (double), size is 360 DATAFILES (16 total): avg-rmsd-all.dat (Standard Data File): AVERAGE_00016[avg] AVERAGE_00016[sd] AVERAGE_00016[ymin] AVERAGE_00016[ymax] AVERAGE_00016[yminidx] AVERAGE_00016[ymaxidx] AVERAGE_00016[names] rmsd-hist.dat (Standard Data File): RMSD COM-hist.dat (Standard Data File): CoM distO2O4-hist.dat (Standard Data File): Hist_00025 distN3-O5p-hist.dat (Standard Data File): Hist_00026 2dhist_rms-dist.gnu (Gnuplot File): Hist_00027 dihed_g1-hist.dat (Standard Data File): gamma_1 dihed_d1-hist.dat (Standard Data File): delta_1 dihed_e1-hist.dat (Standard Data File): epsilon_1 dihed_z1-hist.dat (Standard Data File): zeta_1 dihed_c1-hist.dat (Standard Data File): chi_1 dihed_a2-hist.dat (Standard Data File): alpha_2 dihed_b2-hist.dat (Standard Data File): beta_2 dihed_g2-hist.dat (Standard Data File): gamma_2 dihed_d2-hist.dat (Standard Data File): delta_2 dihed_c2-hist.dat (Standard Data File): chi_2 RUN TIMING: TIME: Init : 0.0000 s ( 0.00%) TIME: Trajectory Process : 0.0759 s ( 7.39%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 0.0017 s ( 0.17%) TIME: Data File Write : 0.9496 s ( 92.44%) TIME: Other : 0.0000 s ( 0.00%) TIME: Run Total 1.0272 s ---------- RUN END --------------------------------------------------- TIME: Total execution time: 1.6036 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.