init
read seq nhsqc.seq
read prot hncoO1.prot

add GFTatom CA CB
add GFTatom CA CA

write prot bbgftI1.prot

simulate 3D N CA+CA   0 1 random
simulate 3D N CA+CA 200 3 random # change the value (200) according the gap between two SRDs
simulate 3D N CA+CB   0 1 random
simulate 3D N CA+CB 200 3 random # change the value (200) according the gap between two SRDs
simulate 3D N CA-CA   0 2 random
simulate 3D N CA-CA 200 4 random # change the value (200) according the gap between two SRDs
simulate 3D N CA-CB   0 2 random
simulate 3D N CA-CB 200 4 random # change the value (200) according the gap between two SRDs

write peaks bbgftI1.peaks