subroutine char SEGMENTS
set vari SUB_PASS ?= 0

set vari TEST_ATOM int = 3666
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

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

 symm sele segm name HELP .a ctab 3 100 end matr 1 tran 0 0 0


 make ct ini sele segm name #* end
 rena atom C0 sele segm name #* end
 make segm sele segm name #* end from atoms
 rena segm NEW sele .not OLD end 
 rena resi SP2 sele segm name NEW end


 set vari WORK_SEGM char glob = NEW
< ?MAIN_UTILS:fill_atoms.com NEW

set weight sele segm name NEW end = 1.0

set clas sele segm name NEW .a .not clas name SP2 SP3 end = HELP

 find residue sele segm name NEW end sele segm name HELP .a ctab 3 1000 \
 end \ 
 map MAP_ACT root C0 side C1 list "SP2 SP3 SP3R" fit 

 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
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