AMBER Archive (2005)

Subject: AMBER: problem in writing

From: Gustavo Pierdominici Sottile (gsottile_at_unq.edu.ar)
Date: Mon Feb 07 2005 - 13:59:49 CST


Hi,
    The problem I have is that when using SANDER the .out file is allways
empty. Perhaps I am doing something wrong but can not realize what. I will
explain from the beggining: I have the pdb, two .in files for non-standard
aminoacids and the frcmod file. Then I use tleap with the following
commands:
addatomtypes {{"CS" "C" "sp3"}}
   addatomtypes {{"HG" "H" "sp3"}}
   loadamberprep trp.in
   loadamberprep trp1.in
   loadamberparams fcrmod.mad
   mad= loadpdb a.pdb
   bond mad.402.CE3 mad.453.CD1
   bond mad.454.SG mad.423.SG
   bond mad.422.SG mad.391.SG
   bond mad.466.SG mad.381.SG
   bond mad.383.SG mad.431.SG
   bond mad.433.SG mad.368.SG
   bond mad.374.SG mad.406.SG
   saveamberparm mad prmtop.top prmcrd.crd
 
Then I do sander -O -i din.in -p prmtop.top -c prmcrd.crd -o din.out, and
the din.in file looks like this:
 &cntrl
   imin=0, irest=0, ntx=1, tempi=0.0,
   ntt=3, temp0=4.0, tautp=0.4,
   ntp=0, taup=2.0, nrespa=1,
   nstlim=40000, dt=0.001,
   ntwe=100, ntwx=100, ntpr=100,
 &end

There is something wrong in my prmtop.top and prmcrd.crd files because when
I use those file taken from the amber/test/val.... directory, there is no
problem.
I do not know what I am missing
Pd= If I try to make a minimization, the same things happend. No file is
written
Thanks in advanve
Gustavo

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu