[gmx-users] REMD
David van der Spoel
spoel at xray.bmc.uu.se
Thu Jul 13 21:35:57 CEST 2006
merc mertens wrote:
>>> remd.xvg was build like described above. output is only one trajectory,
>> containing only box information.
>>> additionally i tried:
>>>
>>> trjcat -f traj0.trr traj1.trr traj2.trr ... -o j0.xtc j1.xtc j2.xtc ...
>> -demux remd.xvg
>>> but from that i got only one trajectory as output, too (also only
>> containing box info).
>> How about the trr files? Are they written with the same frequency as
>> what you have in the remd.xvg file? Otherwise I don't know what will
>> happen...
>>
>> --
>> David.
>
> when trying: gmxdump -f traj0.trr | grep time
>
> sample output is:
> natoms= 14670 step= 0 time= 0 lambda= 0
> natoms= 14670 step= 500 time= 1 lambda= 0
> natoms= 14670 step= 1000 time= 2 lambda= 0
> natoms= 14670 step= 1500 time= 3 lambda= 0
> natoms= 14670 step= 2000 time= 4 lambda= 0
> natoms= 14670 step= 2500 time= 5 lambda= 0
> ...
>
> the "time" stamps correspond to the first column in remd.xvg used for trjcat. further, if *xtc is chosen as output in trjcat, for every single frame the *xtc file contains the line: "not available: x", where there surely is x in the single *trr files!
>
> e.g.: gmxdump_d -f traj0.trr
>
> traj0.trr frame 0:
> natoms= 14670 step= 0 time= 0 lambda= 0
> box (3x3):
> box[ 0]={ 8.48915e+00, 0.00000e+00, 0.00000e+00}
> box[ 1]={ 0.00000e+00, 4.99054e+00, 0.00000e+00}
> box[ 2]={ 0.00000e+00, 0.00000e+00, 3.49561e+00}
> x (14670x3):
> x[ 0]={ 1.98230e+00, 2.07120e+00, 3.06491e+00}
> x[ 1]={ 1.93097e+00, 2.10502e+00, 3.14505e+00}
> x[ 2]={ 1.92314e+00, 2.07599e+00, 2.98319e+00}
> ...
>
> for sure i am doing something wrong here, i just do not get what.
> i am puzzled, thanks for any help, merc
>
have you run
gmxcheck -f file.trr ?
--
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: 46 18 471 4205 fax: 46 18 511 755
spoel at xray.bmc.uu.se spoel at gromacs.org http://folding.bmc.uu.se
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
More information about the gromacs.org_gmx-users
mailing list