init
read seq protein.seq
read prot bbAll01.prot

update atom average HA
update atom average HA1
update atom average HA2
update atom average HB
update atom average HB2
update atom average HB3
update atom average QB

add GFTatom C H

write prot habcabconhI1.prot

simulate 3D N CA+HA  -1 1 random
simulate 3D N CA+HA1 -1 1 random
simulate 3D N CA+HA2 -1 1 random

simulate 3D N CA-HA  -1 2 random
simulate 3D N CA-HA1 -1 2 random
simulate 3D N CA-HA2 -1 2 random

simulate 3D N CB+HB  -1 3 random
simulate 3D N CB+HB2 -1 3 random
simulate 3D N CB+HB3 -1 3 random
simulate 3D N CB+QB  -1 3 random

simulate 3D N CB-HB  -1 4 random
simulate 3D N CB-HB2 -1 4 random
simulate 3D N CB-HB3 -1 4 random
simulate 3D N CB-QB  -1 4 random

write peaks habcabconhI1.peaks