subroutine inte MAP_ORIG real SKEL_CUT
set vari MAP_DENS = nmaps + 1
set vari MAP_SKEL = MAP_DENS + 1

make map MAP_DENS from MAP_ORIG init 9999 around 8 sele active end
make map MAP_DENS from MAP_ORIG copy

delete atom sele SKELETON .or. segm name #* end
key old sele all end

make map MAP_SKEL from MAP_DENS init 0 integer
make point init from map MAP_DENS extreme SKEL_CUT 100.
make map MAP_SKEL from MAP_DENS local SKEL_CUT 100.

make map MAP_SKEL cloud
make point from map MAP_DENS saddle MAP_SKEL

make map MAP_SKEL from MAP_DENS init 0 integer
make map MAP_SKEL from MAP_DENS local SKEL_CUT 100.
make point from map MAP_DENS trace MAP_SKEL

make atom from point reindex
key SKELETON sele .not old end

dele atom sele .not by bond all .and SKELETON end
make segm from atom sele SKELETON end

!dele bond sele bond range 10. 1000. .a SKELETON end

set temp sele SKELETON end = 50.
set weigh sele SKELETON end = 1.

delete map MAP_DENS

image sele SKELETON end col 225  bond
return

set col 
sele SKELETON .a segm range 10 100 end = 201
sele SKELETON .a segm range 100 10000 end = 161
exit

return

 ima sele atom name XE* .a SKELETON end col 110 atom cros
 ima sele atom name XS* .a SKELETON end col 160 atom cros
 ima sele atom name XT* .a SKELETON end col 220 atom cros
 ima sele atom name XE* .a SKELETON end col 110 atom cros
 ima sele atom name XS* .a SKELETON end col 160 atom cros

!image col 220 sele SKELETON .a segm rang 10 100000 end bond

return