How should I change it so it calculates the 40 ps?<br><br><div class="gmail_quote">On Wed, Oct 22, 2008 at 6:34 PM, Xavier Periole <span dir="ltr">&lt;<a href="mailto:X.Periole@rug.nl">X.Periole@rug.nl</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="Ih2E3d">On Wed, 22 Oct 2008 17:38:58 -0400<br>
&nbsp;&quot;Daniel K&quot; &lt;<a href="mailto:chmdan@gmail.com" target="_blank">chmdan@gmail.com</a>&gt; wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I have been trying to use g_rotacf with a simulation run for 40 ps (20000)<br>
with steps of 2fs and the frames where extracted every 50 fs (25 frames).<br>
The trajectory file is correct because it has 801 frames. When I calculate<br>
the g_rotacf I only get 20 ps and the time step is 50 fs. Does any body no<br>
why?<br>
</blockquote></div>
What is wrong with this? By default g_rotacf plots only half of the time<br>
correlation function ...<br>
I do not know what you looking at but 40 ps is pretty short simulation.<div class="Ih2E3d"><br>
&nbsp;<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
This is the header and end of the file that I get<br>
# This file was created Mon Oct 20 11:38:06 2008<br>
# by the following command:<br>
# g_rotacf -P 2 -f md_noT.trr -s md_noT.tpr -n NC.ndx -o md_noT_ACF.xvg -d<br>
#<br>
# g_rotacf is part of G R O M A C S:<br>
#<br>
# Glycine aRginine prOline Methionine Alanine Cystine Serine<br>
#<br>
@ &nbsp; &nbsp;title &quot;Rotational Correlation Function&quot;<br>
@ &nbsp; &nbsp;xaxis &nbsp;label &quot;Time (ps)&quot;<br>
@ &nbsp; &nbsp;yaxis &nbsp;label &quot;C(t)&quot;<br>
@TYPE xy<br>
 &nbsp; &nbsp;0.000 &nbsp; &nbsp; 1.00000<br>
 &nbsp; &nbsp;0.050 &nbsp; &nbsp; 0.95376<br>
 &nbsp; &nbsp;0.100 &nbsp; &nbsp; 0.90907<br>
&nbsp;....<br>
 &nbsp; 19.800 &nbsp; &nbsp;-0.22896<br>
 &nbsp; 19.850 &nbsp; &nbsp;-0.22921<br>
 &nbsp; 19.900 &nbsp; &nbsp;-0.23028<br>
 &nbsp; 19.950 &nbsp; &nbsp;-0.23078<br>
 &nbsp; 20.000 &nbsp; &nbsp;-0.23111<br>
&amp;<br>
</blockquote>
<br></div>
-----------------------------------------------------<br>
XAvier Periole - PhD<br>
<br>
- Molecular Dynamics Group -<br>
NMR and Computation<br>
University of Groningen<br>
The Netherlands<br>
-----------------------------------------------------<br>
_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</blockquote></div><br>