[AMBER] Looping over a trajectory

From: Lee-Ping Wang <leeping.stanford.edu>
Date: Tue, 29 May 2012 00:27:22 -0700

Hi there,

 

I'm relatively new to AMBER, and I'm interested in running a calculation
that loops over the conformations in a trajectory and prints out the
energies and forces for each conformation. I was able to calculate forces
on a single conformation by running the command: "sander -i force.mdin -o
force.mdout -p res.prmtop -c res.inpcrd" with the following input file:

 

#==== File Begins ====#

Single Point energy and Force calculation

&cntrl

imin = 0, ntb = 0, cut = 999, nstlim = 0,

/

&debugf

do_debugf = 1, dumpfrc = 1,

/

#==== File Ends ====#

 

However, when I set imin = 5 and call sander with the extra argument "-y
equilibrate.mdcrd", it still appears to calculate forces for just one
conformation.

 

Is it possible to calculate forces for each conformation in a trajectory?
If not, then I'll create a series of inpcrd files and call sander multiple
times, although it seems most efficient to do it in a single run.

 

Thanks,

 

- Lee-Ping Wang

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue May 29 2012 - 00:30:02 PDT
Custom Search