Dear Ed,
> have to be run out to > 10 ns. To write out an entry to a md.crd
> trajectory file at each time step for analysis will generate
> huge files,
> even if the simulation is broken up into segments. I would
> like to avoid
> this if possible. Is there a way in Amber7 to save the
> coordinates of a
> restricted set of four atoms at every time step, while adding a new
> structure to the md.crd trajectory file at a larger increment in time?
As far as I am aware there is no feature in the code to do this but it
shouldn't be too hard to implement it yourself. See the routine corpac() in
dynlib.f - called from runmd(). This routine handles the writing of
coordinates to an mdcrd file. You could quite easily modify this to check
the step number and either write the 4 atoms you want to a file (say mdcrd2)
or write a standard coordinate frame.
All the best
Ross
/\
\/
|\oss Walker
| Department of Molecular Biology TPC15 |
| The Scripps Research Institute |
| Tel:- +1 858 784 8889 | EMail:- ross.rosswalker.co.uk |
|
http://www.rosswalker.co.uk/ | PGP Key available on request |
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Nov 28 2004 - 01:53:00 PST
