On 3/14/14, 4:28 PM, Nimmy McNimmerson wrote:
I have a .top file looking like this: #include "spc.itp" #include "ffgmx.itp" #include "lipid.itp" #define POSRES_LIPID #include "mylipid.itp" #include "dppc.itp" [ system ] test system [ molecules ] ; Compound #mols DX5 200 DPPC 200 The #define POSRES_LIPID is used to determine whether to include position restraints in mylipid.itp. However, when I try to run grompp on this (with everything the exact same as before I included the position restraints), I get the error "Topology include file "lipid_posre.itp" not found". If I take out all the position restraints, it works fine and grompp does what is expected. What's going on? I never asked to include this file, I don't even know what lipid_posre.itp is and google has no results on it. I tried adding #include "lipid_posre.itp" at the top of the .top file, just to see what happened, but I got the same error. How am I supposed to add position restraints? The GROMACS manual just says to add them and doesn't mention additional includes, and I added them in the same way that the manual said to.
You did ask to #include the position restraints because you're explicitly using #define to set the inclusion conditions to be true.
If you don't know what's calling "lipid_posre.itp," that's what grep is for. It's probably dppc.itp, which includes the option to restrain the lipids under a conditional #ifdef POSRES_LIPID. Your topology says "turn on lipid restraints," and one of the lipid .itp files (again, probably dppc.itp) says "if lipid restraints are turned on, #include "lipid_posre.itp."
Remove the #define statement from the topology; that's not normal use (technically correct in terms of syntax, but then it's always true, which is probably not what you want) and instead use "define -DPOSRES_LIPID" in the .mdp file to conditionally trigger the restraints.
-Justin -- ================================================== Justin A. Lemkul, Ph.D. Ruth L. Kirschstein NRSA Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 http://mackerell.umaryland.edu/~jalemkul ================================================== -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
