! subroutine char SEGMENTS
 delete atom select segment name #* NEW* end
! show vari SEGMENTS
set vari DIST = 40.

 key old select all end
 symm select active end cut  dele on dist DIST neigh 2 around sele active end
 image select .not old  end col 105 bond ! atom cross

return


<amore_shifts SEGMENTS NEW1 349.0    2.5  187.4   49.3   50.9   -5.6
 key old select all end
 symm select segm name NEW1 end around center atom $1 dist DIST neigh 2
 image select .not old  end col 125 bond ! atom cross

<amore_shifts SEGMENTS NEW2 44.0    5.0  129.7   79.4   81.1    8.7
 key old select all end
 symm select segm name NEW2 end around center atom $1 dist DIST neigh 2
 image select .not old  end col 145 bond ! atom cross

return

 symm select all end around center coor 30 30 20 dist 30. neigh 2
 key SYMM_M select .not old end
 image select SYMM_M end col 105 bond ! atom cross
 

return

 key old select all end
 symm select segm name SOL end around center atom $1 dist 60. neigh 2
 key SYMM_WAT select .not old end
 image select SYMM_WAT .and. atom name O* end col 125 atom cros
 return