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 = 150
sele segm name SEGMENTS .a \ 
     atom name H* .a by bond weigh 0.1 1.1 end = 150
 sele segm name MOLA .a weigh 0.1 1.1 end = 150
 sele segm name MOLA .and. \ 
    atom name H* .a by bond weigh 0.1 1.1 end = 150
 sele segm name MOLB .a weigh 0.1 1.1 end = 155
 sele segm name MOLB .and. \ 
    atom name H* .a by bond weigh 0.1 1.1 end = 155
 sele segm name IA .a weigh 0.1 1.1 end = 160
 sele segm name IA .and. \ 
    atom name H* .a by bond weigh 0.1 1.1 end = 160
 sele segm name IB .a weigh 0.1 1.1 end = 165
 sele segm name IB .and. \ 
    atom name H* .a by bond weigh 0.1 1.1 end = 165
sele segm name WAT2 end = 170 
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