[AMBER] How to run LEaP in a shell file?

From: SongLin <songlin.moon.ibp.ac.cn>
Date: Mon, 14 Dec 2009 15:10:00 +0800

Dear professor,

I want to build PDB formatted glycan files using LEaP automatically, so
I create several shell files like this:

cd ~/amber10/bin/
./tleap << !
source leaprc.GLYCAM_06
glycan=sequence{glycan 3VA}
glycan=sequence{glycan 3LB}
glycan=sequence{glycan 0SA}
savepdb glycan glycan.pdb

But when I run this file it returns:
*** Error: tl_getline(): not interactive, use stdio.

But it's ok when I run these commands line by line. Can you tell me how
to run LEaP automatically? (I can create LEaP command lines
automatically, but I can't run them automatically)

Thank you very much!

Yours Sincerely,
Lin Song

Center for Computational and Systems Biology,
Institute of Biophysics, CAS.
15 Datun Road, Chaoyang District
Beijing, 100101

AMBER mailing list
Received on Sun Dec 13 2009 - 23:30:02 PST
Custom Search