subroutine char SEGMENTS
set vari TEST_ATOM int = -3666
! image hist pick 155
if ( nathis .gt. 0 ) then
set vari TEST_ATOM int atoms sele atom numb $1 end
end_if

show vari TEST_ATOM
show vari AUTO_SEL
show vari TRACE_LAST
if ( AUTO_SEL .eq. "ACTIVE" .and. TRACE_LAST .eq. "TRC_HYBR" ) then
 set vari SEGMENTS = HYBR
<?MAIN_CMDS:use_atoms.cmds WORK

echo "SYMM OVERLAPS" ROOT_ATOM SEGMENTS

end_if


if ( AUTO_SEL .eq. "HISTORY" ) key XX_KEY sele history end
if ( AUTO_SEL .eq. "ACTIVE" ) key XX_KEY sele active end
if ( AUTO_SEL .eq. "WORK" ) key XX_KEY sele  segm name SEGMENTS end


set vari IROOT_CO inte atom sele XX_KEY .a atom name C0 end
set vari IROOT_CA inte atom sele XX_KEY .a atom name CA end
set vari ISIDE_C3 inte atom sele XX_KEY .a atom name C3 end
set vari ISIDE_CB inte atom sele XX_KEY .a atom name CB end

dele bond sele atom name IROOT_C0 end

set vari ROOT_ATOM = XXXX
set vari SIDE_ATOM = XXXX

! show vari IROOT_CO
! show vari IROOT_CA

if ( IROOT_C0 .gt. 0 ) set vari ROOT_ATOM = C0
if ( IROOT_CA .gt. 0 ) set vari ROOT_ATOM = CA
if ( ISIDE_C3 .gt. 0 ) set vari SIDE_ATOM = C3
if ( ISIDE_CB .gt. 0 ) set vari SIDE_ATOM = CB

echo " ATOMS: ROOT:" ROOT_ATOM "SIDE:" SIDE_ATOM
key OLD sele all end

set vari MAP_ONE = nmaps + 1
set vari MAP_TWO = MAP_ONE + 1
make map MAP_ONE from 1 init 0 integer  cell
make map MAP_TWO from MAP_ONE init 0 integer 
!  find symm sele active end map MAP_ONE \
!   resid by center calc sele segm name CLUS end
   if ( TEST_ATOM .gt. 0 ) then
     find symm sele XX_KEY end map MAP_ONE \
       resid by sele neigh $1 -1 end
   else
     find symm sele XX_KEY end map MAP_ONE
   end_if


! ! key XX_KEY sele segm name HELP end
! ! find symm sele XX_KEY end map MAP_ONE \
! !        resid by sele segm name A B C D end

set col
sele XX_KEY end col 200
sele XX_KEY .a weigh 0.1 1.1 end col 160
exit

ima sele XX_KEY end set bond

!    center calc sele segm name SYMM_SEGM end 

find residue merge \
     ( second 1.0 alpha 2.0 beta 1.0 whole test TEST_ATOM ) \
     sele XX_KEY end \
     map MAP_ACT work MAP_ONE \
     root ROOT_ATOM side SIDE_ATOM 

 rena segm HY_2 sele XX_KEY end
 rena segm HY_O sele segm name HYBR end
 rena segm HYBR sele segm name SYM end
 set weight sele segm name HYBR end = 1.0

 set vari WORK_SEGM global = HYBR


! <DEF_ALL HYBR
!  mini sele segm name HYBR end  step 80 writ 10
ima era

<?MAIN_CMDS:use_atoms.cmds WORK

! < ?MAIN_CMDS:re_image.cmds WORK_SEGM
set vari TRACE_LAST global = TRC_SYMM

! return

set col
sele weigh -100. 0.1 end = 220
sele weigh 0.1 99 end = 100
sele weigh 1.1 99 end = 105
sele weigh 2.1 99 end = 110
sele weigh 3.1 99 end = 115
sele weigh 4.1 99 end = 120
sele weigh 5.1 99 end = 125

sele weigh 101.1 999 end = 150
sele weigh 102.1 999 end = 155
sele weigh 103.1 999 end = 160
sele weigh 104.1 999 end = 165
sele weigh 99000 99999 end = 160
sele weigh 89000 98999 end = 140
sele weigh 79000 88999 end = 120
sele weigh 19000 78999 end = 90
exit

key image sele .n all end
make bond from resi ROOT_ATOM sele XX_KEY end
dele bond sele atom name ROOT_ATOM .a bond range 4.98 1000. end
image set sele atom name ROOT_ATOM .a XX_KEY end bond


! ima sele atom name ROOT_ATOM end set atom cros
! set rad sele all end = 0.34
! image sele atom name ROOT_ATOM .a .n by bond atom name ROOT_ATOM end \
! col 220 atom ball
! ima col 220  pair stick
key image sele .n all end
ima set  pair stick
dele bond sele atom name ROOT_ATOM end
ima sele segm name HYBR end col 160 bond

dele map MAP_ONE

return

ima sele segm name HY_2 .a atom name C3 .a weig 0.1 1.1 end \
   col 160 seq
ima sele segm name HY_2 .a atom name C1 .a weig 0.1 1.1 end \
   col 160 atom ball
ima sele segm name HYBR .a atom name C2 .a weig 0.1 1.1 end \
   col 220 atom ball


return

set vari II int atom sele protein end
if ( II .gt. 0 ) then
key SYMO_TEST sele segm name SEGMENTS .a .not protein end
dele atom sele TEST .a around key protein dist 1.5 end
key SYMO_TEST drop
end_if
return