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 NEW #* HELP end

set vari MAP_TARG = MAP_ACT
set vari WALK_STEP = 1.42 * 1.1

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

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

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

set vari WORK_SEGM global = NEW 

<DEF_ALL WORK_SEGM
ener vdw off elec off
! mini sele segm name WORK_SEGM end  step 80 writ 10
ener vdw on elec on

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

key image sele segm name NEW HELP SKEL end
show atom
show resi
show segm

trace main sele segm name NEW .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 NEW .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