subroutine char SEGMENTS char SEG_KICKS int MAP_REAL int MAP_IDEAL \
    int MAP_INDEX int i1 int ilast int THE_SLIT int THE_BORDER

! the variable that tells how many residues will be during the
! omit map calculation calculation is THE_SLIT

set vari i2 = i1 + THE_SLIT
set vari ii1 int = i1 - THE_BORDER
if ( ii1 .le. 0 ) set vari ii1 int = 1
key active sele by seq resi numb ii1  i2 + THE_BORDER \
	.a segm name SEG_KICKS end

make map MAP_REAL conv real
make map MAP_REAL zero 
make map MAP_REAL atom symm \
     atom dens \
     sele segm name SEG_KICKS .a .not active end  
 fourier map MAP_REAL 
 reflection select defined end r-values 

make map MAP_REAL zero
make map MAP_REAL conv complex
refl set ampl phase fwork = fcalc * 1.
refl set ampl scale fcalc     fwork = fobs * 1. - fcalc * 1.
refl fill-map MAP_REAL select defined end
four map MAP_REAL back 

make map MAP_INDEX from MAP_IDEAL ini 0 inte around 8 sele active end
make map MAP_INDEX from MAP_IDEAL atom index 0.0 1000. \
                  sele resi numb i1 i2 end

anal dens ( correl base-map = MAP_REAL ideal-map = MAP_IDEAL \
     index = MAP_INDEX ) \
     sele weight 0.1 1.1 .and. segm name SEGMENTS \
	.a resi numb i1  i2 end
set vari i1 = i1 + THE_SLIT
if ( i1 .gt. ilast ) return
rewind file