Dear Amber users,
I would like to compute a vector using ptraj and i would like to save the
vector as the file.
I used the following input for ptraj but it does not produces the output
file. Please help me in this issue.
the following are the input file information
ptraj.in
trajin dyn.nc 1 10
strip :WAT
rms first mass :1-285.CA
vector v0 .547 corr .1699 out 142-236-48.vec order 2
vector v1 .547,.1699,.2839 corrplane order 2
analyze timecorr vec1 v1 vec2 v1 tstep 1 tcorr 10 out 142-236-48.tcorr
--------------------------------------------------------------------------------------------------------------------------------
the following out put i received from ptraj
PTRAJ: rms first mass :1-285.CA
Mask [:1-285.CA] represents 279 atoms
PTRAJ: vector v0 .547 corr .1699 out 142-236-48.vec order 2
Mask [.547] represents 1 atoms
Mask [.1699] represents 1 atoms
PTRAJ: vector v1 .547,.1699,.2839 corrplane order 2
Mask [.547,.1699,.2839] represents 3 atoms
PTRAJ: analyze timecorr vec1 v1 vec2 v1 tstep 1 tcorr 10 out
142-236-48.tcorr
[No output trajectory specified (trajout)]
PTRAJ: Successfully read the input file.
Coordinate processing will occur on 10 frames.
Summary of I/O and actions follows:
INPUT COORDINATE FILES
File (/dyn.nc) is a NetCDF AMBER trajectory with box info with 10 sets
NO OUTPUT COORDINATE FILE WAS SPECIFIED
ACTIONS
1> STRIP: 60717 atoms will be removed from trajectory: :295-20533
2> RMS to first frame using mass weighting
Dumping RMSd vs. time (with time interval 1.00)
Atom selection follows :2-94.CA,:97-189.CA,:192-284.CA
3> VECTOR: storage to array named v0
Calculate the vector between the center of mass of the two atom
selections
which follow (with the origin at the center of mass of the first)
Atom selection 1 is :48.CA
Atom selection 2 is :142.CA
The order of Legendre polynomials is 2
Warning: Output of corr, ired, corrired or corrplane vectors is not
yet supported!
4> VECTOR: storage to array named v1
Calculate the vector perpendicular to the least squares best plane
through the atom selection which follows
Atom selection is :48.CA,:142.CA,:236.CA
The order of Legendre polynomials is 2
ANALYZE
1> ANALYZE TIMECORR: Calculating cross-correlation function of vector(s)
v1 v1,
for a correlation time of 10.000000, using a time step of 1.000000.
Corr. func. are not normalized
Corr. func. are calculated using the FFT approach
Results are written to 142-236-48.tcorr
Set 1 ..........
PTRAJ: Successfully read in 10 sets and processed 10 sets.
Dumping accumulated results (if any)
PTRAJ RMS: dumping RMSd vs time data
PTRAJ: Analyzing accumulated data
_______________________________________________________________________________________________________________-----
with regards
sundar
--
S.Sundar Raman
CSIR-SRF
Chemical Laboratory
Central Leather Research Institute
Adyar, Chennai, India- 600 020
&
DAAD Research Scholar
EML-Research,Villa Bosch,
Schloss-Wolfsbrunnenweg 33
D-69118 Heidelberg
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed May 20 2009 - 17:03:29 PDT
