RE: energy output

From: Deng, Jinxia (Nancy) (jinxianancy.deng_at_zoetis.com)
Date: Wed Dec 09 2015 - 13:12:11 CST

Thank you All for your response. I will start from there ☺

Nancy


From: Jeffrey Potoff [mailto:jpotoffx_at_gmail.com]
Sent: Wednesday, December 09, 2015 12:57 PM
To: namd-l_at_ks.uiuc.edu; Deng, Jinxia (Nancy)
Subject: Re: namd-l: energy output

On 12/9/2015 12:39 PM, Deng, Jinxia (Nancy) wrote:
Dear users,

I am fairly new to use NAMD, and wonder how to get energy profile in a separate output file during time course. Currently, it seems everything just print out on the screen or in a directed output file which has many other information as well…

Thanks and looking forward to hearing from you.

Nancy

Nancy Deng, Ph.D. | Principal Scientist, Computational Chemistry & Molecular Modeling
Zoetis Inc| Global Therapeutics Research | 333 Portage Street, Kalamazoo, MI 49007
Office: 2693599121 | jinxianancy.deng_at_zoetis.com<mailto:jinxianancy.deng_at_zoetis.com>


Hi Nancy,
    You can extract the energies using NAMDplot in VMD, or you can get them out of the log file with an awk script, such as this. The $<number> corresponds to a particular column in the energy table. So you can get whatever energies you want by changing, or adding, the variables in the list.

#!/usr/bin/awk -f
BEGIN{eave_total=0; count=0;nmolecule=500; nequil=200000}
{if($1 = /ENERGY:/ && $2>nequil)
{ print $2, $14;
eave_total += $14;
++count;
}
}
END{print "Potential energy average =",eave_total/nmolecule/count,"kcal/mol\n"}


--

Jeffrey Potoff

Professor & Director of Early Engineering Programs

Department of Chemical Engineering and Materials Science

Wayne State University

http://potoff1.eng.wayne.edu

http://gomc.eng.wayne.edu


This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:22:18 CST