[gmx-users] REMD
merc mertens
merc.mertens at gmx.net
Thu Jul 13 14:50:08 CEST 2006
> > 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
--
"Feel free" – 10 GB Mailbox, 100 FreeSMS/Monat ...
Jetzt GMX TopMail testen: http://www.gmx.net/de/go/topmail
More information about the gromacs.org_gmx-users
mailing list