subroutine char SEGMENTS
  image erase
 delete atom sele segm name #* end
 key image sele .not all end
 
 set color
 sele segm name SEGMENTS end = 200
 sele segm name SEGMENTS .a weigh 0.1 1.1 end = 160
 sele segm name SEGMENTS .a \ 
     atom name H* .a by bond weigh 0.1 1.1 end = 160
sele segm name WAT2 end = 160 
exit


image sele segment name SEGMENTS end set bond
image sele segment name SEGMENTS .a .not by bond all end atom cros
image sele segment name SEGMENTS end col 85 hbond
image sele t_anchor end col 200 atom cros
image sele t_anchor .o by pair t_anchor end col 200 bond
return