FW: [gmx-users] position restraints

[abdullah ahmed]

Hi! 

In your previous mail you mentioned: 

The position restraints must belong to the [moleculetype] of 
the species to be restrained.  Once you #include a new molecule, you start a new 
[moleculetype] entry and the position restraints belong to it.  
 

So I rechecked my .top file and found that [moleculetype] only occurs once. Perhaps I have misunderstood you. So I added the top file below. I did not add the contents of [atoms] [bonds] etc because I felt the mail would become unnecessarily long. 

; Include forcefield parameters
#include "ffoplsaa.itp"

[ moleculetype ]
; Name            nrexcl
Protein_A           3

[ atoms ]

[ bonds ]

[ pairs ]

[ angles ]

[ dihedrals ]

[ dihedrals ]

; Include Position restraint file
#ifdef POSRES
#include "posre.itp"
#endif

; Include water topology
#include "spc.itp"

#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct       fcx        fcy        fcz
   1    1       1000       1000       1000
#endif

; Include generic topology for ions
#include "ions.itp"

[ system ]
; Name
Protein

[ molecules ]
; Compound        #mols
Protein_A           1

Thanks again for your help,
Abdullah 

 		 	   		  
Hotmail: Trusted email with Microsoft’s powerful SPAM protection. Sign up now. 		 	   		  
_________________________________________________________________
Your E-mail and More On-the-Go. Get Windows Live Hotmail Free.
https://signup.live.com/signup.aspx?id=60969
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100603/828d30b0/attachment.html>