<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span>Dear Justin</span></div><div><span>Thanks again. Ok I got it.</span></div><div><span>Only there remains a question:</span></div><div><span>you mean it is possible that the simulation has finished for 1 ns but the .xtc file does not have the information for whole frames or even it is possible for g_gyrate program to crash and give the non-compelete .xvg file?</span></div><div><br><span></span></div><div><span>Regards</span></div><div><span>D.M</span></div><div><br></div>  <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1">  <b><span style="font-weight:bold;">From:</span></b> Justin A. Lemkul &lt;jalemkul@vt.edu&gt;<br> <b><span style="font-weight:
 bold;">To:</span></b> delara aghaie &lt;d_aghaie@yahoo.com&gt;; Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt; <br> <b><span style="font-weight: bold;">Sent:</span></b> Monday, 11 June 2012, 17:46<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [gmx-users] radius of gyration<br> </font> </div> <br><br><br>On 6/11/12 9:09 AM, delara aghaie wrote:<br>&gt; Dear Justin<br>&gt; My question can be related to both. Ok let me talk about the tutorial. The<br>&gt; simulation has finished compeletly and log file and other outputs show the<br>&gt; completion of 1 ns simulation.<br>&gt; <br>&gt; this is part of gyrate.xvg file I see:<br>&gt; -------------<br>&gt; @ s0 legend "Rg"<br>&gt; @ s1 legend "RgX"<br>&gt; @ s2 legend "RgY"<br>&gt; @ s3 legend "RgZ"<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  0&nbsp; &nbsp;  1.40391&nbsp; &nbsp;  1.20269&nbsp; &nbsp;  1.02309&nbsp; &nbsp;  1.20363<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 
 2&nbsp; &nbsp;  1.41596&nbsp; &nbsp;  1.21322&nbsp; &nbsp;  1.03459&nbsp; &nbsp;  1.21146<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  4 1.41283&nbsp; &nbsp;  1.21881&nbsp; &nbsp;  1.02271&nbsp; &nbsp;  1.20861<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  6&nbsp; &nbsp;  1.41882&nbsp; &nbsp;  1.21595&nbsp; &nbsp;  1.03424&nbsp; &nbsp; &nbsp; 1.2157<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  8&nbsp; &nbsp; &nbsp; 1.4139&nbsp; &nbsp;  1.20694&nbsp; &nbsp;  1.03791&nbsp; &nbsp;  1.21006<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 10&nbsp; &nbsp;  1.42295&nbsp; &nbsp;  1.21367&nbsp; &nbsp;  1.03899&nbsp; &nbsp;  1.22356<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 12&nbsp; &nbsp;  1.42184&nbsp; &nbsp;  1.20803&nbsp; &nbsp; &nbsp; 1.0391&nbsp; &nbsp;  1.22646<br>&gt; ----<br>&gt; ---<br>&gt;&nbsp; &nbsp; &nbsp;  518&nbsp; &nbsp;  1.41756 1.20547&nbsp; &nbsp;  1.04411&nbsp; &nbsp;  1.21475<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp;  520&nbsp; &nbsp;  1.42478&nbsp; &nbsp; 
 1.21366&nbsp; &nbsp;  1.04952&nbsp; &nbsp;  1.21882<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp;  522&nbsp; &nbsp;  1.41129&nbsp; &nbsp;  1.20451&nbsp; &nbsp;  1.04001&nbsp; &nbsp;  1.20458<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp;  524&nbsp; &nbsp;  1.42205&nbsp; &nbsp;  1.20918&nbsp; &nbsp;  1.04929&nbsp; &nbsp;  1.21711<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp;  526&nbsp; &nbsp;  1.41778&nbsp; &nbsp;  1.20381&nbsp; &nbsp;  1.04779&nbsp; &nbsp;  1.21375<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp;  528&nbsp; &nbsp;  1.41853&nbsp; &nbsp;  1.20465&nbsp; &nbsp;  1.06006 1.20398<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp;  530&nbsp; &nbsp;  1.42129&nbsp; &nbsp;  1.20401&nbsp; &nbsp;  1.05779&nbsp; &nbsp;  1.21307<br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 288,8&nbsp; &nbsp; &nbsp; &nbsp;  Bot<br>&gt;
 ---------------------------------------<br>&gt; This was the start and end of data reported in gyrate.xvg file.<br>&gt; all the outputs are written in files every 1000 steps (1000*0.002=2ps).<br>&gt; I expect to see in gyrate.xvg file, the last line for 500 which corresponds to<br>&gt; the (500*2=1000ps)<br>&gt; <br>&gt; can you please explain the result for having lines till 530?<br><br>Your output file is not labeled according to frame number, it is labeled according to time and thus 530 corresponds to the time of 530 ps.&nbsp; If the file does indeed end there then either there is a corruption in the .xtc file, your simulation did not go to 1 ns, or g_gyrate exited prematurely.&nbsp; It's up to you to figure out what may have gone wrong.&nbsp; No one else can possibly know that.<br><br>&gt; -----------<br>&gt; About plotting the graph if I know use the two first coloumns, I will not have<br>&gt; graph for Rg versus simulation time which I am
 interested in?<br>&gt; <br><br>That's precisely what you will have.&nbsp; Look at the data labels in the .xvg file and you will see what each data set is.<br><br>-Justin<br><br>-- ========================================<br><br>Justin A. Lemkul, Ph.D.<br>Research Scientist<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]vt.edu | (540) 231-9080<br><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br><br>========================================<br><br><br><br><br> </div> </div>  </div></body></html>