CPPTRAJ: Trajectory Analysis. V16.00b OpenMP ___ ___ ___ ___ | \/ | \/ | \/ | _|_/\_|_/\_|_/\_|_ | Date/time: 06/02/16 19:29:39 | Available memory: 16.002 GB INPUT: Reading input from 'pt-clust.in' [parm ../../tip3p/ff12sb/run1/build/full.topo.hmr [tip3p-ff12sb]] Reading '../../tip3p/ff12sb/run1/build/full.topo.hmr' as Amber Topology [parm ../../tip3p/dacP-ff12sb/run1/build/full.topo.hmr [tip3p-dacP-ff12sb]] Reading '../../tip3p/dacP-ff12sb/run1/build/full.topo.hmr' as Amber Topology [parm ../../tip3p/chen-garcia/run1/build/full.topo.hmr [tip3p-chen-garcia]] Reading '../../tip3p/chen-garcia/run1/build/full.topo.hmr' as Amber Topology [parm ../../tip3p/charmm36/run1/build/full.topo.hmr [tip3p-charmm36]] Reading '../../tip3p/charmm36/run1/build/full.topo.hmr' as Amber Topology CHAMBER topology: 1: CHARMM 361 * \\\\ CHARMM36 All-Hydrogen Nucleic Acid Topology File /// [parm ../../opc/ff12sb/run1/build/full.topo.hmr [opc-ff12sb]] Reading '../../opc/ff12sb/run1/build/full.topo.hmr' as Amber Topology [parm ../../opc/dacP-ff12sb/run1/build/full.topo.hmr [opc-dacP-ff12sb]] Reading '../../opc/dacP-ff12sb/run1/build/full.topo.hmr' as Amber Topology [trajin ../../tip3p/ff12sb/run1/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/ff12sb/run1/traj.1.02,../../tip3p/ff12sb/run1/traj.1.03,../../tip3p/ff12sb/run1/traj.1.04,../../tip3p/ff12sb/run1/traj.1.05,../../tip3p/ff12sb/run1/traj.1.06,../../tip3p/ff12sb/run1/traj.1.07,../../tip3p/ff12sb/run1/traj.1.08,../../tip3p/ff12sb/run1/traj.1.09,../../tip3p/ff12sb/run1/traj.1.10,../../tip3p/ff12sb/run1/traj.1.11,../../tip3p/ff12sb/run1/traj.1.12,../../tip3p/ff12sb/run1/traj.1.13,../../tip3p/ff12sb/run1/traj.1.14,../../tip3p/ff12sb/run1/traj.1.15,../../tip3p/ff12sb/run1/traj.1.16,../../tip3p/ff12sb/run1/traj.1.17,../../tip3p/ff12sb/run1/traj.1.18 parm [tip3p-ff12sb]] Reading '../../tip3p/ff12sb/run1/traj.1.01' as Amber NetCDF [trajin ../../tip3p/ff12sb/run2/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/ff12sb/run2/traj.1.02,../../tip3p/ff12sb/run2/traj.1.03,../../tip3p/ff12sb/run2/traj.1.04,../../tip3p/ff12sb/run2/traj.1.05,../../tip3p/ff12sb/run2/traj.1.06,../../tip3p/ff12sb/run2/traj.1.07,../../tip3p/ff12sb/run2/traj.1.08,../../tip3p/ff12sb/run2/traj.1.09,../../tip3p/ff12sb/run2/traj.1.10,../../tip3p/ff12sb/run2/traj.1.11,../../tip3p/ff12sb/run2/traj.1.12,../../tip3p/ff12sb/run2/traj.1.13,../../tip3p/ff12sb/run2/traj.1.14,../../tip3p/ff12sb/run2/traj.1.15,../../tip3p/ff12sb/run2/traj.1.16,../../tip3p/ff12sb/run2/traj.1.17,../../tip3p/ff12sb/run2/traj.1.18 parm [tip3p-ff12sb]] Reading '../../tip3p/ff12sb/run2/traj.1.01' as Amber NetCDF [trajin ../../tip3p/dacP-ff12sb/run1/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/dacP-ff12sb/run1/traj.1.02,../../tip3p/dacP-ff12sb/run1/traj.1.03,../../tip3p/dacP-ff12sb/run1/traj.1.04,../../tip3p/dacP-ff12sb/run1/traj.1.05,../../tip3p/dacP-ff12sb/run1/traj.1.06,../../tip3p/dacP-ff12sb/run1/traj.1.07,../../tip3p/dacP-ff12sb/run1/traj.1.08,../../tip3p/dacP-ff12sb/run1/traj.1.09,../../tip3p/dacP-ff12sb/run1/traj.1.10,../../tip3p/dacP-ff12sb/run1/traj.1.11,../../tip3p/dacP-ff12sb/run1/traj.1.12,../../tip3p/dacP-ff12sb/run1/traj.1.13,../../tip3p/dacP-ff12sb/run1/traj.1.14,../../tip3p/dacP-ff12sb/run1/traj.1.15,../../tip3p/dacP-ff12sb/run1/traj.1.16,../../tip3p/dacP-ff12sb/run1/traj.1.17,../../tip3p/dacP-ff12sb/run1/traj.1.18 parm [tip3p-dacP-ff12sb]] Reading '../../tip3p/dacP-ff12sb/run1/traj.1.01' as Amber NetCDF [trajin ../../tip3p/dacP-ff12sb/run2/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/dacP-ff12sb/run2/traj.1.02,../../tip3p/dacP-ff12sb/run2/traj.1.03,../../tip3p/dacP-ff12sb/run2/traj.1.04,../../tip3p/dacP-ff12sb/run2/traj.1.05,../../tip3p/dacP-ff12sb/run2/traj.1.06,../../tip3p/dacP-ff12sb/run2/traj.1.07,../../tip3p/dacP-ff12sb/run2/traj.1.08,../../tip3p/dacP-ff12sb/run2/traj.1.09,../../tip3p/dacP-ff12sb/run2/traj.1.10,../../tip3p/dacP-ff12sb/run2/traj.1.11,../../tip3p/dacP-ff12sb/run2/traj.1.12,../../tip3p/dacP-ff12sb/run2/traj.1.13,../../tip3p/dacP-ff12sb/run2/traj.1.14,../../tip3p/dacP-ff12sb/run2/traj.1.15,../../tip3p/dacP-ff12sb/run2/traj.1.16,../../tip3p/dacP-ff12sb/run2/traj.1.17,../../tip3p/dacP-ff12sb/run2/traj.1.18 parm [tip3p-dacP-ff12sb]] Reading '../../tip3p/dacP-ff12sb/run2/traj.1.01' as Amber NetCDF [trajin ../../tip3p/chen-garcia/run1/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/chen-garcia/run1/traj.1.02,../../tip3p/chen-garcia/run1/traj.1.03,../../tip3p/chen-garcia/run1/traj.1.04,../../tip3p/chen-garcia/run1/traj.1.05,../../tip3p/chen-garcia/run1/traj.1.06,../../tip3p/chen-garcia/run1/traj.1.07,../../tip3p/chen-garcia/run1/traj.1.08,../../tip3p/chen-garcia/run1/traj.1.09,../../tip3p/chen-garcia/run1/traj.1.10,../../tip3p/chen-garcia/run1/traj.1.11,../../tip3p/chen-garcia/run1/traj.1.12,../../tip3p/chen-garcia/run1/traj.1.13,../../tip3p/chen-garcia/run1/traj.1.14,../../tip3p/chen-garcia/run1/traj.1.15,../../tip3p/chen-garcia/run1/traj.1.16,../../tip3p/chen-garcia/run1/traj.1.17,../../tip3p/chen-garcia/run1/traj.1.18 parm [tip3p-chen-garcia]] Reading '../../tip3p/chen-garcia/run1/traj.1.01' as Amber NetCDF [trajin ../../tip3p/chen-garcia/run2/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/chen-garcia/run2/traj.1.02,../../tip3p/chen-garcia/run2/traj.1.03,../../tip3p/chen-garcia/run2/traj.1.04,../../tip3p/chen-garcia/run2/traj.1.05,../../tip3p/chen-garcia/run2/traj.1.06,../../tip3p/chen-garcia/run2/traj.1.07,../../tip3p/chen-garcia/run2/traj.1.08,../../tip3p/chen-garcia/run2/traj.1.09,../../tip3p/chen-garcia/run2/traj.1.10,../../tip3p/chen-garcia/run2/traj.1.11,../../tip3p/chen-garcia/run2/traj.1.12,../../tip3p/chen-garcia/run2/traj.1.13,../../tip3p/chen-garcia/run2/traj.1.14,../../tip3p/chen-garcia/run2/traj.1.15,../../tip3p/chen-garcia/run2/traj.1.16,../../tip3p/chen-garcia/run2/traj.1.17,../../tip3p/chen-garcia/run2/traj.1.18 parm [tip3p-chen-garcia]] Reading '../../tip3p/chen-garcia/run2/traj.1.01' as Amber NetCDF [trajin ../../tip3p/charmm36/run1/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/charmm36/run1/traj.1.02,../../tip3p/charmm36/run1/traj.1.03,../../tip3p/charmm36/run1/traj.1.04,../../tip3p/charmm36/run1/traj.1.05,../../tip3p/charmm36/run1/traj.1.06,../../tip3p/charmm36/run1/traj.1.07,../../tip3p/charmm36/run1/traj.1.08,../../tip3p/charmm36/run1/traj.1.09,../../tip3p/charmm36/run1/traj.1.10,../../tip3p/charmm36/run1/traj.1.11,../../tip3p/charmm36/run1/traj.1.12,../../tip3p/charmm36/run1/traj.1.13,../../tip3p/charmm36/run1/traj.1.14,../../tip3p/charmm36/run1/traj.1.15,../../tip3p/charmm36/run1/traj.1.16,../../tip3p/charmm36/run1/traj.1.17,../../tip3p/charmm36/run1/traj.1.18 parm [tip3p-charmm36]] Reading '../../tip3p/charmm36/run1/traj.1.01' as Amber NetCDF [trajin ../../tip3p/charmm36/run2/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../tip3p/charmm36/run2/traj.1.02,../../tip3p/charmm36/run2/traj.1.03,../../tip3p/charmm36/run2/traj.1.04,../../tip3p/charmm36/run2/traj.1.05,../../tip3p/charmm36/run2/traj.1.06,../../tip3p/charmm36/run2/traj.1.07,../../tip3p/charmm36/run2/traj.1.08,../../tip3p/charmm36/run2/traj.1.09,../../tip3p/charmm36/run2/traj.1.10,../../tip3p/charmm36/run2/traj.1.11,../../tip3p/charmm36/run2/traj.1.12,../../tip3p/charmm36/run2/traj.1.13,../../tip3p/charmm36/run2/traj.1.14,../../tip3p/charmm36/run2/traj.1.15,../../tip3p/charmm36/run2/traj.1.16,../../tip3p/charmm36/run2/traj.1.17,../../tip3p/charmm36/run2/traj.1.18 parm [tip3p-charmm36]] Reading '../../tip3p/charmm36/run2/traj.1.01' as Amber NetCDF [trajin ../../opc/ff12sb/run1/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../opc/ff12sb/run1/traj.1.02,../../opc/ff12sb/run1/traj.1.03,../../opc/ff12sb/run1/traj.1.04,../../opc/ff12sb/run1/traj.1.05,../../opc/ff12sb/run1/traj.1.06,../../opc/ff12sb/run1/traj.1.07,../../opc/ff12sb/run1/traj.1.08,../../opc/ff12sb/run1/traj.1.09,../../opc/ff12sb/run1/traj.1.10,../../opc/ff12sb/run1/traj.1.11,../../opc/ff12sb/run1/traj.1.12,../../opc/ff12sb/run1/traj.1.13,../../opc/ff12sb/run1/traj.1.14,../../opc/ff12sb/run1/traj.1.15,../../opc/ff12sb/run1/traj.1.16,../../opc/ff12sb/run1/traj.1.17,../../opc/ff12sb/run1/traj.1.18 parm [opc-ff12sb]] Reading '../../opc/ff12sb/run1/traj.1.01' as Amber NetCDF [trajin ../../opc/ff12sb/run2/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../opc/ff12sb/run2/traj.1.02,../../opc/ff12sb/run2/traj.1.03,../../opc/ff12sb/run2/traj.1.04,../../opc/ff12sb/run2/traj.1.05,../../opc/ff12sb/run2/traj.1.06,../../opc/ff12sb/run2/traj.1.07,../../opc/ff12sb/run2/traj.1.08,../../opc/ff12sb/run2/traj.1.09,../../opc/ff12sb/run2/traj.1.10,../../opc/ff12sb/run2/traj.1.11,../../opc/ff12sb/run2/traj.1.12,../../opc/ff12sb/run2/traj.1.13,../../opc/ff12sb/run2/traj.1.14,../../opc/ff12sb/run2/traj.1.15,../../opc/ff12sb/run2/traj.1.16,../../opc/ff12sb/run2/traj.1.17,../../opc/ff12sb/run2/traj.1.18 parm [opc-ff12sb]] Reading '../../opc/ff12sb/run2/traj.1.01' as Amber NetCDF [trajin ../../opc/dacP-ff12sb/run1/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../opc/dacP-ff12sb/run1/traj.1.02,../../opc/dacP-ff12sb/run1/traj.1.03,../../opc/dacP-ff12sb/run1/traj.1.04,../../opc/dacP-ff12sb/run1/traj.1.05,../../opc/dacP-ff12sb/run1/traj.1.06,../../opc/dacP-ff12sb/run1/traj.1.07,../../opc/dacP-ff12sb/run1/traj.1.08,../../opc/dacP-ff12sb/run1/traj.1.09,../../opc/dacP-ff12sb/run1/traj.1.10,../../opc/dacP-ff12sb/run1/traj.1.11,../../opc/dacP-ff12sb/run1/traj.1.12,../../opc/dacP-ff12sb/run1/traj.1.13,../../opc/dacP-ff12sb/run1/traj.1.14,../../opc/dacP-ff12sb/run1/traj.1.15,../../opc/dacP-ff12sb/run1/traj.1.16,../../opc/dacP-ff12sb/run1/traj.1.17,../../opc/dacP-ff12sb/run1/traj.1.18 parm [opc-dacP-ff12sb]] Reading '../../opc/dacP-ff12sb/run1/traj.1.01' as Amber NetCDF [trajin ../../opc/dacP-ff12sb/run2/traj.1.01 remdtraj remdtrajtemp 298.41 trajnames ../../opc/dacP-ff12sb/run2/traj.1.02,../../opc/dacP-ff12sb/run2/traj.1.03,../../opc/dacP-ff12sb/run2/traj.1.04,../../opc/dacP-ff12sb/run2/traj.1.05,../../opc/dacP-ff12sb/run2/traj.1.06,../../opc/dacP-ff12sb/run2/traj.1.07,../../opc/dacP-ff12sb/run2/traj.1.08,../../opc/dacP-ff12sb/run2/traj.1.09,../../opc/dacP-ff12sb/run2/traj.1.10,../../opc/dacP-ff12sb/run2/traj.1.11,../../opc/dacP-ff12sb/run2/traj.1.12,../../opc/dacP-ff12sb/run2/traj.1.13,../../opc/dacP-ff12sb/run2/traj.1.14,../../opc/dacP-ff12sb/run2/traj.1.15,../../opc/dacP-ff12sb/run2/traj.1.16,../../opc/dacP-ff12sb/run2/traj.1.17,../../opc/dacP-ff12sb/run2/traj.1.18 parm [opc-dacP-ff12sb]] Reading '../../opc/dacP-ff12sb/run2/traj.1.01' as Amber NetCDF [autoimage origin] AUTOIMAGE: To origin based on center of mass, anchor is first molecule. [strip :WAT,Na+,Cl-] STRIP: Stripping atoms in mask [:WAT,Na+,Cl-] [rmsd :1-2&!@H= first mass] RMSD: (:1-2&!@H*), reference is first frame (:1-2&!@H*), mass-weighted. Best-fit RMSD will be calculated, coords will be rotated and translated. [cluster :1,2@C5',C4',C3',O3',P,C2 summary clust-summary.dat repout rep repfmt pdb clusterout ctraj clusterfmt netcdf kmeans clusters 5 sieve 200 random savepairdist pairdist CpptrajPairDist summaryhalf split.dat splitframe 50000,100000,150000,200000,250000,300000,350000,400000,450000,498986,548986,598330] CLUSTER: Using coords dataset _DEFAULTCRD_, clustering using RMSD (mask [:1,2@C5',C4',C3',O3',P,C2]) best-fit K-MEANS: Looking for 5 clusters. Sequentially modify each point. Cluster to cluster distance will be based on cluster centroids. Initial clustering will be randomly sieved (with value 200). Previously calcd pair distances CpptrajPairDist will be used if found. Summary of cluster results will be written to clust-summary.dat Summary comparing parts of trajectory data for clusters will be written to split.dat Frames will be split at: 50000 100000 150000 200000 250000 300000 350000 400000 450000 498986 548986 598330 Cluster trajectories will be written to ctraj, format Amber NetCDF Cluster representatives will be written to separate trajectories, prefix (rep), format PDB Warning: One or more analyses requested creation of default COORDS DataSet. CREATECRD: Saving coordinates from Top full.topo.hmr to "_DEFAULTCRD_" ---------- RUN BEGIN ------------------------------------------------- PARAMETER FILES (6 total): 0: [tip3p-ff12sb] full.topo.hmr, 3302 atoms, 1088 res, box: Orthogonal, 1087 mol, 1075 solvent 1: [tip3p-dacP-ff12sb] full.topo.hmr, 3302 atoms, 1088 res, box: Orthogonal, 1087 mol, 1075 solvent 2: [tip3p-chen-garcia] full.topo.hmr, 3302 atoms, 1088 res, box: Orthogonal, 1087 mol, 1075 solvent 3: [tip3p-charmm36] full.topo.hmr, 3302 atoms, 1088 res, box: Orthogonal, 1087 mol, 1075 solvent 4: [opc-ff12sb] full.topo.hmr, 4377 atoms, 1088 res, box: Orthogonal, 1087 mol, 1075 solvent 5: [opc-dacP-ff12sb] full.topo.hmr, 4377 atoms, 1088 res, box: Orthogonal, 1087 mol, 1075 solvent INPUT TRAJECTORIES (12 total): 0: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 1: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 2: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 3: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 4: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 5: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 6: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 7: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 8: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 9: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 48989 of 48989) Looking for frames at 298.41 K 10: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 50000 of 50000) Looking for frames at 298.41 K 11: REMD trajectories (18 total), lowest replica 'traj.1.01' (reading 49346 of 49346) Looking for frames at 298.41 K Coordinate processing will occur on 598335 frames. BEGIN TRAJECTORY PROCESSING: ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 3236 atoms. Stripped topology: 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) Reference mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] Warning: COORDS data sets do not store temperatures. Warning: COORDS data sets do not store times. Estimated memory usage (100000 frames): 81.600 MB ----- traj.1.01 (1-50000, 1) ----- 0% ----- traj.1.01 (1-50000, 1) ----- 10% ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 3236 atoms. Stripped topology: 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] ----- traj.1.01 (1-50000, 1) ----- 20% ----- traj.1.01 (1-50000, 1) ----- 30% ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 3236 atoms. Stripped topology: 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] ----- traj.1.01 (1-50000, 1) ----- 40% ----- traj.1.01 (1-50000, 1) ----- 50% ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 3225 atoms. Stripped topology: 77 atoms, 13 res, box: Orthogonal, 12 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] WARNING: Could not set up actions for full.topo.hmr: skipping. ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 3225 atoms. Stripped topology: 77 atoms, 13 res, box: Orthogonal, 12 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] WARNING: Could not set up actions for full.topo.hmr: skipping. ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 4311 atoms. Stripped topology: 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] ----- traj.1.01 (1-50000, 1) ----- ----- traj.1.01 (1-48989, 1) ----- 60% ..................................................... ACTION SETUP FOR PARM 'full.topo.hmr' (4 actions): 0: [autoimage origin] Anchor molecule is 1 1086 molecules are mobile. 1: [strip :WAT,Na+,Cl-] Stripping 4311 atoms. Stripped topology: 66 atoms, 2 res, box: Orthogonal, 1 mol 2: [rmsd :1-2&!@H= first mass] Target mask: [:1-2&!@H*](42) 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. 3: [createcrd _DEFAULTCRD_] ----- traj.1.01 (1-50000, 1) ----- 70% ----- traj.1.01 (1-49346, 1) ----- 80% Read 498335 frames and processed 498335 frames. TIME: Avg. throughput= 80.3268 frames / second. ACTION OUTPUT: ANALYSIS: Performing 1 analyses: 0: [cluster :1,2@C5',C4',C3',O3',P,C2 summary clust-summary.dat repout rep repfmt pdb clusterout ctraj clusterfmt netcdf kmeans clusters 5 sieve 200 random savepairdist pairdist CpptrajPairDist summaryhalf split.dat splitframe 50000,100000,150000,200000,250000,300000,350000,400000,450000,498986,548986,598330] Starting clustering. Mask [:1,2@C5',C4',C3',O3',P,C2] corresponds to 11 atoms. Calculating pair-wise distances. Random_Number: seed is <= 0, using wallclock time as seed (1086860000) Estimated pair-wise matrix memory usage: > 12.415 MB Pair-wise matrix set up with sieve, 498335 frames, 2492 sieved frames. 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% Complete. Saving pair-wise distances to CpptrajPairDist Memory used by pair-wise matrix and other cluster data: 14.907 MB Round 0: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 0: 2487 points changed cluster assignment. Round 1: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 1: 232 points changed cluster assignment. Round 2: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 2: 206 points changed cluster assignment. Round 3: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 3: 24 points changed cluster assignment. Round 4: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 4: 4 points changed cluster assignment. Round 5: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 5: 2 points changed cluster assignment. Round 6: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete. K-means round 6: No change. Skipping the rest of the iterations. Restoring sieved frames. Parallelizing calculation with 8 threads 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% Complete. #Clustering: 5 clusters 498335 frames #Cluster 0 has average-distance-to-centroid 1.003100 #Cluster 1 has average-distance-to-centroid 0.779159 #Cluster 2 has average-distance-to-centroid 1.149443 #Cluster 3 has average-distance-to-centroid 1.455181 #Cluster 4 has average-distance-to-centroid 1.099065 #DBI: 1.510950 #pSF: 140593.544724 #Algorithm: Kmeans nclusters 5 maxit 100 #Representative frames: 55396 324547 392777 278909 221653 #Sieve value: 200 (random) #Sieved frames: 309 483 530 614 1065 1195 1367 1523 1632 1893 2124 2193 2529 2631 2770 3391 3403 3418 3934 4121 4169 4282 4681 4956 5037 5086 5131 5328 5576 6391 6469 6698 6789 7125 7160 7336 8050 8086 8633 8717 9092 9098 9510 9554 9613 9643 9798 9903 10114 10221 10237 10351 10392 10491 10552 10669 10750 11175 11222 11405 11482 12339 12423 12957 12965 12967 13610 14004 14090 14461 14528 14865 14994 15048 15414 15830 15838 15936 16405 16412 16474 17026 17037 17303 17366 17532 17622 17724 17905 18318 18388 18680 18881 18922 19043 19164 19747 19928 19997 20414 20415 20722 20738 20907 21018 21359 21518 21591 22049 22349 22414 22481 22659 22760 22880 22927 23029 23122 23224 23281 23407 23481 23505 23616 24046 24252 24570 24592 24652 24754 25006 25258 25616 25810 25814 26036 26127 26304 26943 27549 27746 27783 28295 28380 28430 28670 28931 29065 29643 30003 30014 30016 30115 30299 30452 30605 30927 30942 31771 32795 32916 33295 34338 34360 34367 34896 35886 36349 36396 36435 37274 37290 37572 37797 38189 38644 38715 38721 38815 38845 38891 38981 39257 39279 39847 39865 39900 40132 40647 40906 41149 41433 41649 41662 42002 42124 42242 42391 42516 42586 42765 43232 43442 43480 43636 43807 44135 44244 44266 44474 44743 44821 45327 45515 46171 46269 46478 46579 46749 46862 47364 47417 47443 47448 47709 47715 47717 47892 48157 48226 48305 48332 48572 48658 48738 48885 49008 49036 49254 49351 49606 49639 49695 49814 49861 50074 50530 50754 50864 51033 51097 51352 51817 51837 51879 52254 52277 52336 52617 52686 52741 53399 54581 54777 55210 55236 55315 55396 55895 55922 56211 56470 56811 56883 57054 57239 57265 57354 57375 57889 58217 58543 58628 58976 59248 60117 60890 60988 61019 61253 61468 61678 61696 61853 61918 62062 62174 62246 62437 62464 62475 62484 62817 63052 63217 63612 63625 63758 63840 64111 64182 64247 64383 64726 64768 64850 64894 64994 65131 65401 65451 65887 66476 66750 66772 67358 67618 67649 67694 67976 68440 68964 68981 69604 69608 69660 69801 69915 70152 70413 70798 71000 71084 71170 71223 71452 71455 71640 71664 72005 72160 72693 73045 73176 73188 73414 73755 74012 74216 74388 74390 74448 74538 74595 74600 74702 74872 75700 76068 76108 76149 76153 76214 76580 77078 77277 77720 78031 78267 78341 78417 78456 78589 78592 78836 79335 79374 79440 79624 79856 80324 80372 80632 80781 80970 81163 81315 81393 81498 81506 81744 81929 81984 81987 82165 82210 82441 82463 83048 83168 83479 83801 84033 84689 85195 85275 85374 85436 85538 85669 85811 86001 86015 86206 86389 86596 86737 86840 86886 87043 87484 87620 87724 87981 88587 88802 89046 89171 89277 89673 89745 89847 89878 89905 90333 90454 90481 90908 91168 91973 92263 92688 93617 94256 94552 94653 94656 95112 95399 95631 95675 95791 95823 95825 95863 95898 96011 96060 96498 96871 96902 96997 97054 97185 97231 97565 97583 97630 97632 97652 97813 97964 98258 98281 98490 98669 98684 98743 98752 99669 100041 100178 100508 100758 100810 100924 101217 101246 101284 101686 101739 101818 101886 102234 102548 102557 103016 103043 103645 103915 104473 104492 104816 104851 104944 105100 106567 106597 107056 107348 107548 108034 108227 108248 108525 108766 108773 108812 108813 108943 109012 109252 109281 109311 109735 109949 110253 110358 110497 110672 110675 110760 110828 110904 110932 111292 111493 111573 111614 112491 112518 112560 112640 113103 113155 113385 113539 113926 114177 114328 114404 114432 114636 114836 114905 115023 115438 115553 115746 116086 116661 116885 116951 117024 117146 117163 117443 118106 118178 118434 118461 118566 118709 118756 118895 119131 119163 119167 119172 119502 119755 120204 120328 120339 120774 121175 121238 121487 121671 122042 122207 122376 122391 122448 122621 122718 122721 123204 123369 124082 124154 124179 124392 124467 124533 124535 124745 124882 125175 125540 125961 126000 126153 126157 126161 126176 126609 126647 126787 126811 126814 126840 127112 127576 127960 127992 128207 128323 128511 128709 128774 128864 129445 129922 130431 130613 130634 130963 131054 131070 131789 131927 132054 132094 132138 132363 132411 132494 132733 133170 133474 133566 133851 134376 134394 135394 135969 136789 136924 137194 137637 137643 137723 138277 138351 138571 138743 139273 139390 139798 139963 140190 140211 140807 141286 141300 141616 142013 142040 142166 142257 142472 142523 142708 143021 143054 143166 143239 143389 143529 143608 143791 143876 144195 144280 144629 144904 145074 145113 145265 145291 145390 145836 145883 146083 146611 146912 147070 147456 147615 147819 148053 148436 148469 148580 148686 148727 148742 148767 149120 149264 149487 149652 149804 149924 150252 150752 151116 151198 151396 151558 151903 152337 152511 153094 153326 153364 153618 153708 153743 153879 153969 154094 154137 154273 154377 154406 154542 155145 155152 155400 155570 155929 155930 156105 156200 156502 156644 156846 156914 157128 157524 157529 157615 157646 157661 157668 157811 157985 158084 158469 158516 158687 158789 158827 159038 159187 159267 159278 159292 159326 159732 159763 159945 160152 160727 161131 161190 161251 161476 161513 161762 162046 162277 162816 162861 163241 164016 164230 164258 164658 165099 165102 165139 165414 165598 165602 165803 165908 166078 166111 166299 166329 166456 166620 166712 166823 167224 167380 167497 167653 167738 167754 168064 168123 168795 168927 168940 169387 169498 169622 169805 170016 170088 170435 170528 170564 171175 171246 171356 171633 171710 171947 172215 172407 172492 172530 172582 172687 172990 173023 173261 173530 173580 173695 173822 173993 174103 174327 174432 174516 175128 175345 175390 175429 175979 176173 176271 176328 176410 177030 177041 177424 177715 177925 178016 178272 178291 178735 178963 179131 179329 179363 179409 179431 179576 179577 179605 179868 180145 180217 180305 180397 180444 180499 180548 180709 180781 180921 180969 181212 181766 181871 182398 182466 183424 183682 183689 184014 184102 184159 184567 184678 184765 184970 185140 185440 185918 186044 186259 186764 186799 186930 187101 187138 187206 187260 187264 187296 187377 187417 187595 188377 188407 188764 188811 188858 188884 189097 189337 189584 189634 189652 189849 189850 190098 190257 190571 190611 190654 190681 191103 191227 191340 191596 192159 192176 192402 192454 192855 192877 192903 193561 193950 194590 194673 194712 195120 195399 195544 195546 195627 195713 195722 196224 196285 196379 196504 196512 196515 196550 196847 197095 197211 197277 197357 197415 197430 197540 197565 197743 198443 198552 198770 199149 199157 199158 199267 199273 199366 199541 199733 199983 200163 200307 200367 200378 200557 200641 201162 201487 201764 201940 201981 202887 203305 203378 203593 203711 203715 204475 204606 205183 205205 205283 205292 205808 205815 205991 206625 206668 206846 206874 207519 207816 208108 208420 208454 209162 209404 209797 210242 211818 211866 212078 212359 212494 212725 212844 212861 212942 213008 213025 213077 213239 213732 213965 214383 214742 215167 215445 215447 215507 215566 215861 215899 216134 216188 216369 216633 216674 216747 217027 217222 217297 217394 217526 217634 217962 218179 218182 218360 218501 218595 218633 218652 219136 219221 219254 219257 219262 219302 219386 219638 219831 219916 219958 220191 220315 220492 220606 220615 220812 221434 221506 221653 221730 222606 223081 223109 223163 223507 223635 224420 224533 225322 225696 225794 226194 226372 226503 227039 227194 227280 227663 228050 228635 228924 228936 229026 229188 229288 229386 229408 229505 229893 229895 230316 230381 230485 230595 230606 231193 231206 231260 232508 232700 232722 233003 233414 233939 234070 234250 234419 234765 234954 235000 235119 235191 235462 235579 235852 235937 236260 236296 236356 236736 237156 237159 237273 237766 237816 238078 238154 238522 238629 238822 239727 239757 239852 240388 240732 240765 240798 240977 241048 241121 241301 241420 241799 242210 242318 242383 242384 242457 242772 243032 243587 243696 243945 244074 244220 244246 244247 244544 244639 245257 245407 246202 246526 246918 247750 247841 248331 248599 249621 249700 249800 250138 250399 250532 250630 250781 250801 251093 251326 251358 251524 251652 251933 251960 252217 252352 252446 252545 252902 253083 253295 253331 253388 253429 253732 254654 254824 254875 256381 256776 257084 257684 257806 258083 258226 258392 258547 258809 259002 259049 259059 259663 259735 260052 260297 260620 261280 261363 261471 261869 262020 262110 262222 262284 262583 262755 262763 262764 263297 263464 263471 263527 263701 263711 263764 263860 263887 263982 264067 264245 264339 264614 264684 264822 264865 264898 265161 265719 265975 265978 266032 266140 266636 266875 267340 267471 267496 267535 268138 268434 269048 269177 269312 269387 269455 269681 270098 270891 271255 271266 271359 271595 272489 272598 272957 273061 273246 273701 273782 274003 274037 274105 274193 274837 274855 275020 275393 275643 275931 275952 276234 276278 276306 276504 276627 277600 277733 277859 278205 278242 278757 278801 278861 278909 279009 279119 279313 279321 279635 280011 280324 280470 280520 280792 280795 280910 281024 281077 281460 281799 281911 282216 282245 282284 282524 282533 282619 282983 283449 283985 284029 284626 284668 284795 284904 285282 285782 285936 286010 286221 286776 287265 287669 287850 288714 288899 288904 289075 289209 289351 289382 289534 289571 289704 290087 290379 290462 290664 290986 291195 291196 291277 291309 291314 291674 292141 292168 292207 292463 292481 292514 292548 292631 293046 293314 293586 293667 293686 293687 293851 293948 294280 294390 294525 294670 294762 294908 295199 295289 295325 295759 295808 296096 296113 296123 296163 296345 296507 296648 297126 297320 297342 297462 297533 297546 297801 297946 297967 298975 299179 299515 299658 299685 299954 300012 300168 300368 300809 300845 300940 302189 302553 302623 302854 303037 303125 303157 303181 303371 303585 304415 304590 304742 304743 304861 305134 305692 305785 305930 306026 306057 306078 306195 306327 306654 307093 307398 307511 307607 307931 307980 308043 308081 308379 308511 308636 308762 308816 309154 309344 309413 309905 309981 310003 310170 310217 310290 310331 310351 310933 310981 311109 311149 311586 311707 311761 312011 312087 312122 312365 312540 312876 313032 313251 313491 313875 314049 314069 314363 314461 314627 314686 314695 315349 315367 315637 315941 316063 316166 316204 316448 316727 316739 316923 316982 317041 317612 317725 317836 317847 317897 317909 317954 318012 318085 318166 318457 318468 318478 318611 318811 318908 319011 319115 319376 319638 319666 319759 319832 320161 320223 320249 320394 320523 320792 320856 320989 321003 321027 321089 321188 321279 321715 322244 322523 322717 322757 322883 322954 323586 323991 324023 324067 324135 324547 324886 324908 325278 325558 325645 326299 326495 326658 326719 326862 327724 327848 328050 328103 328216 328477 328815 329100 329121 329267 329270 329673 329731 329799 330332 330537 330811 330812 330972 331329 331740 332573 332668 332954 332995 333122 333177 333258 333278 333357 333363 333601 333907 333933 334036 334088 334611 334618 334744 334825 334894 335186 335202 335442 335743 336091 336372 336385 336462 336475 336532 336779 336921 336922 337087 337126 337453 337792 337803 338039 338095 338097 338276 338583 338603 339195 339225 339284 339511 340020 340247 340649 340869 340880 341519 341905 342008 342346 342472 342740 342996 343139 343412 343566 343576 344153 344170 344401 344492 345085 345250 345264 345293 345736 345975 346027 346487 346560 346803 347194 347200 347222 347316 347351 347527 347575 347615 347728 348074 348231 348272 348287 348384 348464 348525 348727 348926 349106 349267 349400 349676 349715 349765 349994 350110 351188 351346 352378 352412 352463 352715 352900 352962 353236 353260 353587 353849 353958 354320 354406 354493 354789 354847 354869 355278 355345 355444 355681 355776 355828 355886 356061 356368 356444 356466 356591 356881 356969 356975 357481 357648 358050 358834 358844 359231 359267 359490 359561 359904 360277 360343 360456 360568 360706 360872 360884 361061 361285 361620 361876 362024 362037 362136 362303 362376 362581 362762 363159 363601 363786 363817 363904 364211 364223 364249 364619 364722 365131 365278 365601 365770 365973 366130 366458 366543 366634 366732 366850 366973 366983 367320 367826 367954 368828 368951 369091 369193 369219 369221 369241 369353 369614 369698 369937 369963 370283 370336 370553 370565 370873 371021 371058 371074 371233 371580 371723 372075 372159 372713 373081 373390 373724 373946 374089 374458 374852 375607 376225 376245 377089 377593 377630 377821 377876 377930 378342 378353 378537 378593 378636 378796 378858 379499 379692 380263 380417 380570 380734 380890 381077 381815 382011 382064 382509 382865 383022 383936 384223 384234 384843 384854 384989 385236 385625 386232 386256 386376 386665 386705 387544 387773 387937 388254 388350 388490 388551 388863 389127 389723 389924 390154 390338 390422 390444 390515 390796 390882 390968 391160 391211 391295 391354 391521 391602 391793 391971 391984 392484 392521 392761 392777 392885 393103 393235 394223 394308 394459 394528 394541 394763 396184 396321 396336 396445 396898 397200 397233 397306 397308 397696 397789 397911 398205 398675 398765 398882 399390 399502 399619 400306 400467 400845 400980 401283 401510 401908 402259 402297 402326 402405 402935 403019 403025 403140 403277 403356 403720 403894 404149 404368 404562 404657 404762 404789 404943 405103 405322 405595 405620 406125 406135 406350 406370 406914 407085 407238 407442 407455 407582 407667 407688 407813 407814 407859 407861 407898 407920 408273 408470 409027 409504 409635 409984 410095 410102 410206 410363 410972 411093 411535 411776 413001 413201 413350 413505 413593 413686 413734 413878 413907 414395 414658 414846 415078 415265 415501 415542 415696 415851 415903 415913 416179 416784 417024 417243 417582 417613 417639 417675 417710 418078 418451 418650 418770 419063 419327 419725 419940 420221 421452 421563 421896 421936 422006 422021 422203 422308 422460 423021 423024 423027 423041 423311 423327 423440 423458 423481 423628 423632 423762 423828 423963 424177 424434 424772 424776 424786 424976 425187 425309 425570 425607 425639 425695 425825 425829 425872 426493 426778 427116 427157 427232 427767 427800 428304 428332 428373 428619 428815 428897 428985 429185 429651 429843 430248 431195 431795 431889 431969 432059 432440 432532 433289 433310 433727 433733 433989 434155 434406 434786 435318 435423 435501 435746 435825 435952 435994 436051 436466 436582 436806 436831 437114 437184 437212 437296 437716 437823 437851 437908 437909 438891 438894 439167 439334 439366 439403 440091 440124 440204 440249 440707 440737 440759 440817 440843 440937 441023 441096 441123 441269 441660 441713 441943 442018 442215 442306 442499 442524 442569 442705 442877 443083 443092 443555 443922 443942 443966 444449 444535 444556 444577 444783 444891 444949 445177 445299 445448 445484 445705 445715 446412 446550 446683 446781 446889 447338 447484 447713 447804 447833 447846 448119 448243 448670 448688 448842 449061 449481 449497 449557 449693 450038 450052 450227 450229 450238 450277 450428 450595 450597 450681 450789 451125 451661 451963 452072 452877 453024 453309 453468 453676 453709 453785 454215 454258 454871 455012 455013 455208 455209 455281 455390 455658 455874 456009 456024 456334 456708 456985 457100 457107 457217 457596 458304 458567 458695 458760 459347 459530 459640 459748 459885 460010 460042 460386 460482 460670 460761 460912 461192 461291 461346 461686 461739 462142 462392 462715 463004 463138 463336 463432 463481 463744 463749 463877 464067 464131 464250 464259 464404 464511 464632 464667 464854 464914 465100 465188 465284 465298 465514 465756 466262 466363 466773 467065 467079 467271 467811 467948 468011 468059 468304 468690 468952 469091 469161 469639 469794 470067 470310 470352 470402 470495 470496 470835 470872 470977 471318 471488 471713 471730 471918 471950 472064 472213 472666 472731 473305 473516 473682 473824 473876 474009 474129 474441 474459 474462 474580 475067 475161 475298 475517 475597 475993 476178 476223 476674 476862 477146 477396 477409 477507 477516 477819 478284 478419 478521 478804 478845 479097 479118 479403 479429 479627 479668 479774 480322 480349 480486 480885 480932 480966 481145 481164 481259 481580 481585 481723 481825 481866 481868 481928 481980 481984 482489 482592 482664 482678 482872 483100 483365 483607 483713 484059 484116 484299 484345 484466 484644 484656 484917 484963 484987 485243 485549 486186 486824 487293 488194 488432 488668 488771 489005 489077 489282 489395 489531 489564 489741 489755 489971 490182 490297 490592 490705 490720 490949 492339 492421 492479 492823 492954 493179 493431 493486 493678 493679 493728 493735 493816 493842 494032 494600 494755 495579 495837 496241 496259 496260 496357 496368 496589 496710 497527 497604 497712 497965 498018 498172 Warning: split frame 498986 is out of bounds; ignoring. Warning: split frame 548986 is out of bounds; ignoring. Warning: split frame 598330 is out of bounds; ignoring. Writing 'ctraj.c0' as Amber NetCDF Writing 'ctraj.c1' as Amber NetCDF Writing 'ctraj.c2' as Amber NetCDF Writing 'ctraj.c3' as Amber NetCDF Writing 'ctraj.c4' as Amber NetCDF Writing 'rep.c0.pdb' as PDB Warning: No PDB space group specified. Writing 'rep.c1.pdb' as PDB Warning: No PDB space group specified. Writing 'rep.c2.pdb' as PDB Warning: No PDB space group specified. Writing 'rep.c3.pdb' as PDB Warning: No PDB space group specified. Writing 'rep.c4.pdb' as PDB Warning: No PDB space group specified. Cluster timing data: TIME: Cluster Init. : 0.3795 s ( 0.10%) TIME: Pairwise Calc.: 5.1680 s ( 1.30%) TIME: Clustering : 0.3216 s ( 0.08%) TIME: Cluster Post. : 391.0580 s ( 98.52%) TIME: Total: 396.9272 s TIME: Analyses took 396.9369 seconds. DATASETS (3 total): RMSD_00006 "RMSD_00006" (double, rms), size is 498335 _DEFAULTCRD_ "_DEFAULTCRD_" (coordinates), size is 498335 (406.641 MB) Box Coords, 66 atoms Cnum_00008 "Cnum_00008" (integer), size is 498335 RUN TIMING: TIME: Init : 0.0018 s ( 0.00%) TIME: Trajectory Process : 6203.8430 s ( 93.99%) TIME: Action Post : 0.0000 s ( 0.00%) TIME: Analysis : 396.9369 s ( 6.01%) TIME: Data File Write : 0.0001 s ( 0.00%) TIME: Other : 0.0005 s ( 0.00%) TIME: Run Total 6600.7822 s ---------- RUN END --------------------------------------------------- TIME: Total execution time: 6645.1903 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.