AMBER Archive (2006)

Subject: AMBER: the sequence command in leap

From: Shuting Wei (gtg549i_at_mail.gatech.edu)
Date: Wed Sep 20 2006 - 10:32:17 CDT

Dear Amber,

I tried to build a polymer chain in leap with sequence command, but i
got the following error:

> loadamberprep polyamine.prepin
Loading Prep file: ./polyamine.prepin
> loadamberparams polyamine.frcmod
Loading parameters: ./polyamine.frcmod
Reading force field mod type file (frcmod)
> set PAM tail PAM.1.12
> desc PAM
UNIT name: PAM
Head atom: .R<PAM 1>.A<C1 1>
Tail atom: .R<PAM 1>.A<N4 12>
Contents:
R<PAM 1>
> triPAM = sequence { PAM PAM PAM }
ERROR: Comparing atoms N4, H14, H16, H18 to atoms C2, N4, H18, H14
!FATAL ERROR----------------------------------------
!FATAL: In file [chirality.c], line 121
!FATAL: Message: Atom C2 is not in the first list
!
!ABORTING.

What's wrong with my input? Thanks a lot!

Shuting

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu