On 26/07/2012 6:47 PM, reising...@rostlab.informatik.tu-muenchen.de wrote:
Ho,
first I minimize my structure. This is the corresponding mdp file:
define = -DPOSRES
integrator = steep
emtol = 10
nsteps = 1500
nstenergy = 1
energygrps = System
coulombtype = PME
rcoulomb = 0.9
rvdw = 0.9
rlist = 0.9
fourierspacing = 0.12
pme_order = 4
ewald_rtol = 1e-5
pbc = xyz
and then I run a md run. This is the corresponding mdp file:
define = -DPOSRES
integrator = md
dt = 0.001
nsteps = 5000
nstxout = 100
nstvout = 0
nstfout = 0
nstlog = 1000
nstxtcout = 500
nstenergy = 5
energygrps = Protein Non-Protein
nstcalcenergy = 5
nstlist = 10
ns-type = Grid
pbc = xyz
rlist = 0.9
coulombtype = PME
rcoulomb = 0.9
rvdw = 0.9
fourierspacing = 0.12
pme_order = 4
ewald_rtol = 1e-5
gen_vel = yes
gen_temp = 200.0
gen_seed = 9999
constraints = all-bonds
tcoupl = V-rescale
tc-grps = Protein Non-Protein
tau_t = 0.1 0.1
ref_t = 298 298
pcoupl = no
In my topology file I include the restraint files like this:
; Include Position restraint file
#ifdef POSRES
#include "posre.itp"
#endif
#ifdef POSRES
#include "posre_memb.itp"
#endif
This won't work for multiple [moleculetype] entries. See
http://www.gromacs.org/Documentation/How-tos/Position_Restraints
I just recognize that there is a "DPOSRES" in the mdp files and a "POSRES"
in my topology file. Is this the problem. Do I have to write it the same
way in the several files?
http://www.gromacs.org/Documentation/Include_File_Mechanism
Mark
Thank you
Hi,
Did you call your genrestr output in your mdp file using define = * ? (for
example define = -DPOSRES ; where posres was your file name.)
Thanks
Peterson J
--
View this message in context:
http://gromacs.5086.n6.nabble.com/question-about-minimisation-tp4999681p4999685.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists