Mark Abraham wrote:
Vitaly V. Chaban wrote:
Mark,

In "g_energy -h" I see nothing abous viscosity except the option
"-vis". The command line is as I  wrote earlier:
"g_energy -vis -acflen 1001"
The system is liquid water. The simulation time is 2 000 ps. The file
I look into is enecorr.xvg.

What the error of some sort do you mean?

-acflen is measured in trajectory frames. How many relevant frames exist in the trajectory file?

Viscosity in g_energy is computed from the pressure tensor which is stored in the energy file. Unfortunately is seems like the acflen is ignored in this case. Please submit a bugzilla.

On the other hand you need all the data you can get for this quantity, because it converges extremely slowly.


Mark

Vitaly

When using g_energy to calculate viscosity ("g_energy -vis"), in
enecorr.xvg the ACF of stress tensor is output. How can one specify
the "acflen" of this ACF?
I try the key "g_energy -vis -acflen XXX" but the length of the ACF is
still equal to the half of the trajectory length.
If you read g_energy -h, you will need to specify the ACF length in a
particular way. If you'd given us your actual command line and some
information about your trajectory, you might have gotten a definitive
answer... as it is, we can only assume you made an error of some sort! :-)

_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


--
David van der Spoel, Ph.D., Professor of Biology
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
sp...@xray.bmc.uu.se    sp...@gromacs.org   http://folding.bmc.uu.se
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to