subroutine char SEGMENTS
set vari SUB_PASS ?= 0

set vari TEST_ATOM int =   -1
if ( nathis .gt. 0 ) then
set vari TEST_ATOM int atoms sele atom numb $1 end
end_if
show vari TEST_ATOM

dele atom sele segm name HYBR HY_* #* HELP end

set vari MAP_TARG = MAP_ACT

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 MAP_WORK = nmaps + 1
make map MAP_WORK from MAP_ACT integ around 8 sele XX_KEY end init 0
if ( MAP_WORK .ne. nmaps ) return

trace hybrid sele XX_KEY end level 20 residue "SP3 SP3R" \
 map MAP_ACT work_map MAP_WORK ! test TEST_ATOM

rena segm HYBR sele segm name NEW end

set col
sele segm name SKEL end = 220
sele temp 0.1 399 .a segm name SKEL end = 160
sele temp 400 12999 .a segm name SKEL end = 200
sele weigh 0.1 1.2 .a segm name SKEL end = 140
exit

image erase
key image sele .n all end
 ima set sele segm name SKEL end bond

  clo uni 1
  key old sele all end
  read unit 1 coor
  ima sele .n old end col 110 atom cros
  clo uni 2
  key old sele all end
  read unit 2 coor
  ima sele .n old end col 160 atom cros
  clo uni 3
  key old sele all end
  read unit 3 coor
  ima sele .n old end col 140 atom cros
  clo uni 4
  key old sele all end
  read unit 4 coor
  ima sele .n old end col 200 atom cros
  clo uni 9
  key old sele all end
  read unit 9 coor
  ima sele .n old end col 158 atom cros


! image from

set vari WORK_SEGM global = HYBR 
<?MAIN_CMDS:use_atoms.cmds WORK

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

! ima col 200 sele segm name HYBR .a atom name C0 end bond
ima col 110 sele segm name HYBR end col 220 bond stick
! ima col 91 sele segm name HYBR .a atom name C0 end bond dist
! ima col 200 sele segm name HYBR .a .n atom name C0 end atom name
! ima col 1 sele segm name HYBR .a atom name C0 end atom numb

key image sele segm name HYBR end 
show atom
show resi
show segm
set vari TRACE_LAST global = TRC_HYBR

delete map MAP_WORK
return


trace main sele segm name HYBR .a resi name SP3 end \
     in_atoms "C1 C2 C0 C3" \
     out_atoms "N C CA CB"  \
     residue XXX  forward segment HELP test TEST_ATOM

trace main sele segm name HYBR .a resi name SP3R end \
     in_atoms "C2 C1 C0 C3" \
     out_atoms "N C CA CB"  \
     residue XXX  backward segment HELP test TEST_ATOM

show atom
show resi
show segm

key image sele .n all end
set vari WORK_SEGM global = HELP
<DEF_ALL WORK_SEGM
  set geometry sele segm name HELP end peptide
set weight sele segm name HELP end =1.0
mini sele segm name HELP end
find pept sele segm name HELP end \ 
          step 60 \ 
          sec 5.1 \ 
          map MAP_ACT 8.8
mini sele segm name HELP end
find pept sele segm name HELP end \ 
          step 60 \ 
          sec 5.1 \ 
          map MAP_ACT 8.8
mini sele segm name HELP end


ima col 160 sele segm name HELP end bond
key image sele image .o segm name SKEL end
! ima col 1 sele segm name HELP .a atom name CA end atom numb

return