AMBER: Problem with PCSHIFT module

From: Edoardo Saccenti <saccenti.cerm.unifi.it>
Date: Fri, 30 Mar 2007 13:20:37 +0200

Dear All,
I'm trying to include pseudocontact shift as a constraint
for energy minimization. I followed the manual instructions
but I get the following error:

>> Pseudocontact shifts will be read from file: ../../NOE/amberpcs1.in
 namelist reports error reading &pcshf

I think the problem is in the file format but I'm unable to find
differences between my file and the example given in the manual.
any clue??

Many thanks
edoardo


This is my input file:
REM
&cntrl
 imin=1,
 maxcyc=18000, ncyc=200,
 scee=1.2,
 ntpr=100,
 ntwx=100,
 cut=10.0,
 nmropt=2,
 pencut=0.1,
 ntb=0,
/

&wt type='REST', istep1=0, istep2=1, value1=0.1,
           value2=1.0, /

&wt type='END' /
DISANG=../../NOE/NOE
PCSHIFT=../../NOE/amberpcs1.in
LISTOUT=POUT

>>> and this my pcs file
&pcshf
nprot=133,
nme='CE ',
nmpcm=1,
optphi(1)=1.047,
opttet(1)=-0.948,
optomg(1)=-1.071,
opta1(1)=559.189,
opta2(1)=208.369,
optkon=30,
iprot(1)=1, obs(1)=0.000, wt(1)=1.000, tolpro(1)=0.150, mltpro(1)=1,
iprot(2)=2, obs(2)=-0.020, wt(2)=1.000, tolpro(2)=0.150, mltpro(2)=1,
iprot(3)=23, obs(3)=-0.030, wt(3)=1.000, tolpro(3)=0.150, mltpro(3)=1,
iprot(4)=24, obs(4)=-0.020, wt(4)=1.000, tolpro(4)=0.150, mltpro(4)=1,
.
.
.
iprot(133)=1198, obs(133)=-0.030, wt(133)=1.000, tolpro(133)=0.150,
mltpro(133)=1,
&end



-- 
************************
Dott. Edoardo Saccenti
CERM - Magnetic Resonance Center
University of Florence
Via Luigi Sacconi, 6
50019 Sesto Fiorentino (Florence) - Italy
tel. +39055 4574281 -
e-mail: saccenti.cerm.unifi.it
web   : http://www.cerm.unifi.it
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Apr 01 2007 - 06:07:38 PDT
Custom Search