namdenergy and bigdcd

From: Grant Schauer (grantschauer_at_gmail.com)
Date: Thu Jun 02 2016 - 23:58:23 CDT

Hi all,

I tried to run namdenergy on a very large trajectory with bigdcd (see
below), but it seems to be just looping over the same first two frames on
repeat, which maybe makes some sense now because namdenergy assumes we
start at the first frame and I can't pass $frame to it. Is there a
workaround for this?

Or, maybe a better question I should have started with: will namdenergy run
out of memory with very large dcds? Thank you!

source bigdcd.tcl
package require namdenergy

set mol [mol new ./molecule.psf type psf waitfor all]

mol addfile ./molecule.pdb type pdb waitfor all
set sel1 [atomselect top protein]
set sel2 [atomselect top nucleic]

proc energies {frame} {
global mol sel1 sel2 out

namdenergy -nonb -vdw -elec -sel $sel1 $sel2 -ofile out.dat -par
par_all27_prot_na.prm -par
}

bigdcd energies dcd ./giant_dcd.dcd
bigdcd_wait

quit

This archive was generated by hypermail 2.1.6 : Sun Dec 31 2017 - 23:20:32 CST