init
read seq ./myprot.seq
read prot noe.prot

# generate cosy virtual atoms
add GFTatom COSY_C ATTACHED_H

write prot ./hcchAroI1.prot 


simulate 3D CG  COSY_H  0 1 
simulate 3D CG  COSY_C   0 1
simulate 3D CG1 COSY_CH  0 1 
simulate 3D CG1 COSY_C   0 1
simulate 3D CG2 COSY_CH  0 1 
simulate 3D CG2 COSY_C   0 1

simulate 3D CD  COSY_CH  0 2 
simulate 3D CD  COSY_C   0 2
simulate 3D CD1 COSY_CH  0 2 
simulate 3D CD1 COSY_C   0 2
simulate 3D CD2 COSY_CH  0 2 
simulate 3D CD2 COSY_C   0 2

simulate 3D CZ2  COSY_CH  0 3 
simulate 3D CZ2  COSY_C   0 3
simulate 3D CZ3  COSY_CH  0 3 
simulate 3D CZ3  COSY_C   0 3
simulate 3D CE3  COSY_CH  0 4 
simulate 3D CE3  COSY_C   0 4
simulate 3D CH2  COSY_CH  0 4 
simulate 3D CH2  COSY_C   0 4


write peaks ./hcchAroI1.peaks

