subroutine char SEGMENTS
calculate bond neigh select active end
calculate bond distance 2.2 disulfide select active end
return
calculate bond select atom name C N .a segm name SEGMENTS end neighbor
calculate bond select atom name P O3' .a segm name SEGMENTS end neighbor