subroutine char SEGMENTS real KICK int MAP_REAL int MAP_IDEAL

delete atom sele segm name #* end

symm sele segm name SEGMENTS end matr 1 tran 0 0 0
copy weight sele segm name SEGMENTS end sele segm name #* end
copy temp sele segm name SEGMENTS end sele segm name #* end
sh temp sele segm name SEGMENTS .a weight 0.1 1.1 end
!set temp sele segm name #* end = RESULT_0
<?MAIN_CMDS:re_phase.cmds "#*" MAP_IDEAL

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

set coor sele segm name #* end kick KICK

anal resi init

set vari MAP_INDEX = nmaps + 1

make map MAP_REAL conv real
make map MAP_REAL zero 
make map MAP_REAL atom symm \
     atom dens \
     sele segm name #* .a .not by sequ 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 -1.0 1000. \
                  sele active 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 active end

delete atom sele segm name #* end

sh anal sele active end
dele map nmaps
return


set vari MAP_INDEX = 3
set vari MAP_IDEAL = 2
make map 4 from MAP_INDEX real ini 9999 copy
make map set 1 10000 9999.
make map set -1 1 -9999
make map from MAP_IDEAL copy