subroutine int THE_MAP real CUT_OFF
key XX_KEY sele active end
set vari CUSHION int = 10


!show vari CUT_OFF

set vari CUT_UP = 999.
if ( CUT_OFF .lt. 0. ) then
set vari CUT_OFF = CUT_OFF * -1.0
set vari CUT_UP = -999
end_if

show vari CUT_OFF

dele atom sele segm name #CNT end
key old sele all end

key ii_list sele active .a .n atom name H* .a map THE_MAP CUT_OFF CUT_UP end 
make point from atom sele ii_list end 
make atom from point ( use_map THE_MAP )
rena segm #CNT sele .n old end
write over file atom_center_list.pdb sele ii_list end coor pdb
key ii_list drop
dele atom sele segm name #CNT end

return

key old sele all end

make resi sele segm name CNTS end root CEX
dele atom sele segm name CNTS .a old end
rena seq A sele .n old end auto 0
image col 220 sele segm name CNTS end atom cros
set vari IATOM int atom sele segm name CNTS end
show segm 

ima hist pick IATOM
set vari NEXT_RESI global int resi select atom numb $1 end
image center sele resi numb NEXT_RESI end 
image map over