Sir.<br> <br>I am sorry for changing the topic. It was not intentional<br>It just came as I wanted to mention the other part of the<br>problem that data had.<br>I have checked the entire file of the gmxdump output <br>of ener.edr firl for Box dimention as well as Temp.<br>
<br>Thanks<br>With Regards<br>M. Kavyashree<br><div class="gmail_quote"><br><br>On Fri, Jun 24, 2011 at 6:32 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br>
<br>
Kavyashree M wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Sir,<br>
<br><div class="im">
I went through the whole .edr file for one<br>
specific term in which I was getting nan.<br>
Of the three columns (Energy; Av. Energy; Sum Energy)<br>
order (power) for temperature did not change at all in the<br>
first column but in the second column there was a change of<br>
power between +1 and +2 which was +2 in another .edr file from<br>
other simulation.<br>
ven though the second column value is arround 1--e+02 while<br>
those with +1 power is around 9.--e+01 which is quite near to<br>
1.--e+02.<br>
<br>
But I did not find any huge changes in the temperature.<br>
<br>
</div></blockquote>
<br>
Per my previous message, please focus on the issue at hand. I do not recall the specifics of the issue with temperature and I'm not going to go digging through the archive for messages scattered over many days. You started this thread regarding box vectors, for which we have the data below in the thread. Please keep the discussion focused. If you jump around between ideas with somewhat abstract descriptions, no one is going to be interested in troubleshooting. If you want free help, make it easy to help you.<br>
<br>
If you have something relevant to post regarding box vectors, please do so, otherwise, if you want to investigate a different issue please give it proper context (new thread or including information from the old thread).<br>
<br>
Actual output would be useful, although I realize that may be too large of a message to send across the list, so snippets that illustrate anything useful would be appropriate.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Thank you<br>
With Regards<br>
M. Kavyashree<br>
<br></div><div class="im">
On Fri, Jun 24, 2011 at 5:53 PM, Kavyashree M <<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a> <mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a>>> wrote:<br>
<br>
Sir,<br>
<br>
I tried it already but It was huge file. this kind of nan comes only<br>
for RMSD not for any other term. I will try looking into that file<br>
again.<br>
<br>
Thanks<br>
With regards<br>
M. Kavyashree<br>
<br>
<br>
<br>
<br>
On Fri, Jun 24, 2011 at 5:20 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a><br></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
Kavyashree M wrote:<br>
<br>
Sir,<br>
<br>
I am sorry I did not get any mail for that query since I<br>
posted it.<br>
So I had to ask once again..<br>
<br>
1. g_energy -f ener.edr -o box.xvg<br>
selecting 15 16 17 0<br>
<br>
2. gmxcheck -e ener.edr<br>
<br>
Opened ener.edr as single precision energy file<br>
frame: 0 (index 0), t: 0.000 Last<br>
energy frame read 50000 time 100000.000 Found 50001 frames with a timestep of 2 ps.<br>
<br>
3. Gromacs version 4.5.3<br>
<br>
This is the same data where I was mentioning that I<br>
was getting "not a number" (nan) error for all<br>
energy calculations under RMSD section.<br>
<br>
<br>
I've said several times before that this energy file is probably<br>
junk, and again I'd suspect that. The only thing to check would<br>
be to run gmxdump on the .edr file (and redirect into some<br>
output file or pipe the output through "more" to actually see<br>
it) and see where the energy values go screwy. If you find a<br>
bunch of nonsensical information, then you have your answer.<br>
<br>
-Justin<br>
<br>
Thank you fro answering.<br>
<br>
With Regards<br>
M. Kavyashree<br>
<br>
<br>
On Fri, Jun 24, 2011 at 4:57 PM, Justin A. Lemkul<br>
<<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br></div></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>>> wrote:<br>
<br>
<br>
<br>
Kavyashree M wrote:<br>
<br>
Dear users,<br>
<br>
Any suggestions?<br>
<br>
<br>
You haven't provided nearly enough diagnostic information<br>
for anyone<br>
to offer you any useful help (as Mark said yesterday).<br>
For example,<br>
please provide:<br>
<br>
1. Your exact g_energy command line<br>
2. The output of gmxcheck for this .edr file<br>
3. Your Gromacs version<br>
<br>
-Justin<br>
<br>
Thank you<br>
M. Kavyashree<br>
<br>
<br>
On Thu, Jun 23, 2011 at 10:38 AM, Kavyashree M<br>
<<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a> <mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a>><br>
<mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a> <mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a>>><br>
<mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a> <mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a>><br>
<mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a> <mailto:<a href="mailto:hmkvsri@gmail.com" target="_blank">hmkvsri@gmail.com</a>>>>> wrote:<br>
<br>
Dear users,<br>
<br>
In one of the simulations while calculating box<br>
dimensions<br>
using g_energy this output was obtained -<br>
Statistics over 50000001 steps [ 0.0000 through<br>
100000.0000<br>
ps ], 3<br>
data sets<br>
All statistics are over 1978700 points<br>
<br>
Energy Average Err.Est. RMSD<br>
Tot-Drift<br>
------------------------------<u></u>____--------------------------<u></u>--__--__-------------------<br>
Box-X 0.0453229 1.8 0.637535 -9.93022 (nm)<br>
Box-Y 0.0453229 1.8 0.637535 -9.93022 (nm)<br>
Box-Z 0.0320482 1.3 0.450805 -7.02173 (nm)<br>
<br>
but the dimensions of the box is different. plot<br>
attached.<br>
I am not able to figure out why only 1978700 data<br>
point are<br>
considered.<br>
Kindly give some suggestions.<br>
<br>
Thank you<br>
with regards<br>
M.Kavyashree<br>
<br>
<br>
<br>
-- ==============================<u></u>____==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br></div></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> |<div class="im">
<br>
(540) 231-9080<br>
<br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<u></u>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank">vt.edu/Pages/Personal/justin</a><br>
</div>
<<a href="http://vt.edu/Pages/Personal/justin" target="_blank">http://vt.edu/Pages/Personal/<u></u>justin</a>><div class="im"><br>
<<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>_<u></u>_<a href="http://vt.edu/Pages/Personal/justin" target="_blank">vt.edu/Pages/Personal/justin</a><br>
<<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a>>><br>
<br>
==============================<u></u>____==========<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><u></u>><div><div>
</div><div class="h5"><br>
<br>
<a href="http://lists.gromacs.org/____mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/____<u></u>mailman/listinfo/gmx-users</a><br>
<<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a>><br>
<<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a><br>
<<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>>><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/____Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/____<u></u>Support/Mailing_Lists/Search</a><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a>><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>>><br>
before posting!<br>
Please don't post (un)subscribe requests to the list. Use<br>
the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@" target="_blank">gmx-users-request@</a>__<a href="http://gromacs.org" target="_blank">gr<u></u>omacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>>.<br>
<br>
Can't post? Read<br>
<a href="http://www.gromacs.org/____Support/Mailing_Lists" target="_blank">http://www.gromacs.org/____<u></u>Support/Mailing_Lists</a><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a>><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>>><br>
<br>
<br>
<br>
-- ==============================<u></u>__==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</a><br>
<<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a>><br>
<br>
==============================<u></u>__==========<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a><br>
<<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.<br>
Can't post? Read <a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>><br>
<br>
<br>
<br>
</div></div></blockquote><div><div></div><div class="h5">
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don'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't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>