* make psf and build coor from ic in the rtf
*
open unit 1 card read name "top_all36_carb_NHSHB3.rtf"
read RTF card unit 1
close unit 1


open unit 1 card flex read name "par_all36_carb.prm"
read PARA card flex unit 1
close unit 1
read sequ card
* sequence
*
2
BGLCA AGLCNA

gener S setup warn
read coor univ unit 1 offset 0
patch 14ab S 1 S 2 setup
autogen angl dihe
patch amin S 2 setup

AUTOGENereate ANGLes DIHEdrals

ic param
ic seed 1 C1 1 C2 1 C3
ic build

print coor

write psf card name disa.psf unit 1
* PSF
*

write coor pdb name disa.pdb unit 1
* PDB Coords
*

write coor card name disa.crd unit 1
* Coords
*

stop