On Tue, Jul 07, 2009, Ross Walker wrote:
>
> So it looks like nscm is getting set to 1000 both in sander and pmemd
> regardless of the value of ntr.
> This seems to run the nucleosome run nice and stable with the behavior I
> expect. Anyone got any comments?
Why is repositioning necessary at all if you have restraints? Plus, your idea
doesn't seem to take into account restarts -- if you've been changing the
reference coordinates on the fly, you'd have to write them out somewhere in
order to be able to restart. But "reference coords are just that: a reference
that should not change.
My vote is to just enforce nscm=0 when ntr=1. However, I can't reproduce your
report using sander: if I set ntr=1 in a gb calculation, nscm is printed out
as 9999999, not 1000. So, you might want to do some quick debugging here.
...dac
_______________________________________________
AMBER-Developers mailing list
AMBER-Developers.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber-developers
Received on Tue Jul 07 2009 - 18:09:40 PDT
