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

wr sele segm name HYBR end res seq

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

image erase
key image sele .n all end

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

<?MAIN_CMDS:use_atoms.cmds WORK

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

ima col 160 sele segm name HELP end bond


ener vdw off elec off
mini sele segm name HELP end
ener vdw on elec on


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


set weight sele segm name HELP end = 1.0
< ?MAIN_CMDS:re_image.cmds  HELP
! 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
set vari TRACE_LAST global = TRC_ALA

return