subroutine char SEGMENTS char METHOD 

if ( AUTO_SEL .eq. "HISTORY" ) key XX_KEY sele by resi atom numb $1 end
if ( AUTO_SEL .eq. "ACTIVE" ) key XX_KEY sele by resi active end
if ( AUTO_SEL .eq. "WORK" ) key XX_KEY sele  segm name SEGMENTS end

if ( METHOD .eq. ASIS ) then
set geom sele  XX_KEY .a .not atom name C N O end -
    group N dihe dist angle ! extend

end_if

if ( METHOD .eq. EXTEND ) then
set geom sele  XX_KEY .a .not atom name C N O end -
    group N dihe dist angle extend
end_if

if ( METHOD .eq. FIT ) then

set vari THE_MAP = nmaps + 1

make map THE_MAP from MAP_ACT  around 10 sele XX_KEY end init 9999 real
make map THE_MAP from MAP_ACT copy
make map THE_MAP sele segm name WORK_SEGM .a .n XX_KEY end atom ( mask 1.6 )
make map THE_MAP sele XX_KEY .a atom name C N O H HT* OT* OXT end atom ( mask 1.0 )
make map THE_MAP set 9990 99999 -999.0

 find side sele XX_KEY .a .not atom name C N O H HT* OT* OXT end -
     map THE_MAP contour MAP_CONT steps 30 - 
          root CA side CB dihe 10.0 dist 10.0 angle 30.0 poin 10.0

dele map THE_MAP

calc bond sele atom name SG .a XX_KEY end dist 2.5

end_if

 <?MAIN_CMDS:re_image.cmds WORK_SEGM
key XX_KEY drop
return

by resi atom numb $1
set vari IRES local int resi sele by resi atom numb $1 end

set vari IATOM local int atom sele resi numb IRES .a atom name CA end

if ( IATOM .ge. 1 ) then
set vari atom_CA  inte atom sele atom name CA .a resi numb IRES end
set vari atom_C   inte atom sele atom name C  .a resi numb IRES end
set vari atom_N   inte atom sele atom name N  .a resi numb IRES end
else
set vari IATOM local int sele atom numb $1 end
end_if


key FIX_GEOM sele resi name IRES .a .not atom name CA N O C H end

key atom_NEIG sele by bond atom numb IATOM .and. FIX_GEOM end
show key atom_NEIG

if ( IRESULT_4 .eq. 0 ) return

if ( IRESULT_4 .eq. 1 ) then
set vari atom_NEXT inte atom sele atom_NEIG end

set geom dihe atom_NEXT atom_OLD1 atom_OLD2 atom_OLD3 value 180.0 -
   select FIX_GEOM .a cut atom_OLD2 atom_OLD1 end
end_if

if ( IRESULT_4 .eq. 2 ) then
set vari atom_NEXT1 inte atom sele atom_NEIG end
key atom_NEIG sele atom_NEIG .a .not atom numb atom_NEXT1 end
set vari atom_NEXT2 inte atom sele atom_NEIG end
show dihe atom_OLD1 atom_OLD2 atom_NEXT1 atom_NEXT2 
if ( RESULT_0 .gt. 

set geom dihe atom_NEXT atom_OLD1 atom_OLD2 atom_OLD3 value 180.0 -
   select FIX_GEOM .a cut atom_OLD2 atom_OLD1 end
end_if

if ( IRESULT_4 .eq. 3 ) then
end_if

if ( IRESULT_4 .ge. 4 ) then
end_if

return



defin clas sele segm name SEGMENTS end by root CA
<MAIN_UTILS:set_weight_auto_loop.com xx

!subroutine char XX_KEY
set vari IRES local int resi sele XX_KEY end
if ( IRES .le. 0 ) return

set vari IROOT local int atom sele clas numb 1 .and. resi numb IRES \
	.and. XX_KEY end
if ( IROOT .gt. 0 ) set weight sele by resi atom numb IROOT end = 0.0

set vari IROOT local int atom sele clas numb 2 .and. resi numb IRES \
	.and. XX_KEY end
if ( IROOT .gt. 0 ) set weight sele by resi atom numb IROOT .and. \
	.not. atom name C N O .and. clas numb 2 200 end = 0.0


set vari IROOT local int atom sele clas numb 3 .and. resi numb IRES \
	.and. XX_KEY .and. .not atom name C N O H* end
if ( IROOT .gt. 0 ) set weight sele by resi atom numb IROOT .and. \
	.not. atom name C N O .and. clas numb 3 200 end = 0.0

set vari IROOT local int atom sele clas numb 4 .and. resi numb IRES \
	.and. XX_KEY .and. .not atom name C N O H* end
if ( IROOT .gt. 0 ) set weight sele by resi atom numb IROOT .and. \
	.not. atom name C N O .and. clas numb 4 200 end = 0.0

set vari IROOT local int atom sele clas numb 5 .and. resi numb IRES \
	.and. XX_KEY .and. .not atom name C N O H* end
if ( IROOT .gt. 0 ) set weight sele by resi atom numb IROOT .and. \
	.not. atom name C N O .and. clas numb 5 200 end = 0.0

set vari IROOT local int atom sele clas numb 6 .and. resi numb IRES \
	.and. XX_KEY .and. .not atom name C N O H* end
if ( IROOT .gt. 0 ) set weight sele by resi atom numb IROOT .and. \
	.not. atom name C N O .and. clas numb 6 200 end = 0.0

set vari TEST int resid select resid ring IRES .a weight 0.1 1.1 end
if ( TEST .ge. 2 ) set weight sele resi numb IRES end = 1.0

key XX_KEY sele XX_KEY .and. .not. residue number IRES end
rewind file



return