subroutine char SEGMENTS int MAP_FROM real ATOM_RAD int MAP_MASK char FILE_MASK

dele atom sele segm name #* end

key other sele .not segm name SEGMENTS .a segm name WORK_SEGM end

symm sele .not segm name SEGMENTS end around \
     sele segm name SEGMENTS end dist ATOM_RAD + 1. cut
symm sele segm name SEGMENTS end around sele segm name SEGMENTS end \
  dist ATOM_RAD + 1. cut

set rad sele segm name SEGMENTS end sele segm name #* .or other end \
	overlap 1.5 ATOM_RAD

 make map MAP_MASK from MAP_FROM init -9999  \
     around 10 sele segm name SEGMENTS end  real
 make map MAP_MASK atom ( radii mask ) select  segm name SEGMENTS end

 make map MAP_MASK + 1 from MAP_MASK init 9999 copy
make map MAP_MASK from MAP_MASK + 1 fill MAP_MASK

write over file FILE_MASK map MAP_MASK xplor
key other drop
return