subroutine char SEGMENTS
set vari SUB_PASS ?= 0


set vari TEST_ATOM int = 0
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 #* 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

if ( SUB_PASS .gt. 0 ) then
key XX_KEY sele segm name HELP end
end_if

! rena segm HELP sele XX_KEY .o protein end
key OLD sele all end

 find residue sele XX_KEY end map MAP_ACT root C0 side C1 \
    list "SP2 SP3 SP3R" fit test TEST_ATOM level 2

dele atom sele XX_KEY end
! dele atom sele segm name HELP end
! dele atom sele segm name HELP .a .not protein end
! rena segm OLD sele segm name HELP end

 find residue sele segm name NEW end sele segm name NEW end \
 root C0 map MAP_ACT merge step WALK_STEP del
! ima cent coor   14.6510    4.7400   21.3870

 set weight sele segm name NEW end = 1.0

key XX_KEY drop
rena segm HELP sele segm name NEW end
! set vari WORK_SEGM char glob segm sele segm name HELP .o protein end
 set vari WORK_SEGM char glob = HELP


<?MAIN_CMDS:use_atoms.cmds WORK

 dele atom sele segm name HELP .a map MAP_ACT -1000. MAP_CONT - 0.2 end
 delete atom sele segm name WORK_SEGM .a .not by bond all end

 <DEF_ALL WORK_SEGM
 mini sele segm name WORK_SEGM end  step 80 writ 10
 dele atom sele segm name HELP .a map MAP_ACT -1000. MAP_CONT - 0.2 end
< ?MAIN_CMDS:re_image.cmds WORK_SEGM

set vari SUB_PASS = SUB_PASS + 1
ener impr on
if ( SUB_PASS .eq. 2 ) then
ener impr off
end_if
if ( SUB_PASS .le. 3 .and. ener_vdw .gt. 2999.0 ) rewind file

set vari II int atom sele skeleton .o segm name OLD end
 if ( II .le. 0 ) then
 delete atom sele segm name OLD .o skeleton end
 return
 end_if

 return
! 1379 1403 1712
! 2248 2249 2236

make map MAP_TARG + 1 from MAP_TARG integer init 0

make atom walk \
  ( \
    method 4 \
    sele segm name HELP end \
    sele skeleton end \
    map MAP_TARG \
    range MAP_CONT 1000. \ 
    lev 10 \ 
    step 1.42 \
    test TEST_ATOM \
    del \
  )
set col 
sele skeleton .a weigh -0.1 0.1 end = 200
sele skeleton .a weigh 0.1 2.2 end = 120
exit

image set sele skeleton end bond

return

image col 220 sele segm name HELP end bond

image col 160 sele atom name XEXT end atom cros
return

 key OLD sele all end
 clo uni 4
 read file fort.4 coor
 image col 200 sele .not OLD end atom cros
 clo uni 4
 delete atom sele .not OLD end
 
 clo uni 3
 read file fort.3 coor
 image col 150 sele .not OLD end atom cros
 clo uni 3
 delete atom sele .not OLD end
 
 clo uni 2
 read file fort.2 coor
 image col 120 sele .not OLD end atom cros
 clo uni 2
 delete atom sele .not OLD end
 
 clo uni 1
 read file fort.1 coor
 image col 85 sele .not OLD end atom cros
 clo uni 1
 delete atom sele .not OLD end

return

key OLD sele all end
clo uni 4
read file fort.4 coor
image col 200 sele .not OLD end atom cros
clo uni 4
delete atom sele .not OLD end

clo uni 3
read file fort.3 coor
image col 150 sele .not OLD end atom cros
clo uni 3
delete atom sele .not OLD end

clo uni 2
read file fort.2 coor
image col 140 sele .not OLD end atom cros
clo uni 2
delete atom sele .not OLD end

clo uni 1
read file fort.1 coor
image col 100 sele .not OLD end atom cros
clo uni 1
delete atom sele .not OLD end

key OLD sele all end
clo uni 9
read file fort.9 coor
image col 220 sele .not OLD end atom cros
clo uni 9
delete atom sele .not OLD end

key OLD sele all end
clo uni 8
read file fort.8 coor
image col 165 sele .not OLD end atom cros
clo uni 8
delete atom sele .not OLD end

key OLD sele all end
clo uni 11
read file fort.11 coor
image col 114 sele .not OLD end atom cros
clo uni 11
delete atom sele .not OLD end


return

image col 110 sele segm name MOL1 CLUS end bond
key image sele image .a .not segm name MOL1 CLUS end 

return

!image col 220 sele .not OLD .a weight 0.01 1000. end bond
!image col 200 sele .not OLD end bond
image col 220 sele .not OLD end atom name
image col 220 sele .not by bond all .a .not OLD end atom cros