Re: RE: Re: AMBER: A question about Dynamics.

From: David A. Case <case.scripps.edu>
Date: Fri, 22 Sep 2006 17:57:46 -0700

On Fri, Sep 22, 2006, Qingning Shu wrote:
>
> At the beginning of second simulation, even
> the temperature dropped.
>

First run:
> sander -O -i vac_md0.in -p vac.prmtop -c vac_min.rst -x vac_md0.crd -r
> vac_md0.rst -o vac_md0.out -ref vac_min.rst
                                  ^^^^^^^^^^^^

Second run:

>
> sander -O -i vac_md1.in -p vac.prmtop -c vac_md0.rst -x vac_md1.crd -r
> vac_md1.rst -o vac_md1.out -ref vac_md0.rst
                             ^^^^^^^^^^^^^^^^

You are using different restraint files for the two runs, so of course the
energies will differ. Take Ross Walker's advice, and try *short* runs with
ntpr=1. That is always a good idea when you results you don't understand.
Only printing every 20 steps makes debugging or understanding what is going on
very difficult.

...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Sep 24 2006 - 06:07:22 PDT
Custom Search