subroutine char SEGMENTS delete atom select segment name #* end set vari ATOM_SPHERE real = 20.0 key old select all end symm select segment name MOLA end around \ center atom $1 dist ATOM_SPHERE neigh 2 key SYMM_M sele .not old end image select SYMM_M end col 105 bond image select SYMM_M .a .not by bond all end col 105 atom cros key old select all end symm select segm name WAT end around center atom $1 dist ATOM_SPHERE neigh 2 key SYMM_WAT select .not old end image select SYMM_WAT .and. atom name O* end col 125 atom cros return