subroutine char SEGMENTS
calculate bond neigh select segm name SEGMENTS .a active end
calculate bond distance 2.1 disulfide select segm name SEGMENTS .a 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