On Fri, Sep 26, 2014, Parker de Waal wrote:
>
> After reading through some of the old mailing list threads regarding
> the various workarounds for ntr =1 and nmropt=1, I still need some
> clarification.
>
> I have an NVT prod.in as show below:
>
>
> Now I wish to add cartesian restraints, following the snippet here:
> http://ambermd.org/Questions/constraints.html into my rst.dist file as
> seen below:
The group cards for cartesian restraints go into the mdin file, not into the
restraints file.
>
> # N373 D69 VS
> &rst
> ixpk= 0, nxpk= 0, iat= 5235, 6651, r1= 0.00, r2= 0.00, r3= 6, r4= 7,
> rk2=0.0, rk3=25.0, ir6=1, ialtd=0,
> &end
> ... (lines removed)
> # A150 F244 VS
> &rst
> ixpk= 0, nxpk= 0, iat= 2310, 9472, r1= 0.00, r2= 0.00, r3= 6, r4= 7,
> &end
> Hold GPCR fixed
> 20.0
> RES 1 708
> END
> END
>
> However whenever I run this the follow error occurs:
>
> 5. REFERENCE ATOM COORDINATES
>
> default_name
> ----- READING GROUP 1; TITLE:
> DISANG=rst.dist
Note that the program thinks the DISANG card is the title for the cartesian
restraints; so somehow it is not reading the NMR-like restraints correctly.
The first thing you need to do is set ntr=0, and make sure that you are
getting the NMR restraints correctly processed. Follow the example in the
$AMBERHOME/test/nmr directory.
By far the simplest thing to do for cartsian restraints is to set them
in the cntrl namelist:
ntr=1, restraint_wt=20., restraintmask=':1-708',
Then there is no need for any GROUP input.
Beyond this:
You could try using the correct "/" rather than the obsolete (and never
standard "&end" to end namelists. But that is probably not the source of the
problem. You may need to let us know which program you are running
(sander or pmemd?) and which version. Posting the entire mdout file is
probably also required...there is something subtle that is preventing your
file from being parsed correctly.
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Sep 25 2014 - 20:00:03 PDT
