Re: [AMBER] Problem During Gradual Heating

From: Jason Swails <jason.swails.gmail.com>
Date: Fri, 28 Jan 2011 09:42:28 -0500

Carlos raised a good point. All data is written to standard output which
for the master node is the output file. Slave threads will output to stdout
instead of mdout, so the error might be found in the stdout file (specified
by #PBS -o output if you are using PBS scheduling). Can you include the
output of this file?

On Fri, Jan 28, 2011 at 9:23 AM, Carlos Simmerling <
carlos.simmerling.gmail.com> wrote:

> try running it single processor as a test- sometimes if there is an
> error message on one of the non-master nodes it won't go in the output
> file.
>
>
>
> On Fri, Jan 28, 2011 at 9:18 AM, Hirdesh Kumar <hirdesh.iitd.gmail.com>
> wrote:
> > Yes, I checked for job status. and I found that it is still running. I
> think
> > that the amber installation in OK as all other jobs are running OK. I am
> > runnning the job using MPI.
> >
> >
> >
> > On Fri, Jan 28, 2011 at 7:45 PM, Carlos Simmerling <
> > carlos.simmerling.gmail.com> wrote:
> >
> >> it isn't clear that this is an amber problem.
> >> did you check to see if the job is still running?
> >> does your amber installation pass the test cases?
> >>
> >>
> >> On Fri, Jan 28, 2011 at 8:58 AM, Hirdesh Kumar <hirdesh.iitd.gmail.com>
> >> wrote:
> >> > Hi All,
> >> > I am tying to heat my system (protein + bound cofactor) slowly. I am
> >> using
> >> > ntt =3 with gamma_ln = 1.0, Below mentioned is the script which I am
> >> using.
> >> >
> >> >
> >> >
> >> > Heating the system linearly with weak restraints
> >> > &cntrl
> >> > imin = 0,
> >> > irest = 0,
> >> > ntx = 1,
> >> > ntb = 1,
> >> > cut = 12,
> >> > ntr = 1,
> >> > ntc = 2,
> >> > ntf = 2,
> >> > tempi = 0.0,
> >> > temp0 = 300.0,
> >> > ntt = 3,
> >> > ig = 20933
> >> > gamma_ln = 1.0, nmropt=1,
> >> > nstlim = 50000, dt = 0.002
> >> > ntpr = 100, ntwx = 100, ntwr = 1000
> >> > /
> >> > &wt TYPE='TEMP0', istep1=0, istep2=40000,
> >> > value1=0.1, value2=300.0, /
> >> > &wt TYPE='DUMPFREQ',ISTEP1=10, /
> >> > &wt TYPE='END' /
> >> > Keep protein fixed
> >> > 5.0
> >> > RES 1 339
> >> > END
> >> > END
> >> >
> >> > The equilibration runs fine upto 9100 steps. After that there is* no
> >> change
> >> > in the mdout file* but the sander job is running. Previously I didnt
> >> > included Ig value in my script so I thought that it may be due to Ig
> >> value
> >> > as I am using ntt=3. Next time I added the Ig value (a random number)
> but
> >> > this time also the run is upto 7100 steps after that there is no
> change
> >> in
> >> > the mdout file.
> >> >
> >> > Below mentioned is the last line of my output file:
> >> >
> >> > NSTEP = 7000 TIME(PS) = 14.000 TEMP(K) = 49.27 PRESS
> =
> >> > 0.0
> >> > Etot = -211822.1441 EKtot = 5837.1249 EPtot =
> >> > -217659.2690
> >> > BOND = 323.1122 ANGLE = 956.0308 DIHED =
> >> > 3200.9235
> >> > 1-4 NB = 1056.0283 1-4 EEL = 14174.4148 VDWAALS =
> >> > 35009.8809
> >> > EELEC = -272523.6263 EHBOND = 0.0000 RESTRAINT =
> >> > 143.9667
> >> > EAMBER (non-restraint) = -217803.2357
> >> > Ewald error estimate: 0.3030E-04
> >> >
> >>
> ------------------------------------------------------------------------------
> >> >
> >> > NMR restraints: Bond = 0.000 Angle = 0.000 Torsion =
> >> 0.000
> >> >
> >>
> ===============================================================================
> >> >
> >> > NSTEP = 7100 TIME(PS) = 14.200 TEMP(K) = 50.36 PRESS
> =
> >> > 0.0
> >> > Etot = -211664.7202 EKtot = 5966.1111 EPtot =
> >> > -217630.8314
> >> > BOND = 312.2514 ANGLE = 957.5403 DIHED =
> >> > 3192.0819
> >> > 1-4 NB = 1062.4902 1-4 EEL = 14166.2999 VDWAALS =
> >> > 34826.4805
> >> > EELEC = -272299.1537 EHBOND = 0.0000 RESTRAINT =
> >> > 151.1782
> >> > EAMBER (non-restraint) = -217782.0095
> >> > Ewald error estimate: 0.2147E-04
> >> >
> >>
> ------------------------------------------------------------------------------
> >> >
> >> > NMR restraints: Bond = 0.000 Angle = 0.000 Torsion =
> >> 0.000
> >> >
> >>
> ===============================================================================
> >> >
> >> >
> >> >
> >> > Please suggest me where I am wrong..
> >> >
> >> > Thanks
> >> > Hirdesh
> >> > _______________________________________________
> >> > AMBER mailing list
> >> > AMBER.ambermd.org
> >> > http://lists.ambermd.org/mailman/listinfo/amber
> >> >
> >>
> >> _______________________________________________
> >> AMBER mailing list
> >> AMBER.ambermd.org
> >> http://lists.ambermd.org/mailman/listinfo/amber
> >>
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>



-- 
Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Graduate Student
352-392-4032
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jan 28 2011 - 07:00:06 PST
Custom Search